[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