vendredi 15 septembre 2017

Estimating outcome with Proc Mixed involving random slope

Hi I have the following proc mix statement which has both time(fixed eff as categorical i.e. visit) and time (random eff as continuous, i.e. days since baseline). I need to estimate the mean outcome at day 252 for treatment and placebo separately and the difference between the two. I am getting confused how to write the ESTIMATE statement. Is this right?

proc mixed;
  class usubjid treatment visit;
  model Y  = Baseline_Score treatment visit
                            treatment*visit /ddfm=kr;
  random intercept days_since_baseline /  sub = usubjid type = un;
  ESTIMATE ‘Active mean’ intercept 1 days_since_baseline 252 treatment 1 0/CL;
  ESTIMATE ‘placebo mean’ intercept 1 days_since_baseline 252 treatment 0 1/CL;
  ESTIMATE ‘DIFF’ days_since_baseline 252 treatment 1 -1 /CL;

run;

How do you account for the interactions? Thank you in advance.




Aucun commentaire:

Enregistrer un commentaire