Sorry if this is obvious, but my S skills aren't great and I haven't been able to find it documented anywhere. I want to write a new function for use with lme objects; the function will simply calculate an ICC (aka "rho") for each level of a mixed-effects model. What I need for this is pretty simple: (c(var1..varn, residual)) / sum(c(var1..varn, residual)) where var1..varn are the variances of the intercepts of each level's equation, and residual is the residual variance. The problem is getting access to var1..varn (or, as R generally reports it, sd1..sdn) programmatically. I can get the residual standard deviation with: print(model$sigma) and I can view the standard deviations with: summary(model$modelStruct$lmeStruct) But I can't figure out how to get at the standard deviations without the human-readable print format. Thanks- Andy Perrin ---------------------------------------------------------------------- Andrew J Perrin - http://www.unc.edu/~aperrin Assistant Professor of Sociology, U of North Carolina, Chapel Hill clists at perrin.socsci.unc.edu * andrew_perrin (at) unc.edu
Extracting elements from an reStruct
1 message · Andrew Perrin