linguaPhylo
linguaPhylo copied to clipboard
Aged conditional coalescent
Conditional coalescent on root age, it may take the PopulationFunction
as input.
tree ~ AgedConditionalCoalescent( ... , ParametricDistribution);
- simulate 1000 trees from PopulationFunction;
- weight each tree by Pr(tree, age | ParametricDistribution), such as
Pr(tree, age | meanlog, sdlog)
for LogNormal; - resample 1 tree proportional to there weights.