> > package <- factor(package) > > bacteria package logcount 1 1 7.66 2 1 6.98 3 1 7.80 4 2 5.26 5 2 5.44 6 2 5.80 7 3 7.41 8 3 7.33 9 3 7.04 10 4 3.51 11 4 2.91 12 4 3.66 > > tapply(logcount,package,mean) 1 2 3 4 7.48 5.50 7.26 3.36 > > tapply(logcount,package,sd) 1 2 3 4 0.4386342 0.2749545 0.1946792 0.3968627 > > ex_02_01.mod <- aov(logcount ~ package) > > summary(ex_02_01.mod) Df Sum Sq Mean Sq F value Pr(>F) package 3 32.873 10.9576 94.584 1.376e-06 *** Residuals 8 0.927 0.1158 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > predict(ex_02_01.mod) 1 2 3 4 5 6 7 8 9 10 11 12 7.48 7.48 7.48 5.50 5.50 5.50 7.26 7.26 7.26 3.36 3.36 3.36 >