À mon avis, lorsque vous dessinez matrices partir d'une distribution Wishart, vous dessinez vraiment des variables aléatoires univariées spécifiquement liées à chaque pas de temps. C'est-à-dire que seule la partie de tout [= la partie triangulaire supérieure] est aléatoire, et la symétrie vous donne le reste. En d'autres termes, l'autocorrélation est définie par paire pour deux entrées quelconques des vecteurs et pour tout . Bien sûr, cela vous laisse avec un nombre potentiellement infaisablement élevé dep×pS1,S2,…Snp⋅(p+1)/2vech(Si)Sip⋅(p+1)/2vech(Si)vech(Si−h)h>0(p⋅(p+1)/2)2des autocorrélations univariées à suivre, et puisque les entrées de chaque seront corrélées les unes aux autres (voir par exemple la définition donnée via les tirages d'une normale ici: https://en.wikipedia.org/ wiki / Wishart_distribution ), je pourrais bien imaginer que vous jetez des informations en faisant cette analyse univariée. Cela étant dit, les autocorrélations univariées peuvent être calculées en entrée en définissant d'abord
Clairement,vech(Si)
S¯Sh¯=1n∑i=1nvech(Si)=1n−h∑i=h+1nvech(Si)vech(Si−h)T.
S¯est un estimateur naturel pour la partie triangulaire supérieure de l'attente (que vous pouvez remplacer par l'attente réelle de votre distribution de Wishart si elle vous est connue). De même, est un estimateur naturel pour le moment . Enfin, notant que
on arrive aux estimations d'autocorrélation pour via
Comme mentionné précédemment, cela vous donne la
S¯hE(vech(Si)vech(Si−h)T)Cov(vech(Si)vech(Si−h)T)=E(vech(Si)vech(Si−h)T)−E(vech(Si))E(vech(Si))T,
A(h)vech(Si)A(h)=S¯h−S¯S¯T.
(p⋅(p+1)/2)2 autocorrélations de chaque entrée de matrice Wishart avec l'autre entrée de matrice Wishart. Si c'est trop d'informations à afficher, je pense qu'une stratégie que vous pourriez adopter serait de définir la série temporelle univariée
c'est-à-dire que vous prenez simplement la moyenne de la valeur absolue de l'autocorrélation. Si vous ne vous souciez que des autocorrélations positives et ne pensez pas qu'une autocorrélation négative est nuisible, alors épargnez-vous De même, si vous pensez que l'autocorrélation le long de la diagonale est pire que hors de la diagonale ou , vous pouvez ajouter des poids
a(h)=1(p⋅(p+1)/2)2∑i=1(p⋅(p+1)/2∑j=1(p⋅(p+1)/2|A(h)ij|,
wij qui tiennent compte de cette «fonction de perte».