[R-sig-ME] How to assess significance of variance components(Please discard previous e mail... but read this one)

David Duffy David.Duffy at qimr.edu.au
Thu Nov 15 02:39:14 CET 2012


On Thu, 15 Nov 2012, Stephane Chantepie wrote:

> Hi all,
>
> Even if I agree with you that testing everything is not always necessary,in
> some cases, it could be important.
>
> For my case :
>
> I have done random regression animal models to test an increase in genetic
> variance across age.

Did someone suggest the RLRsim package?  The example for exactRLRT runs:

   mA <- lmer(Reaction ~ I(Days-4.5) + (1|Subject) +
              (0 + I(Days-4.5)|Subject), sleepstudy)
   m0 <- update(mA, . ~ . - (0 + I(Days-4.5)|Subject))
   m.slope  <- update(mA, . ~ . - (1|Subject))
   #test for subject specific slopes:
   exactRLRT(m.slope, mA, m0)



More information about the R-sig-mixed-models mailing list