#
# Script demo.glmer.txt.
#

library(lme4)
set.seed(42)
X       <- dget("X.txt")
ind.trn <- sample(1:124,100)
ind.val <- setdiff(1:124,ind.trn)
TS      <- X[ind.trn,]
VS      <- X[ind.val,]
g.trn   <- glmer(cbind(Dead,Alive) ~ (Trt+0)/Dose + (Dose | Rep),
                 data=TS,family=binomial(link="cloglog"),nAGQ=0)
coefs   <- unlist(getME(g.trn,c("theta","beta")))
newdev  <- update(g.trn, data=VS, devFunOnly=TRUE,
                  control=glmerControl(check.nobs.vs.nRE="ignore"))
ll <- -0.5*newdev(coefs)