makeWellBehaved: Make mutation model well behaved

In thoree/mut2: Reversible mutation matrices and pairwise likelihood ratios

Make mutation model well behaved

Description

Replace q_ij by min(q_ij,p_j)

Usage

makeWellBehaved(Q, afreq)

Arguments

Q	a mutation model
afreq	A numeric vector of allele frequencies.

Value

A reversible stepwise mutation model.

Author(s)

Thore Egeland.

Examples

## Not run: 
# Not well behaved:
p = c(0.4, 0.6)
Q = stepwiseReversible(alleles = 1:2, afreq = p, rate = 0.6, range = 0.1)
makeWellBehaved(Q, afreq = p)

## End(Not run)

thoree/mut2 documentation built on May 16, 2023, 7:56 p.m.