Reversible jump

Reversible jump Markov chain Monte Carlo is an extension to standard Markov chain Monte Carlo (MCMC) methodology that allows simulation of the posterior distribution on spaces of varying dimensions. [cite journal
author = Green, P.J. |authorlink=Peter Green (statistician)
year = 1995
title = Reversible Jump Markov Chain Monte Carlo Computation and Bayesian Model Determination
journal = Biometrika
volume = 82
issue = 4
pages = 711–732
url = http://links.jstor.org/sici?sici=0006-3444(199512)82%3A4%3C711%3ARJMCMC%3E2.0.CO%3B2-F
month = 12
doi = 10.1093/biomet/82.4.711] Thus, the simulation is possible even if the number of parameters in the model is not known.

Let :$n\_min\; N\_m=\{1,2,ldots,I\}$ be a model indicator and the parameter space whose number of dimensions $d\_m$ depends on the model $n\_m$. The model indication need not be finite. The stationary distribution is the joint posterior distribution of $(M,N\_m)$ that takes the values $(m,n\_m)$.

The proposal $m\text{'}$ can be constructed with a mapping$g\_\{1mm\text{'}\}$ of $m$ and $u$, where $u$ is drawn from a random component$U$ with density $q$ on $R^\{d\_\{mm\text{'}$. The move to state $(m\text{'},n\_m\text{'})$ can thus be formulated as

:$(m\text{'},n\_m\text{'})=(n\_m\text{'},g\_\{1mm\text{'}\}(m,u))$

The function

:

must be "one to one", differentiable, and have a non-zero support

:$sup(gmm\text{'})\; e\; varnothing$

so that there exists an inverse function

:$g^\{-1\}\_\{mm\text{'}\}=g\_\{m\text{'}m\}$

that is differentiable. Therefore, the $(m,u)$ and $(m\text{'},u\text{'})$ must be of equal dimension, which is the case if the dimension criterion

:$d\_m+d\_\{mm\text{'}\}=d\_\{m\text{'}\}+d\_\{m\text{'}m\}$

is met where $d\_\{mm\text{'}\}$ is the dimension of $u$. This is known as "dimension matching".

If $R^\{d\_m\}subset\; R^\{d\_\{m\text{'}$ then the dimensional matchingcondition can be reduced to

:$d\_m+d\_\{mm\text{'}\}=d\_\{m\text{'}\}$

with

:$(m,u)=g\_\{m\text{'}m\}(m)$.

The acceptance probability will be given by

:$a(m,m\text{'})=minleft(1,\; frac\{p\_\{m\text{'}m\}p\_\{m\text{'}\}f\_\{m\text{'}\}(m\text{'})\}\{p\_\{mm\text{'}\}q\_\{mm\text{'}\}(m,u)p\_\{m\}f\_m(m)\}left|detleft(frac\{partial\; g\_\{mm\text{'}\}(m,u)\}\{partial\; (m,u)\}\; ight)\; ight|\; ight),$where $|cdot\; |$ denotes the absolute value and $p\_mf\_m$ is the joint posterior probability

:$p\_mf\_m=c^\{-1\}p(y|m,n\_m)p(m|n\_m)p(n\_m),$where $c$ is the normalising constant.

References