A flexible version of MASS::glmmPQL
myglmmPQL(formula.glm, formula, offset=NULL, family, data, fixef.init = NULL, 
                 weights=NULL, REML=T, niter=10, verbose=T)An object wiht the same format as lmer.
formula used to fit glm for initial fixed effects
formula used to fit iterative lmer in PQL algorithm
glm offset term
glm family
glm data
initial fixed effects estimates, set to zeros if NULL
glm weights
lmer logical scalar - Should the estimates be chosen to optimize the REML criterion (as opposed to the log-likelihood)?
glmmPQL maximum number of iterations.
glmmPQL logical: print out record of iterations?
Use lme4::lmer instead of nlme::varFixed in PQL iteration to allow REML