MRCA not found

[Xiurong Zhao]

Hi Laurent,

I've noticed in some of my runs to estimate parameters of demographic model I would receive the following error message and a core dump:

Bad parameters leading to non-convergent coalescent are written in file bad.par

Simulations did not converge. MRCA not found! (2 lineages remaining)

Looking in the subdirectory, I would see that output.brent_lhoods have hundreds of entries, suggesting that the parameter estimation did proceed for awhile. But there would be a badParams.txt output, recording the parameter values of the last iteration.  Do you have any insight on why might these runs fail?

Here is my .tpl file

//Parameters for the coalescence simulation program : fastsimcoal.exe

4 samples to simulate :

//Population effective sizes (number of genes)

NPOP1

NPOP2

NPOP3

NPOP4

//Samples sizes and samples age 

132

53

12

21

//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, new growth rate, migration matrix index

3 historical event

TDIV1 3 1 1 RESIZE1 0 0

TDIV2 2 0 1 RESIZE2 0 0

TDIV3 1 0 1 RESIZE3 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

FREQ 1 0 3.13e-7 

########################

.est file

// Priors and rules file

// *********************

[PARAMETERS]

//#isInt? #name #dist.#min #max

//all N are in number of haploid individuals

1 NPOP1 unif 10000 1000000 output

1 NPOP2 unif 10000 1000000 output

1 NPOP3 unif 10000 1000000 output

1 NPOP4 unif 10000 1000000 output

1 TDIV1 unif 6000 12000 output

1 TDIV2 unif 6000 12000 output

1 TDIV3 unif 6000 12000 output

1 NANC1 unif 10000 1000000 output

1 NANC2 unif 10000 1000000 output

1 NANC3 unif 10000 1000000 output

[RULES]

[COMPLEX PARAMETERS]

0       RESIZE1 = NANC1/NPOP2 output

0       RESIZE2 = NANC2/NPOP1 output

0       RESIZE3 = NANC3/NANC2 output

#####

and here is my command line :

fsc2./fsc2709 -t B_hap_4POP2.tpl -n10 -m -e B_hap_4POP2.est -M -L 40 -c6 -q -x

Thank you very much for your kind attention,

Warm greetings,

Xiurong Zhao

[Lu Liu]

Hi, Xiurong. I have the same problem. Do you resolve it?

Sincerely,

Lu Liu

[Xiurong Zhao]

I sovle this problem when I set  that TDIV3 was larger than TDIV1 and TDIV2