[R] plotting muliple CI ellipses for EB estimates

I have empirical Bayes estimates for slopes and intercepts for a number of subjects and I would like to plot the slopes and intercepts with confidence ellipses. These ellipses would be based on the confidence intervals for the slope and intercepts (forming the major and minor axis of each ellips

[R] extracting posterior variances of random effects (lme vs. lmer)

List, I am having trouble locating documentation on how one extracts posterior variances of the random effects from an lme object. I have found equivalent documentation for lmer(using ranef()), but nothing for lme. Any help would be greatly appreciated! Sam -- Note the new contac

[R] glmmML vs. lmer - fitting overdispersed Poisson outcome.

Group, I have count data with one observation per subject. I would like to fit a glmm to these data in order to account for overdispersion in the outcome. The lmer() function does not appear to be able to handle data that have only one observation per-cluster id, even though separate variance