Description Author(s) References See Also Examples
Exam3.5 fixed location, factorial treatment structure, Gaussian response
Muhammad Yaseen (myaseen208@gmail.com)
Adeela Munawar (adeela.uaf@gmail.com)
Stroup, W. W. (2012). Generalized Linear Mixed Models: Modern Concepts, Methods and Applications. CRC Press.
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
# , ...
)
)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.