#
# Demo script; two random effects.
#

library(lme4)
Dat   <- dget("data.txt")
fit   <- glmer(cbind(Good, Bad) ~ (trtmnt+0)/x + (x | batch),
               family = binomial, data = Dat)
ccc   <- fixef(fit)
Sigma <- VarCorr(fit)[[1]]

# Roll-your-own:
set.seed(101)
nrep  <- length(levels(Dat$batch))
Z     <- MASS::mvrnorm(nrep,c(0,0),Sigma)
beta0 <- ccc[1:4][as.numeric(Dat$trtmnt)]
beta1 <- ccc[5:8][as.numeric(Dat$trtmnt)]
lnpr  <- with(Dat,beta0 + beta1*x + Z[batch,1] + Z[batch,2]*x)
p     <- 1/(1+exp(-lnpr)) # Inverse of logit link.
size  <- with(Dat,Good+Bad)
s.ryo <- rbinom(nrow(Dat),size,p)

# Using simulate.merMod:
set.seed(101)
s.mer <- simulate(fit)
s.mer <- s.mer[,1][,1]

# The results are not equal.