valeurs initiales modèle logistique

Bonsoir,

Si quelqu'un peut m'aider à trouver les valeurs initiales de mon modèle de régression logisitique non linéaire. J'ai toujours ce type d'erreur :

Code:

model_sapin_arbre<-saturation_lumen_pct~a+(b-a)/(1+exp(c*(d-ry_Relatif_sapin)))
start<-list(a=0.003,b=0.005,c=0.007,d=0.009)
names(start)=c("a","b","c","d")
model_non_lin_sapin<-nls(model_sapin_arbre,algorithm="default",data=tab_cernes_sapin_tot,start=start,control=nls.control(maxiter=100,warnOnly=T))
Warning message:
In nls(model_sapin_arbre, algorithm = "default", data = tab_cernes_sapin_tot,  :
 gradient singulier

Bonjour,

Regarde le package nls2. Il y a quelques fonctions pour t'aider à trouver des valeurs de départ et obtenir une convergence.

Nik