Estoy viendo un problema cuando callibrating un MP de distribución. Asumir un registro de regreso de la serie para el SP500, con las siguientes dimensiones
dim(xts.sp500.ret.stocksonly)
==> [1] 1133 478
sp500.cor <- cor.empirical(xts.sp500.ret.stocksonly)
sp500.eigens <- eigen(sp500.cor)$values
sp500.eigen.density <- density(sp500.eigens,n=5000)
plot(sp500.eigen.density,xlim=c(0,4),main="sp500 returns eigenvalue density")
Supongo que mi 'Q' valor es 1133/478 =
Problema: Aunque la 'forma' y atajos que parecen OK -- el valor de la densidad (eje y) parecen muy APAGADO. Pico de 1.5 para la serie real - 5 o así que para el teórico, (nota estoy truncar la trama así que el mercado de valores no se muestran, son enormes alrededor de 200).
Pregunta: 1) Es esta espera? 2) ¿Cómo afecta esto a callibration? Debo confiar en los resultados y basta con mirar los puntos de corte? 3) También cuando 'limpiar' la matriz veo la mayoría de código (por ejemplo, filtro.RMT en tawny) simplemente se reemplaza por debajo de valores de Lambda+ con el promedio, lo que acerca de Lambda - aunque?
gracias!