Entering edit mode
8.0 years ago
prostoesh
▴
20
Anyone knows how to run MSMS? I am trying to simulate a scenario, where a population splits 2 times (separately), resulting in 3 subpopulations. I used a command
msms -ms 5 3 -t 5 -N 5000000 -I 3 5 0 0 -es 1.5 1 0.5 -es 7.5 1 0.5
but when i parse the output in R's PopGenome with
GENOME.class@populations
i got a list of 1 population, not 3!
is there a way to make this work? i've read the manuals for both ms and msms, but i still got a feeling i'm not 100% sure what exactly i am doing