Paired t-test and proc mianalyze using Rubin's method

[Liamb]

Hello,

After using this syntax on imputed data (10 imputations), I can not retrieve parameterEstimates for the 10 imputations and then use Rubin's method for the global analysis. Is it possible? If yes, how to do it?

proc ttest data=ttest_cd;
paired cd1*cd2;
by _imputation_;
ods output parameterEstimates=diff ;
run;

Thanks

[SAS_Rob · Accepted solution]

You would need the Statistics ODS table from TTEST to do this.  Below is a full example you should be able to follow.

/*Create a Sample dataset*/

data a;
do rep=1 to 40;
if ranuni(3211)>.6 then y1=.;
else y1=1.6+rannor(3211);
if ranuni(3211)>.8 then y2=.;
else y2=2.3+rannor(3211);
output;
end;
run;

proc mi data=a seed=10 out=outmi;
var y1 y2;
run;

proc ttest data=outmi;
by _imputation_;
paired y1*y2;
ods output statistics=t_stats;
run;

proc mianalyze data=t_stats;
modeleffects mean;
stderr stderr;
run;

[Liamb]

Thanks a lot Rob