考虑一类长期方差估计量
JT^≡γ^0+2∑j=1T−1k(jℓT)γ^j
k是核函数或加权函数,γ^j是样本自协方差。k, 其中必须是对称的并且有k(0)=1.ℓT是带宽参数。
Newey & West (Econometrica 1987)提出了 Bartlett 核
k(jℓT)={(1−jℓT)0for0⩽j⩽ℓT−1forj>ℓT−1
Hansen & Hodrick (Journal of Political Economy 1980) 的估计相当于取一个截断的核,即对于一些,对于 ,否则正如 Newey & West 所讨论的,这个估计量是一致的,但不能保证是半正定的(在估计矩阵时),而 Newey & West 的核估计量是。k=1j≤MMk=0
对具有强负系数的 MA(1) 过程尝试。人口数量已知为,但 Hansen-Hodrick 估计量可能不是: M=1θJ=σ2(1+θ)2>0
set.seed(2)
y <- arima.sim(model = list(ma = -0.95), n = 10)
acf.MA1 <- acf(y, type = "covariance", plot = FALSE)$acf
acf.MA1[1] + 2 * acf.MA1[2]
## [1] -0.4056092
对于长期方差,这不是一个令人信服的估计。
使用 Newey-West 估计器可以避免这种情况:
acf.MA1[1] + acf.MA1[2]
## [1] 0.8634806
使用sandwich
包,这也可以计算为:
library("sandwich")
m <- lm(y ~ 1)
kernHAC(m, kernel = "Bartlett", bw = 2,
prewhite = FALSE, adjust = FALSE, sandwich = FALSE)
## (Intercept)
## (Intercept) 0.8634806
Hansen-Hodrick 估计可以得到:
kernHAC(m, kernel = "Truncated", bw = 1,
prewhite = FALSE, adjust = FALSE, sandwich = FALSE)
## (Intercept)
## (Intercept) -0.4056092
另请参见NeweyWest()
和lrvar()
fromsandwich
以获得方便的接口,以分别获得线性模型的 Newey-West 估计量和时间序列的长期方差。
Andrews (Econometrica 1991)提供了更一般条件下的分析。
至于您关于重叠数据的子问题,我不知道主题原因。我怀疑传统是这种普遍做法的根源。