Page 87 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:\cama3\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 88 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 99 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 101 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 116 Figure 6.11 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