Настройка степеней свободы для коррекции сферичности при повторном измерении Anova из автомобильного пакета

Я выполняю повторные измерения ANOVA с использованием пакета car. Что отлично работает и возвращает вывод, подобный этому:

Univariate Type III Repeated-Measures ANOVA Assuming Sphericity

                SS num Df Error SS den Df        F    Pr(>F)    
(Intercept) 7260.0      1   603.33     15 180.4972 9.100e-10 ***
phase        167.5      2   169.17     30  14.8522 3.286e-05 ***
hour         106.3      4    73.71     60  21.6309 4.360e-11 ***
phase:hour    11.1      8   122.92    120   1.3525    0.2245    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 


Mauchly Tests for Sphericity

        Test statistic  p-value
phase             0.70470 0.086304
hour              0.11516 0.000718
phase:hour        0.01139 0.027376


Greenhouse-Geisser and Huynh-Feldt Corrections
for Departure from Sphericity

            GG eps Pr(>F[GG])    
phase      0.77202  0.0001891 ***
hour       0.49842  1.578e-06 ***
phase:hour 0.51297  0.2602357    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Теперь он показывает мне, что в некоторых случаях необходимо применять поправку на сферичность. Насколько я понимаю, эта поправка влияет не только на p-значения, но и на степени свободы (df). Однако вывод этого не показывает. Итак, как я могу отобразить скорректированный df?


r anova r-car
person JereB    schedule 25.07.2014    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Возможно, вы уже разобрались с этим, но чтобы скорректировать степени свободы в случаях поправок Гринхауза-Гейссера или Хюина-Фельдта, вы просто умножаете каждую степень свободы на соответствующее значение эпсилон. Вот пример, основанный на вашем результате:

# degrees of freedom for phase:hour
num_Df <- 8
den_Df <- 120

# Greenhouse-Geisser epsilon
gg_e <- 0.51297

# adjusted degrees of freedom
num_Df_adj <- gg_e * num_Df 
den_Df_adj <- gg_e * den_Df
person xilliam    schedule 07.02.2020