Graph_Droite de régression_Modèles mixtes_lme()

Bonjour,
J'ai recours aux modèles mixtes pour traiter mes données. J'utilise la fonction "lme" du package "nlme".
Mon code est le suivant:

Code:

---

1 2	a1=read.table("Zenith12_16deltaTTIMETC.txt",header=TRUE) m522=lme(deltaT~1+latitude+ftype, random=~0+ftype|fID2, na.action=na.omit, method="ML", data=a1)

---

Je souhaite obtenir un graph avec mes données et la droite de régression. Pour cela je procède comme suit:

Code:

---

1 2	plot(a1$deltaT~a1$latitude, data=a1) abline(m522)

---

Mais j'obtiens le message d'erreur suivant:

Code:

---

1 2 3 4 5 6

> abline(m522) Erreur dans int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : l'objet (list) ne peut être couverti automatiquement en un type 'double' De plus : Warning message: In abline(m522) : utilisation de seulement les deux premiers de 4coefficients de régression

---

Lorsque j'utilise d'autres type de régression (glm par exemple) je ne rencontre pas ce pb. Je suis donc conscient que le problème doit être lié au fait que j'utilise les modèles mixtes mais je ne comprends pas bien pourquoi!

Merci pour votre aide!

Bonjour,

D'après la doc de la fonction abline, en cas d'utilisation d'un objet de régression, il faut que celui-ci dispose d'une méthode coef... Pour les modèles mixtes, coef retourne les coefficients par groupe (terme à droite du "|"). Vous pourriez utiliser fixef pour extraire les effets fixes du modèle mixte et donner ça comme argument à abline :

Code:

---

1 2 3 4

library(nlme) example(lme) plot(distance ~ age, data=Orthodont) abline(a=fixef(fm2)[1], b=fixef(fm2)[2])

---

HTH

Vincent