pqlmer function

pqlmer

Fit a GLMM model with multivariate normal random effects using Penalized Quasi-Likelihood for mermod objects.

pqlmer(formula, family, data, niter = 40, verbose = T)

Arguments

• formula: The lme4 model formula.
• family: a family function of the error distribution and link function to be used in the model.
• data: the dataframe containing the variables in the model.
• niter: Maximum number of iterations to perform.
• verbose: if TRUE, iterations are printed to console.

Returns

A pseudo linear mixed model of class "lme" .

Examples

# Compare lmer PQL with lme PQL

library(MASS)

lmePQL = glmmPQL(y ~ trt + week + I(week > 2), random = ~ 1 | ID,
                  family = binomial, data = bacteria,
                  verbose = FALSE)

merPQL= pqlmer(y ~ trt + week + I(week > 2) + (1 | ID),
               family = binomial, data = bacteria,
               verbose = FALSE)

summary(lmePQL)
summary(merPQL)

See Also

glmmPQL