Simulation Negative Selection

Question

[Rubal7]

Hello,

Does anyone know a genetics simulation software that can simulate regions under long-term negative-selection? Coalescent simulators such as msms are only capable of simulating selection at specific sites. We are interested in simulating negative selection over an entire region (so multiple sites simultaneously). I have a suspicion that some forward simulators are capable of doing this but I have no experience with forward simulators. If anyone thinks they know a simulator that can do what we are looking for please let me know. Thanks in advance!

Best,

Rubal

Answer 1

[tiagoantao]

Simupop can do this. It is robust, well documented and with many examples. But it is a programmable simulator: you have to know (or learn) Python. I cannot think of a non-programmable simulator that can do this (easypop comes to mind, but will not do it). Another possibility might be Rmetasim (in R), but I really do not know if it does what you want. simupop is a very safe bet if you know python or are willing to learn.

[Rubal7]

Thanks this looks promising, I'll look into it, appreciate your help!

Answer 2

[Adam]

SFSCODE can simulate data under a huge range of models: http://sfscode.sourceforge.net/SFS_CODE/index/index.html

[Rubal7]

Ah yes I've heard of SFSCODE, this sounds good thanks