[R-sig-ME] VarCorr vs ranef

David Afshartous dafshartous at med.miami.edu
Mon Sep 8 16:05:39 CEST 2008


For a nice discussion of this topic see p.459 of Gelman & Hill (2007).
(see http://www.stat.columbia.edu/~gelman/)


On 8/31/08 8:52 AM, "D Chaws" <cat.dev.urandom at gmail.com> wrote:

> Thanks so much for the reply.  This still seems very strange.  Even if the
> differences between population and subject effects is the issue, wouldn't
> one expect a bit more similarity between the actual effects for the subjects
> and the population effects inferred from those effects?  Dr. Bates or anyone
> else, can you resolve this mystery?  Alternatively, is there a way to get
> population estimates of the random effects for subjects (contradiction in
> terms?), like fitted.lme with the level = 0 argument?
> 
> All this is in service of an attempt to gain a simple scatterplot between
> two random effects that closely reflect the estimates from VarCorr or
> summary.  I'm sure someone must have a method for this already worked out.
> pairs.lme plots the raw data from ranef, so the discrepancy is still a
> problem there.
> 
> Thanks so much for your help.
> 
> - DC
> 
> On Sun, Aug 31, 2008 at 6:13 AM, Daniel Ezra Johnson <
> danielezrajohnson at gmail.com> wrote:
> 
>> On Sun, Aug 31, 2008 at 6:53 AM, D Chaws <cat.dev.urandom at gmail.com>
>> wrote:
>>> Can someone tell me why correlations between raw random effects are
>>> different from that provided in VarCorr for lme models?
>>> For example:
>>> 
>>> fm1 = lme(distance ~ I(age-8), random = ~ 1 + I(age-8) | Subject, data =
>>> Orthodont)
>>> R# VarCorr(fm1)
>>> Subject = pdLogChol(1 + I(age - 8))
>>>            Variance StdDev Corr
>>> (Intercept) 3.55937  1.8866 (Intr)
>>> I(age - 8)  0.05127  0.2264 0.209
>>> Residual    1.71620  1.3100
>>> 
>>> and
>>> 
>>> R# cor(ranef(fm1))
>>>            (Intercept) I(age - 8)
>>> (Intercept)      1.0000     0.5764
>>> I(age - 8)       0.5764     1.0000
>>> 
>> 
>> This isn't a complete answer, but the figures in VarCorr and the model
>> summary are the population estimates for the random effects (the
>> parameters) while everything derived from ranef() refers to the actual
>> Subjects in the data (the BLUPs).
>> 
>> Look at:
>> 
>>> sd(ranef(fm1))
>> (Intercept)  I(age - 8)
>>  1.7359554   0.1557322
>> 
>> Those figures don't match the VarCorr standard deviations either,
>> especially the second.
>> 
>> I don't know why the BLUPs pattern differently, exactly, but I did
>> look at plot(coefs(fm1)) which suggested Sex should be added as a
>> fixed effect. Once I did that, the correlation between the random
>> effects changed quite a lot (but was still different between VarCorr
>> and ranef; the population correlation was actually negative...)
>> 
>> D
>> 
> 
> [[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models




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