自带CO2数据集。
> wb<-subset(CO2,Treatment=="chilled")
> fit<-aov(uptake~conc*Type+Error(Plant/conc),wb)
> summary(fit)
Error: Plant
Df Sum Sq Mean Sq F valuePr(>F)
Type 1 2667.22667.2 60.41 0.00148 **
Residuals4176.6 44.1
---
Signif. codes:0 ?**?0.001 ?*?0.01 ??0.05 ??0.1 ??1

Error: Plant:conc
Df Sum Sq Mean Sq F value Pr(>F)
conc 1888.6 888.6215.46 0.000125 ***
conc:Type1239.2 239.2 58.01 0.001595 **
Residuals4 16.5 4.1
---
Signif. codes:0 ?**?0.001 ?*?0.01 ??0.05 ??0.1 ??1

Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 30 869 28.97
> Anova(fit)
Error in vcov.default(mod) :
there is no vcov() method for models of class aovlist, listof