Difference between revisions of "Probability of fixation"
From Genetics Wiki
| Line 1: | Line 1: | ||
| + | This was derived in [[Kimura 1962]]. | ||
| + | |||
<math>u(p)=\frac{1-e^{4N_esp}}{1-e^{4N_es}}</math> | <math>u(p)=\frac{1-e^{4N_esp}}{1-e^{4N_es}}</math> | ||
| Line 12: | Line 14: | ||
<math>u(p)_2 \approx 1-e^{2s} \approx 1-1+2s = 2s</math>. | <math>u(p)_2 \approx 1-e^{2s} \approx 1-1+2s = 2s</math>. | ||
| + | |||
| + | This agrees with the results of [[Fisher 1930]] and [[Wright 1931]]. | ||
Revision as of 21:16, 22 September 2018
This was derived in Kimura 1962.
[math]u(p)=\frac{1-e^{4N_esp}}{1-e^{4N_es}}[/math]
If we are considering the initial frequency of a single new mutation in the population p=1/(2Ne),
[math]u(p)_1=\frac{1-e^{4N_es\frac{1}{2N_e}}}{1-e^{4N_es}}=\frac{1-e^{2s}}{1-e^{4N_es}}[/math].
And if 4Nes is large
[math]u(p)_2\approx\frac{1-e^{2s}}{1}=1-e^{2s}[/math].
[math]e^{2s}\approx 1+2s[/math]
[math]u(p)_2 \approx 1-e^{2s} \approx 1-1+2s = 2s[/math].
This agrees with the results of Fisher 1930 and Wright 1931.
