Practical Multivariate Analysis, Fifth Edition, by Afifi, May and Clark Chapter 6: Simple regression and correlation | SAS Textbook Examples

Page 83 Figure 6.1  Scatter diagram and regression line of FEV1 versus height for fathers in lung function data.

NOTE:  Throughout this chapter and the next, the variables FFEV1 and MFEV1 have been divided by 100 before use in analyses.  We have named these modified variables FFEV1a and MFEV1a, respectively.

data lung;
set "c:\pma5\lung";
ffev1a = ffev1/100;
run;
goptions reset = all;
symbol1 i=r v=circle;
axis1 order=(58 to 78 by 4) label=('fheight');
axis2 order=(2 to 6.5 by .5);
proc gplot data = lung;
plot ffev1a*fheight / haxis=axis1 vaxis=axis2;
run;
quit;

Page 84 The descriptive statistics at the top of the page.

proc means data = lung mean std;
var ffev1a fheight;
run;

The MEANS Procedure

Variable    Label                                 Mean         Std Dev
----------------------------------------------------------------------
ffev1a                                       4.0932667       0.6507523
FHEIGHT     height of father in inches      69.2600000       2.7791892
----------------------------------------------------------------------

Page 94 The correlation in the middle of the page

proc corr data = lung;
var ffev1a fheight;
run;
<some output omitted>
      Pearson Correlation Coefficients, N = 150
              Prob > |r| under H0: Rho=0

                                  ffev1a       FHEIGHT

ffev1a                           1.00000       0.50440
                                                <.0001

FHEIGHT                          0.50440       1.00000
height of father in inches        <.0001

Page 96 Table 6.2  ANOVA example from Figure 6.1.

NOTE:  We need to generate the residuals for use in the next figure.

proc reg data = lung;
model ffev1a = fheight;
output out = resid r=resid;
run;
quit;

<some output omitted>
The REG Procedure
Model: MODEL1
Dependent Variable: ffev1a

                             Analysis of Variance

                                    Sum of           Mean
Source                   DF        Squares         Square    F Value    Pr > F

Model                     1       16.05317       16.05317      50.50    <.0001
Error                   148       47.04513        0.31787
Corrected Total         149       63.09830

Page 109 Figure 6.10  This has been skipped for now.

Page 113 Figure 6.12  Normal probability plot of the residuals of the regression of FEV1 on height for fathers in the lung function data.

NOTE:  The graph produced by proc univariate has the axes reversed from that shown in the text.

proc univariate data = resid plots;
var resid;
run;

The UNIVARIATE Procedure
Variable:  resid  (Residual)

                       Normal Probability Plot
     1.5+                                                * *
        |                                              *+++
        |                                          ****+
     0.9+                                        **+
        |                                    ****
        |                                 ****
     0.3+                             *****
        |                        *****
        |                     ****
    -0.3+                  ****
        |               ****
        |           *****
    -0.9+        +**
        |    +****
        | +*+*
    -1.5+*
         +----+----+----+----+----+----+----+----+----+----+
             -2        -1         0        +1        +2