## using a new data set to predict the probability of oncholysis at certain times

treatment2 <- rep(c("itraconazole", "terbinafine"), 3)
time <- c(0, 0, 4, 4, 12, 12)
newdata <- data.frame(treatment2,
                      time)
predicted.probability <- predict(fit.glmm.bi, newdata = newdata,
        re.form = NA,
        type = "response", allow.new.levels = FALSE,
        na.action = na.pass)

tabel <- cbind(newdata, round(predicted.probability, 4))
colnames(tabel) <- c("Treatment", "Time in weeks", "Predicted probability of oncholysis")
pandoc.table(tabel)