Hi !
If just run the simplest model (no gene flow), there won't be a problem.
But with the addition of the migration matrix, there is an error:
Error: deme 1 going extinct while still having remaining lineages
I've looked at some issues with the same error;
1. obs files can be found correctly by fsc27
2. The RESIZE parameter seems to have been simplified in newer versions of fsc.
So I didn't write the relevant parameters.
(Iike: 0 RESIZE3 = ANCPOPOSIZE/MIDPOPSIZE hide; If the first number is set to 1, I may also encounter this error)
3. I also wrote an additional migration matrix
my model is:
If I remove the migration matrix, there is no error.
my tpl file is:
//Parameters for the coalescence simulation program : fastsimcoal.exe
2 samples to stimulate
//Population effective sizes (number of genes)
CAPOPSIZE
JUPOPSIZE
//Haploid samples sizes
21
23
//Growth rates: negative growth implies population expansion
0
0
//Number of migration matrices : 0 implies no migration between demes
2
//Migration matrix 0
0 M0
M1 0
//Migration matrix 1
0 0
0 0
//historical event: time, source, sink, migrants, new deme size, new growth rate, migration matrix index
1 historical event
TDIV 0 1 1 ANCPOPOSIZE 0 1 absoluteResize
//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
FREQ 1 0 3.03e-9 OUTEXP
my est file is:
// Priors and rules file
// *********************
[PARAMETERS]
//#isInt? #name #dist.#min #max
//all Ns are in number of haploid individuals
1 CAPOPSIZE logunif 1e2 1e8 output
1 JUPOPSIZE logunif 1e2 1e8 output
1 ANCPOPSIZE logunif 1e2 1e8 output
1 TDIV unif 1 1340000 output bounded
0 M0 logunif 1e-7 1 output
0 M1 logunif 1e-7 1 output
[COMPLEX PARAMETERS]
in my tpl file, I wote:
TDIV 0 1 1 ANCPOPOSIZE 0 1 absoluteResize
It seems to be undefined where deme1 comes from and goes to?
So I changed:
"
TDIV 0 1 1 ANCPOPOSIZE 0 1 absoluteResize
TDIV 1 1 1 ANCPOPOSIZE 0 1 absoluteResize
"
It can not work and met : "Segmentation fault (core dumped)"