σμ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.4·~50 minIntermediário

Diagnóstico de resíduos VIF e Durbin–Watson

Posição no roteiroMódulo 1 de 3 · Tópico 4 de 8
1.4 Diagnóstico de resíduos VIF e Durbin–Watson
—Nota do lab— / 100XP total—Sequência—

Verificar os pressupostos antes de confiar nas estrelas

Os erros-padrão de um resumo de regressão só estão certos se três coisas valerem: os resíduos não são correlacionados entre si, sua dispersão não muda no tempo e os regressores não são quase cópias uns dos outros. O statsmodels dá um teste para cada uma. A estatística de Durbin–Watson olha os resíduos vizinhos. O fator de inflação da variância, VIF, olha quão bem cada regressor pode ser reconstruído a partir dos outros. O teste de Breusch–Pagan pergunta se a dispersão dos resíduos depende dos regressores, e o de Jarque–Bera se a assimetria e as caudas pesadas são grandes o bastante para importar. Um mau resultado não torna os coeficientes inúteis, diz quais erros-padrão pôr em dúvida.

DW=∑t=2n(et−et−1)2∑t=1net2≈2(1−ρ^1),VIFj=11−Rj2,sej∝VIFj\mathrm{DW}=\frac{\sum_{t=2}^{n}(e_t-e_{t-1})^{2}}{\sum_{t=1}^{n}e_t^{2}}\approx2(1-\hat\rho_1),\qquad \mathrm{VIF}_j=\frac{1}{1-R_j^{2}},\qquad \mathrm{se}_j\propto\sqrt{\mathrm{VIF}_j}DW=∑t=1n​et2​∑t=2n​(et​−et−1​)2​≈2(1−ρ^​1​),VIFj​=1−Rj2​1​,sej​∝VIFj​​(1)
Resíduos e dois regressores correlacionadosTomemos os resíduos −0,4;0,8;0;−0,8;0,4-0{,}4;0{,}8;0;-0{,}8;0{,}4−0,4;0,8;0;−0,8;0,4 do ajuste de cinco meses. As diferenças entre vizinhos são 1,2;−0,8;−0,8;1,21{,}2;-0{,}8;-0{,}8;1{,}21,2;−0,8;−0,8;1,2, seus quadrados somam 4,164{,}164,16 e ∑et2=1,6\sum e_t^2=1{,}6∑et2​=1,6, então DW=2,6\mathrm{DW}=2{,}6DW=2,6 e ρ^1≈1−2,6/2=−0,3\hat\rho_1\approx1-2{,}6/2=-0{,}3ρ^​1​≈1−2,6/2=−0,3: os resíduos tendem a trocar de sinal. Um valor perto de 222 significa sem correlação de primeira ordem, perto de 000 correlação positiva e perto de 444 negativa. Se dois regressores têm correlação 0,90{,}90,9, cada um é explicado pelo outro com Rj2=0,81R_j^2=0{,}81Rj2​=0,81 e VIF=1/0,19=5,26\mathrm{VIF}=1/0{,}19=5{,}26VIF=1/0,19=5,26. Com correlação 0,990{,}990,99 o VIF é cerca de 505050, e com Rj2=0,99R_j^2=0{,}99Rj2​=0,99 é 100100100, o que estica os erros-padrão 100=10\sqrt{100}=10100​=10 vezes. Em dados reais, a regressão de três fatores da Apple (258258258 meses desde março de 2005) tem DW=1,94\mathrm{DW}=1{,}94DW=1,94 e valor p de Breusch–Pagan de 0,210{,}210,21, então há pouco a temer, embora Jarque–Bera rejeite a normalidade (p=0,015p=0{,}015p=0,015, curtose 3,83{,}83,8). Agora acrescente o retorno excedente do fundo SPY ao lado do fator de mercado: a correlação é 0,9940{,}9940,994, os dois VIF saltam para 82,982{,}982,9, o erro-padrão da exposição ao mercado vai de 0,100{,}100,10 a 0,910{,}910,91 e sua estimativa passa a −0,23-0{,}23−0,23 enquanto o R2R^2R2 quase não se move (0,3470{,}3470,347). Por fim, regrida log HD em log LOW em 6 6966\,6966696 dias: R2=0,93R^2=0{,}93R2=0,93 mas DW=0,002\mathrm{DW}=0{,}002DW=0,002. Um DW muito abaixo do R2R^2R2 avisa que dois preços com tendência estão sendo pareados por acaso.
Um alerta não é um vereditoCada diagnóstico aponta uma hipótese. Corrija o problema correspondente: erros-padrão robustos para resíduos correlacionados ou desiguais (próximo tópico), remover ou combinar regressores colineares, ou diferenciar séries com tendência.
  • Durbin–Watson perto de 2 é saudável, perto de 0 ou 4 indica que resíduos vizinhos estão ligados.
  • VIF acima de uns 10 indica que os regressores se emaranham e as estatísticas t individuais não são confiáveis.
  • Um R2R^2R2 alto com um DW minúsculo é o sinal clássico de uma regressão entre séries com tendência.

Lab: o que é avaliado

  • Uma estatística de Durbin–Watson20 pts
  • A autocorrelação implícita10 pts
  • Um fator de inflação da variância20 pts
  • Quanto um erro-padrão se estica20 pts
  • Associar os diagnósticos10 pts
  • Regressores colineares20 pts
0/ 6 respondidas
Tarefas
1Resposta numérica20 ptsSem resposta

Os resíduos são −0,4;0,8;0;−0,8;0,4-0{,}4;0{,}8;0;-0{,}8;0{,}4−0,4;0,8;0;−0,8;0,4. As diferenças ao quadrado entre vizinhos somam 4,164{,}164,16 e ∑et2=1,6\sum e_t^2=1{,}6∑et2​=1,6. Quanto vale DW\mathrm{DW}DW? Dê uma casa decimal.

dw=dw =dw=aceito: ± 0.05
2Resposta numérica10 ptsSem resposta

Use ρ^1≈1−DW/2\hat\rho_1\approx1-\mathrm{DW}/2ρ^​1​≈1−DW/2 com DW=2,6\mathrm{DW}=2{,}6DW=2,6. Quanto vale ρ^1\hat\rho_1ρ^​1​?

rho=rho =rho=aceito: ± 0.01
3Resposta numérica20 ptsSem resposta

Dois regressores têm correlação 0,90{,}90,9. Regredir um no outro dá Rj2=0,81R_j^2=0{,}81Rj2​=0,81. Quanto vale VIF=1/(1−Rj2)\mathrm{VIF}=1/(1-R_j^2)VIF=1/(1−Rj2​)? Dê duas casas decimais.

v=v =v=aceito: ± 0.03
4Resposta numérica20 ptsSem resposta

Os dois VIF valem 82,982{,}982,9. Por que fator cada erro-padrão é maior do que seria com regressores não correlacionados, VIF\sqrt{\mathrm{VIF}}VIF​? Dê uma casa decimal.

s=s =s=aceito: ± 0.1
5Associar10 ptsSem resposta

Associe cada diagnóstico ao problema que ele procura.

Escolha o par de cada item da esquerda.

Durbin–Watson
VIF
Breusch–Pagan
Jarque–Bera
6Explicar com palavras20 ptsSem resposta

Acrescentar o retorno do fundo SPY ao lado do fator de mercado deixou o R2R^2R2 em 0,3470{,}3470,347 mas multiplicou por nove o erro-padrão da exposição ao mercado. Em uma ou duas frases, explique por quê.

Modo matemáticaEntrada por teclado ou LaTeXTarefa 1 de 6 · 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