Testing for simple population branching - No output produced
Maria
Dear users,
I am having problems running a simple analysis perhaps someone can help here.
I am using fsc26 and SFS stats that I obtained with easySFS. I have 4 populations and want depict a simple branching topology with no migration. This is my command:
/cluster/home/s_menb/software/fsc26_linux64/fsc26 -t model_1_braya.tpl -e model_1_braya.est -E 1 -m -0 -C 10 -n 10000 -L 40 -s 0 -M --multiSFS -x -d -c 64 -B 64
See below the scenario topology, .est and .tpl file.
model_1_braya.est
// Search ranges and rules file
// ****************************
//$NPOP0$
//$NPOP1$
//$NPOP2$
//$NPOP3$
[PARAMETERS]
//#isInt? #name #dist.#min #max
//all Ns are in number of haploid individuals
1 $TDIV$ unif 300000 400000 output
1 $TIME1$ unif 20000 40000 output
1 $TIME2$ unif 4000 3000 output
1 $NANC$ unif 10 1e7 output
1 $NPOP0$ unif 10 1e7 output
1 $NPOP1$ unif 10 1e7 output
1 $NPOP2$ unif 10 1e7 output
1 $NPOP3$ unif 10 1e7 output
[RULES]
1 $TDIV$ > $TIME1$
1 $TIME1$ > $TIME2$
model_1_braya.tpl
//Parameters for the coalescence simulation program : simcoal.exe
4 samples to simulate :
//Population effective sizes (number of genes)
$NPOP0$
$NPOP1$
$NPOP2$
$NPOP3$
//Samples sizes and samples age
35
14
27
14
//Growth rates: negative growth implies population expansion
0
0
0
0
//Number of migration matrices : 0 implies no migration between demes
0
//historical event: time, source, sink, migrants, new deme size, growth rate, migr mat index
3 historical event
$TIME2$ 1 2 0 1 0 0
$TIME1$ 3 1 0 1 0 0
$TDIV$ 0 3 0 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 gen recomb and mut rates
SNP 1448 0 5.77e-9
The program has been running for 8 days. It has written the sim and parm files but no output produced. For such a small dataset and based on previous experience with fsc26 I assume the programme is just hanging and stuck.
Does anyone have a suggestion where to look for a mistake possibly in the files?
Best regards,
Maria