R/birth_death_evolution_rates.R

In BrunoVilela/FARM: A package that tells you how farm spread

# BirthDeath function
bd <- function(tree) {
  tree$edge.length <- tree$edge.length / max(tree$edge.length) 
  x <- birthdeath(tree)  
  b <- x$para[2] / (1 - x$para[1])
  d <- b - x$para[2]
  c(setNames(c(b, d), c("b", "d")), x$para)
}