.par file starts, but no genealogies are built

[bekah...@gmail.com]

Hello,

I am trying to simulate a simple admixture scenario with two populations. Fastsimcoal will start, however, it stops at this step:

"fastsimcoal2 is building 100 genealogies ...

Simulating 1 independent chromosomes using 12 batches and 1 threads."

No genealogies are built and the results file only contains a blank .arp file, a blank .arb file and the .simparam file. The simparam file is correctly identifying my input parameters.

Previous input files I have work as do the example input files that came with the program. I have tried re-making my input file several times and using the example input files as a template, but with no luck. Anyone have any other ideas as to what might be happening?

I am running the latest version of fastsimcoal25 (version 2.5.2.8) on Windows 10, 64-bit machine.

My .par file is:

//Number of population samples (demes)
2
//Population effective sizes (number of genes)
1300
82000
//Samples sizes and samples age
300
300
//Growth rates : negative growth implies population expansion
0
0
//Number of migration matrices : 0 implies no migration between demes
0
//historical event: time, source, sink, migrants, new deme size, new growth rate, migration matrix index
1 historical event
2000 1 0 0.35 1 0 0
//Number of independent loci [chromosome]
1 0
//Per chromosome: Number of contiguous linkage Block: a block is a set of contiguous loci
1
//per Block:data type, number of loci, per generation recombination and mutation rates and optional parameters
DNA 402 0 0.000004123 0.8967

The code used for running the program:
fsc252 -i Grp17_Admix_Simulation.par -n 100 -S

Thanks,
Bekah

[Laurent Excoffier]

Hi

You indeed have a problem with

1 historical event
2000 1 0 0.35 1 0 0

With this historical event you are sending 35% of the lineages of deme 1 into deme 0 but then these lineages will remain isolated for ever hence preventing convergence of the coalescent process

Laurent