> digest.mod1 <- lm(digest ~ sequence + sequence/steer + period + trt + carryover) > anova(digest.mod1) Analysis of Variance Table Response: digest Df Sum Sq Mean Sq F value Pr(>F) sequence 5 331.67 66.333 7.3723 0.0006434 *** period 2 288.17 144.083 16.0134 0.0001011 *** trt 2 559.50 279.750 31.0913 1.448e-06 *** carryover 2 18.37 9.187 1.0211 0.3801418 sequence:steer 6 114.33 19.056 2.1178 0.1015777 Residuals 18 161.96 8.998 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > digest.mod2 <- lm(digest ~ sequence + sequence/steer + period + carryover + trt ) > anova(digest.mod2) Analysis of Variance Table Response: digest Df Sum Sq Mean Sq F value Pr(>F) sequence 5 331.67 66.333 7.3723 0.0006434 *** period 2 288.17 144.083 16.0134 0.0001011 *** carryover 2 129.60 64.800 7.2019 0.0050362 ** trt 2 448.27 224.137 24.9106 6.534e-06 *** sequence:steer 6 114.33 19.056 2.1178 0.1015777 Residuals 18 161.96 8.998 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > digest.mod3 <- lm(digest ~ sequence + sequence/steer + period + trt1 + trt2 + carry1 + carry2) > summary(digest.mod3) Call: lm(formula = digest ~ sequence + sequence/steer + period + trt1 + trt2 + carry1 + carry2) Residuals: Min 1Q Median 3Q Max -5.417 -1.422 0.000 1.625 3.604 Coefficients: (60 not defined because of singularities) Estimate Std. Error t value Pr(>|t|) (Intercept) 50.8542 1.9037 26.713 6.20e-16 *** sequence2 -2.8750 2.5246 -1.139 0.269710 sequence3 -6.3542 2.5246 -2.517 0.021535 * sequence4 -1.0208 2.5246 -0.404 0.690711 sequence5 1.6667 2.4492 0.681 0.504847 sequence6 0.4583 2.5246 0.182 0.857964 period2 5.9167 1.2246 4.832 0.000134 *** period3 6.0833 1.2246 4.968 9.96e-05 *** trt1 4.5625 0.7905 5.772 1.80e-05 *** trt2 0.5000 0.7905 0.633 0.534989 carry1 0.9375 1.0605 0.884 0.388353 carry2 -1.5000 1.0605 -1.414 0.174311 sequence1:steer2 3.6667 2.4492 1.497 0.151697 sequence2:steer2 NA NA NA NA sequence3:steer2 NA NA NA NA sequence4:steer2 NA NA NA NA sequence5:steer2 NA NA NA NA sequence6:steer2 NA NA NA NA sequence1:steer3 NA NA NA NA sequence2:steer3 -4.0000 2.4492 -1.633 0.119796 sequence3:steer3 NA NA NA NA sequence4:steer3 NA NA NA NA sequence5:steer3 NA NA NA NA sequence6:steer3 NA NA NA NA sequence1:steer4 NA NA NA NA sequence2:steer4 NA NA NA NA sequence3:steer4 NA NA NA NA sequence4:steer4 NA NA NA NA sequence5:steer4 NA NA NA NA sequence6:steer4 NA NA NA NA sequence1:steer5 NA NA NA NA sequence2:steer5 NA NA NA NA sequence3:steer5 -3.3333 2.4492 -1.361 0.190310 sequence4:steer5 NA NA NA NA sequence5:steer5 NA NA NA NA sequence6:steer5 NA NA NA NA sequence1:steer6 NA NA NA NA sequence2:steer6 NA NA NA NA sequence3:steer6 NA NA NA NA sequence4:steer6 NA NA NA NA sequence5:steer6 NA NA NA NA sequence6:steer6 NA NA NA NA sequence1:steer7 NA NA NA NA sequence2:steer7 NA NA NA NA sequence3:steer7 NA NA NA NA sequence4:steer7 -3.3333 2.4492 -1.361 0.190310 sequence5:steer7 NA NA NA NA sequence6:steer7 NA NA NA NA sequence1:steer8 NA NA NA NA sequence2:steer8 NA NA NA NA sequence3:steer8 NA NA NA NA sequence4:steer8 NA NA NA NA sequence5:steer8 NA NA NA NA sequence6:steer8 NA NA NA NA sequence1:steer9 NA NA NA NA sequence2:steer9 NA NA NA NA sequence3:steer9 NA NA NA NA sequence4:steer9 NA NA NA NA sequence5:steer9 -3.6667 2.4492 -1.497 0.151697 sequence6:steer9 NA NA NA NA sequence1:steer10 NA NA NA NA sequence2:steer10 NA NA NA NA sequence3:steer10 NA NA NA NA sequence4:steer10 NA NA NA NA sequence5:steer10 NA NA NA NA sequence6:steer10 NA NA NA NA sequence1:steer11 NA NA NA NA sequence2:steer11 NA NA NA NA sequence3:steer11 NA NA NA NA sequence4:steer11 NA NA NA NA sequence5:steer11 NA NA NA NA sequence6:steer11 -3.3333 2.4492 -1.361 0.190310 sequence1:steer12 NA NA NA NA sequence2:steer12 NA NA NA NA sequence3:steer12 NA NA NA NA sequence4:steer12 NA NA NA NA sequence5:steer12 NA NA NA NA sequence6:steer12 NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 3 on 18 degrees of freedom (1 observation deleted due to missingness) Multiple R-squared: 0.8901, Adjusted R-squared: 0.7864 F-statistic: 8.578 on 17 and 18 DF, p-value: 1.803e-05 > anova(digest.mod3) Analysis of Variance Table Response: digest Df Sum Sq Mean Sq F value Pr(>F) sequence 5 331.67 66.33 7.3723 0.0006434 *** period 2 288.17 144.08 16.0134 0.0001011 *** trt1 1 541.50 541.50 60.1821 3.787e-07 *** trt2 1 18.00 18.00 2.0005 0.1743109 carry1 1 0.38 0.38 0.0417 0.8405264 carry2 1 18.00 18.00 2.0005 0.1743109 sequence:steer 6 114.33 19.06 2.1178 0.1015777 Residuals 18 161.96 9.00 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > digest.mod4 <- lm(digest ~ sequence + sequence/steer + period + carry1 + carry2 + trt1 + trt2) > summary(digest.mod4) Call: lm(formula = digest ~ sequence + sequence/steer + period + carry1 + carry2 + trt1 + trt2) Residuals: Min 1Q Median 3Q Max -5.417 -1.422 0.000 1.625 3.604 Coefficients: (60 not defined because of singularities) Estimate Std. Error t value Pr(>|t|) (Intercept) 50.8542 1.9037 26.713 6.20e-16 *** sequence2 -2.8750 2.5246 -1.139 0.269710 sequence3 -6.3542 2.5246 -2.517 0.021535 * sequence4 -1.0208 2.5246 -0.404 0.690711 sequence5 1.6667 2.4492 0.681 0.504847 sequence6 0.4583 2.5246 0.182 0.857964 period2 5.9167 1.2246 4.832 0.000134 *** period3 6.0833 1.2246 4.968 9.96e-05 *** carry1 0.9375 1.0605 0.884 0.388353 carry2 -1.5000 1.0605 -1.414 0.174311 trt1 4.5625 0.7905 5.772 1.80e-05 *** trt2 0.5000 0.7905 0.633 0.534989 sequence1:steer2 3.6667 2.4492 1.497 0.151697 sequence2:steer2 NA NA NA NA sequence3:steer2 NA NA NA NA sequence4:steer2 NA NA NA NA sequence5:steer2 NA NA NA NA sequence6:steer2 NA NA NA NA sequence1:steer3 NA NA NA NA sequence2:steer3 -4.0000 2.4492 -1.633 0.119796 sequence3:steer3 NA NA NA NA sequence4:steer3 NA NA NA NA sequence5:steer3 NA NA NA NA sequence6:steer3 NA NA NA NA sequence1:steer4 NA NA NA NA sequence2:steer4 NA NA NA NA sequence3:steer4 NA NA NA NA sequence4:steer4 NA NA NA NA sequence5:steer4 NA NA NA NA sequence6:steer4 NA NA NA NA sequence1:steer5 NA NA NA NA sequence2:steer5 NA NA NA NA sequence3:steer5 -3.3333 2.4492 -1.361 0.190310 sequence4:steer5 NA NA NA NA sequence5:steer5 NA NA NA NA sequence6:steer5 NA NA NA NA sequence1:steer6 NA NA NA NA sequence2:steer6 NA NA NA NA sequence3:steer6 NA NA NA NA sequence4:steer6 NA NA NA NA sequence5:steer6 NA NA NA NA sequence6:steer6 NA NA NA NA sequence1:steer7 NA NA NA NA sequence2:steer7 NA NA NA NA sequence3:steer7 NA NA NA NA sequence4:steer7 -3.3333 2.4492 -1.361 0.190310 sequence5:steer7 NA NA NA NA sequence6:steer7 NA NA NA NA sequence1:steer8 NA NA NA NA sequence2:steer8 NA NA NA NA sequence3:steer8 NA NA NA NA sequence4:steer8 NA NA NA NA sequence5:steer8 NA NA NA NA sequence6:steer8 NA NA NA NA sequence1:steer9 NA NA NA NA sequence2:steer9 NA NA NA NA sequence3:steer9 NA NA NA NA sequence4:steer9 NA NA NA NA sequence5:steer9 -3.6667 2.4492 -1.497 0.151697 sequence6:steer9 NA NA NA NA sequence1:steer10 NA NA NA NA sequence2:steer10 NA NA NA NA sequence3:steer10 NA NA NA NA sequence4:steer10 NA NA NA NA sequence5:steer10 NA NA NA NA sequence6:steer10 NA NA NA NA sequence1:steer11 NA NA NA NA sequence2:steer11 NA NA NA NA sequence3:steer11 NA NA NA NA sequence4:steer11 NA NA NA NA sequence5:steer11 NA NA NA NA sequence6:steer11 -3.3333 2.4492 -1.361 0.190310 sequence1:steer12 NA NA NA NA sequence2:steer12 NA NA NA NA sequence3:steer12 NA NA NA NA sequence4:steer12 NA NA NA NA sequence5:steer12 NA NA NA NA sequence6:steer12 NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 3 on 18 degrees of freedom (1 observation deleted due to missingness) Multiple R-squared: 0.8901, Adjusted R-squared: 0.7864 F-statistic: 8.578 on 17 and 18 DF, p-value: 1.803e-05 > anova(digest.mod4) Analysis of Variance Table Response: digest Df Sum Sq Mean Sq F value Pr(>F) sequence 5 331.67 66.33 7.3723 0.0006434 *** period 2 288.17 144.08 16.0134 0.0001011 *** carry1 1 97.20 97.20 10.8028 0.0040994 ** carry2 1 32.40 32.40 3.6009 0.0739062 . trt1 1 444.67 444.67 49.4210 1.467e-06 *** trt2 1 3.60 3.60 0.4001 0.5349895 sequence:steer 6 114.33 19.06 2.1178 0.1015777 Residuals 18 161.96 9.00 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > dev.off() null device 1 >