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))