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

  1. #1
    Futur Membre du Club
    PROC MIXED : afficher les estimates et les erreurs standards
    Bonjour,
    j'utilise actuellement un modèle mixte pour analyser des données répétées. La variable "album" est quantitative, il s'agit de mon critère de jugement. La variable "visite" qui correspond au temps est en classe (1-INCLUSION / 2-M6 / 3-M12). J'ai deux groupes de sujets (variable BRAS) qui est en classe aussi (Traitement / Placebo). Le code tel quel fonctionne mais je n'arrive pas à afficher les estimates et leurs erreurs standard à chaque temps. Pour les estimations il est facile de les calculer manuellement mais pour les erreurs je suis coincé. J'ai laissé dans le code, en commentaire, ce qui j'ai essayé mais qui ne fonctionne pas.

    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    PROC MIXED DATA=mytable NOCLPRINT=4 method=reml;
    	CLASS idSujet visite (ref="1-INCLUSION") bras;
    	MODEL album=visite bras visite*bras /solution DDFM=sat ;
    	RANDOM intercept /subject=idSujet type=un    ;
    	*estimate 'Traitement Inclusion' intercept 1  bras 1 0  ;
    	*estimate 'Placebo Inclusion' intercept 1  bras 0 1  ;
    	*estimate 'Traitement Inclusion' visite 1 0 0  bras 1 0  visite*bras 1 0 0 0 0 0;
    	*estimate 'Traitement M6' visite 1  bras 1 visite*bras 1 1;
    	*estimate 'Traitement  M12' visite 2  bras 1 visite*bras 2 1;
    RUN;


    Merci par avance !!

  2. #2
    Expert confirmé
    Bonjour.
    Il est possible que l'instruction
    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    SLICE visite*bras / SLICEBY=visite MEANS CL ;
    te fournisse les estimations recherchées (moyennes ajustées type LS-Means par visite sur chaque bras).
    Bon courage.
    Olivier

  3. #3
    Futur Membre du Club
    Merci cela fonctionne !

###raw>template_hook.ano_emploi###