Simulated data

2016-11-25T14:15:57Z (GMT) by Michael Blum
<div>##################################################</div><div>Introduction</div><div>##################################################</div><div><br></div><div>The data was simulated with ms (Hudson, Bioinformatics, 2002). The command that invoked ms (including the seed) and the simulated data are provided. For each simulation, the ms output was parsed into two files that are directly used for the analysis.</div><div>100 unlinked loci were simulated, where each of them contains 100 SNPs. For each of the 100 demes, 20 chromosomes are simulated.</div><div><br></div><div>##################################################</div><div>Folders and parameters</div><div>##################################################</div><div><br></div><div>- Secondary contact models -</div><div><br></div><div>Parameters and files:</div><div>mu mu_B t_C d t_s bottleneck_intensity archivename</div><div>0.005 0.005 0 0.125 19 0.1 seccon.sim0.tgz </div><div>0.005 0.005 1 0.125 19 0.1 seccon.sim1.tgz </div><div>0.005 0.005 2 0.125 19 0.1 seccon.sim2.tgz </div><div>0.005 0.005 3 0.125 19 0.1 seccon.sim3.tgz </div><div>0.005 0.005 4 0.125 19 0.1 seccon.sim4.tgz </div><div>0.005 0.005 5 0.125 19 0.1 seccon.sim5.tgz </div><div>0.005 0.0005 1 0.125 19 0.1 seccon.sim31.tgz </div><div>0.005 0.0005 2 0.125 19 0.1 seccon.sim32.tgz </div><div>0.005 0.0005 3 0.125 19 0.1 seccon.sim33.tgz </div><div>0.005 0.0005 4 0.125 19 0.1 seccon.sim34.tgz </div><div>0.005 0.0005 5 0.125 19 0.1 seccon.sim35.tgz </div><div>0.005 0.005 1 0.125 19 0.25 seccon.sim51.tgz </div><div>0.005 0.005 1 0.125 19 0.5 seccon.sim52.tgz </div><div>0.005 0.005 1 0.25 19 0.1 seccon.sim61.tgz </div><div>0.005 0.005 1 0.375 19 0.1 seccon.sim62.tgz </div><div>0.01 0.01 1 0.125 19 0.1 seccon.sim71.tgz </div><div>0.05 0.05 1 0.125 19 0.1 seccon.sim72.tgz </div><div><br></div><div>main text (fig 3): sim0-sim5</div><div>SI:</div><div>- bottleneck intensity (fig S3): sim1, sim51-52</div><div>- expansion speed (fig S4): sim1, sim61-62</div><div>- migration rate (fig S5): sim1, sim71-72</div><div>- secondary contact model with barrier (fig S6): sim0, sim31-35</div><div><br></div><div><br></div><div>- Stepping stone models -</div><div><br></div><div>Parameters and files:</div><div>mu mu_B archivename</div><div>0.005 0.005 step.sim10.tgz</div><div>0.005 0.0005 step.sim11.tgz</div><div>0.005 0.00005 step.sim12.tgz</div><div>0.005 0.001 step.sim13.tgz</div><div>0.005 0.0001 step.sim14.tgz</div><div>0.005 0.00001 step.sim15.tgz</div><div>0.01 0.00002 step.sim71.tgz</div><div>0.05 0.0001 step.sim72.tgz</div><div><br></div><div>main text (fig 3): sim10-sim15</div><div>SI: </div><div>- migration rate (fig S5): sim15, sim71-72</div><div><br></div><div><br></div><div>############################################</div><div>Files: </div><div>############################################</div><div><br></div><div>- Filenames -</div><div><br></div><div>sim* (see archivename) + number of the replicate (1-100)</div><div><br></div><div>- Filetypes - </div><div><br></div><div>- *.mscmd: Command that invoked ms.</div><div>- *.posout: Positions of the segregating sites.</div><div> One line per unlinked locus. For each locus, the positions of the 100 SNPs on a scale of (0, 1) are given, separated by blanks. </div><div>- *.parseout: Segregating sites.</div><div> One line per chromosome. Each line consists of the number of the deme (1-100, i. e. there are 20 lines per deme) followed by a sequence of 0s, and 1s, 0 for the ancestral and 1 for the derived allele at each polymorphic site, separated by blanks.</div><div><br></div>