initialize() { initializeMutationRate(u); initializeMutationType("m1",0.5,"f",0.0); initializeMutationType("m2",0.5,"f",0.05); 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.000236,0.000236)); p2.setMigrationRates(c(p1,p5),c(0.000236,0.000236)); p5.setMigrationRates(c(p1,p2),c(0.000236,0.000236));} 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");} 40403 { sim.simulationFinished(); } fitness(m2, p1) {return 1.5 - sim.mutationFrequencies(p1, mut);} fitness(m2, p2) {return 1.5 - sim.mutationFrequencies(p2, mut);} fitness(m2, p3) {return 1.5 - sim.mutationFrequencies(p3, mut);} fitness(m2, p4) {return 1.5 - sim.mutationFrequencies(p4, mut);} fitness(m2, p5) {return 1.5 - sim.mutationFrequencies(p5, mut);} fitness(m2, p6) {return 1.5 - sim.mutationFrequencies(p6, mut);}