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