σμ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.2·~1 hAvanzado

Recorte de valores propios y Marchenko-Pastur

Posición en la rutaMódulo 4 de 5 · Tema 2 de 3
4.2 Recorte de valores propios y Marchenko-Pastur
—Puntuación del lab— / 100XP total—Racha—

Distinguir los valores propios de señal del ruido

Si los rendimientos fueran ruido puro sin estructura, los valores propios de su matriz de correlación igualmente se repartirían en un intervalo conocido, el intervalo de Marchenko-Pastur [λ−,λ+][\lambda_-,\lambda_+][λ−​,λ+​] con λ±=(1±q)2\lambda_\pm=(1\pm\sqrt q)^2λ±​=(1±q​)2 y q=N/Tq=N/Tq=N/T. Ese intervalo es una regla. Un valor propio de la matriz de correlación muestral por encima de λ+\lambda_+λ+​ es mayor de lo que el ruido puede explicar, así que refleja un factor común real. Los valores propios dentro del intervalo son indistinguibles del ruido: sus vectores propios son direcciones aleatorias que no informan sobre el futuro.

λi′={λi,λi>λ+λˉnoise,λi≤λ+,λˉnoise=1k∑λi≤λ+λi,Cclean=QΛ′QT\lambda_i'=\begin{cases}\lambda_i,&\lambda_i>\lambda_+\\[2pt]\bar\lambda_{\mathrm{noise}},&\lambda_i\le\lambda_+\end{cases},\qquad \bar\lambda_{\mathrm{noise}}=\frac{1}{k}\sum_{\lambda_i\le\lambda_+}\lambda_i,\qquad C_{\mathrm{clean}}=Q\Lambda'Q^{T}λi′​={λi​,λˉnoise​,​λi​>λ+​λi​≤λ+​​,λˉnoise​=k1​λi​≤λ+​∑​λi​,Cclean​=QΛ′QT(1)

El recorte de valores propios actúa sobre este veredicto. Se conservan los valores propios por encima de λ+\lambda_+λ+​ y se sustituyen todos los demás por su propia media λˉnoise\bar\lambda_{\mathrm{noise}}λˉnoise​, lo que aplana la parte ruidosa del espectro. Promediar en lugar de poner cero conserva la traza, que en una matriz de correlación es NNN, y mantiene la matriz definida positiva y por tanto invertible, lo que necesita un optimizador. Después se reconstruye Cclean=QΛ′QTC_{\mathrm{clean}}=Q\Lambda'Q^{T}Cclean​=QΛ′QT y se reescala para que su diagonal sea exactamente 1, porque la sustitución cambia ligeramente la diagonal. Hay dos refinamientos habituales. El modo de mercado absorbe varianza, así que el nivel de ruido es en realidad σ2≈1−λ1/N\sigma^2\approx1-\lambda_1/Nσ2≈1−λ1​/N y no 1. Y la ley de Marchenko-Pastur supone observaciones independientes y colas finas, así que el borde es una guía, no una prueba precisa; compruebe la matriz depurada fuera de muestra.

Cuatro activos, dieciséis observacionesSupongamos N=4N=4N=4, T=16T=16T=16, así que q=0,25q=0{,}25q=0,25 y λ+=(1+0,5)2=2,25\lambda_+=(1+0{,}5)^2=2{,}25λ+​=(1+0,5)2=2,25. La matriz de correlación muestral tiene valores propios 2,60; 0,80; 0,40; 0,202{,}60;\,0{,}80;\,0{,}40;\,0{,}202,60;0,80;0,40;0,20 (suman N=4N=4N=4). Solo 2,602{,}602,60 supera 2,252{,}252,25: un valor propio de señal, el factor común de mercado. Los otros tres son ruido, con media (0,80+0,40+0,20)/3=0,4667(0{,}80+0{,}40+0{,}20)/3=0{,}4667(0,80+0,40+0,20)/3=0,4667. El espectro depurado es 2,60; 0,4667; 0,4667; 0,46672{,}60;\,0{,}4667;\,0{,}4667;\,0{,}46672,60;0,4667;0,4667;0,4667, que sigue sumando 4 y tiene un valor propio mínimo de 0,470{,}470,47 en lugar de 0,200{,}200,20, de modo que su número de condición baja de 131313 a unos 5,65{,}65,6. Para un problema mayor con q=0,4q=0{,}4q=0,4 el borde superior es (1+0,4)2=2,665(1+\sqrt{0{,}4})^2=2{,}665(1+0,4​)2=2,665.
  • Borde: los valores propios por encima de (1+q)2(1+\sqrt q)^2(1+q​)2 son señal; los del interior del intervalo son ruido.
  • Recorte: sustituir los valores propios de ruido por su media, conservando la traza y el carácter definido positivo.
  • Cierre: reconstruir QΛ′QTQ\Lambda'Q^{T}QΛ′QT y reescalar la diagonal a 1; validar siempre fuera de muestra.
Una referencia, no un veredictoLos rendimientos reales tienen colas gruesas y volatilidad que cambia con el tiempo, lo que dispersa los valores propios de ruido más allá del borde teórico. Trate λ+\lambda_+λ+​ como un umbral inicial razonable para cuántas componentes conservar y deje que una prueba fuera de muestra decida si la depuración realmente ayuda. Eso es exactamente lo que hace el proyecto final.

Lab: qué se califica

  • El borde superior15 pts
  • Contar valores propios de señal10 pts
  • El valor de sustitución15 pts
  • La traza se conserva10 pts
  • El algoritmo de recorte15 pts
  • Por qué la media15 pts
  • Qué significa un valor propio por encima del borde20 pts
0/ 7 respondidos
Ejercicios
1Respuesta numérica15 ptsSin responder

Para N = 200 activos y T = 500 observaciones (q = 0,4), ¿cuál es el borde superior de Marchenko-Pastur (1+q)2(1+\sqrt q)^2(1+q​)2? Dé tres decimales.

lambda+=lambda_+ =lambda+​=aceptado: ± 0.005
2Respuesta numérica10 ptsSin responder

Una matriz de correlación con q = 0,25 (borde superior 2,25) tiene valores propios 2,60; 0,80; 0,40 y 0,20. ¿Cuántos superan el borde?

signal=signal =signal=aceptado: ± 0.01
3Respuesta numérica15 ptsSin responder

Para la misma matriz, los tres valores propios de ruido 0,80; 0,40 y 0,20 se sustituyen por su media. ¿Cuál es el valor común? Dé cuatro decimales.

average=average =average=aceptado: ± 0.001
4Completar los huecos10 ptsSin responder

Una matriz de correlación de N=4N=4N=4 activos tiene traza 4. Tras el recorte, los valores propios siguen sumando lo mismo. Complete el espacio.

∑iλi′=\displaystyle \sum_i\lambda_i'=\htmlData{blank=a}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}}i∑​λi′​=
5Ordenar los pasos15 ptsSin responder

Ponga en orden los pasos del recorte de valores propios.

Arrastra para reordenar, o enfoca un elemento, pulsa Espacio para levantarlo, las flechas para moverlo y Espacio para soltarlo.

  1. 1Reescalar para que la diagonal sea exactamente 1
  2. 2Sustituir los valores propios iguales o inferiores al borde por su media
  3. 3Calcular la matriz de correlación muestral y su descomposición propia
  4. 4Reconstruir la matriz como QΛ′QTQ\Lambda'Q^{T}QΛ′QT
  5. 5Hallar el borde (1+N/T)2(1+\sqrt{N/T})^2(1+N/T​)2
6Elección15 ptsSin responder

¿Por qué se sustituyen los valores propios de ruido por su media y no por cero?

7Explicar con palabras20 ptsSin responder

En una o dos frases, diga qué significa que un valor propio de la matriz de correlación muestral esté por encima del borde superior de Marchenko-Pastur y qué significa que esté dentro.

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