Комментарии 9
Похоже, остатки в вашей линейной модели гетероскедастичны:
plot(hubble$D, resid(fit), pch=21, col="black", bg="lightblue", cex=.9)
![](https://habrastorage.org/getpro/habr/upload_files/5a2/6b9/c9a/5a26b9c9a534e30fb2c03a597f4173cd.png)
последние оценки постоянной Хаббла на основе реликтового излучения показывают 67,4, а на основе сверхновых в удаленных галактиках - 7
кстати а знаете, что Постоянная Хаббла - не является постоянной? Вот думаю этим и будут объясняться отличия, в одном из методов оценки не учитывается какой-то фактор, о котором нам пока не известно
Свободный член уравнения приравнен к 0, поскольку состояние Вселенной в начальный момент Большого Взрыва предполагает, что галактик еще не существовало и они не могли удаляться друг от друга
Рекомендую обратиться к более ранней публикации на эту тему: https://perfdynamics.blogspot.com/2010/06/linear-modeling-in-r-and-hubble-bubble.html
Hubble took this angular effect into account with the three additional X, Y, Z terms in his linear model:
r K + X cos(α) cos(δ) + Y sin (α) cos(δ) + Z sin(δ) = vwhich is unnumbered and appears about half way between Table 1 and Table 2 in his paper
Угловые эффекты - вот суть проблемы!
В таких расчетах очень желательно еще погрешность величины вычислять...
Вычисляем возраст Вселенной в R