#use the data set chickwts from base R
head(chickwts)
N <- table(chickwts$feed)
stdev <- tapply(chickwts$weight, chickwts$feed, sd)
mn <- tapply(chickwts$weight, chickwts$feed, mean)
anovaSummarized(N, mn, stdev)
Run the code above in your browser using DataLab