Тема: Диета...
Показать сообщение отдельно
Старый 17.07.2013, 21:48   #554
Hogfather
Platinum Member
 
Аватар для Hogfather
 
Регистрация: 22.07.2010
Адрес: Санкт-Петербург
Сообщений: 3,304
По умолчанию

Цитата:
Сообщение от 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: У женщин фаза луны может играть роль. Но это требует весовых замеров женской тушки.
---------
DNF is not an option
Hogfather вне форума   Ответить с цитированием
Реклама