data(isric)
isric1 <- isric[1:32,]
library(nlme)
omuto <- nlsList(y ~ SSomuto(x,ths1,alp1,ths2,alp2) | Sample, isric1)
omuto.nlme <- nlme(omuto)## for fitting mixed-effects models
summary(omuto.nlme)
Run the code above in your browser using DataLab