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
|
seed used = 1
BASEML (in paml version 4.3, August 2009) ../Tests/PAML/alignment.phylip REV dGamma (ncatG=5)
Frequencies..
T C A G
Homo_sapie 0.2027 0.1622 0.3694 0.2658
Pan_troglo 0.1982 0.1667 0.3694 0.2658
Gorilla_go 0.2072 0.1577 0.3694 0.2658
Pongo_pygm 0.2027 0.1622 0.3694 0.2658
Macaca_mul 0.1937 0.1667 0.3739 0.2658
Homogeneity statistic: X2 = 0.00092 G = 0.00092
Average 0.200901 0.163063 0.370270 0.265766
# constant sites: 218 (98.20%)
ln Lmax (unconstrained) = -316.049385
Distances: TN93 (kappa) (alpha set at 0.50)
This matrix is not used in later m.l. analysis.
Homo_sapie
Pan_troglo 0.0047(999.0000)
Gorilla_go 0.0047(999.0000) 0.0097(999.0000)
Pongo_pygm 0.0000(999.0000) 0.0047(999.0000) 0.0047(999.0000)
Macaca_mul 0.0093( 3.6925) 0.0145( 7.6730) 0.0145( 7.7118) 0.0093( 3.6925)
TREE # 1: (((1, 2), 3), 4, 5); MP score: 4.00
lnL(ntime: 7 np: 13): -319.774326 +0.000000
6..7 7..8 8..1 8..2 7..3 6..4 6..5
0.000004 0.000004 0.000004 0.004533 0.004529 0.000004 0.009174 999.000000 160.593137 0.000010 0.000010 0.000010 200.111998
SEs for parameters:
-1.00000 -1.00000 -1.00000 801727.63247 730462.67590 -1.00000 1610106.73811 -1.00000 -1.00000 3749008252.60285 16592006264.44326 4678796535.94680 212206538264286.93750
tree length = 0.01825
(((Homo_sapie, Pan_troglo), Gorilla_go), Pongo_pygm, Macaca_mul);
(((Homo_sapie: 0.000004, Pan_troglo: 0.004533): 0.000004, Gorilla_go: 0.004529): 0.000004, Pongo_pygm: 0.000004, Macaca_mul: 0.009174);
Detailed output identifying parameters
Parameters in the rate matrix (REV) (Yang 1994 J Mol Evol 39:105-111):
Rate parameters: 999.00000 160.59314 0.00001 0.00001 0.00001
Base frequencies: 0.20090 0.16306 0.37027 0.26577
Rate matrix Q, Average Ts/Tv = 2.7478
-2.483319 1.819246 0.664073 0.000000
2.241391 -2.241391 0.000000 0.000000
0.360312 0.000000 -0.363280 0.002968
0.000000 0.000000 0.004135 -0.004135
alpha (gamma, K=5) = 200.11200
rate: 0.90307 0.96129 0.99837 1.03642 1.10086
freq: 0.20000 0.20000 0.20000 0.20000 0.20000
|
