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

SAS STAT Discussion :

MMRM (mixed model repeated measures)


Sujet :

SAS STAT

  1. #1
    Nouveau membre du Club
    Femme Profil pro
    Inscrit en
    Novembre 2011
    Messages
    46
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France, Val de Marne (Île de France)

    Informations forums :
    Inscription : Novembre 2011
    Messages : 46
    Points : 31
    Points
    31
    Par défaut MMRM (mixed model repeated measures)
    Bonjour,

    Pour une analyse de sensibilité du critère primaire je vais utiliser la méthode MMRM avec ( score as a covariate, visit , treatment groups and interaction treatment*visit as factors).
    Et le code sera
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    proc mixed data=table1
    class treatment visit patient;
    model y=  treatment visit covariable treatment*visit;
    repeated visit/subject=patient type=un;
    run;
    Maintenant on me demande de réaliser le même model "assuming a constant treatment effect over visits"

    Faut-il juste changer la matrice de variance covariance dans le code au dessus par type=cs?

    Merci pour votre aide

  2. #2
    Expert confirmé
    Avatar de olivier.decourt
    Homme Profil pro
    Formateur R/SAS/statistiques
    Inscrit en
    Avril 2008
    Messages
    2 064
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 46
    Localisation : France

    Informations professionnelles :
    Activité : Formateur R/SAS/statistiques
    Secteur : Conseil

    Informations forums :
    Inscription : Avril 2008
    Messages : 2 064
    Points : 4 478
    Points
    4 478
    Par défaut
    Non, le TYPE dans REPEATED indique la nature des corrélations entre répétitions d'un même sujet. Rien à voir (directement) avec l'effet du traitement, même si ça a évidemment ensuite une influence sur les tests effectués.

    J'aurais plutôt pensé que ce 2e modèle n'a pas d'interaction traitement*visite : on suppose que l'effet du traitement ne varie pas au fil du temps.
    Bon courage.
    Olivier

  3. #3
    Nouveau membre du Club
    Femme Profil pro
    Inscrit en
    Novembre 2011
    Messages
    46
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France, Val de Marne (Île de France)

    Informations forums :
    Inscription : Novembre 2011
    Messages : 46
    Points : 31
    Points
    31
    Par défaut
    J'avoue que je ne comprend pas :
    model "assuming a constant treatment effect over visits"

    je ne voyais que modifier la matrice de covariance.

    Mais supposons qu'on supprime le facteur "treatment*visit", en utilisant ce même code:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
     
    proc mixed data=table1
    class treatment visit patient;
    model y=  treatment visit covariable ;
    repeated visit/subject=patient type=un;
    run;
    comment je pourrais chercher l’effet du traitement à une visite donnée?
    car d'habitude si l’interaction est significative j'utilise:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    proc mixed data=table1
    class treatment visit patient;
    model y=  treatment visit covariable treatment*visit  ;
    repeated visit/subject=patient type=un
    lsmeans treatment*visit/slice=visit;
    run;
    J'espère que je suis claire

  4. #4
    Membre chevronné
    Homme Profil pro
    Biostatisticien
    Inscrit en
    Juin 2009
    Messages
    1 206
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Irlande

    Informations professionnelles :
    Activité : Biostatisticien
    Secteur : Industrie Pharmaceutique

    Informations forums :
    Inscription : Juin 2009
    Messages : 1 206
    Points : 1 868
    Points
    1 868
    Par défaut
    +1 Olivier

    Petite question pour clarification: pourquoi veux tu chercher l'effet du traitement a une visite donnee si tu supposes un effet traitement constant dans le temps?

  5. #5
    Expert confirmé
    Avatar de olivier.decourt
    Homme Profil pro
    Formateur R/SAS/statistiques
    Inscrit en
    Avril 2008
    Messages
    2 064
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 46
    Localisation : France

    Informations professionnelles :
    Activité : Formateur R/SAS/statistiques
    Secteur : Conseil

    Informations forums :
    Inscription : Avril 2008
    Messages : 2 064
    Points : 4 478
    Points
    4 478
    Par défaut
    Si tu enlèves l'interaction en supposant que l'effet traitement est constant à toutes les visites, tu vois l'effet du traitement à n'importe quelle visite avec
    Bon courage.
    Olivier

  6. #6
    Nouveau membre du Club
    Femme Profil pro
    Inscrit en
    Novembre 2011
    Messages
    46
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France, Val de Marne (Île de France)

    Informations forums :
    Inscription : Novembre 2011
    Messages : 46
    Points : 31
    Points
    31
    Par défaut
    Merci pour vos réponses pertinentes.

    En effet , j'ai les idées mélangées:

    pour récapituler:

    j'utilise ce code
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
     
    proc mixed data=table1
    class treatment visit patient;
    model y=  treatment visit covariable treatment*visit  ;
    repeated visit/subject=patient type=un
    lsmeans treatment*visit/slice=visit;
    run;
    pour voir si l'effet traitement varie au fil du temps

    et ce code

    j'utilise ce code
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
     
    proc mixed data=table1
    class treatment visit patient;
    model y=  treatment visit covariable  ;
    repeated visit/subject=patient type=un
    lsmeans treatment/pdiff;
    run;
    Si on suppose un effet traitement constant dans le temps

  7. #7
    Expert confirmé
    Avatar de olivier.decourt
    Homme Profil pro
    Formateur R/SAS/statistiques
    Inscrit en
    Avril 2008
    Messages
    2 064
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 46
    Localisation : France

    Informations professionnelles :
    Activité : Formateur R/SAS/statistiques
    Secteur : Conseil

    Informations forums :
    Inscription : Avril 2008
    Messages : 2 064
    Points : 4 478
    Points
    4 478
    Par défaut
    Je ne sais pas si c'est exactement ce qui était demandé mais en tout cas c'est la manière dont j'interprète la question, oui.
    Bon courage.
    Olivier
    Bon courage.
    Olivier

  8. #8
    Nouveau membre du Club
    Femme Profil pro
    Inscrit en
    Novembre 2011
    Messages
    46
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France, Val de Marne (Île de France)

    Informations forums :
    Inscription : Novembre 2011
    Messages : 46
    Points : 31
    Points
    31
    Par défaut
    Merci beaucoup Olivier.

  9. #9
    Membre chevronné
    Homme Profil pro
    Biostatisticien
    Inscrit en
    Juin 2009
    Messages
    1 206
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Irlande

    Informations professionnelles :
    Activité : Biostatisticien
    Secteur : Industrie Pharmaceutique

    Informations forums :
    Inscription : Juin 2009
    Messages : 1 206
    Points : 1 868
    Points
    1 868
    Par défaut
    Par contre, je ne suis encore qu'un debutant avec la methode MMRM (mixed model repeated measure), j'en viens donc la question: en quoi est-ce du MMRM?

Discussions similaires

  1. mixed models - Time serie
    Par CupOfAppleTea dans le forum R
    Réponses: 0
    Dernier message: 18/04/2013, 17h59
  2. calcul de puissance avec mixed repeated
    Par jnjnono dans le forum SAS STAT
    Réponses: 4
    Dernier message: 30/09/2010, 16h35
  3. Modele GEE : mixed ?
    Par mariebozo dans le forum SAS STAT
    Réponses: 1
    Dernier message: 20/10/2009, 11h11
  4. Repeater, Button et databind; mix ?
    Par darkendorf dans le forum Silverlight
    Réponses: 4
    Dernier message: 08/10/2009, 09h01
  5. MIXED Model ?
    Par jfdreyfus dans le forum R
    Réponses: 0
    Dernier message: 07/09/2009, 05h16

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