[R-sig-ME] 2 correlated random effects with quadrature?

[Ross Boylan]

Is there a way to fit generalized linear mixed model with 2 correlated
random effects in R, using quadrature?  At the moment, I'm only
concerned with binary outcomes.

When I try glmer from lme4 with the quadrature argument I get
Error: AGQ only defined for a single scalar random-effects term

Yes, I know 2 dimensional quadrature is slow.

Ross Boylaln