# UPDATED 7 MAR 2006 TO INCLUDE INTERACTION IN MODEL 3 Record of R session after displaying the 6 model fits from the script Rfit.R: > display (M1) lmer(formula = y ~ x + (1 | group)) coef.est coef.se (Intercept) 5.21 1.82 x -3.29 0.70 Error terms: Groups Name Std.Dev. group (Intercept) 3.25 Residual 6.27 # of obs: 120, groups: group, 15 deviance = 795.1 > display (M2) lmer(formula = y ~ x + u.full + (1 | group)) coef.est coef.se (Intercept) 1.25 2.62 x -3.21 0.70 u.full 2.65 1.28 Error terms: Groups Name Std.Dev. group (Intercept) 2.57 Residual 6.30 # of obs: 120, groups: group, 15 deviance = 791.0 > display (M3) lmer(formula = y ~ x + u.full + x:u.full + (1 + x | group)) coef.est coef.se (Intercept) -2.40 5.24 x -1.49 2.03 u.full 5.02 3.17 x:u.full -1.13 1.26 Error terms: Groups Name Std.Dev. Corr group (Intercept) 4.94 x 1.53 -0.86 Residual 6.14 # of obs: 120, groups: group, 15 deviance = 788.7 > display (M4) lmer(formula = y ~ x + (1 | group), family = binomial(link = "logit")) coef.est coef.se (Intercept) 6.13 1.51 x -3.25 0.67 Error terms: Groups Name Std.Dev. group (Intercept) 2.54 No residual sd # of obs: 120, groups: group, 15 deviance = 101.2 overdispersion parameter = 0.7 > display (M5) lmer(formula = y ~ x + (1 | group), family = quasipoisson(link = "log"), offset = log.z) coef.est coef.se (Intercept) 6.36 0.49 x -3.57 0.06 Error terms: Groups Name Std.Dev. group (Intercept) 1.71 Residual 2.34 # of obs: 120, groups: group, 15 deviance = 686.0 > display (M6) lmer(formula = y ~ 1 + (1 | state) + (1 | occupation) + (1 | state.occupation)) coef.est coef.se (Intercept) 4.54 0.94 Error terms: Groups Name Std.Dev. state.occupation (Intercept) 5.70 state (Intercept) 2.12 occupation (Intercept) 2.34 Residual 6.18 # of obs: 1200, groups: state.occupation, 400; state, 50; occupation, 8 deviance = 8326.1