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