Calculate distribution of allele difference — relationship_allele_diff

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.

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