# STAT_515_Lec_17

bacteria <- c(7.66,6.98,7.80,5.26,5.44,5.80,
              7.41,7.33,7.04,
              3.51,2.91,3.66)

packaging <- as.factor(c(rep("Commercial",3),
                         rep("Vacuum",3),
                         rep("Mixed Gas",3),
                         rep("C02",3)))

boxplot(bacteria ~ packaging)
stripchart(bacteria~ packaging, vertical=T, pch = 19)

lm_out <- lm(bacteria ~ packaging) # linear model
anova(lm_out)
plot(lm_out,which=2)
plot(lm_out,which=1, add.smooth= F)