IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

R Discussion :

Graph_Droite de régression_Modèles mixtes_lme()


Sujet :

R

  1. #1
    Nouveau membre du Club
    Profil pro
    Inscrit en
    Mars 2012
    Messages
    69
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Mars 2012
    Messages : 69
    Points : 38
    Points
    38
    Par défaut 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 : 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!

  2. #2
    Membre confirmé
    Homme Profil pro
    MCU
    Inscrit en
    Juillet 2010
    Messages
    185
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Nord (Nord Pas de Calais)

    Informations professionnelles :
    Activité : MCU
    Secteur : Santé

    Informations forums :
    Inscription : Juillet 2010
    Messages : 185
    Points : 486
    Points
    486
    Par défaut
    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 : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    library(nlme)
    example(lme)
    plot(distance ~ age, data=Orthodont)
    abline(a=fixef(fm2)[1], b=fixef(fm2)[2])
    HTH

    Vincent

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo