How to write about historical events in .tpl

118 views
Skip to first unread message

yulu ye

unread,
Aug 27, 2024, 10:31:39 PM8/27/24
to fastsimcoal2
Hi everyone, I want to describe a topology like this and I don't know if this is the right way to write the .tpl

屏幕截图 2024-08-28 103039.jpg


the tpl file
//Number of population samples
5
//Population effective sizes (number of genes)
NPOP_0
NPOP_1
NPOP_2
NPOP_3
NPOP_4
//Sample sizes
20
20
20
20
20
//Growth rates : negative growth implies population expansion
0
0
0
0
0
//Number of migration matrices : 0 implies no migration between demes
0
//historical event: time, source, sink, migrants, new size, new growth rate, migr. matrix
4 historical events
TDIV_34 3 4 1 NANC_CHG_34 0 0
TDIV_24 2 4 1 NANC_CHG_24 0 0
TDIV_14 1 4 1 NANC_CHG_14 0 0
TDIV_40 4 0 1 NANC_CHG_40 0 0
//Number of independent loci [chromosome]
1 0
//Per chromosome: Number of linkage blocks
1
//per Block: data type, num loci, rec. rate and mut rate + optional parameters
FREQ 1 0 1e-8 OUTEXP

the est file
// Search ranges and rules file
// ****************************

[PARAMETERS]
//#isInt? #name   #dist.#min  #max
//all Ns are in number of haploid individuals
1    NPOP_0     unif    1e4  1e7     output
1    NPOP_1     unif    1e4  1e7     output
1    NPOP_2     unif    1e4  1e7     output
1    NPOP_3     unif    1e4  1e7     output
1    NPOP_4     unif    1e4  1e7     output
1    TDIV_34      unif    1e4  1e6     output    bounded
1    TDIV_24      unif    TDIV_34  1e6     output   paramInRange
1    TDIV_14      unif    TDIV_24  2.5e6     output paramInRange
1    TDIV_40      unif    TDIV_14  2.5e6     output paramInRange
0    NANC_CHG_34    unif    0.1  5       hide   bounded
0    NANC_CHG_24    unif    0.1  5       hide   bounded
0    NANC_CHG_14    unif    0.1  5       hide   bounded
0    NANC_CHG_40    unif    0.1  5       hide   bounded

[COMPLEX PARAMETERS]
1    NANC_34 = NPOP_4*NANC_CHG_34    output
1    NANC_24 = NANC_34*NANC_CHG_24    output
1    NANC_14 = NANC_24*NANC_CHG_14    output
1    NANC_40 = NANC_14*NANC_CHG_40    output
Reply all
Reply to author
Forward
0 new messages