Description Usage Arguments Value Author(s) References See Also Examples
LTT.plot.gen calculates the number of lineages through time for each input tree, as well as the average number of lineages over all trees. The trees may be simulated using any function in TreeSim, or may be empirical trees. The method works for ultrametric and non-ultrametric trees which are binary or have polytomies.
1 | LTT.plot.gen(trees,bound=10^(-12))
|
trees |
List with one or two entries. First component: list of phylogenetic trees; second component: vector with time of origins (can be empty). |
bound |
Determines the value by which leaf times may differ in an ultrametric tree. If two tips are further apart, they are considered as sequentially sampled tips. |
out |
out[[1]]: First column are the branching times of ALL input trees. Second column is the number of lineages after the branching time. out[[i]]: Equivalent vector as out [[1]], but for tree i-1. |
Tanja Stadler
T. Stadler. Simulating trees on a fixed number of extant species. Syst. Biol., 60: 676-684, 2011.
LTT.plot,sim.bd.taxa, sim.bd.age, sim.rateshift.taxa, sim.gsa.taxa, birthdeath.tree
1 2 3 4 5 6 7 8 9 10 | # Simulation of a tree with a mrca at time 10 in the past,
# under the density-dependent model
numbsim<-10
age<-10
lambda<-0.3
mu<-0.2
K<-40
# You can produce LTT plots as follows.
# (for now this is un-commented, as some combinations of geiger / TreeSim on certain platforms
# produce problems. If this is the case for you, please report to tanjs.stadler@bsse.ethz.ch).
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.