-
df=ddply(pH,c("Treatment","Week"),summarise, mean=mean(pH),sd=sd(pH),n=length(pH),se=sd/sqrt(n)) colnames(df)[3]="pH" df install.packages("plyr") library(plyr) p