# Exam3.5: Example 3.5 from Generalized Linear Mixed Models: Modern... In StroupGLMM: R Codes and Datasets for Generalized Linear Mixed Models: Modern Concepts, Methods and Applications by Walter W. Stroup

## Description

Exam3.5 fixed location, factorial treatment structure, Gaussian response

## References

1. Stroup, W. W. (2012). Generalized Linear Mixed Models: Modern Concepts, Methods and Applications. CRC Press.

`DataSet3.2`
 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82``` ```data(DataSet3.2) DataSet3.2\$A <- factor(x = DataSet3.2\$A) DataSet3.2\$B <- factor(x = DataSet3.2\$B) DataSet3.2\$loc <- factor(x = DataSet3.2\$loc, level = c(8, 1, 2, 3, 4, 5, 6, 7)) Exam3.5.lm <- lm( formula = Y~ A + B +loc , data = DataSet3.2 # , subset # , weights # , na.action , method = "qr" , model = TRUE # , x = FALSE # , y = FALSE , qr = TRUE , singular.ok = TRUE , contrasts = NULL # , offset # , ... ) ##---a0 marginal mean list3.5.a0 <- list(B = c("0" = 1,"1" = 0) ) library(phia) Test3.5.a0 <- summary(testFactors( model = Exam3.5.lm , levels = list3.5.a0) ) ##---b0 marginal mean list3.5.b0 <- list(B = c("0" = 1,"1" = 0) ) Test3.5.b0 <- summary(testFactors( model = Exam3.5.lm , levels = list3.5.b0) ) ##---Simple effect of A on B0 Test3.5.AB0 <- summary(testInteractions( model = Exam3.5.lm , custom = list3.5.b0 , across = "B") ) ##---Simple effect of B on A0 Test3.5.BA0 <- summary(testInteractions( model = Exam3.5.lm , custom = list3.5.a0 , across = "A") ) ##---Simple Effect of A over B (SimpleEffect3.5.AB <- summary(testInteractions( model = Exam3.5.lm , fixed = "A" , across = "B") ) ) ##---Simple Effect of B over A (SimpleEffect3.5.BA <- summary(testInteractions( model = Exam3.5.lm , fixed = "B" , across = "A") ) ) #------------------------------------------------------------- ## Individula least squares treatment means #------------------------------------------------------------- (Lsm3.5 <- lsmeans::lsmeans( object = Exam3.5.lm , specs = ~A*B # , ... ) ) ```