The REG Procedure

OUTPUT Statement

OUTPUT <OUT=SAS-data-set> <keyword=names> <…keyword=names> ;

The OUTPUT statement creates a new SAS data set that saves diagnostic measures calculated after fitting the model. The OUTPUT statement refers to the most recent MODEL statement. At least one keyword=names specification is required.

All the variables in the original data set are included in the new data set, along with variables created in the OUTPUT statement. These new variables contain the values of a variety of statistics and diagnostic measures that are calculated for each observation in the data set. If you want to create a SAS data set in a permanent library, you must specify a two-level name. For more information about permanent libraries and SAS data sets, see SAS Language Reference: Concepts.

The OUTPUT statement cannot be used when a TYPE=CORR, TYPE=COV, or TYPE=SSCP data set is used as the input data set for PROC REG. See the section Input Data Sets for more details.

The statistics created in the OUTPUT statement are described in this section. More details are given in the section Predicted and Residual Values and the section Influence Statistics. Also see Chapter 4: Introduction to Regression Procedures, for definitions of the statistics available from the REG procedure.

You can specify the following options in the OUTPUT statement:

OUT=SAS data set

gives the name of the new data set. By default, the procedure uses the DATAn convention to name the new data set.

keyword=names

specifies the statistics to include in the output data set and names the new variables that contain the statistics. Specify a keyword for each desired statistic (see the following list of keywords), an equal sign, and the variable or variables to contain the statistic.

In the output data set, the first variable listed after a keyword in the OUTPUT statement contains that statistic for the first dependent variable listed in the MODEL statement; the second variable contains the statistic for the second dependent variable in the MODEL statement, and so on. The list of variables following the equal sign can be shorter than the list of dependent variables in the MODEL statement. In this case, the procedure creates the new names in order of the dependent variables in the MODEL statement.

For example, the following SAS statements create an output data set named b:

proc reg data=a;
   model y z=x1 x2;
   output out=b
      p=yhat zhat
      r=yresid zresid;
run;

In addition to the variables in the input data set, b contains the following variables:

  • yhat, with values that are predicted values of the dependent variable y

  • zhat, with values that are predicted values of the dependent variable z

  • yresid, with values that are the residual values of y

  • zresid, with values that are the residual values of z

You can specify the following keywords in the OUTPUT statement. See the section Model Fit and Diagnostic Statistics for computational formulas.

Table 83.6: Keywords for OUTPUT Statement

Keyword

Description

COOKD=names

Cook’s D influence statistic

COVRATIO=names

standard influence of observation on covariance of betas, as
discussed in the section Influence Statistics

DFFITS=names

standard influence of observation on predicted value

H=names

leverage, $\mb {x}_ i(\mb {X}’\mb {X})^{-1}\mb {x}_ i^{\prime }$

LCL=names

lower bound of a $100(1-\alpha )$% confidence interval for an
individual prediction. This includes the variance of the
error, as well as the variance of the parameter estimates.

LCLM=names

lower bound of a $100(1-\alpha )$% confidence interval for the
expected value (mean) of the dependent variable

PREDICTED | P=names

predicted values

PRESS=names

ith residual divided by $(1-h)$, where h is the leverage,
and where the model has been refit without the ith
observation

RESIDUAL | R=names

residuals, calculated as ACTUAL minus PREDICTED

RSTUDENT=names

a studentized residual with the current observation deleted

STDI=names

standard error of the individual predicted value

STDP=names

standard error of the mean predicted value

STDR=names

standard error of the residual

STUDENT=names

studentized residuals, which are the residuals divided by their
standard errors

UCL=names

upper bound of a $100(1-\alpha )$% confidence interval for an
individual prediction

UCLM=names

upper bound of a $100(1-\alpha )$% confidence interval for the
expected value (mean) of the dependent variable