1 votos

Normal multivariante cuando la descomposición de Cholesky falla en Sigma

Estoy tratando de hacer distribuciones multivariadas de rendimientos en cubos donde todos los rendimientos están correlacionados al menos 0,6 con un nivel de confianza del 95%. Tengo los cubos, pero sus Sigmas no se pueden descomponer y mi generador de números aleatorios falla.

¿Qué alternativas tengo? Si divido los cubos en otros más pequeños, pierdo la correlación en la simulación. ¿O estoy completamente equivocado?

0 votos

¿Es su matriz $\Sigma$ ¿Singular?

0 votos

Sólo puedo suponer que lo es. Estoy ejecutando cubos más pequeños y el uso de medios de rodadura para tratar de añadir en algunos de la correlación histórica.

0 votos

Creo que no entiendo tu configuración pero puedes comprobar si la matriz es singular y si esa es la causa del problema.

1voto

Tienes que ajustar tu matriz de correlación para que se convierta en positiva definida.

Hay una rutina de R que hará esto por ti - enlace .

O, si quieres hacerlo tú mismo, creo que el método general es hacer una descomposición de valores propios, poner a cero cualquier valor propio negativo y luego reconstruir la matriz original.

Sin embargo, si vas a seguir este camino, es mejor que obtengas tu $\Sigma^\frac{1}{2}$ matriz utilizando la EVD en lugar de la cholesky.

Finanhelp.com

FinanHelp es una comunidad para personas con conocimientos de economía y finanzas, o quiere aprender. Puedes hacer tus propias preguntas o resolver las de los demás.

Powered by:

X