dot-f0: Probability of allele loss

In lbozhilova/kimura: Fitting Kimura Distirbutions to Heteroplasmy Data

Probability of allele loss

Description

Estimate the probability that an allele is lost due to random genetic drift.

Usage

.f0(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 lost 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

.f0(0.3, 0.5)

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