Bonjour,
Je fais des graphes avec des courbes de tendances logistiques. 21 variétés pour 21 graphes.
Voici la fonction que j'utilise pour la variété 13 par exemple:
plot(Stade[Var==15]~Jour[Var==15]
para13<-getInitial(Stade[Var==13]~SSlogis(Jour[Var==13], Asym, xmid, scal))
lines(Jour[Var==13], SSlogis(Jour[Var==13], para13[1], para13[2], para13[3]))
Tout marche bien sauf pour 2 variétés. Voici le premier message d'erreur pour la variété 16:
Erreur dans nls(y~1/(1+exp((xmid-x)/scal)), data=xy, start=list(xmid=aux[1],:
le pas 0.000488281 est devenu inférieur à 'minFactor' de 0.000976562
et celui pour la variété 18:
Erreur dans nls(y~1/(1+exp((xmid-x)/scal)), data=xy, start=list(xmid=aux[1],:
gradient singulier
Que dois-je faire pour qu'il me trace quand même les courbes???
Merci pour tout!
Antoine