1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138 | 100 1 4 2
1 1
m s
0
1 0 100
1 0 100
1 0 100
1 0 100
1 0 100
1 0 100
1 0 100
3 0 1000
1 0 0.5
1 1 1
0
4 1 10 1 1
3 0 10
2 1 0.1 1 1
1 -2 0 0 0
1 -2 0 0 0
0
0
0
----------------------------------------------------------------------------------------
PopABC - Mark Beaumont & Joao Lopes 01/05/09
>no_iterations, generation_time, no_populations, no_loci
>escalar per locus (autosome - 1; X-linked - 0.75; Y-linked or mitDNA - 0.25)
>type of DNA data (s - sequence; m - microssatelites)
>topology: 0 - uniform distribution;
1 - choose topology from a list;
2 - specify topology manually [e.g. ((Pop1,Pop2)Pop3) -> 1 2 2 3];
3 - uniform distribution (and choose a Model marker);
4 - choose topology from a list (and choose a Model marker);
5 - specify topology manually (and choose a Model marker).
>ne1 params: 1 - uniform distribtuion;
2 - generalized gamma distribution;
>ne2 params
>ne3 params
>ne4 params
>neanc1 params
>neanc2 params
>neanc3 params
>t1 params: 1 - uniform distribtuion;
2 - generalized gamma distribution;
3 - uniform distribtuion (for all time events);
4 - generalized gamma distribution (for all time events);
[for 1 and 2 t(n) is added to t(n+1)]
[for 3 and 4 set only one priors for all t(n)]
>t2 params
>t3 params
>mig1 params: 0 - no migration;
1 - uniform distribtuion;
2 - generalized gamma distribution;
3 - uniform distribtuion (on number of migrations);
4 - generalized gamma distribution (on number of migrations).
[for 3 and 4 real mig rate is calculated as nmig/Ne]
>mig2 params
>mig3 params
>mig4 params
>miganc1 params
>miganc2 params
>mutM params: 0 - zero mutation;
1 - lognormal distribution: (mean of mean(log10); stdev of mean(log10);
mean of Sdev(log10); stdev of stdev(log10). Stdev truncated at 0.
2 - normal distribution: (mean of mean; stdev of mean; mean of Sdev;
stdev of stdev. Stdev truncated at 0.
>mutS params:
>recM params: 0 - zero recombination;
1 - lognormal distribution: (mean of mean(log10); stdev of mean(log10);
mean of Sdev(log10); stdev of stdev(log10). Stdev truncated at 0.
2 - normal distribution: (mean of mean; stdev of mean; mean of Sdev;
stdev of stdev. Stdev truncated at 0.
>recS params
>migweight: 0 - do not use migweights matrix;
1 - use migweights matrix as following:
0 mw112 mw113 mw114
0 mw122 mw123 mw124
0 mw132 mw133 mw134
mw211 0 mw213 mw214
mw221 0 mw223 mw224
mw231 0 mw233 mw234
mw311 mw312 0 mw314
mw321 mw322 0 mw324
mw331 mw332 0 mw334
mw411 mw412 mw413 0
mw421 mw422 mw423 0
mw431 mw432 mw433 0
, where mwitj is the prob that the fraction of migrantes in pop i comes
from pop j at a period of time before time event t. Sum of prob should
be equal to 1.
[only use migweight if the topology is specified (option 1,2,4 or 5)]
----------------------------------------------------------------------------------------
Tree topology:
|| Neanc3 Neanc3
|| | |
t3|| ------------- ----------------
|| | | | |
|| Neanc2 | Neanc2 |
|| | | | |
t2|| ----------- | OR -------- |
|| | | | | | |
|| Neanc1 | | | | Neanc1
|| | | | | | |
t1|| -------- | | | | --------
|| | | | | | | | |
\/ Ne Ne Ne Ne Ne Ne Ne Ne
|