#
# Script try.newdata.txt.
#

library(lme4)
fit   <- glmer(cbind(incidence, size - incidence) ~ 0 + period + (1 | herd),
               family = binomial, data = cbpp)
ccc   <- getME(fit,"beta")
sigma <- getME(fit,"theta")
X     <- data.frame(herd=factor(rep(1:30,each=4)),size=rep(40,120),
                    period=factor(rep(1:4,30)))

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

# Using simulate.merMod, first try.
s.mer1 <- simulate(fit,newdata=X,allow.new.levels=TRUE)

# Using simulate.merMod, second try.
newpar <- list(theta=getME(fit,"theta"),beta=getME(fit,"beta"))
set.seed(101)
s.mer2 <- simulate(~0+period +(1|herd),newdata=X,newparams=newpar)