Calculate distribution of allele difference after m meioses.

relationship_allele_diff_dist(meioses, mu_dw, mu_up, method = "explicit")

Arguments

meioses

number of meioses separating the two individuals

mu_dw

mutation rate for 1-step down-mutation

mu_up

mutation rate for 1-step up-mutation

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)