> fp.mod1 <- aov(tnl ~ fee_trt*expblock + fee_trt:dcc + week:expblock) > summary(fp.mod1) Df Sum Sq Mean Sq F value Pr(>F) fee_trt 1 0.21985 0.219854 130.3805 < 2.2e-16 *** expblock 4 0.15384 0.038460 22.8078 4.769e-15 *** fee_trt:expblock 4 0.13995 0.034988 20.7487 6.727e-14 *** fee_trt:dcc 8 0.37906 0.047383 28.0996 < 2.2e-16 *** expblock:week 15 0.07145 0.004763 2.8247 0.0006002 *** Residuals 167 0.28160 0.001686 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > > library(lme4) Loading required package: Matrix Loading required package: lattice Attaching package: 'Matrix' The following object(s) are masked from 'package:base': det Attaching package: 'lme4' The following object(s) are masked from 'package:stats': AIC Warning message: package 'lme4' was built under R version 2.12.1 > > fp.mod2 <- lmer(tnl ~ fee_trt*expblock + (1|fee_trt:dcc) + (1|expblock:week)) > summary(fp.mod2) Linear mixed model fit by REML Formula: tnl ~ fee_trt * expblock + (1 | fee_trt:dcc) + (1 | expblock:week) AIC BIC logLik deviance REMLdev -576 -533.1 301 -666.6 -602 Random effects: Groups Name Variance Std.Dev. expblock:week (Intercept) 0.00030769 0.017541 fee_trt:dcc (Intercept) 0.00228484 0.047800 Residual 0.00168625 0.041064 Number of obs: 200, groups: expblock:week, 20; fee_trt:dcc, 10 Fixed effects: Estimate Std. Error t value (Intercept) 0.21665 0.02298 9.429 fee_trtcontrol 0.02316 0.03358 0.690 expblockin1 0.04730 0.01716 2.757 expblockin2 0.10970 0.01716 6.394 expblockin3 0.12727 0.01716 7.418 expblockpost 0.11849 0.01716 6.906 fee_trtcontrol:expblockin1 -0.05822 0.01874 -3.106 fee_trtcontrol:expblockin2 -0.12172 0.01874 -6.494 fee_trtcontrol:expblockin3 -0.13371 0.01874 -7.134 fee_trtcontrol:expblockpost -0.14055 0.01874 -7.499 Correlation of Fixed Effects: (Intr) f_trtc expbl1 expbl2 expbl3 expblc f_tr:1 f_tr:2 f_tr:3 fe_trtcntrl -0.585 expblockin1 -0.373 0.122 expblockin2 -0.373 0.122 0.500 expblockin3 -0.373 0.122 0.500 0.500 expblockpst -0.373 0.122 0.500 0.500 0.500 f_trtcntr:1 0.163 -0.279 -0.437 -0.218 -0.218 -0.218 f_trtcntr:2 0.163 -0.279 -0.218 -0.437 -0.218 -0.218 0.500 f_trtcntr:3 0.163 -0.279 -0.218 -0.218 -0.437 -0.218 0.500 0.500 f_trtcntrl: 0.163 -0.279 -0.218 -0.218 -0.218 -0.437 0.500 0.500 0.500 > anova(fp.mod2) Analysis of Variance Table Df Sum Sq Mean Sq F value fee_trt 1 0.007796 0.007796 4.6235 expblock 4 0.054462 0.013616 8.0745 fee_trt:expblock 4 0.139950 0.034988 20.7487 > vcov(fp.mod2) 10 x 10 Matrix of class "dpoMatrix" [,1] [,2] [,3] [,4] [,5] [1,] 0.0005279367 -0.0004510150 -0.0001471822 -0.0001471822 -0.0001471822 [2,] -0.0004510150 0.0011275374 0.0000702605 0.0000702605 0.0000702605 [3,] -0.0001471822 0.0000702605 0.0002943644 0.0001471822 0.0001471822 [4,] -0.0001471822 0.0000702605 0.0001471822 0.0002943644 0.0001471822 [5,] -0.0001471822 0.0000702605 0.0001471822 0.0001471822 0.0002943644 [6,] -0.0001471822 0.0000702605 0.0001471822 0.0001471822 0.0001471822 [7,] 0.0000702605 -0.0001756512 -0.0001405210 -0.0000702605 -0.0000702605 [8,] 0.0000702605 -0.0001756512 -0.0000702605 -0.0001405210 -0.0000702605 [9,] 0.0000702605 -0.0001756512 -0.0000702605 -0.0000702605 -0.0001405210 [10,] 0.0000702605 -0.0001756512 -0.0000702605 -0.0000702605 -0.0000702605 [,6] [,7] [,8] [,9] [,10] [1,] -0.0001471822 0.0000702605 0.0000702605 0.0000702605 0.0000702605 [2,] 0.0000702605 -0.0001756512 -0.0001756512 -0.0001756512 -0.0001756512 [3,] 0.0001471822 -0.0001405210 -0.0000702605 -0.0000702605 -0.0000702605 [4,] 0.0001471822 -0.0000702605 -0.0001405210 -0.0000702605 -0.0000702605 [5,] 0.0001471822 -0.0000702605 -0.0000702605 -0.0001405210 -0.0000702605 [6,] 0.0002943644 -0.0000702605 -0.0000702605 -0.0000702605 -0.0001405210 [7,] -0.0000702605 0.0003513025 0.0001756512 0.0001756512 0.0001756512 [8,] -0.0000702605 0.0001756512 0.0003513025 0.0001756512 0.0001756512 [9,] -0.0000702605 0.0001756512 0.0001756512 0.0003513025 0.0001756512 [10,] -0.0001405210 0.0001756512 0.0001756512 0.0001756512 0.0003513025 >