Exam5.2: Example 5.2 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

Example 5.2 from Generalized Linear Mixed Models: Modern Concepts, Methods and Applications by Walter W. Stroup(p-164)

Description

Exam5.2 three factor main effects only design

Author(s)

1. Muhammad Yaseen ([e-mail])

2. Adeela Munawar ([e-mail])

References

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

See Also

DataSet5.2

Examples

DataSet5.2$a <- factor( x = DataSet5.2$a)
DataSet5.2$b <- factor( x = DataSet5.2$b)
DataSet5.2$c <- factor(x  = DataSet5.2$c)

##---first adding factor a in model
Exam5.2.lm1 <- lm(formula = y~ a, data  = DataSet5.2)
summary(Exam5.2.lm1)
library(parameters)
model_parameters(Exam5.2.lm1)

library(emmeans)
##---A first
emmeans(object  = Exam5.2.lm1, specs = ~a)
contrast(emmeans(object  = Exam5.2.lm1, specs = ~a), method = "pairwise")
anova(object = Exam5.2.lm1)

##---then adding factor b in model
Exam5.2.lm2 <- lm(formula = y~ a + b, data  = DataSet5.2)
summary(Exam5.2.lm2)
model_parameters(Exam5.2.lm2)

emmeans(object  = Exam5.2.lm2, specs = ~b)
contrast(emmeans(object  = Exam5.2.lm2, specs = ~b), method = "pairwise")
anova(object = Exam5.2.lm2)

##---then adding factor c in model
Exam5.2.lm3 <- lm(formula = y~ a + b + c, data = DataSet5.2)

summary(Exam5.2.lm3)
model_parameters(Exam5.2.lm3)

emmeans(object  = Exam5.2.lm3, specs = ~c)
contrast(emmeans(object  = Exam5.2.lm3, specs = ~c), method = "pairwise")
anova(object = Exam5.2.lm3)

##---Now Change the order and add b first in model
Exam5.2.lm4 <- lm(formula = y ~  b, data = DataSet5.2)
summary(Exam5.2.lm4)
model_parameters(Exam5.2.lm4)

emmeans(object  = Exam5.2.lm4, specs = ~b)
contrast(emmeans(object  = Exam5.2.lm4, specs = ~b), method = "pairwise")
anova(object = Exam5.2.lm4)


##---then adding factor a in model
Exam5.2.lm5 <- lm(formula = y ~ b + a, data  = DataSet5.2)
summary(Exam5.2.lm5)
model_parameters(Exam5.2.lm5)

emmeans(object  = Exam5.2.lm5, specs = ~a)
contrast(emmeans(object  = Exam5.2.lm5, specs = ~a), method = "pairwise")
anova(object = Exam5.2.lm5)

StroupGLMM documentation built on Oct. 2, 2024, 1:07 a.m.