title: Quantifying MCMC exploration of phylogenetic tree space.
abstract: In order to gain an understanding of the effectiveness of phylogenetic Markov chain Monte Carlo (MCMC), it is important to understand how quickly the empirical distribution of the MCMC converges to the posterior distribution. In this article, we investigate this problem on phylogenetic tree topologies with a metric that is especially well suited to the task: the subtree prune-and-regraft (SPR) metric. This metric directly corresponds to the minimum number of MCMC rearrangements required to move between trees in common phylogenetic MCMC implementations. We develop a novel graph-based approach to analyze tree posteriors and find that the SPR metric is much more informative than simpler metrics that are unrelated to MCMC moves. In doing so, we show conclusively that topological peaks do occur in Bayesian phylogenetic posteriors from real data sets as sampled with standard MCMC approaches, investigate the efficiency of Metropolis-coupled MCMC (MCMCMC) in traversing the valleys between peaks, and show that conditional clade distribution (CCD) can have systematic problems when there are multiple peaks.
authors: Program in Computational Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA cwhidden@fhcrc.org.
authors: Chris
authors: C
authors: Whidden
authors: Program in Computational Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA.
authors: Frederick A
authors: FA
authors: Matsen
mesh_terms: D001105:Archaea / Q000145:classification / Q000235:genetics; D001419:Bacteria / Q000145:classification / Q000235:genetics; D001499:Bayes Theorem; D002965:Classification / Q000379:methods*; D056890:Eukaryota / Q000145:classification / Q000235:genetics; D008390:Markov Chains; D008957:Models, Genetic; D009010:Monte Carlo Method; D010802:Phylogeny
keywords: Markov chain Monte Carlo; phylogenetic methods; subtree prune-and-regraft; topological peaks; tree space
references: Syst Biol. 2012 May;61(3):539-42
references: 22357727
references: J Mol Evol. 1996 Sep;43(3):287-92
references: 8703095
references: Syst Biol. 2014 Jul;63(4):566-81
references: 24695589
references: Syst Biol. 2008 Feb;57(1):86-103
references: 18278678
references: Genome Res. 2003 Nov;13(11):2498-504
references: 14597658
references: Syst Biol. 2013 Nov;62(6):901-12
references: 23925510
references: J Mol Evol. 1985;22(2):160-74
references: 3934395
references: Syst Biol. 2005 Jun;54(3):471-82
references: 16012112
references: Mol Biol Evol. 2012 Aug;29(8):1969-73
references: 22367748
references: Mycologia. 2003 Jul-Aug;95(4):561-4
references: 21148964
references: Bioinformatics. 2006 Aug 15;22(16):2047-8
references: 16679334
references: Syst Biol. 2005 Apr;54(2):241-53
references: 16012095
references: Syst Biol. 2006 Aug;55(4):553-65
references: 16857650
references: Mol Biol Evol. 1990 Nov;7(6):607-33
references: 2283953
references: Evol Bioinform Online. 2008 Feb 09;4:17-27
references: 19204804
references: Syst Biol. 2006 Aug;55(4):652-61
references: 16969941
references: Syst Biol. 2008 Jun;57(3):406-19
references: 18570035
references: Mol Ecol. 2004 Apr;13(4):757-73
references: 15012754
references: Syst Biol. 2012 Jan;61(1):1-11
references: 21828081
references: Science. 2006 Apr 21;312(5772):367; author reply 367
references: 16627724
references: PLoS Comput Biol. 2014 Apr;10(4):e1003537
references: 24722319
references: Syst Biol. 2013 Jul;62(4):501-11
references: 23479066
references: Bioinformatics. 2001 Aug;17(8):754-5
references: 11524383
references: Science. 2005 Sep 30;309(5744):2207-9
references: 16195459
references: Syst Biol. 2003 Oct;52(5):705-16
references: 14530137
references: BMC Evol Biol. 2007;7:214
references: 17996036