長期分散推定量のクラスを検討する
kは、カーネルまたは重み付け関数であり、 γの jはサンプル自己共分散です。kは、とりわけ対称的で、k(0)=1でなければなりません。ℓTは、帯域幅パラメータです。
JT^≡γ^0+2∑j=1T−1k(jℓT)γ^j
kγ^jkk(0)=1ℓT
Newey&West(Econometrica 1987)は、バートレットカーネル
k (jℓT) = { ( 1 − jℓT)0にとって0 ⩽ J ⩽ ℓT− 1にとってj > ℓT− 1
ハンセン&Hodrickの(政治経済学1980年のジャーナル)推定量切り捨てカーネル、すなわち、服用するのためのJ ≤ MいくつかのためのM、およびK = 0それ以外を。Newey&Westで説明されているように、この推定量は一貫していますが、(行列を推定する場合)正の半正定値であるとは保証されていませんが、Newey&Westのカーネル推定量はそうです。k = 1J ≤ MMk = 0
強く負の係数θを持つMA(1)プロセスに対してを試してください。人口量があることが知られているJ = σ 2(1 + θ )2 > 0が、ハンセン-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
また、ハンセン-ホドリックの推定値は次のように取得できます。
kernHAC(m, kernel = "Truncated", bw = 1,
prewhite = FALSE, adjust = FALSE, sandwich = FALSE)
## (Intercept)
## (Intercept) -0.4056092
線形モデルのNewey-West推定量と時系列の長期分散をそれぞれ取得するための便利なインターフェイスについてもNeweyWest()
およびlrvar()
from sandwich
を参照してください。
Andrews(Econometrica 1991)は、より一般的な条件下での分析を提供しています。
重複データに関するサブ質問については、主題の理由を知りません。私は、伝統がこの一般的な慣行のルーツにあると思います。