> ex_12_02.mod1 <- aov(y ~ a+s+m+c+t+a:s+a:m+a:c+a:t+s:m+s:c+s:t+m:c+m:t+c:t) > summary.lm(ex_12_02.mod1) Call: aov(formula = y ~ a + s + m + c + t + a:s + a:m + a:c + a:t + s:m + s:c + s:t + m:c + m:t + c:t) Residuals: ALL 16 residuals are 0: no residual degrees of freedom! Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 45.31875 NA NA NA a -1.96875 NA NA NA s -0.38125 NA NA NA m 4.85625 NA NA NA c -0.11875 NA NA NA t 1.80625 NA NA NA a:s -1.84375 NA NA NA a:m -0.20625 NA NA NA a:c 0.06875 NA NA NA a:t 0.36875 NA NA NA s:m -0.04375 NA NA NA s:c -0.36875 NA NA NA s:t 0.08125 NA NA NA m:c 0.36875 NA NA NA m:t 0.54375 NA NA NA c:t 0.21875 NA NA NA Residual standard error: NaN on 0 degrees of freedom Multiple R-squared: 1, Adjusted R-squared: NaN F-statistic: NaN on 15 and 0 DF, p-value: NA > > effects <- ex_12_02.mod1$coef[2:16] > effects a s m c t a:s a:m a:c -1.96875 -0.38125 4.85625 -0.11875 1.80625 -1.84375 -0.20625 0.06875 a:t s:m s:c s:t m:c m:t c:t 0.36875 -0.04375 -0.36875 0.08125 0.36875 0.54375 0.21875 >