[R] Using lme (nlme) to find the conditional variance of therandom effects

Rick Bilonick rab at nauticom.net
Tue Nov 13 07:03:32 CET 2007


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.



More information about the R-help mailing list