setwd("...")

##import data
ps <- read.table("ps.txt")
colnames(ps) <- c("Solution", "Growth")
head(ps, 16)

##ANOVA test
##Please note the factorization step
ps$Solution <- as.factor(ps$Solution)
anova(lm(Growth ~ Solution, ps))

##Following up steps
pairwise.t.test(ps$Growth, ps$Solution)
tapply(ps$Growth, ps$Solution, mean)