[R] multiple nls - next fit even after convergence problem

[Gabor Grothendieck]

See:

https://stat.ethz.ch/pipermail/r-help/2005-May/072035.html

On 8/4/07, GOUACHE David <D.GOUACHE at arvalisinstitutduvegetal.fr> wrote:
> Hello R-gurus,
> I'm trying to adjust different growth curves to a rather extensive dataset.
> I wrote up a function to go through all of them, but am encountering a problem :
> among the more than 1000 curves I have, obviously for some of them I encounter conversion problems.
> I'd like for my function to keep going to the next curve and store the fact that for curve number X I had a convergence problem.
> This is my original function :
>
> comp.fit.2<-function(tab)
> {
>  fit.log<-nls(surf.vert.tot ~ 100/(1+exp(((log(81))/a)*(sum.T.levee-b))), start=list( a=ifelse(sum(tab$surf.vert.tot>76)<1 | sum(tab$surf.vert.tot<15)<1,400,-max(tab$sum.T.levee[tab$surf.vert.tot>76],na.rm=T)+min(tab$sum.T.levee[tab$surf.vert.tot<15],na.rm=T)), b=tab$sum.T.levee[abs(tab$surf.vert.tot-50)==min(abs(tab$surf.vert.tot-50),na.rm=T)]),data=tab,control=list(maxiter=100))
>  fit.exp<-nls(surf.vert.tot ~ ifelse(100-10*exp((log(9)/a)*(sum.T.levee-b))<0,0,100-10*exp((log(9)/a)*(sum.T.levee-b))), start=list( a=ifelse(sum(tab$surf.vert.tot>76)<1 | sum(tab$surf.vert.tot<15)<1,400,-max(tab$sum.T.levee[tab$surf.vert.tot>76],na.rm=T)+min(tab$sum.T.levee[tab$surf.vert.tot<15],na.rm=T)), b=tab$sum.T.levee[abs(tab$surf.vert.tot-90)==min(abs(tab$surf.vert.tot-90),na.rm=T)]),data=tab,control=list(maxiter=100))
>  fit.gomp<-nls(surf.vert.tot ~ 100*exp(-(log(10/9))*exp(-( ( log(log(10/9))-log(log(10)) )/a )*(sum.T.levee-b))), start=list( a=ifelse(sum(tab$surf.vert.tot>76)<1 | sum(tab$surf.vert.tot<15)<1,400,-max(tab$sum.T.levee[tab$surf.vert.tot>76],na.rm=T)+min(tab$sum.T.levee[tab$surf.vert.tot<15],na.rm=T)), b=tab$sum.T.levee[abs(tab$surf.vert.tot-90)==min(abs(tab$surf.vert.tot-90),na.rm=T)]),data=tab,control=list(maxiter=100))
>  rmse.log<-sqrt(mean(residuals(fit.log)^2))
>  rmse.exp<-sqrt(mean(residuals(fit.exp)^2))
>  rmse.gomp<-sqrt(mean(residuals(fit.gomp)^2))
>  data.frame(rmse.log=rmse.log,rmse.gomp=rmse.gomp,rmse.exp=rmse.exp,semis=unique(tab$semis),densite=unique(tab$densite),traitement=unique(tab$traitement),bloc=unique(tab$bloc),num.feuille.def=unique(tab$num.feuille.def))
> }
>
> I've thought of just storing the 3 model results in a list, and then going through each model object,
> but even then, when I have a convergence problem, the function breaks off.
> What I'd like is to find a way for the function to keep running despite the convergence problem, and store in the model object NA or something like that...
> Does anybody have an idea for this ?
> Below is a sample of my data set, called tab.ex, with the points for 3 curves. The 1st and 3rd converge, ans the second fails and breaks off the function.
> Thanks in advance for your suggestions !!!
>
>
> tab.ex$semis<-as.factor(tab.ex$semis)
> tab.ex$densite<-as.factor(tab.ex$densite)
> tab.ex$bloc<-as.factor(tab.ex$bloc)
> tab.ex$num.feuille.def<-as.factor(tab.ex$num.feuille.def)
> tab.ex$num.not<-as.factor(tab.ex$num.not)
> tab.ex$sum.T.levee<-tab.ex$sum.T.semis-117
>
> semis densite traitement bloc num.feuille.def surf.vert.tot sum.T.semis
> 1 1 NT 1 1 100 1764.95
> 1 1 NT 1 1 100 1867.3
> 1 1 NT 1 1 98.50833333 1983.25
> 1 1 NT 1 1 37.91416667 2200
> 1 1 NT 1 1 0 2308.45
> 1 1 NT 1 1 0 2436.05
> 1 1 NT 1 1 0 2549.95
> 1 1 NT 1 1 0 2678.05
> 1 2 NT 1 1 100 1641.05
> 1 2 NT 1 1 100 1764.95
> 1 2 NT 1 1 100 1867.3
> 1 2 NT 1 1 99.58333333 1983.25
> 1 2 NT 1 1 0 2200
> 1 2 NT 1 1 0 2308.45
> 1 2 NT 1 1 0 2436.05
> 1 2 NT 1 1 0 2549.95
> 1 2 NT 1 1 0 2678.05
> 1 1 T 1 3 100 1440.1
> 1 1 T 1 3 99.45545455 1525.7
> 1 1 T 1 3 98.11727273 1641.05
> 1 1 T 1 3 97.65545455 1764.95
> 1 1 T 1 3 95.113 1867.3
> 1 1 T 1 3 93.04545455 1983.25
> 1 1 T 1 3 87.00666667 2200
> 1 1 T 1 3 50.895 2308.45
> 1 1 T 1 3 21.67416667 2436.05
> 1 1 T 1 3 4.333333333 2549.95
> 1 1 T 1 3 0 2678.05
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>