Quin.func(Subject, time, conc, dose, interval, lV, lKa, lKe)
Quin.nlme <-
nlme(conc ~ Quin.func(Subject, time, conc, dose, interval,
lV, lKa, lKe),
fixed = list(lV ~ ., lKa ~ ., lKe ~ .),
random = list(lV ~ ., lKe ~ .),
cluster = ~ Subject, data = Quinidine,
start = list(fixed = c(5, -0.3, -3)),
re.structure = "d", na.pattern = ~ !is.na(conc))