birthdeath.simulation {PhySim}R Documentation

Simulate batches of birth death Trees

Description

simulates batches of birth death trees for a given time period.

Usage

birthdeath.simulation(b, d, T, N)

Arguments

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

Details

To simulate pure birth (Yule) trees set extinction to 0.

Value

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.

Author(s)

Jason T. Weir

References

Weir, J. T. and Schluter, D. 2007. Science in press

See Also

birthdeath.tree to simulate a single tree


[Package PhySim version 1.0 Index]