#
# Script demo1.R
#
# Using library lme4.

X <- dget("X.txt")
library(lme4)

# 1. cloglog link:
fit <- glmer(cbind(Dead,Alive) ~ (0+Trt)/Dose + (Dose | Rep),
             family=binomial(link="cloglog"),data=X)
predResp  <- predict(fit,type="response")
predLink  <- predict(fit,type="link")
g   <- binomial(link="cloglog")$linkfun
ginv <- binomial(link="cloglog")$linkinv
cat("cloglog link:\n")
cat("Range of linkfun(predict(fit,type=\"response\")) minus\n")
cat("         predict(fit,type=\"link\"):\n")
print(range(g(predResp) - predLink))
cat("\n")
cat("Range of linkinv(predict(fit,type=\"link\")) minus\n")
cat("         predict(fit,type=\"response\"):\n")
print(range(ginv(predLink) - predResp))
cat("\n")

# 2. logit link:
fit <- glmer(cbind(Dead,Alive) ~ (0+Trt)/Dose + (Dose | Rep),
             family=binomial(link="logit"),data=X)
predResp  <- predict(fit,type="response")
predLink  <- predict(fit,type="link")
g   <- binomial(link="logit")$linkfun
ginv <- binomial(link="logit")$linkinv
cat("logit link:\n")
cat("Range of linkfun(predict(fit,type=\"response\")) minus\n")
cat("         predict(fit,type=\"link\"):\n")
print(range(g(predResp) - predLink))
cat("\n")
cat("Range of linkinv(predict(fit,type=\"link\")) minus\n")
cat("         predict(fit,type=\"response\"):\n")
print(range(ginv(predLink) - predResp))
cat("\n")