[R] Extracting Variance Components

John Wilkinson (pipex) wilks at dial.pipex.com
Thu Oct 27 18:01:25 CEST 2005


Mike,

use ---

VarCorr(lme.object)

or for a user friendly output use varcomp from the 'ape' package--

require(ape)
varcomp(lme.object)

varcomp also allows scaling of components to unity (*100 gives %)
and also allows for cumulative sum of components.

Note. varcomp doesn't work for lmer objects.

HTH,

John
--------------------------------------------------

Michel Friesenhahn wrote---------
 
Dear List,

Is there a way to extract variance components from lmeObjects or 
summary.lme objects without using intervals()?  For my purposes I don't 
need the confidence intervals which I'm obtaining using parametric 
bootstrap.

Thanks,

Mike




More information about the R-help mailing list