pr_double_hit: Calculate probability of two mutations on three-species tree

In guerreror/pepo: Estimate Hemiplasy and Homoplasy on a Phylogeny

Description

One of two types of double-mutation events will be returned, depending on the value of where. Mutation and back-mutation are assumed to have the same rate.

"ac"–> product of the following probabilities: mutation on the [a] and [c] branches, no mutations on [b], ancestral, or internal branches. "banc" –> product of the following probabilities: mutation on the [abc] ancestor and on the [b] branch, no mutation on [a],[c] or internal branches.

Usage

pr_double_hit(mu_list, where)

Arguments

mu_list	List with "anc", "a", "b", "c", and "bc" entries, each of which is a mutation probability.
where	String. "ac" or "banc".

Value

Numeric. Probability

See Also

homoplasy_mutations

guerreror/pepo documentation built on May 23, 2019, 9:01 p.m.