[problème utilisation] GNLS : moindres carrés généralisés

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:

---

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 :mrgreen:
Merci ^^

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

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

Code:

---

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:

---

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

---