Hello,
I make multiple imputation with mixed model and I wish combine results for p-value of global interaction. I saw that there might be a macro, is it still relevant? how can I proceed ?
Here, my code :
proc mixed data=v1_v2_im method=reml ;
class prod (ref="Placebo") TV_VISIT_TMP(ref="V1") id_bis;
model changeVn_V0 = prod TV_VISIT_TMP prod*TV_VISIT_TMP / s ddfm=kenwardroger;
repeated TV_VISIT_TMP / subject=id_bis type=CS;
lsmeans prod*TV_VISIT_TMP / slice=TV_VISIT_TMP cl;
slice prod*TV_VISIT_TMP / sliceby(TV_VISIT_TMP="V1") cl pdiff=control("Placebo" "V1");
slice prod*TV_VISIT_TMP / sliceby(TV_VISIT_TMP="V2") cl pdiff=control("Placebo" "V2");
by simul _Imputation_;
ods output SolutionF=mixparms ;
ods output SliceDiffs=slicev1v2;
run;
proc mianalyze parms(classvar=full)=mixparms ;
class prod TV_VISIT_TMP;
modeleffects Intercept prod TV_VISIT_TMP prod*TV_VISIT_TMP ;
by simul;
ods output ParameterEstimates=param_simul;
run;
Thanks,
Clemence