A beta-splitting model for evolutionary trees
2014, 2015, Raazesh Sainudiin and Amandine Veber
This is code accompanying the technical report 'A beta-splitting model for evolutionary trees' by Raazesh Sainudiin and Amandine Veber. The code is for clarity and simplicity of implementation (it is not efficient for large scale simulations).
Probability of tree under the generating-organizing DTDS Markov chain at various resolutions.
Demonstration with a=b=0 (Yule model)
((1, 3, 2), 0.1667, 0.1700)
((3, 2, 1), 0.1667, 0.1666)
((2, 1, 3), 0.1667, 0.1625)
((3, 1, 2), 0.1667, 0.1664)
((1, 2, 3), 0.1667, 0.1683)
((2, 3, 1), 0.1667, 0.1662)
([., [[., .], .]], 0.16667, 0.17000)
([., [., [., .]]], 0.16667, 0.16830)
([[[., .], .], .], 0.16667, 0.16660)
([[., [., .]], .], 0.16667, 0.16620)
([[., .], [., .]], 0.33333, 0.32890)
(0.66667, 0.67110)
(0.33333, 0.32890)
Demonstration for more balanced trees with
((1, 3, 2), 0.1251, 0.1229)
((3, 2, 1), 0.1251, 0.1242)
((2, 3, 1), 0.1251, 0.1230)
((3, 1, 2), 0.2499, 0.2571)
((1, 2, 3), 0.1251, 0.1216)
((2, 1, 3), 0.2499, 0.2512)
([[., [., .]], .], 0.1251, 0.1230)
([., [., [., .]]], 0.1251, 0.1216)
([[[., .], .], .], 0.1251, 0.1242)
([., [[., .], .]], 0.1251, 0.1229)
([[., .], [., .]], 0.4998, 0.5083)
(0.50025, 0.49170)
(0.49975, 0.50830)
(0.093984, 0.098800)
(0.37463, 0.37940)
(0.015703, 0.014600)
(0.015687, 0.016600)
(0.24988, 0.24850)
(0.25012, 0.24210)