σμ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 Matrix Decompositions›Módulo 4 · Estimar a covariância›Tópico 4.3·~55 minAvançado

O encolhimento de Ledoit-Wolf

Posição no roteiroMódulo 4 de 5 · Tópico 3 de 3
4.3 O encolhimento de Ledoit-Wolf
—Nota do lab— / 100XP total—Sequência—

Puxar uma estimativa em direção a algo simples

O corte decide autovalor por autovalor se é sinal. O encolhimento toma um caminho mais suave: misturar a ruidosa matriz amostral SSS com um alvo simples e bem condicionado FFF, Σ∗=(1−δ)S+δF\Sigma^{*}=(1-\delta)S+\delta FΣ∗=(1−δ)S+δF, em que a intensidade δ∈[0,1]\delta\in[0,1]δ∈[0,1] diz quanto confiar no alvo. O alvo clássico é uma identidade escalada F=μIF=\mu IF=μI com μ=tr⁡(S)/N\mu=\operatorname{tr}(S)/Nμ=tr(S)/N, a variância média: diz que todos os ativos têm a mesma variância e são não correlacionados. Certamente é falso, mas é estável, e um pouco dele cancela muito ruído.

Σ∗=(1−δ)S+δμI,μ=tr⁡SN,δ=min⁡ ⁣(1,b2d2),d2=∥S−μI∥F2N\Sigma^{*}=(1-\delta)S+\delta\mu I,\qquad \mu=\frac{\operatorname{tr}S}{N},\qquad \delta=\min\!\Big(1,\frac{b^{2}}{d^{2}}\Big),\qquad d^{2}=\frac{\lVert S-\mu I\rVert_F^{2}}{N}Σ∗=(1−δ)S+δμI,μ=NtrS​,δ=min(1,d2b2​),d2=N∥S−μI∥F2​​(1)

A intensidade não é um parâmetro livre. Ledoit e Wolf mostraram que o melhor δ\deltaδ para esse alvo tem forma fechada: o numerador b2b^2b2 mede o quão ruidosa é SSS (a distância quadrática média das matrizes de um único dia xtxtT\mathbf x_t\mathbf x_t^{T}xt​xtT​ até SSS, dividida por TTT e limitada por d2d^2d2) e o denominador d2d^2d2 mede o quão longe SSS está do alvo. Uma SSS ruidosa próxima do alvo recebe um δ\deltaδ grande. O encolhimento tem três propriedades limpas: preserva o traço, aproxima cada autovalor da média (λ↦(1−δ)λ+δμ\lambda\mapsto(1-\delta)\lambda+\delta\muλ↦(1−δ)λ+δμ), de modo que os grandes diminuem e os pequenos aumentam, e torna Σ∗\Sigma^{*}Σ∗ definida positiva para qualquer δ>0\delta>0δ>0, mesmo que SSS seja singular. À medida que TTT cresce, o termo de ruído cai e δ→0\delta\to0δ→0.

Um exemplo 2 por 2Tome S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​), então μ=(4+2)/2=3\mu=(4+2)/2=3μ=(4+2)/2=3 e S−μI=(111−1)S-\mu I=\begin{pmatrix}1&1\\1&-1\end{pmatrix}S−μI=(11​1−1​) com norma de Frobenius ao quadrado 4, logo d2=4/2=2d^2=4/2=2d2=4/2=2. Se o termo de ruído é b2=0,5b^2=0{,}5b2=0,5, a intensidade é δ=0,5/2=0,25\delta=0{,}5/2=0{,}25δ=0,5/2=0,25. Com δ=0,2\delta=0{,}2δ=0,2 em vez disso, Σ∗=0,8S+0,2⋅3I=(3,80,80,82,2)\Sigma^{*}=0{,}8S+0{,}2\cdot3I=\begin{pmatrix}3{,}8&0{,}8\\0{,}8&2{,}2\end{pmatrix}Σ∗=0,8S+0,2⋅3I=(3,80,8​0,82,2​). Os autovalores de SSS são 3±2=4,4143\pm\sqrt2=4{,}4143±2​=4,414 e 1,5861{,}5861,586; os de Σ∗\Sigma^{*}Σ∗ são 0,8(4,414)+0,6=4,1310{,}8(4{,}414)+0{,}6=4{,}1310,8(4,414)+0,6=4,131 e 0,8(1,586)+0,6=1,8690{,}8(1{,}586)+0{,}6=1{,}8690,8(1,586)+0,6=1,869. O número de condição melhora de 2,782{,}782,78 para 2,212{,}212,21.
  • Mistura: Σ∗=(1−δ)S+δF\Sigma^{*}=(1-\delta)S+\delta FΣ∗=(1−δ)S+δF com o alvo F=μIF=\mu IF=μI, μ=tr⁡S/N\mu=\operatorname{tr}S/Nμ=trS/N.
  • Intensidade: δ=min⁡(1,b2/d2)\delta=\min(1,b^2/d^2)δ=min(1,b2/d2); matrizes ruidosas próximas do alvo são mais encolhidas.
  • Efeitos: traço preservado, autovalores aproximados da média, definida positiva para todo δ>0\delta>0δ>0; δ→0\delta\to0δ→0 quando TTT cresce.
Um viés de propósitoO encolhimento adiciona viés deliberadamente, porque o alvo é errado, em troca de uma grande redução de variância. O erro quadrático esperado da estimativa cai e, na estimação de covariância com muitos ativos, a troca costuma ser muito favorável. É a mesma barganha viés-variância da regressão ridge, e é o motivo de o Ledoit-Wolf ser o estimador de covariância padrão em muitas bibliotecas.

Lab: o que é avaliado

  • O nível do alvo10 pts
  • Uma matriz encolhida20 pts
  • A intensidade15 pts
  • Um autovalor encolhido15 pts
  • O novo número de condição15 pts
  • Mais dados, menos encolhimento15 pts
  • A troca viés-variância10 pts
0/ 7 respondidas
Tarefas
1Resposta numérica10 ptsSem resposta

Para S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​) o alvo é μI\mu IμI com μ=tr⁡S/N\mu=\operatorname{tr}S/Nμ=trS/N. Quanto vale μ\muμ?

mu=mu =mu=aceito: ± 0.01
2Matriz · vetor20 ptsSem resposta

Com δ=0,2\delta=0{,}2δ=0,2 e o alvo 3I3I3I, calcule Σ∗=0,8S+0,2⋅3I\Sigma^{*}=0{,}8S+0{,}2\cdot3IΣ∗=0,8S+0,2⋅3I para S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​).

Σ∗=\Sigma^{*} =Σ∗=
3Resposta numérica15 ptsSem resposta

Com um termo de ruído b2=0,5b^2=0{,}5b2=0,5 e uma distância d2=2d^2=2d2=2, qual é a intensidade de Ledoit-Wolf δ=min⁡(1,b2/d2)\delta=\min(1,b^2/d^2)δ=min(1,b2/d2)?

delta=delta =delta=aceito: ± 0.005
4Resposta numérica15 ptsSem resposta

O maior autovalor de SSS é 3+2≈4,4143+\sqrt2\approx4{,}4143+2​≈4,414. Com δ=0,2\delta=0{,}2δ=0,2 e μ=3\mu=3μ=3, qual é o maior autovalor de Σ∗\Sigma^{*}Σ∗, por (1−δ)λ+δμ(1-\delta)\lambda+\delta\mu(1−δ)λ+δμ? Dê com três casas decimais.

lambdamax=lambda_max =lambdam​ax=aceito: ± 0.005
5Resposta numérica15 ptsSem resposta

A matriz encolhida tem autovalores 4,131 e 1,869. Qual é o seu número de condição (maior sobre menor)? Dê com duas casas decimais.

kappa=kappa =kappa=aceito: ± 0.01
6Escolha15 ptsSem resposta

Com N fixo, o que acontece com a intensidade ótima de Ledoit-Wolf quando T, o número de observações, cresce?

7Explicar com palavras10 ptsSem resposta

Em uma ou duas frases, diga por que encolher em direção a um alvo errado ainda pode melhorar uma estimativa de covariância.

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