# simulate tree
t = ape::rtree(6)
# simulate fossils
f = sim.fossils.poisson(rate = 2, tree = t)
# transform format
t2 = SAtree.from.fossils(t,f)$tree
# prune fossils
t4 = prune.fossils(t2)
# or transform to sampled tree first
t3 = sampled.tree.from.combined(t2)
t4 = prune.fossils(t3)
plot(t4)
Run the code above in your browser using DataLab