powered by
Finds the stationary distribution of allele frequencies, if it exists, w.r.t. a given mutation matrix.
findStationary(mutmat)
A vector of length ncol(mutmat), or NULL.
ncol(mutmat)
A mutation matrix.
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