Human Evolution

Suppose we have a population of 50 million. Suppose that we have a gene that gives you a 10% boost to your worth as a member of a community, in the sense that it increases the number of your children who survive to reproductive ager by 10%.

In [10]:
# log odds ratio for 1/10 million

import numpy as np

np.log(1/(10*10**6))

Out[10]:
-16.11809565095832
In [6]:
# log odds ratio for 1/10

np.log(1/10)

Out[6]:
-2.3025850929940455
In [7]:
# log odds ratio for 9/10

np.log(10)

Out[7]:
2.302585092994046
In [11]:
# time to go from 10% to 90%

25*(4.6/.1)

Out[11]:
1149.9999999999998
In [12]:
25*(13.8/.1)

Out[12]:
3450.0
In [ ]: