This is the companion web page for the paper

Olivier Cappé, Christian Robert and Tobias Rydén. Reversible jump, birth-and-death and more general continuous time Markov chain Monte Carlo samplers. Journal of the Royal Statistical Society: Series B (Statistical Methodology), Volume 65, Issue 3, pages 679-700, 2003.

The purpose of the software available from this page is to implement transdimensional Markov Chain Monte Carlo (MCMC) for inference in (scalar) Gaussian mixture models, with unknown number of components. Two methods are implemented: Reversible Jump (RJ) MCMC for rj_mix and Continuous Time (CT) for ct_mix

This is the C code that was used to produce the figures in Section 4 of paper. This sofware implements transdimensional MCMC (markov Chain Monte Carlo) for (scalar) Gaussian mixture using either Reversible Jump of Continuous Time simulation.

The software is C code which is distributed in source form (please read the documentation for further details):

The current version of the software is 3.6 (the current version and time stamp are printed on each page of the manual, each individual file has its own version number which may differ).

You can find a preprint version of the JRSSB paper here (gzip-compressed ps file) (note: this is not the published version and is provided for your personal use only).

If you use this code for anything useful (other than finding bugs!), we would appreciate that you acknowledge this page through something like


@misc{ctrj_mix,
  author    = {O. Capp\'{e} and C. Robert and T. Ryd\'{e}n},
  title     = {{CT/RJ-Mix}: {T}ransdimensional {MCMC}
                  for {G}aussian mixtures (available as {C} source code)},
  url       = {http://www.tsi.enst.fr/~cappe/Code/CTRJ_mix},
  month     = dec,
  year      = {2003},
}