dot-f1: Probability of fixing allele

In lbozhilova/kimura: Fitting Kimura Distirbutions to Heteroplasmy Data

Probability of fixing allele

Description

Estimate the probability that an allele is fixed (i.e. the wildtype is lost) due to random genetic drift.

Usage

.f1(p, b)

Arguments

p	number in (0; 1): initial heteroplasmy.
b	number in (0; 1): drift parameter.

Details

The parameter p corresponds to the initial heteroplasmy, and therefore is also equal to the mean heteroplasmy after genetic drift.

The parameter b corresponds to the amount of drift, and is equal to b = exp(-t/N_eff), where t denotes time and N_eff denotes effective sample size.

Value

The probability the allele is fixed due to genetic drift.

References

Wonnapinij, Passorn, Patrick F. Chinnery, and David C. Samuels. "The distribution of mitochondrial DNA heteroplasmy due to random genetic drift." The American Journal of Human Genetics 83.5 (2008): 582-593.

Kimura, Motoo. "Solution of a process of random genetic drift with a continuous model." Proceedings of the National Academy of Sciences of the United States of America 41.3 (1955): 144.

Examples

.f1(0.3, 0.5)

lbozhilova/kimura documentation built on May 7, 2023, 12:07 p.m.