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: [EMAIL PROTECTED] > > [mailto:[EMAIL PROTECTED] 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. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.