Combinaison de plusieurs chaînes MCMC parallèles en une seule chaîne plus longue

11

Disons que l'on a couru mchaînes MCMC parallèles où chaque chaîne a subi un rodage. Soit les chaînes résultantes désignées par

X1(je),,XN(je) pour je=1,,m,
N est la longueur de chaque chaîne après rodage.

Si l'on veut combiner ces chaînes en une longue chaîne, est-ce aussi simple que de les concaténer

X1(1),,XN(1),,X1(m),,XN(m)?

Dans mon cas, chacun Xje est un paramétreur θje. Mon objectif est de prélever sur la partie postérieure

p(θy),
ysont les données. La raison pour laquelle je m'intéresse aux chaînes parallèles est qu'elles sont nécessaires au calcul du facteur de réduction d'échelle potentiel (PSRF).
Lotus3000
la source

Réponses:

10

Notez que les chaînes individuelles ont une dépendance série; les valeurs de chaînes distinctes ne le font pas, donc si vous vouliez qu'elle ressemble à une longue chaîne, les concaténer ne serait pas juste.

Cependant, si vous êtes uniquement intéressé par la distribution, l'ordre dans la chaîne n'est pas pertinent. Vous ne cherchez pas réellement à concaténer les chaînes pour cela, vous voulez simplement regrouper toutes les informations de distribution (les traiter comme un grand échantillon). Certes, si les chaînes convergent toutes vers leur distribution stationnaire, elles seront toutes des échantillons de la même distribution - vous pouvez les combiner.

En effet, certaines personnes connaissent une période de rodage et ne tirent qu'une seule valeur de plusieurs chaînes distinctes.

(Garder les pistes séparées pourrait aider à juger si elles ont réellement convergé.)

Si vous calculez la variance qui tient compte de la structure de dépendance, cependant, vous la baserez sur le fait que les différentes exécutions sont indépendantes, mais les valeurs de la même exécution dépendent.

Glen_b -Reinstate Monica
la source
Ainsi, par exemple, si je souhaite trouver une estimation de mon vecteur de paramètres, je pourrais simplement le calculer comme
θ^=je=1mj=1Nθj(je)nm?
Lotus3000,
1
Oui. C'est juste la variance qui est (légèrement) délicate.
Glen_b -Reinstate Monica
Est-ce Gibbs ou autre chose?
Glen_b -Reinstate Monica
Metropolis-Hastings
Lotus3000
Oh, d'accord, pas de problème. J'allais juste suggérer un estimateur différent si c'était Gibbs.
Glen_b -Reinstate Monica