Artículos de revistas
Comparing the performance of a reversible jump Markov chain Monte Carlo algorithm for DNA sequences alignment
Registro en:
Journal Of Statistical Computation And Simulation. Taylor & Francis Ltd, v. 76, n. 7, n. 567, n. 584, 2006.
0094-9655
WOS:000237679800001
10.1080/10629360500109226
Autor
Alvarez, LJ
Garcia, NL
Rodrigues, ER
Institución
Resumen
Assume that K independent copies are made from a common prototype DNA sequence whose length is a random variable. In this paper, the problem of aligning those copies and therefore the problem of estimating the prototype sequence that produced the copies is addressed. A hidden Markov chain is used to model the copying procedure, and a reversible jump Markov chain Monte Carlo algorithm is used to sample the parameters of the model from their posterior distribution. Using the sample obtained, the Bayesian model and the prototype sequence may be selected using the maximum a posteriori estimate. A prior distribution for the prototype DNA sequence that incorporates a correlation among neighbouring bases is also considered. In addition, an analysis of the performance of the algorithm is presented when different scenarios are taken into account. 76 7 567 584