自带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 value  Pr(>F)   
Type       1 2667.2  2667.2   60.41 0.00148 **
Residuals  4  176.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       1  888.6   888.6  215.46 0.000125 ***
conc:Type  1  239.2   239.2   58.01 0.001595 **
Residuals  4   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