Error: deme 1 going extinct while still having remaining lineages

327 views
Skip to first unread message

Michael Larson

unread,
Dec 20, 2023, 12:39:20 AM12/20/23
to fastsimcoal2
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:
ae8c1397655658b855ad5d5a0390e13.png
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)"

Laurent Excoffier

unread,
Dec 20, 2023, 3:35:06 AM12/20/23
to fastsimcoal2
 Hi

From what I can see, ANPOPOSIZE is not defined in your est file, and fsc reads it as zero, hence the error message.

It’s always good to check the par file generated by fsc to see if there are undefined keywords remaining 

Hope it helps

Laurent 

Michael Larson

unread,
Dec 20, 2023, 4:24:53 AM12/20/23
to fastsimcoal2
Thanks for your reply. 
I was careless. I have noticed It was just a typo. 
Im so sorry
Reply all
Reply to author
Forward
0 new messages