Description Author(s) References See Also Examples
Examp3.2 is used for inspecting probability distribution and to define a plausible process through linear models and generalized linear models.
Muhammad Yaseen (myaseen208@gmail.com)
Duchateau, L. and Janssen, P. and Rowlands, G. J. (1998).Linear Mixed Models. An Introduction with applications in Veterinary Research. International Livestock Research Institute.
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 | #-------------------------------------------------------------
## 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
, control = lmerControl()
, start = NULL
, verbose = 0L
# , subset
# , weights
# , na.action
# , offset
, contrasts = list(sex = "contr.SAS", breed1 = "contr.SAS")
, devFunOnly = FALSE
# , ...
)
summary(fm3.4)
anova(object = fm3.4, ddf = "Satterthwaite")
lsmeansLT(model = fm3.4)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.