Discussion :

[deltamoins]

Bonjour les amis,

C'est la première fois que je viens du côté de "R", j'ai un soucis avec une fonction. Je veux ajuster par moindres carrés un modèle donné et pour cela j'utilise une fonction de régression avec une variance hétérogène via gnls. Mais j'ai toujours une erreur qui apparaît.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
# les declarations
# le fichier dat contient les températures dans la première colonne et les poids dans la seconde colonne.
ble=read.table("ble.dat",header=TRUE)
attach(ble)
library(nlme)
 
# la fonction de régression avec deux paramètres : theta 1 et theta 2 
Freg = function(x,theta) {
  theta[1]*exp((theta[2]*x))
}
 
# ajustement par moindres carrés de la fonction 
ajust3 = gnls(poids~Freg(temperature,theta),data=ble,start=list(theta = c(75,0.002)),weights=varPower())

En lançant ajust3 j’obtiens cette erreur :

Erreur dans model.frame.default(formula = ~. + poids + temperature + theta, :
les longueurs des variables diffèrent (trouvé pour 'theta')

Voilà, j'espère que ça vous inspire et que vous avez une petite idée
Merci ^^

[deltamoins]

Alors le problème est RESOLU, il fallait faire ceci :

#changer la forme de la déclaration de la fonction :

Code :

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

1
2
3
Freg2 = function(x,theta1,theta2) {
  theta1*exp((theta2*x))
}

#et faire l'appel de la fonction dans la méthode des moindres carres comme ça :

Code :

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

ajustGNLS2 = gnls(poids~Freg2(temperature,theta1,theta2),data=ble,start=c(theta1=75,theta2=0.002),weights=varPower())