r - 方差分析重复测量

Question

我计算了以下线性模型的 ANOVA 的 Greenhouse-Geisser-Epsilon（基于二次形式分布的 Box-Approximation）。

Y<-data.frame(response=c(10,19,27,28,9,13,25,29,4,10,20,18,5,6,12,17),
               treatment=factor(rep(1:4,4)),
               subject=factor(rep(1:4,each=4)))
fit<-lm(response~treatment,Y)
anova(fit)

有没有办法使用像 car 或 ez 这样的 R 包来获得 GG-Epsilon 以获得上述 anova 使用的测试统计量的近似分布。

我只知道如何获得 GG-Epsilon 用于以下 ANOVA 中的测试统计分布：

ezANOVA(pigeon, dv=.(response), wid=.(subject), within=.(treatment), detailed=T)

但是这个检验统计量由不同的二次形式组成，而不是在第一个例子中由 anova 函数使用的统计量。