ANOVA重复测量

Question

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

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-package来获取GG-Epsilon，以获得上面anova使用的测试统计量的近似分布。

我只知道如何获得GG-Epsilon用于以下ANOVA中使用的检验统计量的分布：

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

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