### nlmer() convergence testing / monitoring / ...
## -------------------
### The output of tests here are *not* 'diff'ed (<==> no *.Rout.save file)
library(lme4)
## 'Theoph' Data modeling
if (lme4:::testLevel() > 1) {
Th.start <- c(lKe=-2.5, lKa=0.5, lCl=-3)
(nm2 <- nlmer(conc ~ SSfol(Dose, Time,lKe, lKa, lCl) ~ lKe + lKa + lCl|Subject,
Theoph, start = Th.start))
(nm3 <- nlmer(conc ~ SSfol(Dose, Time,lKe, lKa, lCl) ~ (lKe|Subject) +
(lKa|Subject) + (lCl|Subject),
Theoph, start = Th.start))
## dropping lKe from random effects:
(nm4 <- nlmer(conc ~ SSfol(Dose, Time,lKe, lKa, lCl) ~ lKa + lCl|Subject,
Theoph, start = Th.start,
control = nlmerControl(tolPwrss=1e-8)))
(nm5 <- nlmer(conc ~ SSfol(Dose, Time,lKe, lKa, lCl) ~ (lKa|Subject) + (lCl|Subject),
Theoph, start = Th.start))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.