three_level: A dataset simulated from a three-level model
In glmmsr: Fit a Generalized Linear Mixed Model

A dataset simulated from a three-level model

1	three_level

An object of class list of length 4.

# Fit a three-level model with the Laplace approximation to the likelihood
(mod_Laplace <- glmm(response ~ covariate + (1 | cluster) + (1 | group),
                     data = three_level, family = binomial,
                     method = "Laplace"))

# if we try to fit with adaptive Gaussian quadrature, we get an error
## Not run: 
  (mod_AGQ <- glmm(response ~ covariate + (1 | cluster) + (1 | group),
                   data = three_level, family = binomial, method = "AGQ",
                   control = list(nAGQ = 15)))

## End(Not run)

# We can fit with the Sequential Reduction approximation
## Not run: 
  (mod_SR <- glmm(response ~ covariate + (1 | cluster) + (1 | group),
                  data = three_level, family = binomial, method = "SR",
                  control = list(nSL = 3)))

## End(Not run)
# the estimates of the random effects standard deviations
# are larger than those using the Laplace approximation