© 2008 Society of Systematic Biologists
A Bayesian Perspective on a Non-parsimonious Parsimony Model
1 Department of Integrative Biology, University of California Berkeley 3060 VLSB No. 3140, Berkeley, CA 94720-3140, USA; E-mail: johnh{at}berkeley.edu
2 Department of Botany, University of Wisconsin 430 Lincoln Drive, Madison, WI 53706, USA
3 Department of Statistics, University of Wisconsin 1300 University Avenue, Madison, WI 53706, USA
4 Swedish Museum of Natural History Box 50007, SE-104 05 Stockholm, Sweden
Edited by Marc Suchard
 |
Abstract |
|---|
Several stochastic models of character change, when implemented in a maximum likelihood framework, are known to give a correspondence between the maximum parsimony method and the method of maximum likelihood. One such model has an independently estimated branch-length parameter for each site and each branch of the phylogenetic tree. This model—the no-common-mechanism model—has many parameters, and, in fact, the number of parameters increases as fast as the alignment is extended. We take a Bayesian approach to the no-common-mechanism model and place independent gamma prior probability distributions on the branch-length parameters. We are able to analytically integrate over the branch lengths, and this allowed us to implement an efficient Markov chain Monte Carlo method for exploring the space of phylogenetic trees. We were able to reliably estimate the posterior probabilities of clades for phylogenetic trees of up to 500 sequences. However, the Bayesian approach to the problem, at least as implemented here with an independent prior on the length of each branch, does not tame the behavior of the branch-length parameters. The integrated likelihood appears to be a simple rescaling of the parsimony score for a tree, and the marginal posterior probability distribution of the length of a branch is dependent upon how the maximum parsimony method reconstructs the characters at the interior nodes of the tree. The method we describe, however, is of potential importance in the analysis of morphological character data and also for improving the behavior of Markov chain Monte Carlo methods implemented for models in which sites share a common branch-length parameter.
Keywords: Bayesian phylogenetic inference; Markov chain Monte Carlo; maximum likelihood; parsimony model
Received May 11, 2007; Revised August 8, 2007; Accepted January 14, 2008
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  J. Wu and E. Susko General Heterotachy and Distance Method Adjustments Mol. Biol. Evol., December 1, 2009; 26(12): 2689 - 2697. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. Kim and M. J. Sanderson Penalized Likelihood Phylogenetic Inference: Bridging the Parsimony-Likelihood Gap Syst Biol, October 1, 2008; 57(5): 665 - 674. [Abstract] [Full Text] [PDF]
|
|