simulated_tree <- ape::rphylo(n = 10, birth = 1, death = 0)
gamma_statistic(simulated_tree) # should be around 0.
if (requireNamespace("DDD")) {
ddd_tree <- DDD::dd_sim(pars = c(1, 0, 10), age = 7)$tes
gamma_statistic(ddd_tree) # because of diversity dependence, should be < 0
}
Run the code above in your browser using DataLab