σμ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 3 · Cortes transversais e caudas›Tópico 3.1·~45 minAvançado

Regressão quantílica para o risco de cauda

Posição no roteiroMódulo 3 de 3 · Tópico 1 de 3
3.1 Regressão quantílica para o risco de cauda
—Nota do lab— / 100XP total—Sequência—

Modelar as bordas, não o meio

Os mínimos quadrados modelam a média. O risco vive nas caudas: os piores dias de uma ação podem reagir a uma variável que mal mexe no seu dia médio. A regressão quantílica ajusta diretamente um quantil escolhido τ\tauτ da resposta, por exemplo o quantil de 5 % em que se constrói um limite de valor em risco. Ela minimiza a perda de verificação (pinball) ρτ(u)=u (τ−1[u<0])\rho_\tau(u)=u\,(\tau-\mathbf 1[u<0])ρτ​(u)=u(τ−1[u<0]) sobre os resíduos u=y−qu=y-qu=y−q. Resíduos positivos são ponderados por τ\tauτ e os negativos por 1−τ1-\tau1−τ, então para τ=0,9\tau=0{,}9τ=0,9 um palpite baixo demais custa nove vezes mais que um alto demais, o que empurra o palpite para o percentil 90. Para τ=0,5\tau=0{,}5τ=0,5 a perda é metade do erro absoluto e o minimizador é a mediana. No statsmodels escreve-se smf.quantreg('y ~ vol', data).fit(q=0.05).

ρτ(u)=u (τ−1[u<0]),β^τ=arg⁡min⁡β∑iρτ(yi−xi⊤β),Qτ(y∣x)=x⊤βτ\rho_\tau(u)=u\,(\tau-\mathbf 1[u<0]),\qquad \hat\beta_\tau=\arg\min_\beta\sum_i\rho_\tau(y_i-x_i^{\top}\beta),\qquad Q_\tau(y\mid x)=x^{\top}\beta_\tauρτ​(u)=u(τ−1[u<0]),β^​τ​=argβmin​i∑​ρτ​(yi​−xi⊤​β),Qτ​(y∣x)=x⊤βτ​(1)
A perda à mão e uma reta de valor em riscoTomemos y=1,2,3,10,20y=1,2,3,10,20y=1,2,3,10,20. A mediana 333 tem perda absoluta 2+1+0+7+17=272+1+0+7+17=272+1+0+7+17=27, enquanto a média 7,27{,}27,2 tem 31,231{,}231,2. Para τ=0,9\tau=0{,}9τ=0,9 e o palpite q=10q=10q=10 os resíduos são −9,−8,−7,0,10-9,-8,-7,0,10−9,−8,−7,0,10 e a perda é 0,1(9+8+7)+0,9⋅10=11,40{,}1(9+8+7)+0{,}9\cdot10=11{,}40,1(9+8+7)+0,9⋅10=11,4, enquanto o palpite q=20q=20q=20 custa só 0,1(19+18+17+10)=6,40{,}1(19+18+17+10)=6{,}40,1(19+18+17+10)=6,4, porque quatro de cinco pontos ficam abaixo, que é como um quantil de 90 % deve parecer. Agora dados reais: o retorno do SPY do dia seguinte em percentual contra a volatilidade dos últimos 20 dias (também em percentual), 5 3845\,3845384 dias. Os mínimos quadrados quase não acham efeito, com inclinação 0,0110{,}0110,011. A regressão quantílica acha inclinação −1,50-1{,}50−1,50 em τ=0,05\tau=0{,}05τ=0,05 (erro-padrão 0,050{,}050,05) e +1,37+1{,}37+1,37 em τ=0,95\tau=0{,}95τ=0,95, e 0,030{,}030,03 na mediana. Com volatilidade 1%1\%1% os quantis ajustados são −1,77%-1{,}77\%−1,77%, 0,07%0{,}07\%0,07% e 1,64%1{,}64\%1,64%, e com 2%2\%2% são −3,27%-3{,}27\%−3,27%, 0,10%0{,}10\%0,10% e 3,01%3{,}01\%3,01%: uma faixa de largura 6,286{,}286,28 pontos. As caudas se abrem com a volatilidade enquanto o centro mal se move. Como limite de valor em risco, a reta ajustada de τ=0,05\tau=0{,}05τ=0,05 é ultrapassada em 5,0%5{,}0\%5,0% dos dias no total, e em 7,2%7{,}2\%7,2% dos 263263263 dias com volatilidade acima de 2%2\%2%, onde o limite normal constante de −1,92%-1{,}92\%−1,92% é ultrapassado em 22,4%22{,}4\%22,4% e em apenas 1,3%1{,}3\%1,3% dos 2 6002\,6002600 dias calmos abaixo de 0,8%0{,}8\%0,8%.

Duas cautelas. Retas quantílicas ajustadas separadamente podem se cruzar, de modo que a de 5 % termina acima da de 50 % nos extremos de xxx, o que um analista cuidadoso verifica. E os erros-padrão da regressão quantílica são por bootstrap ou assintóticos, com as mesmas preocupações de autocorrelação de antes, então estimativas de cauda em amostras curtas são ruidosas.

  • Escolha o quantil que importa: 5 % para o risco de queda, 50 % para o caso típico, 95 % para a alta.
  • Verifique a cobertura por regime: um limite certo no total pode estar errado em cada período calmo e em cada período agitado.

Lab: o que é avaliado

  • Uma perda absoluta15 pts
  • Uma perda pinball15 pts
  • Um quantil de 5 % previsto15 pts
  • A largura da faixa10 pts
  • Violações esperadas15 pts
  • Por que o limite constante falha10 pts
  • A média esconde a cauda20 pts
0/ 7 respondidas
Tarefas
1Resposta numérica15 ptsSem resposta

Para y=1,2,3,10,20y=1,2,3,10,20y=1,2,3,10,20 e o palpite q=3q=3q=3 (a mediana), quanto vale ∑∣yi−q∣\sum|y_i-q|∑∣yi​−q∣?

l1=l1 =l1=aceito: ± 0.1
2Resposta numérica15 ptsSem resposta

Para τ=0,9\tau=0{,}9τ=0,9 e o palpite q=10q=10q=10 os resíduos y−qy-qy−q são −9,−8,−7,0,10-9,-8,-7,0,10−9,−8,−7,0,10. Resíduos negativos custam 0,10{,}10,1 por unidade e positivos 0,90{,}90,9 por unidade. Qual é a perda total? Dê uma casa decimal.

pin=pin =pin=aceito: ± 0.05
3Resposta numérica15 ptsSem resposta

A reta de 5 % ajustada é −0,277−1,495⋅v-0{,}277-1{,}495\cdot v−0,277−1,495⋅v, com a volatilidade de 20 dias em percentual. Quanto vale o quantil de 5 % do retorno de amanhã quando a volatilidade é 1,51{,}51,5? Dê duas casas decimais.

v=v =v=aceito: ± 0.02
4Resposta numérica10 ptsSem resposta

Com volatilidade 2%2\%2% os quantis ajustados de 5 % e 95 % são −3,27-3{,}27−3,27 e 3,013{,}013,01. Qual é a largura da faixa entre eles? Dê duas casas decimais.

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

Um limite de 5 % correto deveria ser ultrapassado em 5%5\%5% dos dias. Entre os 263263263 dias com volatilidade acima de 2%2\%2%, quantas violações se esperam? Dê uma casa decimal.

n=n =n=aceito: ± 0.1
6Escolha10 ptsSem resposta

Um limite normal constante de -1,92 % é ultrapassado em 22,4 % dos dias de alta volatilidade e 1,3 % dos calmos. Por quê?

7Explicar com palavras20 ptsSem resposta

Os mínimos quadrados dão uma inclinação de 0,0110{,}0110,011 para o retorno do dia seguinte sobre a volatilidade, enquanto a inclinação do quantil de 5 % é −1,50-1{,}50−1,50. Em uma ou duas frases, explique como as duas podem ser verdadeiras.

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