Ancient hybrid speciation model

[张德志]

Dear Laurent and all users,

I am using fsc26 to estimate an ancient hybrid speciation model for 4 species (attached figure). My MAF data was generated by easySFS.py by using about 500000 SNPs. I used the following command: 

./fsc26 -t 4species_nomissing.tpl -e 4species_nomissing.est -m -n 500000 -L 50 -M -c 30 -B 30

fsc26 is runing, but also got message like this:

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

Num. steps in coalescent process: 1000001

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

Num. steps in coalescent process: 1000001

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

Num. steps in coalescent process: 1000001

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

Brent minimization result: Tdiver = 1212947.0000000 (lhood = -3019629.3646568)

Here is my tpl and est file:

###TPL FILE//Parameters for the coalescence simulation program : fastsimcoal.exe

4 samples to simulate :

//Population effective sizes (number of genes)

NPOP0

NPOP1

NPOP2

NPOP3

//Samples sizes and samples age 

10

10

10

7

//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

4 historical event

T10 1 0 1 RESIZE10 0 0

Tadmix 2 0 pror 1 0 0

Tadmix 2 3 1 1 0 0

Tdiver 0 3 1 RESIZE03 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 generation recombination and mutation rates and optional parameters

FREQ 1 0 5.61e-9

###EST FILE

// Priors and rules file

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

[PARAMETERS]

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

//all Ns are in number of haploid individuals

1  NPOP0     unif     10  100000000   output

1  NPOP1        unif     10  100000000     output

1  NPOP2      unif     10  100000000   output

1  NPOP3     unif     10  100000000      output

1  NPOPA10     unif     10  100000000      hide

1  NPOPA03     unif     10  100000000      output

1  T10        unif     10000   500000   output

1  Tadmix        unif     10000   1500000   output

1  Tdiver        unif     10000   1500000   output

0  pror        unif     0   1   output    bounded

[RULES]

[COMPLEX PARAMETERS]

0  RESIZE10   = NPOPA10/NPOP0      hide

0  RESIZE03   = NPOPA03/NPOP3      hide

Any suggestions will help a lot,

My best wishes,

Dezhi