Fonction _ nls _ courbe

Bonjour,

Je suis confronté à un blocage, je veux faire une courbe d’exponentielle de décroissance. J’utilise le code suivant :

Code:

---

1 2 3 4 5 6 7

tab<-read.table("C:\\Users\\Momadou\\Desktop\\TEST R\\regression.csv",header=TRUE,sep=";",dec=",") tab plot (tab$ROS~tab$FMC) test <- seq(0,200,1) c <- -0.05; b<- 0; a <- 5 lines(b+a*exp(c*test),col=4) nlsreg <- nls(ROS ~ b+a*exp(c*test),data=Table,start=list(b=0,a=5,c=-0.05),trace=T)

---

Mais comme résultat j’ai ceci en dessous avec le message erreur à la fin:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

550.1471 : 0.00 5.00 -0.05 266.0114 : 1.13846736 0.10412358 -0.05087402 265.9792 : 1.13828708 0.10449806 -0.09271866 265.9161 : 1.1378039 0.1476209 -0.1870380 265.6392 : 1.1456370 0.4906249 -1.1208783 265.5984 : 1.1359695 0.6315808 -1.2541616 265.4990 : 1.1370413 0.6784515 -2.0325196 265.2979 : 1.1405639 0.8240748 -10.3728936 Erreur dans nls (ROS ~ b + a * exp(c * test), data = tab, start = list(b = 0, : gradient singulier De plus : Il y a eu 27 avis (utilisez warnings() pour les visionner)

---

Si je fais warnings(), j’obtiens :

Code:

---

1 2 3 4 5 6

Messages d'avis : 1: In lhs - rhs : la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court 2: In .swts * (lhs - rhs) : la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court 3: In .swts * attr(rhs, "gradient") :

---

………… …………………. …………………… ……………………
Merci pour votre aide

Bonjour et bienvenue,

Pour commencer, il y a une incohérence dans le nom de la data.frame utilisée

Citation:

---

Envoyé par Komine

Code:

---

1 2 3 4

tab<-read.table("C:\\Users\\Momadou\\Desktop\\TEST R\\regression.csv",header=TRUE,sep=";",dec=",") tab ... nlsreg <- nls(ROS ~ b+a*exp(c*test),data=Table,start=list(b=0,a=5,c=-0.05),trace=T)

---

---

Ensuite, si tu as toujours des problèmes de gradient singulier, joue sur les valeurs de départ : e.g. b=0.001, a=16/3,...

Merci Yves pour ta réponse, j'ai réglé le problème de la cohérance du dataframe mais le problème avec le message d'erreur reste le meme malgré le changement apporté sur les valeurs de a, b, et c.
Encore en attente d'un coup de pousse.
Merci

Bonjour,
Je suis un nouvel utilisateur de nls R. Je travaille sur des données. En fait, c'est sur l'évolution d'une mesure -(pc, pj, py) en fonction du temps.
Au tracé, on voit bien que les mesures baissent en fonction du temps. Les courbes ont des allures exponentielles decroissantes f(pc)<f(pj)<f(py).
J'ai choisi le modèle nls(pc~a*exp(-b*(pi/0.01)^2*temps)+c, start=list(a=1, b=10e-10,c=1). J’ai ainsi déterminé les paramètres a,b, et c des différents modèles.
Contrairement à ce qu'on observe, b de pc> à b de pj et b de py> à b de pj et pc.
Comment cela est-il possible?