data(isric)
isric1 <- isric[1:32,]
library(nlme)
library(minpack.lm)
gardner <- nlsLM(y ~ SSgard(x,thr,ths,alp,nscal), isric1)
gardner
gard <- nlsList(y ~ SSgard(x,thr,ths,alp,scal) | Sample, isric1)
coef(gard)
Run the code above in your browser using DataLab