#
# Script demo.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,]
f.trn   <- glmer(cbind(Dead,Alive) ~ (Trt+0)/Dose + (Dose | Rep),
                 data=TS,family=binomial(link="cloglog"),nAGQ=0)
coefs   <- unlist(getME(f.trn,c("theta","beta")))
newdev  <- update(f.trn, data=VS, devFunOnly=TRUE)
newdev(coefs)