Интерпретация кодов значимости R для таблицы ANOVA?

кадр тестовых данных:

> foo
      x     y     z
1 0.191 0.324 0.620
2 0.229 0.302 0.648
3 0.191 0.351 0.626
4 0.229 0.324 0.630
5 0.152 0.374 0.656
6 0.191 0.295 0.609
7 0.229 0.267 0.665
8 0.152 0.353 0.657
9 0.152 0.355 0.655

Две линейные модели:

model1 <- lm(z~polym(x,y,degree = 1),data=foo)
model2 <- lm(z~polym(x,y,degree = 2),data=foo)

ANOVA для двух моделей возвращает:

> anova(model1,model2)
Analysis of Variance Table

Model 1: z ~ polym(x, y, degree = 1)
Model 2: z ~ polym(x, y, degree = 2)
  Res.Df      RSS Df Sum of Sq    F Pr(>F)  
1      6 0.002988                           
2      3 0.000169  3   0.00282 16.6  0.023 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Почему сингл *? 0,05 > 0,023 > 0,01, поэтому не должен ли он печатать символ .?


r regression anova
person DeltaIV    schedule 02.08.2016    source источник
comment
Табличные коды таблицы могут быть немного трудными для расшифровки, но: 0,023 является значимым на уровне 0,05 (*), но не на уровне 0,01 (**). Символ (.) соответствует уровню 0,1 или десяти процентов.   -  person lmo    schedule 02.08.2016
comment
0 ‘***’/0.001 ‘**’/0.01 ‘*’/0.05 ‘.’/0.1 ‘ ’/ и за 1 ничего не печатается. Ваш 0,023 находится между 0,01 и 0,05, поэтому он отмечен одной звездочкой.   -  person Roman    schedule 02.08.2016

Ответы (1)


arrow_upward
13
arrow_downward

Нет ничего плохого.

0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

означает:

annotation          p-value          significance level
   ***            [0, 0.001]                0.001
    **         (0.001, 0.01]                0.01
     *          (0.01, 0.05]                0.05
     .           (0.05, 0.1]                0.1
                    (0.1, 1]                1

0,023 находится в пределах (0.01, 0.05], поэтому он должен быть аннотирован *.

person Zheyuan Li    schedule 02.08.2016