#
# Script scr.riz
#

library(GLMMadaptive)
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   <- mixed_model(cbind(Dead,Alive) ~ (Trt+0)/Dose,
                       random= ~ Dose | Rep,
                       data=TS,family=binomial(link="cloglog"))
beta    <- fixef(f.trn)
D       <- f.trn$D
f.val   <- mixed_model(cbind(Dead,Alive) ~ (Trt+0)/Dose,
                       random= ~ Dose | Rep,
                       data=VS,family=binomial(link="cloglog"),
                       iter_EM=0,iter_qN_outer=0,
                       initial_values=list(beta=beta,D=D))
print(logLik(f.val))