setwd("C:/Users/cbt12/Documents/Sean switching expt") CFU <- read.csv("2019-03-22 wspA and rpfR survival at pH 4.csv", header=TRUE) summary(CFU) #change strain and pH to factors CFU$Strain <- as.factor(CFU$Strain) CFU$pH <- as.factor(CFU$pH) summary(CFU) #ANOVA at pH 4 ph4 <- subset(CFU, pH==4) ph4.aov <- aov(CFU.mL~Strain, ph4) summary(ph4.aov) TukeyHSD(ph4.aov) #ANOVA at pH 7 ph7 <- subset(CFU, pH==7) ph7.aov <- aov(CFU.mL~Strain, ph7) summary(ph7.aov)