Error when using paramInRange

13 views
Skip to first unread message

Vincent Kendlbacher

unread,
Jul 30, 2025, 7:51:56 AMJul 30
to fastsimcoal2
Dear all, 
I am having a very hard time getting my tpl & est files to run.
The tpl. file is the following: 
-----------------------------------------------
//Parameters for the coalescence simulation program : fsimcoal2.exe
4 samples to simulate :
//Population effective sizes (number of genes): P=Peru, H=Hig, Q=Queule, A=Arg
NP
NH
NQ
NA
//Samples sizes and samples age
10
60
28
50
//Growth rates : negative growth implies population expansion
0
0
0
0
//Number of migration matrices : 0 implies no migration between demes
3
//current Migration matrix 0
0.0000 MIGHP 0.0000 0.0000
MIGHP 0.0000 MIGQH 0.0000
0.0000 MIGQH 0.0000 MIGAQ
0.0000 0.0000 MIGAQ 0.0000
//current Migration matrix 1
0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 MIGAQ
0.0000 0.0000 MIGAQ 0.0000
//current Migration matrix 2
0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000
//historical event: time, source, sink, migrants, new deme size, new growth rate, migration matrix index
4   historical event
TDIVSEC 0 0 0 1 0 0
TDIVHIG 1 2 1 1 0 1
TDIVQUE 2 3 1 1 0 2
TDIVOCEAN 3 0 1 1 0 2
//Number of independent loci [chromosome]
1 0
//Per chromosome: Number of linkage Blocks
1
//Per Block:data type, number of loci, per gen recomb and mut rates and optional parameters
FREQ 1 0 7e-10 OUTEXP
-------------------------------------------------------------
and the est file the following: 
-------------------------------------------------------------
// Priors and rules file
// *********************

[PARAMETERS]
//#isInt? #name #dist.#min #max
//all N are in number of diploid individuals
0  MIGHP       unif     0.001    0.1       output
0  MIGQH       unif     0.001    0.1       output
0  MIGAQ       unif     0.001    0.1       output
1  NP          logunif  100      1000000   output
1  NH          logunif  100      1000000   output
1  NQ          logunif  100      1000000   output
1  NA          logunif  100      1000000   output
1  TDIVSEC     logunif  100      5500      output bounded
1  TDIVHIG     logunif  TDIVSEC  5500      output paramInRange
1  TDIVQUE     logunif  TDIVHIG  5500      output paramInRange
1  TDIVOCEAN   logunif  TDIVQUE  7000      output paramInRange
------------------------------------------------------------------
When submitting to my hpc cluster, I always get a "(core dumped") segmantation fault error. However, I dont get this error when specifying normal numbers as "min" & "max" for my "TDIVs" -> Hence, this has to have something to do with you I am using the paramInRange command...
Can anyone help?
I am already quite desperate...
Best,
Vincent
Reply all
Reply to author
Forward
0 new messages