σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyRoadmapsRoteiros
Ctrl K
AO
+1,2 mil alunos
EntrarParticipar ↗
RoadmapsRoteiros
AO
+1,2 mil alunos
Entrar
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyRoadmapsRoteiros
Ctrl K
AO
+1,2 mil alunos
EntrarParticipar ↗
RoadmapsRoteiros
AO
+1,2 mil alunos
Entrar
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Roteiro de statsmodels›Módulo 1 · A regressão e suas verificações›Tópico 1.6·~50 minIntermediário

Erros-padrão robustos de Newey–West

Posição no roteiroMódulo 1 de 3 · Tópico 6 de 8
1.6 Erros-padrão robustos de Newey–West
—Nota do lab— / 100XP total—Sequência—

Quando erros vizinhos estão ligados

Os erros-padrão comuns supõem que cada observação traz informação nova. Dados financeiros muitas vezes quebram isso. Retornos medidos em três meses e amostrados todo mês se sobrepõem: cada observação divide dois meses com a vizinha e os erros se movem juntos. O mesmo acontece quando a volatilidade se agrupa. A correção mantém os coeficientes e muda só os erros-padrão. O estimador de Newey–West, também chamado HAC por consistente com heterocedasticidade e autocorrelação, estima a variância de longo prazo somando as autocovariâncias defasadas γℓ\gamma_\ellγℓ​ dos resíduos, cada uma multiplicada por um peso de Bartlett 1−ℓ/(L+1)1-\ell/(L+1)1−ℓ/(L+1) que cai a zero na defasagem de corte LLL. No statsmodels basta acrescentar cov_type='HAC' e cov_kwds={'maxlags': L} a .fit().

ω^=γ^0+2∑ℓ=1L(1−ℓL+1)γ^ℓ,γ^ℓ=1n∑t=ℓ+1netet−ℓ,se(zˉ)=ω^/n,L≈⌊4(n100)2/9⌋\hat\omega=\hat\gamma_0+2\sum_{\ell=1}^{L}\Bigl(1-\frac{\ell}{L+1}\Bigr)\hat\gamma_\ell,\qquad \hat\gamma_\ell=\frac1n\sum_{t=\ell+1}^{n}e_te_{t-\ell},\qquad \mathrm{se}(\bar z)=\sqrt{\hat\omega/n},\qquad L\approx\Bigl\lfloor4\Bigl(\frac{n}{100}\Bigr)^{2/9}\Bigr\rfloorω^=γ^​0​+2ℓ=1∑L​(1−L+1ℓ​)γ^​ℓ​,γ^​ℓ​=n1​t=ℓ+1∑n​et​et−ℓ​,se(zˉ)=ω^/n​,L≈⌊4(100n​)2/9⌋(1)
Uma série de cinco númerosTomemos z=3,3,2,1,1z=3,3,2,1,1z=3,3,2,1,1 de média 222 e desvios 1,1,0,−1,−11,1,0,-1,-11,1,0,−1,−1. Então γ^0=4/5=0,8\hat\gamma_0=4/5=0{,}8γ^​0​=4/5=0,8 e γ^1=(1+0+0+1)/5=0,4\hat\gamma_1=(1+0+0+1)/5=0{,}4γ^​1​=(1+0+0+1)/5=0,4, positivo porque os vizinhos se parecem. Com L=1L=1L=1 o peso de Bartlett é 1−1/2=0,51-1/2=0{,}51−1/2=0,5, logo ω^=0,8+2⋅0,5⋅0,4=1,2\hat\omega=0{,}8+2\cdot0{,}5\cdot0{,}4=1{,}2ω^=0,8+2⋅0,5⋅0,4=1,2 e se=1,2/5=0,490\mathrm{se}=\sqrt{1{,}2/5}=0{,}490se=1,2/5​=0,490. O erro-padrão comum é 1/5=0,447\sqrt{1/5}=0{,}4471/5​=0,447 (com s2=SSE/(n−1)=1s^2=\mathrm{SSE}/(n-1)=1s2=SSE/(n−1)=1), então a estatística t da média cai de 2/0,447=4,472/0{,}447=4{,}472/0,447=4,47 para 2/0,490=4,082/0{,}490=4{,}082/0,490=4,08. Quando a autocorrelação é negativa, o erro-padrão diminui em vez disso.

Os dados reais mostram o tamanho do efeito. O retorno médio de 3 meses do SPY, amostrado todo mês, é 2,78%2{,}78\%2,78% em 256256256 observações sobrepostas. Os erros-padrão comuns dão 0,456%0{,}456\%0,456% e t=6,10t=6{,}10t=6,10. As autocorrelações da série centrada são 0,680{,}680,68, 0,350{,}350,35 e 0,080{,}080,08 nas defasagens 111, 222 e 333, o triângulo que a sobreposição de dois meses em três cria. Com Newey–West em L=2L=2L=2 o erro-padrão é 0,664%0{,}664\%0,664% (t=4,18t=4{,}18t=4,18), em L=4L=4L=4 é 0,732%0{,}732\%0,732% (t=3,80t=3{,}80t=3,80) e em L=6L=6L=6 é 0,762%0{,}762\%0,762% (t=3,65t=3{,}65t=3,65), e os 868686 trimestres sem sobreposição dão t=3,55t=3{,}55t=3,55. Numa regressão mensal de três fatores para a Apple (319319319 meses) a mudança é menor mas real: o erro-padrão da exposição a valor sobe de 0,1440{,}1440,144 para 0,1700{,}1700,170, então sua estatística t vai de −4,98-4{,}98−4,98 a −4,22-4{,}22−4,22.

  • Mesmos coeficientes, novos erros-padrão: o HAC muda só a incerteza, nunca a estimativa.
  • A sobreposição define a defasagem: com retornos de hhh meses amostrados todo mês, use pelo menos L=h−1L=h-1L=h−1 e muitas vezes alguns a mais.
  • Informe o t honesto: um t de 6 pode virar um t de 4 quando os erros vizinhos são considerados.

Lab: o que é avaliado

  • O termo de variância10 pts
  • A primeira autocovariância10 pts
  • A variância de longo prazo15 pts
  • O erro-padrão HAC15 pts
  • A estatística t robusta10 pts
  • A chamada do statsmodels10 pts
  • Associar o vocabulário10 pts
  • Por que o t comum é grande demais20 pts
0/ 8 respondidas
Tarefas
1Resposta numérica10 ptsSem resposta

Os desvios da média são 1,1,0,−1,−11,1,0,-1,-11,1,0,−1,−1. Quanto vale γ^0=1n∑et2\hat\gamma_0=\frac1n\sum e_t^2γ^​0​=n1​∑et2​? Dê uma casa decimal.

g0=g0 =g0=aceito: ± 0.01
2Resposta numérica10 ptsSem resposta

Com os mesmos desvios, quanto vale γ^1=1n∑t=2netet−1\hat\gamma_1=\frac1n\sum_{t=2}^{n}e_te_{t-1}γ^​1​=n1​∑t=2n​et​et−1​? Dê uma casa decimal.

g1=g1 =g1=aceito: ± 0.01
3Resposta numérica15 ptsSem resposta

Com L=1L=1L=1 o peso de Bartlett é 0,50{,}50,5. Com γ^0=0,8\hat\gamma_0=0{,}8γ^​0​=0,8 e γ^1=0,4\hat\gamma_1=0{,}4γ^​1​=0,4, quanto vale ω^=γ^0+2⋅0,5⋅γ^1\hat\omega=\hat\gamma_0+2\cdot0{,}5\cdot\hat\gamma_1ω^=γ^​0​+2⋅0,5⋅γ^​1​? Dê uma casa decimal.

w=w =w=aceito: ± 0.02
4Resposta numérica15 ptsSem resposta

Com ω^=1,2\hat\omega=1{,}2ω^=1,2 e n=5n=5n=5, quanto vale se(zˉ)=ω^/n\mathrm{se}(\bar z)=\sqrt{\hat\omega/n}se(zˉ)=ω^/n​? Use 0,24=0,4899\sqrt{0{,}24}=0{,}48990,24​=0,4899 e dê três casas decimais.

se=se =se=aceito: ± 0.004
5Resposta numérica10 ptsSem resposta

A média da série é 222 e seu erro-padrão HAC é 0,4900{,}4900,490. Quanto vale a estatística t? Dê duas casas decimais.

t=t =t=aceito: ± 0.03
6Escolha10 ptsSem resposta

Qual chamada ajusta o modelo com erros-padrão de Newey–West com 3 defasagens?

7Associar10 ptsSem resposta

Associe cada termo ao seu significado.

Escolha o par de cada item da esquerda.

HAC
Peso de Bartlett
maxlags
Retornos sobrepostos
8Explicar com palavras20 ptsSem resposta

O retorno médio de 3 meses tem t=6,10t=6{,}10t=6,10 com erros comuns e t=3,80t=3{,}80t=3,80 com Newey–West em L=4L=4L=4. Em uma ou duas frases, explique por que retornos sobrepostos tornam o ttt comum grande demais.

Modo matemáticaEntrada por teclado ou LaTeXTarefa 1 de 8 · Ctrl ⏎ para verificar

Entre para executar o seu código.

Todas as aulas são gratuitas para ler. Crie uma conta ou entre para executar código, rodar os testes e enviar este lab.

  • IDE no navegador com terminal bash — nada para instalar
  • Correção instantânea, XP e o seu lugar no roteiro
Criar uma contaEntrar
OU
Continuar com o GitHub