Using gret l for Principles of Econometrics, 4th Edition

Bandwidth and Kernel

HAC is not quite as automatic as the heteroskedasticity consistent (HCCME) estimator in chapter 8. To be robust with respect to autocorrelation you have to specify how far away in time the autocorrelation is likely to be significant. Essentially, the autocorrelated errors over the chosen time window are averaged in the computation of the HAC standard errors; you have to specify how many periods over which to average and how much weight to assign each residual in that average. The language of time-series analysis can often be opaque. This is the case here. The weighted average is called a kernel and the number of errors to average in this respect is called bandwidth. Just think of the kernel as another name for weighted average and bandwidth as the term for number of terms to average.

Now, what this has to do with gretl is fairly simple. You get to pick a method of averaging (Bartlett kernel or Parzen kernel) and a bandwidth (nwl, nw2 or some integer). Gretl defaults to the Bartlett kernel and the bandwidth nw 1 = 0.75 x N1/3. As you can see, the bandwidth nw1 is computed based on the sample size, N. The nw2 bandwidth is nw2 = 4 x (N/100)2/9. This one appears to be the default in other programs like EViews.

Implicity there is a trade-off to consider. Larger bandwidths reduce bias (good) as well as precision (bad). Smaller bandwidths exclude more relevant autocorrelations (and hence have more bias), but use more observations to compute the overall covariance and hence increase precision (smaller variance). The general principle is to choose a bandwidth that is large enough to contain the largest autocorrelations. The choice will ultimately depend on the frequency of observation and the length of time it takes for your system to adjust to shocks.

The bandwidth or kernel can be changed using the set command from the console or in a script. The set command is used to change various defaults in gretl and the relevant switches for our use are hac_lag and hac_kernel. The use of these is demonstrated below. The following script changes the kernel to bartlett and the bandwidth to nw2. Then the differences of the unemployment rate are generated. The Phillips curve is estimated by OLS using the ordinary covariance estimator and then by the HAC estimator. The results are collected in a model table.

1 open "@gretldirdatapoephillips_aus. gdt"

2 set hac_kernel bartlett

3 set hac_lag nw2

4 diff u

5 ols inf const d_u

6 modeltab add

7 ols inf const d_u —robust

8 modeltab add

9 modeltab show

The results from the model table are

OLS estimates
Dependent variable: inf






















Standard errors in parentheses
* indicates significance at the 10 percent level
** indicates significance at the 5 percent level
HAC: bandwidth 3 - Bartlett kernel

You can see that the HAC standard errors are quite a bit larger than the usual (and inconsistent) ones. Once gretl recognizes that your data are time-series, then the —robust option will auto­matically apply the HAC estimator of standard errors with the default values of the kernel and bandwidth (or the ones you have set with the set command).

Notice that the standard errors computed using HAC are a little different from those in Hill et al. (2011). No worries, though. They are statistically valid and suggest that EViews and gretl are doing the computations a bit differently.

Добавить комментарий

Using gret l for Principles of Econometrics, 4th Edition


In appendix 10F of POE4, the authors conduct a Monte Carlo experiment comparing the performance of OLS and TSLS. The basic simulation is based on the model y = x …

Hausman Test

The Hausman test probes the consistency of the random effects estimator. The null hypothesis is that these estimates are consistent-that is, that the requirement of orthogonality of the model’s errors …

Time-Varying Volatility and ARCH Models: Introduction to Financial Econometrics

In this chapter we’ll estimate several models in which the variance of the dependent variable changes over time. These are broadly referred to as ARCH (autoregressive conditional heteroskedas - ticity) …

Как с нами связаться:

тел./факс +38 05235  77193 Бухгалтерия
+38 050 512 11 94 — гл. инженер-менеджер (продажи всего оборудования)

+38 050 457 13 30 — Рашид - продажи новинок
Схема проезда к производственному офису:
Схема проезда к МСД

Партнеры МСД

Контакты для заказов шлакоблочного оборудования:

+38 096 992 9559 Инна (вайбер, вацап, телеграм)
Эл. почта: