R/calc_expected_loctree_values.R

In wadedismukes/treeduckenTools: treeduckenTools

# based on these being exponentially distributed times to event we just take the
# rate and multiply by the tree length
expected_number_transfers <- function(tr_rate, tree_length) {
    return(tr_rate * tree_length)
}

expected_number_births <- function(gene_brate, tree_length) {
    return(gene_brate * tree_length)
}

expected_number_births <- function(gene_drate, tree_length) {
    return(gene_drate * tree_length)
}