σμ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 Particle filtering and non-stationary time series›Módulo 3 · Um mercado que não para de mudar›Tópico 3.1·~55 minAvançado

Detecção de deriva e pontos de mudança

Posição no roteiroMódulo 3 de 3 · Tópico 1 de 3
3.1 Detecção de deriva e pontos de mudança
—Nota do lab— / 100XP total—Sequência—

Quando o mercado muda de ideia

Uma série é estacionária se suas estatísticas, como a média e a dispersão, não mudam com o tempo. Os mercados não são. Uma estratégia pode ter média estável por anos e depois perdê-la, a volatilidade pode ficar perto de 7%7\%7% e depois saltar, e uma correlação pode quebrar numa crise. Essas mudanças têm três formas: um salto repentino (mudança de nível), um deslizamento lento (deriva) e uma troca entre regimes. O retorno do SPY tem volatilidade anual de 6,8%6{,}8\%6,8% em 2017, 41,1%41{,}1\%41,1% em 2008 e 33,7%33{,}7\%33,7% em 2020. Uma única estimativa com todo o histórico não descreve bem nenhum desses anos, então é preciso notar as mudanças e depois se adaptar.

xt=rt2σ02,St=max⁡(0, St−1+xt−1−k),St>hx_t=\frac{r_t^{2}}{\sigma_0^{2}},\qquad S_t=\max\left(0,\ S_{t-1}+x_t-1-k\right),\qquad S_t>hxt​=σ02​rt2​​,St​=max(0, St−1​+xt​−1−k),St​>h(1)

O detector CUSUM (soma acumulada) é a ferramenta clássica. Fixa-se uma volatilidade de base σ0\sigma_0σ0​ de um período calmo e mede-se cada dia por xt=rt2/σ02x_t=r_t^{2}/\sigma_0^{2}xt​=rt2​/σ02​, que tem média 111 enquanto nada mudou. Soma-se xt−1−kx_t-1-kxt​−1−k a um total corrente SSS a cada dia, em que kkk é uma folga que ignora pequenos excessos, e nunca se deixa o total cair abaixo de zero. Quando o total passa de um limiar hhh, dispara-se o alarme. Voltar a zero é a chave: um longo período calmo não pode acumular crédito que depois esconda uma mudança real. Um kkk ou hhh pequeno reage mais rápido mas dá mais alarmes falsos.

Uma passada à mão e o SPY realTomemos x=0,5, 3, 4, 2, 0,2x=0{,}5,\ 3,\ 4,\ 2,\ 0{,}2x=0,5, 3, 4, 2, 0,2 com k=0,5k=0{,}5k=0,5 e h=4h=4h=4. O total vai S=0, 1,5, 4,0, 4,5, 3,2S=0,\ 1{,}5,\ 4{,}0,\ 4{,}5,\ 3{,}2S=0, 1,5, 4,0, 4,5, 3,2. Passa de 444 pela primeira vez no dia 444. Um dia de −2%-2\%−2%, com σ0=0,431%\sigma_0=0{,}431\%σ0​=0,431% diário, dá x=(2/0,431)2=21,5x=(2/0{,}431)^{2}=21{,}5x=(2/0,431)2=21,5, um passo enorme. No SPY, com a base tirada de 2017 (6,85%6{,}85\%6,85% ao ano), k=0,5k=0{,}5k=0,5 e h=5h=5h=5, o primeiro alarme é em 26 de janeiro de 2018, dias antes da queda de fevereiro de 2018, e em 16 de março de 2020 o total passa de 300030003000. A razão de variâncias entre 2008 e 2017 é (41,1/6,75)2=37(41{,}1/6{,}75)^{2}=37(41,1/6,75)2=37 vezes.

Detectar é só metade do trabalho. Um filtro de partículas também pode levantar uma bandeira sozinho: quando o retorno de hoje é muito improvável sob a nuvem atual, a verossimilhança média é minúscula e o tamanho amostral efetivo desaba, como em 10 de outubro de 2018. Após um alarme pode-se alargar o ruído do filtro para que siga mais rápido o novo nível, reiniciá-lo de uma crença mais plana ou passar a um estimador de memória mais curta, que é o próximo tópico.

  • Três formas de mudança: um salto, uma deriva e uma troca de regime.
  • CUSUM: um total corrente de excesso de evidência que volta a zero e dispara um alarme acima de um limiar.
  • Velocidade contra alarmes falsos: um limiar mais baixo acha as mudanças mais cedo mas dá mais alarmes falsos.

Lab: o que é avaliado

  • O total no dia dois10 pts
  • O total no dia quatro10 pts
  • O primeiro dia de alarme10 pts
  • Um retorno ao quadrado padronizado15 pts
  • Uma razão de variâncias10 pts
  • Baixar o limiar10 pts
  • Um laço de monitoramento10 pts
  • Associar o vocabulário10 pts
  • Por que voltar a zero15 pts
0/ 9 respondidas
Tarefas
1Resposta numérica10 ptsSem resposta

Parta de S=0S=0S=0. Os valores são x=0,5, 3, 4, 2, 0,2x=0{,}5,\ 3,\ 4,\ 2,\ 0{,}2x=0,5, 3, 4, 2, 0,2 com k=0,5k=0{,}5k=0,5. Use St=max⁡(0, St−1+xt−1−k)S_t=\max(0,\ S_{t-1}+x_t-1-k)St​=max(0, St−1​+xt​−1−k). Quanto vale SSS no dia dois? Dê uma casa decimal.

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

Continue o mesmo cálculo. No dia três S=4,0S=4{,}0S=4,0 e x4=2x_4=2x4​=2. Quanto vale SSS no dia quatro? Dê uma casa decimal.

S=S =S=aceito: ± 0.05
3Resposta numérica10 ptsSem resposta

Os totais são S=0, 1,5, 4,0, 4,5, 3,2S=0,\ 1{,}5,\ 4{,}0,\ 4{,}5,\ 3{,}2S=0, 1,5, 4,0, 4,5, 3,2 e o alarme dispara quando S>4S>4S>4. Em que dia o alarme dispara pela primeira vez?

day=day =day=aceito: ± 0.01
4Resposta numérica15 ptsSem resposta

A volatilidade diária de base é σ0=0,431%\sigma_0=0{,}431\%σ0​=0,431%. Um dia com retorno −2%-2\%−2% dá x=r2/σ02x=r^{2}/\sigma_0^{2}x=r2/σ02​. Quanto vale xxx? Dê uma casa decimal.

x=x =x=aceito: ± 0.3
5Resposta numérica10 ptsSem resposta

A volatilidade anual do SPY foi 41,1%41{,}1\%41,1% em 2008 e 6,75%6{,}75\%6,75% em 2017. Por qual fator a variância de 2008 é maior? Dê uma casa decimal.

r=r =r=aceito: ± 0.5
6Escolha10 ptsSem resposta

Você baixa o limiar hhh do CUSUM. Qual é o efeito habitual?

7Ordenar os passos10 ptsSem resposta

Coloque em ordem os passos diários de um monitor CUSUM.

Arraste para reordenar, ou foque um item, Espaço para levantá-lo, setas para movê-lo e Espaço para soltá-lo.

  1. 1Nunca deixar o total cair abaixo de zero
  2. 2Padronizar o retorno ao quadrado de hoje pela variância de base
  3. 3Comparar o total com o limiar e disparar o alarme se o passar
  4. 4Somar ao total o excesso sobre um mais a folga k
8Associar10 ptsSem resposta

Associe cada termo ao seu significado.

Escolha o par de cada item da esquerda.

Mudança de nível
Deriva
Troca de regime
CUSUM
9Explicar com palavras15 ptsSem resposta

Em uma ou duas frases, explique por que nunca se deixa o total CUSUM cair abaixo de zero.

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