tests/test-glmernbref.R

In lme4/lme4: Linear Mixed-Effects Models using 'Eigen' and S4

## DON'T load lme4; test is to see if glmer.nb works when
## lme4 is not loaded

set.seed(101)
dd <- data.frame(x=runif(200), f= rep(1:20, each=10))
b <- rnorm(20)
dd <- transform(dd, y = rnbinom(200, mu  = exp(1 + 2*x + b[f]), size = 2))
g <- lme4::glmer.nb(y~x + (1|f), data = dd)
stopifnot(inherits(g, "glmerMod"))