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