Many phylogenetic inference software packages, such as BEAST2 and RevBayes, can handle trees containing fossil samples. In the zero-edge format, tip samples are represented normally in the tree, but sampled ancestors, i.e fossil samples which have sampled descendants, are represented as tips at the end of edges with length 0.
This format is included in FossilSim
as the
SAtree
objects.
SAtree
objects can be built from a fossils
object and the associated tree. Tips (both extinct and extant) will be
labelled with the species they were sampled from followed by an index:
the oldest sample for a given species will get index 1 and all other
samples will be ordered from oldest to youngest. The SAtree
format also includes a field complete
which indicates
whether the tree should be considered as including all lineages of the
process or only sampled lineages. Note that if the tree is complete, the
youngest tip of a given extinct species represents the extinction event
for this species and not a fossil sample.
t = ape::rtree(6)
f = sim.fossils.poisson(rate = 2, tree = t)
SAt = SAtree.from.fossils(tree = t, fossils = f)
print(SAt$tree)
##
## Phylogenetic tree with 18 tips and 17 internal nodes.
##
## Tip labels:
## t2_2, t1_5, t3_1, t5_3, t4_1, t6_2, ...
##
## Rooted; includes branch length(s).
## sp edge hmin hmax h tip.label
## 1 1 1 0.6305393 0.6305393 0.6305393 t2_1
## 5 2 2 0.5127315 0.5127315 0.5127315 t1_1
## 2 2 2 0.4771845 0.4771845 0.4771845 t1_2
## 3 2 2 0.1484455 0.1484455 0.1484455 t1_3
## 4 2 2 0.1336602 0.1336602 0.1336602 t1_4
## 7 4 4 0.1122017 0.1122017 0.1122017 t5_1
## 6 4 4 0.1034095 0.1034095 0.1034095 t5_2
## 8 6 6 0.5934814 0.5934814 0.5934814 t6_1
## 9 10 10 1.0961561 1.0961561 1.0961561 t10_1
## 10 10 10 0.7913941 0.7913941 0.7913941 t10_2
## ...
## Fossil record with 12 occurrences representing 6 species
## Fossil record not simulated using taxonomy: all speciation events are assumed to be symmetric
## [1] TRUE
Other functions are present in FossilSim
to modify
trees. prune.fossils
will remove all intermediate fossil
samples and keep only the first and last occurrences of any species.
sampled.tree.from.combined
removes all unsampled
lineages from the tree, and can optionally apply a sampling probability
to extant samples, or keep only a specified list of extant samples.
These two functions can be applied to any object of type
phylo
, but are designed to work on trees which include
fossil samples.