Discussion :

[Komine]

Bonjour,

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[pitipoisson]

Bonjour et bienvenue,

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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,...

[Komine]

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

[kbmad]

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?