reasonableMutationMatrix: Checks if pj < mij

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

Checks if pj < mij

Description

Entries where pj < mij are unreasonable and are detected

Usage

reasonableMutationMatrix(mutmat, p, verbose = FALSE)

Arguments

mutmat	A mutation matrix.
p	A vector with allele frequencies of the same length as the size of mutmat.
verbose	Logical.

Value

Number of cases where pj < mij

Author(s)

Thore Egeland

Examples

n = 4
p = 1:n/sum(1:n)
names(p) = 1:n
mutmat = mutationMatrix("equal", p, rate = 0.5)
reasonableMutationMatrix(mutmat, p, verbose = FALSE)

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