data(chickwts)
chickwts[c(1:3,69:71),]


a=aov(weight~feed,data=chickwts)       #  analysis of variance solution
summary(a)
 
l=lm(weight~feed,data=chickwts)          # linear regression solution
summary(l)

anova(l)