glmmPQL Question

John Michel <jwmichel at ...> writes:

 I am working on some analyses in which I have a two level random
effects model and a dichotomous dependent variable. I am using
glmmPQL with family=binomial in my model command. I have calculated
my models successfully; however, I am now trying to calculate my
variance components. I have done a lot using the lme function and
have used the following formula to estimate variance components:

VarCorr(Model.1)
GRPpdSymm(1)

What is the syntax to get the variance components in glmmPQL? If
this is not simple, I am assuming I could use lmer; however, I also
cannot figure out how to get the variance components in lmer. In any
even, I would appreciate any suggestions here.

Did you try VarCorr() ??

 From ?glmmPQL:

x <- glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,
  family = binomial, data = bacteria)

 VarCorr(x)
ID = pdLogChol(1) 
            Variance  StdDev   
(Intercept) 1.9898963 1.4106368
Residual    0.6084797 0.7800511

  You probably already know that the object returned by VarCorr()
has a slightly weird format.
 
 I don't know where your second line, GRPpdSymm(1) comes from??

  VarCorr() should work just fine on lmer objects too.

  Perhaps you can give a reproducible example of what you've tried/
what fails for you?
 
  Follow-ups should probably go to r-sig-mixed-models <at> r-project.org