initialize() { initializeMutationRate(u); initializeMutationType("m1",0.5,"f",0.0); initializeMutationType("m2",2000,"f",0.00001); initializeGenomicElementType("g1",c(m1,m2),c(671,1)); initializeGenomicElement(g1, 1, 672); initializeRecombinationRate(6e-7); } 1 { sim.addSubpop("p1",2550);} 25000 { sim.addSubpopSplit("p2",1630,p1); p1.setSubpopulationSize(1535);} 26967 { sim.addSubpopSplit("p3",1491,p2); sim.addSubpopSplit("p4",3376,p2); p2.setSubpopulationSize(2800);} 33369 { sim.addSubpopSplit("p5",4754,p2); sim.addSubpopSplit("p6",1447,p2); p2.setSubpopulationSize(2683);} 40285 { p1.setMigrationRates(c(p2,p5),c(0.00005,0.00005)); p2.setMigrationRates(c(p1,p5),c(0.00005,0.00005)); p5.setMigrationRates(c(p1,p2),c(0.00005,0.00005));} 40403 late() { s1 = sample(p2.individuals,20); s2 = sample(p5.individuals,20); s3 = sample(p4.individuals,20); s4 = sample(p3.individuals,20); s5 = sample(p1.individuals,40); s6 = sample(p6.individuals,18); tot = c(s1,s2,s3,s4,s5,s6); tot.genomes.outputMS(filePath="output_"+simnum+"_"+locus+".ms");}