The Metropolized Partial Importance Sampling MCMC Mixes Slowly on Minimum Reversal Rearrangement Paths

Markov chain Monte Carlo has been the standard technique for inferring the posterior distribution of genome rearrangement scenarios under a Bayesian approach. We present here a negative result on the rate of convergence of the generally used Markov chains. We prove that the relaxation time of the Markov chains walking on the optimal reversal sorting scenarios might grow exponentially with the size of the signed permutations, namely, with the number of syntheny blocks.

• Type: Journal
• Year: 2010
• Where: TCBB
• Authors: István Miklós, Bence Melykuti, Krister Swenson
• Keywords: Genome Rearrangement Scenarios, Markov Chain Monte Carlo, Markov Chains, Software Engineering, TCBB 2010