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 : Sélectionner tout - Visualiser dans une fenêtre à part
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 : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
plot(a1$deltaT~a1$latitude, data=a1)
abline(m522)
Mais j'obtiens le message d'erreur suivant:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
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!