Discussion :

[Vladoz]

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 !!

[olivier.decourt]

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).

[Vladoz]

Merci cela fonctionne !