Message-ID: <1194933812.3256.20.camel@localhost.localdomain>
Date: 2007-11-13T06:03:32Z
From: Rick Bilonick
Subject: Using lme (nlme) to find the conditional variance of therandom effects
In-Reply-To: <2323A6D37908A847A7C32F1E3662C80E012DE67F@dc1ex01.air.org>
On Mon, 2007-11-12 at 16:45 -0500, Doran, Harold wrote:
> No, don't reach into the bVar slot. Use the proper extractor function
> ranef() with postVar=T. There is no similar function for lme()
>
> > -----Original Message-----
> > From: r-help-bounces at r-project.org
> > [mailto:r-help-bounces at r-project.org] On Behalf Of Rick Bilonick
> > Sent: Monday, November 12, 2007 4:40 PM
> > To: R Help
> > Subject: [R] Using lme (nlme) to find the conditional
> > variance of therandom effects
> >
> > Using lmer in the lme4 package, you can compute the
> > conditional variance-covariance matrix of the random effects
> > using the bVar slot:
> >
> > bVar: A list of the diagonal inner blocks (upper triangles
> > only) of the positive-definite matrices on the diagonal of
> > the inverse of ZtZ+Omega.
> > With the appropriate scale factor (and conversion to a
> > symmetric matrix) these are the conditional
> > variance-covariance matrices of the random effects.
> >
> > Is there anything similar in the nlme package using the lme function?
> >
> > Rick B.
> >
> > ______________________________________________
Is there some way to get ranef with postVar=TRUE to show what the
variances are, or what the lower and upper bounds are? qqmath makes nice
plots but I need to obtain the numerical values.
Rick B.