Selection model - KamilSJaron/reproductive_mode_TE_dynamics GitHub Wiki

I used the same function as Dolgin and Charlesworth 2006 (implemented in get_fitness.R, the function uses two parameters that correspond to decrease of fitness due to individual TE insertions and multiplicative effect given number of TEs present already in genome (epistasis). These parameters were selected so they show ~2% fitness decrease between 50 and 51 TE copies, a standard yeast TE copy number that was used in the experiment for estimation of fitness effects in yeast. Evaluation is in script calculate_sa_and_sb.R and the resulting fitness function is shown in following figure, where expected simulation space is marked with vertical lines.