Data Missing;
                Infile 'cancerdot.dat';
                Input SexP DeptP AnxtP GSItP DeptS AnxtS GSItS SexChild Totbpt;
run;
proc means data = missing;
        var SexP DeptP AnxtP GSItP DeptS AnxtS GSItS Sexchild Totbpt;
    run;

Proc Reg data = Missing;
        Model Totbpt = SexP DeptP AnxtP DeptS AnxtS;
        Run;

Proc MI Data = Missing out = miout seed = 35399 nimpute = 5;
        Var SexP DeptP AnxtP  DeptS AnxtS   Totbpt;/*GSItP GSItS SexChild*/
        EM;
        MCMC NBITER = 500 NITER = 500;

Run;
Proc Reg data = miout outest = outreg covout ;
        Model Totbpt = SexP DeptP AnxtP DeptS AnxtS;
        by _Imputation_;
Run;
Proc print data = outreg;
        var _Imputation_ _Type_ _Name_
        intercept sexp deptp anxtp depts anxts;
        Title 'Parameter estimates from imputed data sets';
run;

proc MIANALYZE data = outreg;
        ModelEffects  SexP DeptP AnxtP DeptS AnxtS intercept;
Run;