#
# Demo script; single random effect.
#

library(lme4)
fit   <- glmer(cbind(incidence, size - incidence) ~ 0 + period + (1 | herd),
             family = binomial, data = cbpp)
ccc   <- getME(fit,"beta")
sigma <- getME(fit,"theta")

# Roll-your-own:
set.seed(101)
Z     <- rnorm(length(levels(cbpp$herd)),0,sigma)
lnpr  <- with(cbpp,ccc[period] +  Z[herd]) 
p     <- 1/(1+exp(-lnpr))
s.ryo <- rbinom(nrow(cbpp),cbpp$size,p)

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

# Check for equality:
print(all.equal(s.ryo,s.mer))