DNase.lis <- nlsList(density ~ fpl(conc, A, B, ld50, scal), data = DNase, cluster = ~ Run) DNase.nlme <- nlme(DNase.lis)