# Boxplot of chick weights for different feed types
plot(weight~feed,data = chickwts)
# Compute ANOVA table where groups are based on
# feed
anova(lm(weight~feed,data=chickwts))

# Another way to get the same ANOVA table:
# First create an ANOVA object.
weight.aov<-aov(weight~feed, data=chickwts)
# Then display the summary.
summary(weight.aov)