重复测量方差分析 - 自由度

Question

I'm studying repeated measures ANOVAs, found this example online and can't understand how are the df calculated. What would be the mathematical expression for calculating both df values in this example?

data(obk.long, package = "afex")

# estimate mixed ANOVA on the full design:
aov_car(value ~ treatment * gender + Error(id/(phase*hour)), 
        data = obk.long, observed = "gender")


# the three calls return the same ANOVA table:
# Anova Table (Type 3 tests)
# 
# Response: value
#                         Effect          df   MSE         F  ges p.value
# 1                    treatment       2, 10 22.81    3.94 + .198    .055
# 2                       gender       1, 10 22.81    3.66 + .115    .085
# 3             treatment:gender       2, 10 22.81      2.86 .179    .104
# 4                        phase 1.60, 15.99  5.02 16.13 *** .151   <.001
# 5              treatment:phase 3.20, 15.99  5.02    4.85 * .097    .013
# 6                 gender:phase 1.60, 15.99  5.02      0.28 .003    .709
# 7       treatment:gender:phase 3.20, 15.99  5.02      0.64 .014    .612
# 8                         hour 1.84, 18.41  3.39 16.69 *** .125   <.001
# 9               treatment:hour 3.68, 18.41  3.39      0.09 .002    .979
# 10                 gender:hour 1.84, 18.41  3.39      0.45 .004    .628
# 11       treatment:gender:hour 3.68, 18.41  3.39      0.62 .011    .641
# 12                  phase:hour 3.60, 35.96  2.67      1.18 .015    .335
# 13        treatment:phase:hour 7.19, 35.96  2.67      0.35 .009    .930
# 14           gender:phase:hour 3.60, 35.96  2.67      0.93 .012    .449
# 15 treatment:gender:phase:hour 7.19, 35.96  2.67      0.74 .019    .646
# ---
# Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1 

Answer 1

Those degrees of freedom are the result of applying the Greenhouse–Geisser correction, which is the default correction applied by aov_car, as far as I know.

Salkind (2010), states that:

In addition to the usual assumptions of normality of the error and homogeneity of variance, the F test for repeated-measurement designs assumes a condition called sphericity. Intuitively, this condition indicates that the ranking of the subjects does not change across experimental treatments. This is equivalent to stating that the population correlation (computed from the subjects' scores) between two treatments is the same for all pairs of treatments. This condition implies that there is no interaction between the subject factor and the treatment. If the sphericity assumption is not valid, then the F test becomes too liberal (ie, the proportion of rejections of the null hypothesis is larger than the α level when the null hypothesis is true).

As per this link :

The degree to which sphericity is present, or not, is represented by a statistic called epsilon (ε). An epsilon of 1 (ie, ε = 1) indicates that the condition of sphericity is exactly met. The further epsilon decreases below 1 (ie, ε < 1), the greater the violation of sphericity. Therefore, you can think of epsilon as a statistic that describes the degree to which sphericity has been violated. The lowest value that epsilon (ε) can take is called the lower-bound estimate. Both the Greenhouse-Geisser and the Huynd-Feldt procedures attempt to estimate epsilon (ε), albeit in different ways (it is an estimate because we are dealing with samples, not populations). For this reason, the estimates of sphericity (ε) tend to always be different depending on which procedure is used. By estimating epsilon (ε), all these procedures then use their sphericity estimate (ε) to correct the degrees of freedom for the F-distribution.

The calculation of ε is rather long, and one would need to paste here a series of formulas and notations to fully reproduce the logic, so I just refer you to the quoted sources.

You can see what ε values R is using to perform the correction through the summary command:

b<-aov_car(value ~ treatment * gender + Error(id/(phase*hour)), 
        data = obk.long, observed = "gender")
b
summary(b)

The degrees of freedom you see in your question are the result of multiplying the original degrees of freedom by ε.

If you want to see those original degrees of freedom, you can specify that you don't want any correction:

c<-aov_car(value ~ treatment * gender + Error(id/(phase*hour)), 
           data = obk.long, observed = "gender", anova_table = list(correction = "none"))
c

Sources: Salkind, NJ (Ed.). (2010). Encyclopedia of research design. SAGE Publications.