¿Se puede calcular la autocorrelación de matrices de covarianza muestreadas por MCMC?


8

Imagine que muestreamos una matriz de covarianza de una distribución Wishart por MCMC.

En cada iteración, obtenemos una nueva matriz de muestra Si de la distribución Wishart.

Pregunta : Dada la traza que contiene todas las muestrasS1,...Sn, ¿puedo trazar la autocorrelación de estas muestras?

He visto a alguien usando la autocorrelación de log(det(S)), pero no encontré justificación.

Respuestas:


1

A mi modo de ver, cuando dibuja matrices de una distribución Wishart, realmente está dibujando variables aleatorias univariadas específicamente relacionadas en cada paso de tiempo. Es decir, solo la parte de cualquier [= la parte triangular superior] es aleatoria, y la simetría te da el resto. En otras palabras, la autocorrelación se define en parejas para cualquiera de las dos entradas de -dimensional vectors y para cualquier . Por supuesto, esto te deja con un número posiblemente inviable dep×pS1,S2,Snp(p+1)/2vech(Si)Sip(p+1)/2vech(Si)vech(Sih)h>0(p(p+1)/2)2autocorrelaciones univariadas para rastrear, y dado que las entradas de cada estarán significativamente relacionadas entre sí (ver, por ejemplo, la definición dada a través de dibujos de una normal aquí: https://en.wikipedia.org/ wiki / Wishart_distribution ), me imagino que descartas información haciendo este análisis univariante. Dicho esto, las autocorrelaciones univariadas se pueden calcular por entrada definiendo primero Claramente,vech(Si)

S¯=1ni=1nvech(Si)Sh¯=1nhi=h+1nvech(Si)vech(Sih)T.
S¯es un estimador natural de la parte triangular superior de la expectativa (que puede reemplazar con la verdadera expectativa de su distribución Wishart si la conoce). Del mismo modo, es un estimador natural por el momento . Por último, observando que se llega a las estimaciones de autocorrelación para través de Como se mencionó anteriormente, esto le da laS¯hE(vech(Si)vech(Sih)T)
Cov(vech(Si)vech(Sih)T)=E(vech(Si)vech(Sih)T)E(vech(Si))E(vech(Si))T,
A(h)vech(Si)
A(h)=S¯hS¯S¯T.
(p(p+1)/2)2 autocorrelaciones de cada entrada de matriz de Wishart entre sí la entrada de matriz de Wishart. Si esa es demasiada información para mostrar, creo que una estrategia que podría tomar sería definir la serie de tiempo univariada es decir, simplemente toma el promedio del valor absoluto de la autocorrelación. Si solo le interesan las autocorrelaciones positivas y no cree que la autocorrelación negativa sea perjudicial, entonces evite el valor absoluto. Del mismo modo, si cree que la autocorrelación a lo largo de la diagonal es peor que fuera de la diagonal o al revés, podría agregar pesos
a(h)=1(p(p+1)/2)2i=1(p(p+1)/2j=1(p(p+1)/2|A(h)ij|,
wij que toman en cuenta esta 'función de pérdida'.
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.