Learn R Programming

pedmut (version 0.8.0)

findStationary: Find the stationary frequency distribution

Description

Finds the stationary distribution of allele frequencies, if it exists, w.r.t. a given mutation matrix.

Usage

findStationary(mutmat)

Value

A vector of length ncol(mutmat), or NULL.

Arguments

mutmat

A mutation matrix.

Examples

Run this code

m1 = mutationMatrix("equal", alleles = 1:4, rate = 0.1)
findStationary(m1)

m2 = mutationMatrix("random", alleles = 1:3, seed = 123)
a = findStationary(m2)

a %*% m2 - a  # check

Run the code above in your browser using DataLab