> ex11_2.mod1 <- aov(y ~ replicate + replicate/block + a*b*c) > anova(ex11_2.mod1) Analysis of Variance Table Response: y Df Sum Sq Mean Sq F value Pr(>F) replicate 2 111.00 55.50 3.7569 0.057073 . a 1 600.00 600.00 40.6154 5.274e-05 *** b 1 253.50 253.50 17.1600 0.001637 ** c 1 54.00 54.00 3.6554 0.082274 . replicate:block 3 108.00 36.00 2.4369 0.119640 a:b 1 6.25 6.25 0.4231 0.528757 a:c 1 1.00 1.00 0.0677 0.799533 b:c 1 6.25 6.25 0.4231 0.528757 a:b:c 1 13.50 13.50 0.9138 0.359634 Residuals 11 162.50 14.77 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > summary.lm(ex11_2.mod1) Call: aov(formula = y ~ replicate + replicate/block + a * b * c) Residuals: Min 1Q Median 3Q Max -5.6250 -1.6875 0.0625 2.1562 5.6250 Coefficients: (12 not defined because of singularities) Estimate Std. Error t value Pr(>|t|) (Intercept) 32.1250 2.1486 14.952 1.18e-08 *** replicate2 6.3750 3.0386 2.098 0.05980 . replicate3 6.7500 3.0386 2.221 0.04824 * a 5.0000 0.7846 6.373 5.27e-05 *** b 3.2500 0.7846 4.142 0.00164 ** c 1.5000 0.7846 1.912 0.08227 . replicate1:block2 3.7500 3.3286 1.127 0.28389 replicate2:block3 -4.5000 3.3286 -1.352 0.20354 replicate3:block5 0.7500 3.3286 0.225 0.82586 a:b 0.6250 0.9609 0.650 0.52876 a:c -0.2500 0.9609 -0.260 0.79953 b:c -0.6250 0.9609 -0.650 0.52876 a:b:c -0.7500 0.7846 -0.956 0.35963 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 3.844 on 11 degrees of freedom Multiple R-squared: 0.8765, Adjusted R-squared: 0.7418 F-statistic: 6.507 on 12 and 11 DF, p-value: 0.002008 >