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