Modèle non linéaire nls.

Bonjour,

Quelqu'un parmi vous peux m'expliquer comment se fait les chois des paramètres (a,b,c) de l'option start d'un modèle non linéaire nls :
NL65<-nls(DV_max~a*(Ratio^c)*exp(-b*Ratio),start=c(a=1,b=1,c=1),data=AST_Zwi5) ,!c'est la première fois que j'ai utilisé le modèle non linéaire et parfois quand j'ai mal fais le choix des paramètres il m'affiche un message d'erreur suivant " ERREUR:
le nombre d'itérations a dépassé le maximum de 50 " !

Merci beaucoup d'avance.

A priori tu dois avoir une idée des valeurs de départ. Tu peux les estimer par exemple via un modèle linéaire ou glm et ensuite intégrer les paramètres ainsi obtenus en paramètre de départ. Par défault le nombre d'itération est de 50, mais tu peux augmenter ce nombre au besoin via l'argument control.

Code:

x <- rnorm(100)
y <- 0.25*exp(1.2*x)

g1 <- glm(y~x, family=gaussian(link="log"))
nl1 <- nls(y~a*exp(b*x), start=c(a=0.2,b=0.8), control=nls.control(maxiter=500))

Regarde également la fonction nls2 (package nls2) qui propose des algo pour définir les valeurs de départ à partir d'une grille de valeurs. ça permet d'éviter en partie le problème des minima locaux.

Merci beaucoup pour les réponse .