The freely available R package glmmADMB can do Adaptive
Gaussian Quadrature for this type of model,
since it is built using AD Model Builder's random effects
module which incorporates this feature.
There is now a beta version of the software for
people using R on the Mac intel platform.
http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html
    Cheers,
    Dave
-- 
David A. Fournier
P.O. Box 2040,
Sidney, B.C. V8l 3S3
Canada
Phone/FAX 250-655-3364
http://otter-rsch.com