data(CO2)

CO2.aov <- aov(uptake ~ Type + Treatment,data=CO2)
print(summary(CO2.aov))
##             Df Sum Sq Mean Sq F value   Pr(>F)    
## Type         1   3366    3366   50.92 3.68e-10 ***
## Treatment    1    988     988   14.95 0.000222 ***
## Residuals   81   5353      66                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1