In richelbilderbeek/ribir: ribir, basic phylogenetics page

set.seed(42)
trees <- TreeSim::sim.bd.age(
  age = 10, numbsim = 2, lambda = 0.2, mu = 0.0,
  complete = FALSE
)
ape::plot.phylo(trees[[1]])


gammas <- rep(NA, length(trees))

for (i in seq_along(trees)) {
  gammas[i] <- ape::gammaStat(trees[[i]])
}

gammas

richelbilderbeek/ribir documentation built on March 19, 2021, 3:55 a.m.