SAS/ETS procedures PROC MODEL and PROC AUTOREG both provide the functionality to specify Newey-West heteroscedasticity and autocorrelation standard error correction.
PROC MODEL provides the KERNEL= option with the GMM estimator in the FIT statement to correct standard errors for heteroscedasticity and autocorrelation. The procedure supports three different kernels — Bartlett, Parzen, and Quadratic.
The general form of the KERNEL= option is
KERNEL=( PARZEN | QS | BART, c, e )
where c and e are nonnegative values used to compute the bandwidth parameter
l(n) = cne
The Newey-West standard error correction is a commonly used heteroscedasticity and autocorrelation correction. The formula for the Newey-West covariance matrix estimator can be found in Greene (2000). The Newey-West estimator corresponds to the Bartlett kernel with bandwidth parameter L+1, where L is the maximum lag length. To specify the Newey-West kernel with lag length L using PROC MODEL, specify KERNEL=(BART, L+1, 0), which produces bandwidth parameter
l(n) = (L+1)n0 = L+1
For more details on the three kernels supported in PROC MODEL, see "Estimation Methods" in the Details section of the PROC MODEL documentation.
PROC AUTOREG provides the COVEST=HAC and COVEST=NEWEYWEST options in the MODEL statement to specify heteroscedasticity and autocorrelation correction. The COVEST=HAC option supports the following kernels with the KERNEL= option: BARTLETT, PARZEN, QUADRATICSPECTRAL, TRUNCATED, and TUKEYHANNING. To obtain the Newey-West covariance matrix estimator with maximum lag length L, specify the BARTLETT kernel with fixed value L+1 for the bandwidth parameter using COVEST=HAC(KERNEL=BARTLETT, BANDWIDTH=L+1). Alternatively, specify COVEST=NEWEYWEST(GAMMA=0, RATE=0, CONSTANT=L+1), where the bandwidth parameter is calculated as
b = [GAMMA*TRATE+CONSTANT], where [x] denotes the largest integer less than or equal to x, and T is the sample size.
For more details on the COVEST=HAC and COVEST=NEWEYWEST options, see "MODEL Statement" in the Syntax section in the PROC AUTOREG documentation.
These data are presented in Example 13.1 in Greene (2004). These statements create data set ONE with variables containing data on GNP, investment, price index, and nominal interest rate.
These steps create and display the real variables used in the regression analysis. Real GNP and real investment are obtained by dividing the nominal figures by the price index. An approximation to the real interest rate is obtained by subtracting the rate of change in the price index from the interest rate.
|
These statements regress real investment on real interest rate and real GNP using ordinary least squares (OLS) using PROC MODEL.
The parameter estimates and the uncorrected standard errors from the model are shown below.
|
To perform Newey-West standard error correction, PROC MODEL is run again specifying the GMM estimation method in the FIT statement. KERNEL=(BART, 5, 0) is also specified which requests the Bartlett kernel with a lag length of 4. The VARDEF=n option is specified to be consistent with the original Newey-West formula.
The parameter estimates and the Newey-West corrected standard errors from the above analysis appear below. Note that with GMM estimation, when the instruments used are all the exogenous variables in the OLS regression model, then the parameter estimates are identical to the OLS estimates. The results are consistent with table 13.2 in Greene (2000).
|
In the above example, the lag length of 4 is chosen to be the same as in Greene (2000). Neither PROC MODEL nor PROC AUTOREG provides an option to determine the maximum lag length in Newey West correction. Greene (2000) contains this discussion about choosing the maximum lag length:
"The maximum lag L must be determined in advance to be large enough that autocorrelations at lags longer than L are small enough to ignore. For a moving-average process, this can be expected to be a relatively small number. For autoregressive process or mixtures, however, the autocorrelations are never zero, and the researcher must make a judgment as to how far back it is necessary to go."
To specify the Newey-West standard error correction using PROC AUTOREG, either one of the following steps can be used. The same results are obtained as those produced in PROC MODEL:
or
To obtain the Newey-West standard error correction for the sample mean of a series, fit an intercept-only model to the series. The estimated intercept is the sample mean, and the Newey-West standard error correction can be obtained as in Example 1.
The following PROC MEANS step computes the sample mean of the variable R_INVEST.
|
These statements fit an intercept-only model to the real investment series and perform the Newey-West standard error correction on the sample mean. The intercept-only model is specified by fitting the equation R_INVEST=B0. The INTONLY option in the INSTRUMENTS statement indicates that the only instrument used in the estimation is the intercept.
The parameter estimates and the standard errors appear below. The B0 parameter estimate is the sample mean of the R_INVEST series. The reported standard errors are Newey-West corrected standard errors for the B0 estimate.
|
The same results can be obtained using the AUTOREG procedure with either the COVEST=HAC option or the COVEST=NEWEYWEST option:
or
____
Greene, William H. (2000), Econometric Analysis, Fourth Edition, Upper Saddle River, NJ: Prentice-Hall.
Product Family | Product | System | SAS Release | |
Reported | Fixed* | |||
SAS System | SAS/ETS | z/OS | ||
OpenVMS VAX | ||||
Microsoft® Windows® for 64-Bit Itanium-based Systems | ||||
Microsoft Windows Server 2003 Datacenter 64-bit Edition | ||||
Microsoft Windows Server 2003 Enterprise 64-bit Edition | ||||
Microsoft Windows XP 64-bit Edition | ||||
Microsoft® Windows® for x64 | ||||
OS/2 | ||||
Microsoft Windows 95/98 | ||||
Microsoft Windows 2000 Advanced Server | ||||
Microsoft Windows 2000 Datacenter Server | ||||
Microsoft Windows 2000 Server | ||||
Microsoft Windows 2000 Professional | ||||
Microsoft Windows NT Workstation | ||||
Microsoft Windows Server 2003 Datacenter Edition | ||||
Microsoft Windows Server 2003 Enterprise Edition | ||||
Microsoft Windows Server 2003 Standard Edition | ||||
Microsoft Windows Server 2003 for x64 | ||||
Microsoft Windows Server 2008 | ||||
Microsoft Windows Server 2008 for x64 | ||||
Microsoft Windows XP Professional | ||||
Windows 7 Enterprise 32 bit | ||||
Windows 7 Enterprise x64 | ||||
Windows 7 Home Premium 32 bit | ||||
Windows 7 Home Premium x64 | ||||
Windows 7 Professional 32 bit | ||||
Windows 7 Professional x64 | ||||
Windows 7 Ultimate 32 bit | ||||
Windows 7 Ultimate x64 | ||||
Windows Millennium Edition (Me) | ||||
Windows Vista | ||||
Windows Vista for x64 | ||||
64-bit Enabled AIX | ||||
64-bit Enabled HP-UX | ||||
64-bit Enabled Solaris | ||||
ABI+ for Intel Architecture | ||||
AIX | ||||
HP-UX | ||||
HP-UX IPF | ||||
IRIX | ||||
Linux | ||||
Linux for x64 | ||||
Linux on Itanium | ||||
OpenVMS Alpha | ||||
OpenVMS on HP Integrity | ||||
Solaris | ||||
Solaris for x64 | ||||
Tru64 UNIX |
Type: | Usage Note |
Priority: | |
Topic: | SAS Reference ==> Procedures ==> MODEL Analytics ==> Econometrics Analytics ==> Time Series Analysis SAS Reference ==> Procedures ==> AUTOREG |
Date Modified: | 2019-11-05 15:28:25 |
Date Created: | 2010-06-23 16:25:40 |