##Used with MesocosmBGIS1.csv for primary experiment
data<-read.csv("MesocosmBGIS1.csv",header=T,colClasses = c("factor","factor","numeric"))
attach(data)


ANOVA=aov(CFU ~ KNO3+Day, data=data)
summary(ANOVA)
# Tukey test to study each family of treatment :
tukey.out <- TukeyHSD(ANOVA,conf.level=0.95)
tukey.out



##Used with MO157O26bijoumulti2.csv for preliminary experiment
data<-read.csv("O157O26bijoumulti2.csv",header=T,colClasses = c("factor","factor","factor","factor","numeric"))
attach(data)


ANOVA=aov(CFU ~ KNO3+Day+Strain+Water_type, data=data)
summary(ANOVA)
# Tukey test to study each family of treatment :
tukey.out <- TukeyHSD(ANOVA,conf.level=0.95)
tukey.out