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