Цитата:
Сообщение от Hogfather
Писать много лень, да и некогда.
|
Вот. На автокорреляцию не проверил, за что и был наказан. Если добавить в качестве регрессора предыдущее значение веса тушки, то полином становится не 4, а 2 степени и фаза луны также оказывается лишней. В общем, чудеса разоблачены, астрология -- лженаука.
Еще игры в R
Код:
> summary(lm1<-lm(weight~poly(day,2)+weigth1,data=MyW))
Call:
lm(formula = weight ~ poly(day, 2) + weigth1, data = MyW)
Residuals:
Min 1Q Median 3Q Max
-1.8507 -0.4794 0.0044 0.4653 3.3179
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 46.25645 7.21014 6.415 1.84e-09 ***
poly(day, 2)1 -42.27425 6.75988 -6.254 4.19e-09 ***
poly(day, 2)2 2.56471 0.90183 2.844 0.0051 **
weigth1 0.59624 0.06264 9.518 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.7593 on 146 degrees of freedom
(1 observation deleted due to missingness)
Multiple R-squared: 0.9925, Adjusted R-squared: 0.9923
F-statistic: 6428 on 3 and 146 DF, p-value: < 2.2e-16
> AIC(lm1)
[1] 349.0238
> AIC(lm2<-lm(weight~poly(day,4)+mds,data=MyW))
[1] 401.9843
2.5 % 97.5 %
(Intercept) 32.0067225 60.5061783
poly(day, 2)1 -55.6341040 -28.9143911
poly(day, 2)2 0.7823834 4.3470294
weigth1 0.4724352 0.7200359
Как видим, вторая модель лучше. Гипотезы о зависимости от дня недели и фазы луны опровергнуты, все зависит только от человека и случайной составляющей.
UPD: У женщин фаза луны может играть роль. Но это требует весовых замеров женской тушки.