DATA HOSP; INPUT X1 X2 X3 X4 X5 Y; /* Defines variables */ CARDS; 15.57 2463 472.92 18.0 4.45 566.52 44.02 2048 1339.75 9.5 6.92 696.82 20.42 3940 620.25 12.8 4.28 1033.15 18.74 6505 568.33 36.7 3.90 1603.62 49.20 5723 1497.60 35.7 5.50 1611.37 44.92 11520 1365.83 24.0 4.60 1613.27 55.48 5779 1687.00 43.3 5.62 1854.17 59.28 5969 1639.92 46.7 5.15 2160.55 94.39 8461 2872.33 78.7 6.18 2305.58 128.02 20106 3655.08 180.5 6.15 3503.93 96.00 13313 2912.00 60.9 5.88 3571.89 131.42 10771 3921.00 103.7 4.88 3741.40 127.21 15543 3865.67 126.8 5.50 4026.52 252.90 36194 7684.10 157.7 7.00 10343.81 409.20 34703 12446.33 169.4 10.78 11732.17 463.70 39204 14098.40 331.4 7.05 15414.94 510.22 86533 15524.00 371.6 6.35 18854.45 ; PROC PRINT; RUN; /* define symbol characteristics */ symbol interpol=NONE /* regression analysis with */ value=diamond /* plot symbol */ height=3 /* plot symbol height */ cv=red /* plot symbol color */ ci=blue /* regression line color */ co=green /* confidence limits color */ width=2; /* line width */ /* produce plot */ PROC GPLOT DATA=HOSP; PLOT Y * (X1 X2 X3 X4 X5); /* Plots y against X1, X2, X3, X4, and X5. */ RUN; PROC CORR; /* Prints correlation matrix */ RUN; PROC REG; MODEL Y = X1 X2 X3 X4 X5/PARTIAL; /*Gives partial leverage residual plots */ RUN; /* VIF calculates variance*/ PROC REG; /*inflation factors for model*/ MODEL Y = X2 X3 X5/P VIF; /*Y = b0 + b1x2 + b2x3 + b3x5 + e*/ RUN;