Hitung kesalahan standar Newey-West tanpa objek lm di R

Saya menanyakan pertanyaan ini kemarin di StackOverflow, dan mendapat jawaban, tetapi kami sepakat bahwa itu agak sedikit meretas dan mungkin ada cara yang lebih baik untuk melihatnya.

Pertanyaannya: Saya ingin menghitung kesalahan standar Newey-West (HAC) untuk vektor (dalam hal ini vektor pengembalian saham). Fungsi NeweyWest()dalam sandwichpaket melakukan ini, tetapi mengambil lmobjek sebagai input. Solusi yang ditawarkan Joris Meys adalah memproyeksikan vektor ke 1, yang mengubah vektor saya menjadi residu untuk dimasukkan NeweyWest(). Itu adalah:

as.numeric(NeweyWest(lm(rnorm(100) ~ 1)))

untuk varian mean.

Haruskah saya melakukannya seperti ini? Atau ada cara untuk lebih langsung melakukan apa yang saya inginkan? Terima kasih!

r standard-error autocorrelation heteroscedasticity Richard Herron
sumber

Pertanyaan tidak jelas. Apa yang Anda maksud dengan "standard error for a vector"? Biasanya kami menginginkan kesalahan standar estimasi parameter. Parameter apa yang Anda perkirakan? Kode yang Anda berikan menghasilkan estimasi Newey West dari kesalahan standar kuadrat dari rata-rata. Itukah yang kamu inginkan?

Cyrus S

@Cyrus - Dengan "vektor" Maksud saya bukan lmobjek. Saya sering memiliki vektor (katakanlah serangkaian pengembalian saham) yang saya tidak ingin terlibat dalam regresi (karena saya tidak peduli tentang proyeksi itu, selain pada 1), tetapi saya masih menginginkan HAC kesalahan standar. Dalam hal ini estimasi parameter adalah return saham. Jawaban di atas melakukan itu, tetapi membutuhkan penghitungan lmobjek, yang saya benar-benar tidak perlu. Jadi saya bertanya-tanya apakah ada rutin di R yang melakukan ini tanpa membuat lmobjek.

Richard Herron

Maaf, masih belum jelas: "Dalam hal ini estimasi parameter adalah pengembalian saham." Maksud Anda, "rata-rata pengembalian saham dalam seri"? Jika ya, maka apa yang Anda miliki baik-baik saja.

Cyrus S

@Cyrus - Saya tahu apa yang saya punya, tapi saya berharap ada cara untuk menghitung UK tanpa melewati lmobjek untuk kasus satu vektor. Saya rasa tidak. Terima kasih telah membantu saya menjelaskan pertanyaan saya!

Richard Herron

Jawaban:

Misalkan kita memiliki regresi

\begin{aligned} y = X β + kamu \end{aligned}

$\begin{align*} y=X\beta+u \end{align*}$

Kemudian OLS memperkirakan adalah dan dengan asumsi bahwa adalah berisi perkiraan kita memiliki $\hat{\beta}$

\begin{aligned} \hat{β} - β = (X^{'} X)^{- 1} X^{'} kamu \end{aligned}

$\begin{align*} \widehat{\beta}-\beta=(X'X)^{-1}X'u \end{align*}$

\hat{β}

$\hat{\beta}$

\begin{aligned} V Sebuah r (\hat{β}) = E [(X^{'} X)^{- 1} X^{'} kamu {kamu}^{'} X (X^{'} X)^{- 1}] \end{aligned}

$\begin{align*} Var(\widehat{\beta})=E\left[(X'X)^{-1}X'uu'X(X'X)^{-1}\right] \end{align*}$

Asumsi OLS yang biasa adalah bahwa dan $E(u|X)=0$ yang memberi kita ini matriks kovarians biasanya dilaporkan dalam paket statistik. $E(uu'|X)=\sigma^2I_n$

\begin{aligned} V Sebuah r (\hat{β}) = σ^{2} E (X^{'} X)^{- 1} \end{aligned}

$\begin{align*} Var(\widehat{\beta})=\sigma^2E(X'X)^{-1} \end{align*}$

Jika adalah heteroscedastic dan (atau) autocorellated, maka dan output yang biasa memberikan hasil yang menyesatkan. Untuk mendapatkan hasil yang benar, kesalahan standar HAC dihitung. Semua metode untuk kesalahan HAC menghitung $u_i$ $E(uu'|X)\neq\sigma^2I_n$

\begin{aligned} d saya Sebuah g (E (X^{'} X)^{- 1} X^{'} kamu {kamu}^{'} X (X^{'} X)^{- 1}) . \end{aligned}

$\begin{align*} diag(E(X'X)^{-1}X'uu'X(X'X)^{-1}). \end{align*}$ Mereka berbeda pada asumsi mereka seperti apa

E (u u^{'} | X)

$E(uu'|X)$

Jadi wajar kalau fungsi itu NeweyWestmeminta model linier. Metode Newey-West menghitung kesalahan standar yang benar dari penduga model linier. Jadi solusi Anda adalah sempurna benar jika Anda menganggap bahwa return saham Anda mengikuti model dan Anda ingin memperkirakan menjaga terhadap penyimpangan dalam .

\begin{aligned} r_{t} = μ + {kamu}_{t} \end{aligned}

$\begin{align} r_t=\mu+u_t \end{align}$

V a r (μ)

$Var(\mu)$

u_{t}

$u_t$

$Var(r_t)$

\begin{aligned} r_{t} = σ_{t} ε_{t} \end{aligned}

$\begin{align*} r_t=\sigma_t\varepsilon_t \end{align*}$

ε_{t}

$\varepsilon_t$

σ_{t}

$\sigma_t$

V a r (r_{t}) = V a r (σ_{t})

$Var(r_t)=Var(\sigma_t)$ dan Anda memiliki perkiraan "benar" dari varians Anda, menjaga terhadap kekhasan pengembalian saham biasa seperti pengelompokan volatilitas, kemiringan dan lain-lain.

mpiktas
sumber

Terima kasih! Mungkin tidak ada cara yang lebih efisien bagi saya untuk kode ini daripada membentuk lmobjek.

Richard Herron

Saya kira lmobjeknya adalah cara untuk pergi! Terima kasih untuk ringkasan yang bagus ... kadang-kadang dalam aplikasi saya terlalu jauh dari teori.

Richard Herron