# potato.R.txt
spud <- read.table("potato2.data")
rot <- spud$Rot ; bact <- spud$Bact ; temp <- spud$Temp
# Means and standard deviations
Nz <- aggregate(rot,by=list(Temp=temp,Bact=bact),length)
Meanz <- aggregate(rot,by=list(temp,bact),mean)
Varz <- aggregate(rot,by=list(temp,bact),var)
SummaryStats <- cbind(Nz,Meanz[,3],sqrt(Varz[,3]))
dimnames(SummaryStats)[[2]][3] <- "N"
dimnames(SummaryStats)[[2]][4] <- "Mean"
dimnames(SummaryStats)[[2]][5] <- "St Dev"
print(SummaryStats)
print(anova(lm(Rot~factor(Bact)*factor(Temp),data=spud)))