Calculate distribution of allele difference for symmetric mutation rates

Calculate distribution of allele difference after m meioses.

relationship_allele_diff_dist_sym(meioses, mu_updw, method = "explicit")

Arguments

meioses

number of meioses separating the two individuals

mu_updw

mutation rate for 1-step down- and up-mutations, i.e. total mutation rate is 2*mu_updw

method

"explicit" (default): use known formulas for eigenvalues and eigenvectors. Can cause numerical problems. "matmult": do matrix multiplication instead of diagonalisation. "matmult_mpfr": as "matmult" but with the Rmpfr library (note that this returns list instead of data.frame). "r_eigen": use R's eigen() function to find eigen values. Mostly for debugging.

Value

data.frame with columns d (allele difference) and p (prob)