Replica-exchange Monte Carlo scheme for Bayesian data analysis

We develop a sampling algorithm to explore the probability densities arising in Bayesian data analysis problems. Our algorithm is a multiparameter generalization of a replica-exchange Monte Carlo scheme. The strategy relies on gradual weighing of experimental data and on Tsallis generalized statistics. We demonstrate the effectiveness of the method on nuclear magnetic resonance data for a folded protein.