r - ANOVA 反復測定

Question

次の線形モデルの ANOVA の Greenhouse-Geisser-Epsilon (二次形式の分布のボックス近似に基づく) を計算します。

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 パッケージを使用して、上記の anova で使用される検定統計量のおおよその分布の GG-Epsilon を取得する方法はありますか。

次の ANOVA で使用される検定統計量の分布の GG-Epsilon を取得する方法しか知りません。

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

しかし、この検定統計量は、最初の例で anova 関数によって使用された統計量とは異なる 2 次形式で構成されています。