Examp2.4.2.2: Examp2.4.2.2 from Duchateau, L. and Janssen, P. and Rowlands,...

In VetResearchLMM: Linear Mixed Models - An Introduction with Applications in Veterinary Research

Description

Examp2.4.2.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 ([e-mail])

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

#-------------------------------------------------------------
## Example 2.4.2.2 p-64
#-------------------------------------------------------------
# PROC MIXED DATA=ex125 METHOD=ML;
# CLASS drug dose region;
# MODEL pcv=drug dose drug*dose;
# RANDOM region drug*region;
# RUN;
# 
# PROC MIXED DATA=ex125 METHOD=REML;
# CLASS drug dose region;
# MODEL pcv=drug dose drug*dose;
# RANDOM region drug*region;
# RUN;
 
 
library(lme4)
str(ex125)

fm2.4 <- 
  lme4::lmer(
         formula    = Pcv ~ dose*Drug + (1|Region/Drug)
       , data       = ex125
       , REML       = FALSE
       , control    = lmerControl()
       , start      = NULL
       , verbose    = 0L
    #  , subset
    #  , weights
    #  , na.action
    #  , offset
       , contrasts  = NULL
       , devFunOnly = FALSE
    #  , ...
       ) 
summary(fm2.4)
anova(fm2.4)

fm2.5 <- 
  lme4::lmer(
         formula    = Pcv ~ dose*Drug + (1|Region/Drug)
       , data       = ex125
       , REML       = TRUE
       , control    = lmerControl()
       , start      = NULL
       , verbose    = 0L
    #  , subset
    #  , weights
    #  , na.action
    #  , offset
       , contrasts  = NULL
       , devFunOnly = FALSE
    #  , ...
       )  
summary(fm2.5)
anova(fm2.5)

library(lmerTest)

fm2.6 <- 
    lmerTest::lmer(
         formula    = Pcv ~ dose*Drug + (1|Region/Drug)
      	, data       = ex125
      	, REML       = FALSE
      	, control    = lmerControl()
      	, start      = NULL
      	, verbose    = 0L
      #	, subset
      #	, weights
      #	, na.action
      #	, offset
      	, contrasts  = NULL
      	, devFunOnly = FALSE
      #	, ...
      	)
summary(fm2.6)
anova(fm2.6)

fm2.7 <- 
    lmerTest::lmer(
         formula    = Pcv ~ dose*Drug + (1|Region/Drug)
      	, data       = ex125
      	, REML       = TRUE
      	, control    = lmerControl()
      	, start      = NULL
      	, verbose    = 0L
      #	, subset
      #	, weights
      #	, na.action
      #	, offset
      	, contrasts  = NULL
      	, devFunOnly = FALSE
      #	, ...
      	) 
summary(fm2.7)
anova(fm2.7)

VetResearchLMM documentation built on May 1, 2019, 8:39 p.m.