birthdeath.simulation {PhySim} | R Documentation |
simulates batches of birth death trees for a given time period.
birthdeath.simulation(b, d, T, N)
b |
birth or speciation rate |
d |
death or extinction rate |
T |
number of time units to simulate each tree for |
N |
number of trees to simulate |
To simulate pure birth (Yule) trees set extinction to 0.
Writes trees to two files with extinct species included (.treeExt) or excluded (.treeNoE). Returns the number of trees that went extinct, left one lineage and two or more lineages at time = T.
Jason T. Weir
Weir, J. T. and Schluter, D. 2007. Science in press
birthdeath.tree
to simulate a single tree