boundsGamma: Upper limit for overall mutation rate

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

Upper limit for overall mutation rate

Description

Upper limit for overall mutation rate

Usage

boundsGamma(alleles, afreq, range)

Arguments

alleles	A character vector with allele labels.
afreq	A numeric vector of allele frequencies.
range	A positive number.

Value

Bounds on gamma

Author(s)

Thore Egeland.

Examples

alleles = 1:3; afreq = c(0.2, 0.3,  0.5);  range = 0.1
bounds = boundsGamma(alleles, afreq ,range)
R1 = stepwiseReversible(alleles, afreq, rate = bounds[1], range)
R2 = stepwiseReversible(alleles, afreq, rate = bounds[2] range)
isRegular(R2, afreq)
                        

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