Examp3.2 from Duchateau, L. and Janssen, P. and Rowlands, G. J. (1998).Linear Mixed Models. An Introduction with applications in Veterinary Research. International Livestock Research Institute.

Share:

Description

Examp3.2 is used for inspecting probability distribution and to define a plausible process through linear models and generalized linear models.

Author(s)

  1. Muhammad Yaseen (myaseen208@gmail.com)

References

  1. Duchateau, L. and Janssen, P. and Rowlands, G. J. (1998).Linear Mixed Models. An Introduction with applications in Veterinary Research. International Livestock Research Institute.

See Also

ex124

Examples

 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
#-------------------------------------------------------------
## Example 3.3 p-88
#-------------------------------------------------------------
# PROC MIXED DATA=ex32;
# CLASS sex sire_id breed;
# MODEL ww = sex agew breed/SOLUTION DDFM=SATTERTH;
# RANDOM sire_id(breed)/SOLUTION;
# LSMEANS breed/ADJUST = TUKEY;
# RUN;
 
 library(lmerTest)
 str(ex32)
 ex32$sire_id1 <- factor(ex32$sire_id)
 ex32$breed1   <- factor(ex32$breed)
 
 fm3.4 <- 
          lmerTest::lmer(
                          formula   = Ww ~ sex + agew + breed1 + (1|sire_id1:breed1)
                        , data      = ex32
                        , REML      = TRUE
                        , contrasts = list(sex = "contr.SAS", breed1 = "contr.SAS")
                     )
 summary(fm3.4)
 lmerTest::anova(object = fm3.4, ddf = "Satterthwaite")
 lmerTest::lsmeans(model = fm3.4, test.effs = "breed1")

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.