treeave <- read.table("treeave.dat",header=T)
# soil is automatically a factor() since it has characters

treeave.aov <- aov(score~soil,treeave)
summary(treeave.aov)
plot(treeave.aov)
boxplot(split(treeave$score,treeave$soil))