Documentation Example 10 for PROC MIANALYZE
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: MIANAX10 */
/* TITLE: Documentation Example 10 for PROC MIANALYZE */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: multiple imputation */
/* PROCS: MI, MIANALYZE, CORR */
/* DATA: */
/* */
/* SUPPORT: Yang Yuan UPDATE: Aug 21, 2007 */
/* REF: PROC MIANALYZE, EXAMPLE 10 */
/* MISC: */
/****************************************************************/
*----------------- Data on Physical Fitness -----------------*
| These measurements were made on men involved in a physical |
| fitness course at N.C. State University. |
| Only selected variables of |
| Oxygen (oxygen intake, ml per kg body weight per minute), |
| Runtime (time to run 1.5 miles in minutes), and |
| RunPulse (heart rate while running) are used. |
| Certain values were changed to missing for the analysis. |
*------------------------------------------------------------*;
data Fitness1;
input Oxygen RunTime RunPulse @@;
datalines;
44.609 11.37 178 45.313 10.07 185
54.297 8.65 156 59.571 . .
49.874 9.22 . 44.811 11.63 176
. 11.95 176 . 10.85 .
39.442 13.08 174 60.055 8.63 170
50.541 . . 37.388 14.03 186
44.754 11.12 176 47.273 . .
51.855 10.33 166 49.156 8.95 180
40.836 10.95 168 46.672 10.00 .
46.774 10.25 . 50.388 10.08 168
39.407 12.63 174 46.080 11.17 156
45.441 9.63 164 . 8.92 .
45.118 11.08 . 39.203 12.88 168
45.790 10.47 186 50.545 9.93 148
48.673 9.40 186 47.920 11.50 170
47.467 10.50 170
;
proc mi data=Fitness1 seed=3237851 noprint out=outmi;
var Oxygen RunTime RunPulse;
run;
proc corr data=outmi fisher(biasadj=no);
var Oxygen RunTime;
by _Imputation_;
ods output FisherPearsonCorr= outz;
run;
proc print data=outz;
title 'Fisher''s Correlation Statistics';
var _Imputation_ NObs ZVal;
run;
data outz;
set outz;
StdZ= 1. / sqrt(NObs-3);
run;
proc mianalyze data=outz;
ods output ParameterEstimates=parms;
modeleffects ZVal;
stderr StdZ;
run;
proc print data=parms;
title 'Parameter Estimates with 95% Confidence Limits';
var Estimate LCLMean UCLMean;
run;
data corr_ci;
set parms;
r= tanh( Estimate);
r_lower= tanh( LCLMean);
r_upper= tanh( UCLMean);
run;
proc print data=corr_ci;
title 'Estimated Correlation Coefficient'
' with 95% Confidence Limits';
var r r_lower r_upper;
run;