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
|
#!/usr/bin/python
# Author: Francois-Jose Serra
# Creation Date: 2010/04/26 17:17:06
from ete3.evol import EvolTree
import sys, re
typ = 'S'
#while typ != 'L' and typ != 'S':
# typ = raw_input (\
# "choose kind of example [L]ong or [S]hort, hit [L] or [S]:\n")
TREE_PATH = "data/S_example/measuring_%s_tree.nw" % (typ)
ALG_PATH = "data/S_example/alignment_%s_measuring_evol.fasta" % (typ)
WORKING_PATH = "data/S_example/paml/"
MY_PATH = ''
TREE_PATH = MY_PATH + re.sub('\./', '', TREE_PATH)
ALG_PATH = MY_PATH + re.sub('\./', '', ALG_PATH )
T = EvolTree (TREE_PATH)
T.link_to_alignment (ALG_PATH)
T.workdir = (WORKING_PATH)
T.link_to_evol_model(T.workdir + '/fb/out','fb')
T.link_to_evol_model(T.workdir + '/M1/out','M1')
T.link_to_evol_model(T.workdir + '/M2/out','M2')
T.show(histfaces=['M2'])
sys.stderr.write('\n\nThe End.\n\n')
|
