search for: contraintuit

...e fitting a model to a (part) of a dataset, just for illustration. I can extract the random interaction and the residual variance in group meth==1 using VarCorr, but how do I get the other residual variance? Is there any way to get the other variances in numerical form directly - it seems a litte contraintuitive to use "as.numeric" when extracting estimates, it's a bit like good old days writing SAS-programs that reads the SAS output files... Bendix ------------------------------------------------------------------------ ------- library( nlme ) dfr <- structure(list(meth = structure...