Bonjour,

Je débute tout juste R.
Jusqu'à maintenant j'avais juste fait tourner des programmes qu'on m'avait passés...

Je dois réaliser un fit de la fonction DL0 = Deq*(1-alpha*exp(-4*L/Leq))
J'ai un jeu de valeurs L/DL0 et je cherche alpha, Deq et Leq.

J'ai donc voulu utiliser la fonction nls mais je ne dois pas réussir à rentrer les infos comme il faut, ça me retourne l'erreur "erreur gradient singulier".

Voici le mini code que j'ai fait :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
setwd("C:/Users/pierre/Desktop/R")
a1<-read.table('toshiba.txt',col.names=c('L','DL0'))
attach(a1)
result_fit <- nls(DL0 ~ Deq*(1-alpha*exp(-4*L/Leq)), data = a1, start = list(Deq = 1, alpha = 0.5, Leq = 40))
Par ailleurs, j'ai un avertissement comme quoi la dernière ligne du fichier toshiba.txt est incomplète, mais elle me semble normal (après ce problème ne l'empêche pas de tourner).

Quel est le problème ?