PROC MIXED : Calculer les pentes quand l'interaction avec le temps est significative

[oualina]

Bonjour à tous,
Suis nouvelle au forum et excuse de commencer par une question.

J'ai écrit un modèle avec proc mixed à mesures répétées.

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
ods Output SolutionF=SolF01;
proc mixed DATA=aaar method=ml empirical ORDER=internal;
   class Sbjid VisNbr ;
   repeated VisNbr /  subject=SbjId type=VC r;
 
    lsmeans VisNbr;
 
    model rTotVolr= Visnbr|ictp_i Visnbr|p3np_old  Visnbr|MMP3_i  Visnbr|MMP_i  /* modele sans fMMP9_i*/
                   Visnbr|timp2_i /s  residual outp=r01(keep=Sbjid Visnbr Pred--PearsonResid) ;
 
  *estimate 'Effet ICTP at Incl' intercept 1 
                                Visnbr 1 0
                                ictp_i 1
                                Visnbr*ictp_i ;
 
  WHERE VisNbr IN (0 1);
  title ' ';
run;
 
symbol1 v=dot c=red;
symbol2 v=star c=blue;
 
goptions reset=ALL;
proc gplot DATA=r01;
  plot resid*pred=1 StudentResid*pred=2 ;
run;

J'ai un effet significatif entre le temps (variable VisNbr) et les variables biomarkeurs. Voudrait écrire mes equations de modèle car voudrait calculer un score, ne sait pas si je dois utiliser l'estimate mais par exemple voudrait calculer les deux pentes associées à l'effet de ICTP_i aux temps Visnbr=0, 1.
Comment faire?
Merci de votre aide...

[Manoutz]

Tu peux obtenir les estimates de tes différents paramètres dans la partie "solution pour effets fixes". Je te conseilles également d'utiliser l'option noint.

[oualina]

Merci de ta réponse,

J'ai mes estimates mais il considere la reference 0 a temps = 1an, est-ce que c'est legitime d'enlever l'intercept comme ca sans justification?

[Manoutz]

C'est justement pour cette raison de "référence" associée à un groupe que je préfère retirer les interceptes lorsque le contexte et l'objectif le permettent. Cela ne reflète que mon opinion et je n'ai par ailleurs pas tous les prérequis du problème. Je suppose que le choix de retenir cette écriture (en particulier pour les effets fixes) vient de ton objectif (?).