The REG Procedure

Influence Statistics

This section discusses the INFLUENCE option, which produces several influence statistics, and the PARTIAL option, which produces partial regression leverage plots.

The INFLUENCE Option

The INFLUENCE option (in the MODEL statement) requests the statistics proposed by Belsley, Kuh, and Welsch (1980) to measure the influence of each observation on the estimates. Influential observations are those that, according to various criteria, appear to have a large influence on the parameter estimates.

Let $\mb{b}(i)$ be the parameter estimates after deleting the ith observation; let $s(i)^2$ be the variance estimate after deleting the ith observation; let $\mb{X}(i)$ be the $\mb{X}$ matrix without the ith observation; let $\hat{y}(i)$ be the ith value predicted without using the ith observation; let $r_ i = y_ i - \hat{y}_ i$ be the ith residual; and let $h_ i$ be the ith diagonal of the projection matrix for the predictor space, also called the hat matrix:

\[ h_ i = \mb{x}_ i(\mb{X}’\mb{X})^{-1}\mb{x}_ i^{\prime } \]

Belsley, Kuh, and Welsch (1980) propose a cutoff of $2p/n$, where n is the number of observations used to fit the model and p is the number of parameters in the model. Observations with $h_ i$ values above this cutoff should be investigated.

For each observation, PROC REG first displays the residual, the studentized residual (RSTUDENT), and the $h_ i$. The studentized residual RSTUDENT differs slightly from STUDENT since the error variance is estimated by $s_{(i)}^2$ without the ith observation, not by $s^2$. For example,

\[ \mbox{RSTUDENT} = \frac{r_ i}{s_{(i)} \sqrt {(1 - h_ i)}} \]

Observations with RSTUDENT larger than 2 in absolute value might need some attention.

The COVRATIO statistic measures the change in the determinant of the covariance matrix of the estimates by deleting the ith observation:

\[ \mbox{COVRATIO} = \frac{ \mbox{det} \left( s^2(i) (\mb{X}_{(i)}^{\prime }\mb{X}_{(i)})^{-1} \right) }{ \mbox{det} \left( s^2 (\mb{X}'\mb{X})^{-1} \right) } \]

Belsley, Kuh, and Welsch (1980) suggest that observations with

\[ |\mbox{COVRATIO} - 1| \ge \frac{3p}{n} \]

where p is the number of parameters in the model and n is the number of observations used to fit the model, are worth investigation.

The DFFITS statistic is a scaled measure of the change in the predicted value for the ith observation and is calculated by deleting the ith observation. A large value indicates that the observation is very influential in its neighborhood of the $\mb{X}$ space.

\[ \mbox{DFFITS} = \frac{\hat{y}_ i - \hat{y}_{(i)}}{s_{(i)} \sqrt {h_{(i)}}} \]

Large values of DFFITS indicate influential observations. A general cutoff to consider is 2; a size-adjusted cutoff recommended by Belsley, Kuh, and Welsch (1980) is $2\sqrt {p/n}$, where n and p are as defined previously.

The DFFITS statistic is very similar to Cook’s D, defined in the section Predicted and Residual Values.

The DFBETAS statistics are the scaled measures of the change in each parameter estimate and are calculated by deleting the ith observation:

\[ \mbox{DFBETAS}_ j = \frac{b_ j - b_{(i)j}}{s_{(i)} \sqrt {(\mb{X}’\mb{X})_{jj}}} \]

where $(\mb{X}’\mb{X})_{jj}$ is the $(j,j)$ element of $(\mb{X}’\mb{X})^{-1}$.

In general, large values of DFBETAS indicate observations that are influential in estimating a given parameter. Belsley, Kuh, and Welsch (1980) recommend 2 as a general cutoff value to indicate influential observations and $2/\sqrt {n}$ as a size-adjusted cutoff.

The following statements use a subset of the data in the population example in the section Polynomial Regression. The INFLUENCE option produces the tables shown in Figure 97.38 and Figure 97.39.

proc reg data=USPopulation;
   where Year <= 1970;
   model Population=Year YearSq / influence;
run;

Figure 97.38: Regression Using the INFLUENCE Option

The REG Procedure
Model: MODEL1
Dependent Variable: Population

Output Statistics
Obs Residual RStudent Hat Diag
H
Cov
Ratio
DFFITS DFBETAS
Intercept Year YearSq
1 -1.1094 -0.4972 0.3865 1.8834 -0.3946 -0.2842 0.2810 -0.2779
2 0.2691 0.1082 0.2501 1.6147 0.0625 0.0376 -0.0370 0.0365
3 0.9305 0.3561 0.1652 1.4176 0.1584 0.0666 -0.0651 0.0636
4 0.7908 0.2941 0.1184 1.3531 0.1078 0.0182 -0.0172 0.0161
5 0.2110 0.0774 0.0983 1.3444 0.0256 -0.0030 0.0033 -0.0035
6 -0.6629 -0.2431 0.0951 1.3255 -0.0788 0.0296 -0.0302 0.0307
7 -0.8869 -0.3268 0.1009 1.3214 -0.1095 0.0609 -0.0616 0.0621
8 -0.2501 -0.0923 0.1095 1.3605 -0.0324 0.0216 -0.0217 0.0218
9 -0.7593 -0.2820 0.1164 1.3519 -0.1023 0.0743 -0.0745 0.0747
10 -0.5757 -0.2139 0.1190 1.3650 -0.0786 0.0586 -0.0587 0.0587
11 0.7938 0.2949 0.1164 1.3499 0.1070 -0.0784 0.0783 -0.0781
12 1.1492 0.4265 0.1095 1.3144 0.1496 -0.1018 0.1014 -0.1009
13 3.1664 1.2189 0.1009 1.0168 0.4084 -0.2357 0.2338 -0.2318
14 1.6746 0.6207 0.0951 1.2430 0.2013 -0.0811 0.0798 -0.0784
15 2.2406 0.8407 0.0983 1.1724 0.2776 -0.0427 0.0404 -0.0380
16 -6.6335 -3.1845 0.1184 0.2924 -1.1673 -0.1531 0.1636 -0.1747
17 -6.0147 -2.8433 0.1652 0.3989 -1.2649 -0.4843 0.4958 -0.5076
18 1.6770 0.6847 0.2501 1.4757 0.3954 0.2240 -0.2274 0.2308
19 3.9895 1.9947 0.3865 0.9766 1.5831 1.0902 -1.1025 1.1151



Figure 97.39: Residual Statistics

Sum of Residuals -5.8175E-11
Sum of Squared Residuals 123.74557
Predicted Residual SS (PRESS) 188.54924



In Figure 97.38, observations 16, 17, and 19 exceed or are near the cutoff value of 2 for RSTUDENT. None of the observations exceeds the general cutoff of 2 for DFFITS or the DFBETAS, but observations 16, 17, and 19 exceed at least one of the size-adjusted cutoffs for these statistics. Observations 1 and 19 exceed the cutoff for the hat diagonals, and observations 1, 2, 16, 17, and 18 exceed the cutoffs for COVRATIO. Taken together, these statistics indicate that you should look first at observations 16, 17, and 19 and then perhaps investigate the other observations that exceeded a cutoff.

When ODS Graphics is enabled, you can request influence diagnostic plots by using the PLOTS= option in the PROC REG statement as shown in the following statements:

ods graphics on;

proc reg data=USPopulation
      plots(label)=(CooksD RStudentByLeverage DFFITS DFBETAS);
   where Year <= 1970;
   id Year;
   model Population=Year YearSq;
run;

The LABEL suboption specified in the PLOTS(LABEL)= option requests that observations that exceed the relevant cutoffs for the statistics being plotted are labeled. Since Year has been named in an ID statement, the value of Year is used for the labels. The requested plots are shown in Figure 97.40.

Figure 97.40: Influence Diagnostics

Influence Diagnostics
External File:images/dtInfluencePlots1.png
External File:images/dtInfluencePlots2.png
External File:images/dtInfluencePlots3.png


When ODS Graphics is enabled, you can request a chart of the studentized residuals and Cook’s D by using the PLOTS=RESIDUALCHART option as shown in the following statements:

proc reg data=USPopulation plots(only)=residualchart;
   where Year <= 1970;
   id Year;
   model Population=Year YearSq;
run;

The chart is shown in Figure 97.41.

Figure 97.41: Studentized Residuals and Cook’s D Chart

Studentized Residuals and Cook’s  Chart


The PARTIAL and PARTIALDATA Options

The PARTIAL option in the MODEL statement produces partial regression leverage plots. If ODS Graphics is not enabled, this option requires the use of the LINEPRINTER option in the PROC REG statement. One plot is created for each regressor in the current full model. For example, plots are produced for regressors included by using ADD statements; plots are not produced for interim models in the various model-selection methods but only for the full model. If you use a model-selection method and the final model contains only a subset of the original regressors, the PARTIAL option still produces plots for all regressors in the full model. If ODS Graphics is enabled, these plots are produced as high-resolution graphics, in panels with a maximum of six partial regression leverage plots per panel. Multiple panels are displayed for models with more than six regressors.

For a given regressor, the partial regression leverage plot is the plot of the dependent variable and the regressor after they have been made orthogonal to the other regressors in the model. These can be obtained by plotting the residuals for the dependent variable against the residuals for the selected regressor, where the residuals for the dependent variable are calculated with the selected regressor omitted, and the residuals for the selected regressor are calculated from a model where the selected regressor is regressed on the remaining regressors. A line fit to the points has a slope equal to the parameter estimate in the full model.

When ODS Graphics is not enabled, points in the plot are marked by the number of replicates appearing at one position. The symbol '*' is used if there are 10 or more replicates. If an ID statement is specified, the leftmost nonblank character in the value of the ID variable is used as the plotting symbol.

The PARTIALDATA option in the MODEL statement produces a table that contains the partial regression data that are displayed in the partial regression leverage plots. You can request partial regression data even if you do not requests plots with the PARTIAL option.

The following statements use the fitness data in Example 97.2 with the PARTIAL option and ODS Graphics to produce the partial regression leverage plots. The plots are shown in Figure 97.42.

ods graphics on;

proc reg data=fitness;
   model Oxygen=RunTime Weight Age / partial;
run;

ods graphics off;

Figure 97.42: Partial Regression Leverage Plots

 Partial Regression Leverage Plots