σμ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 AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • 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 AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Ruta de Matrix Decompositions›Módulo 4 · Estimar la covarianza›Tema 4.3·~55 minAvanzado

La contracción de Ledoit-Wolf

Posición en la rutaMódulo 4 de 5 · Tema 3 de 3
4.3 La contracción de Ledoit-Wolf
—Puntuación del lab— / 100XP total—Racha—

Acercar una estimación a algo sencillo

El recorte decide valor propio por valor propio si es señal. La contracción toma un camino más suave: mezclar la ruidosa matriz muestral SSS con un objetivo simple y bien condicionado FFF, Σ∗=(1−δ)S+δF\Sigma^{*}=(1-\delta)S+\delta FΣ∗=(1−δ)S+δF, donde la intensidad δ∈[0,1]\delta\in[0,1]δ∈[0,1] dice cuánto confiar en el objetivo. El objetivo clásico es una identidad escalada F=μIF=\mu IF=μI con μ=tr⁡(S)/N\mu=\operatorname{tr}(S)/Nμ=tr(S)/N, la varianza media: dice que todos los activos tienen la misma varianza y no están correlacionados. Seguro que es falso, pero es estable, y un poco de él cancela mucho ruido.

Σ∗=(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)

La intensidad no es un parámetro libre. Ledoit y Wolf mostraron que el mejor δ\deltaδ para este objetivo tiene forma cerrada: el numerador b2b^2b2 mide cuán ruidosa es SSS (la distancia cuadrática media de las matrices de un solo día xtxtT\mathbf x_t\mathbf x_t^{T}xt​xtT​ a SSS, dividida por TTT y acotada por d2d^2d2) y el denominador d2d^2d2 mide cuán lejos está SSS del objetivo. Una SSS ruidosa cercana a su objetivo recibe un δ\deltaδ grande. La contracción tiene tres propiedades limpias: conserva la traza, acerca cada valor propio a la media (λ↦(1−δ)λ+δμ\lambda\mapsto(1-\delta)\lambda+\delta\muλ↦(1−δ)λ+δμ), de modo que los grandes bajan y los pequeños suben, y hace Σ∗\Sigma^{*}Σ∗ definida positiva para cualquier δ>0\delta>0δ>0, incluso si SSS es singular. A medida que TTT crece, el término de ruido baja y δ→0\delta\to0δ→0.

Un ejemplo de 2 por 2Tomemos S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​), con lo que μ=(4+2)/2=3\mu=(4+2)/2=3μ=(4+2)/2=3 y S−μI=(111−1)S-\mu I=\begin{pmatrix}1&1\\1&-1\end{pmatrix}S−μI=(11​1−1​) con norma de Frobenius al cuadrado 4, de ahí d2=4/2=2d^2=4/2=2d2=4/2=2. Si el término de ruido es b2=0,5b^2=0{,}5b2=0,5, la intensidad es δ=0,5/2=0,25\delta=0{,}5/2=0{,}25δ=0,5/2=0,25. Con δ=0,2\delta=0{,}2δ=0,2 en su lugar, Σ∗=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​). Los valores propios de SSS son 3±2=4,4143\pm\sqrt2=4{,}4143±2​=4,414 y 1,5861{,}5861,586; los de Σ∗\Sigma^{*}Σ∗ son 0,8(4,414)+0,6=4,1310{,}8(4{,}414)+0{,}6=4{,}1310,8(4,414)+0,6=4,131 y 0,8(1,586)+0,6=1,8690{,}8(1{,}586)+0{,}6=1{,}8690,8(1,586)+0,6=1,869. El número de condición mejora de 2,782{,}782,78 a 2,212{,}212,21.
  • Mezcla: Σ∗=(1−δ)S+δF\Sigma^{*}=(1-\delta)S+\delta FΣ∗=(1−δ)S+δF con el objetivo F=μIF=\mu IF=μI, μ=tr⁡S/N\mu=\operatorname{tr}S/Nμ=trS/N.
  • Intensidad: δ=min⁡(1,b2/d2)\delta=\min(1,b^2/d^2)δ=min(1,b2/d2); las matrices ruidosas cercanas al objetivo se contraen más.
  • Efectos: traza conservada, valores propios acercados a la media, definida positiva para todo δ>0\delta>0δ>0; δ→0\delta\to0δ→0 al crecer TTT.
Un sesgo a propósitoLa contracción añade sesgo deliberadamente, porque el objetivo es erróneo, a cambio de una gran reducción de varianza. El error cuadrático esperado de la estimación baja, y en la estimación de covarianza con muchos activos el trato suele ser muy favorable. Es el mismo trato sesgo-varianza que la regresión ridge y es la razón por la que Ledoit-Wolf es el estimador de covarianza por defecto en muchas bibliotecas.

Lab: qué se califica

  • El nivel del objetivo10 pts
  • Una matriz contraída20 pts
  • La intensidad15 pts
  • Un valor propio contraído15 pts
  • El nuevo número de condición15 pts
  • Más datos, menos contracción15 pts
  • El equilibrio sesgo-varianza10 pts
0/ 7 respondidos
Ejercicios
1Respuesta numérica10 ptsSin responder

Para S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​) el objetivo es μI\mu IμI con μ=tr⁡S/N\mu=\operatorname{tr}S/Nμ=trS/N. ¿Cuánto vale μ\muμ?

mu=mu =mu=aceptado: ± 0.01
2Matriz · vector20 ptsSin responder

Con δ=0,2\delta=0{,}2δ=0,2 y el objetivo 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^{*} =Σ∗=
3Respuesta numérica15 ptsSin responder

Con un término de ruido b2=0,5b^2=0{,}5b2=0,5 y una distancia d2=2d^2=2d2=2, ¿cuál es la intensidad de Ledoit-Wolf δ=min⁡(1,b2/d2)\delta=\min(1,b^2/d^2)δ=min(1,b2/d2)?

delta=delta =delta=aceptado: ± 0.005
4Respuesta numérica15 ptsSin responder

El mayor valor propio de SSS es 3+2≈4,4143+\sqrt2\approx4{,}4143+2​≈4,414. Con δ=0,2\delta=0{,}2δ=0,2 y μ=3\mu=3μ=3, ¿cuál es el mayor valor propio de Σ∗\Sigma^{*}Σ∗, mediante (1−δ)λ+δμ(1-\delta)\lambda+\delta\mu(1−δ)λ+δμ? Dé tres decimales.

lambdamax=lambda_max =lambdam​ax=aceptado: ± 0.005
5Respuesta numérica15 ptsSin responder

La matriz contraída tiene valores propios 4,131 y 1,869. ¿Cuál es su número de condición (mayor entre menor)? Dé dos decimales.

kappa=kappa =kappa=aceptado: ± 0.01
6Elección15 ptsSin responder

Con N fijo, ¿qué ocurre con la intensidad óptima de Ledoit-Wolf cuando T, el número de observaciones, crece?

7Explicar con palabras10 ptsSin responder

En una o dos frases, diga por qué acercar a un objetivo erróneo puede aun así mejorar una estimación de covarianza.

Modo matemáticasEntrada por teclado o LaTeXEjercicio 1 de 7 · Ctrl ⏎ para comprobar

Inicia sesión para ejecutar tu código.

Todas las clases se pueden leer gratis. Crea una cuenta o inicia sesión para ejecutar código, lanzar las pruebas y enviar este lab.

  • IDE en el navegador con terminal bash — sin instalar nada
  • Corrección instantánea, XP y tu lugar en la ruta
Crear una cuentaIniciar sesión
O
Continuar con GitHub