summaryrefslogtreecommitdiff
path: root/examples/evol/5_branchsite_cladetest.py
diff options
context:
space:
mode:
Diffstat (limited to 'examples/evol/5_branchsite_cladetest.py')
-rw-r--r--examples/evol/5_branchsite_cladetest.py6
1 files changed, 3 insertions, 3 deletions
diff --git a/examples/evol/5_branchsite_cladetest.py b/examples/evol/5_branchsite_cladetest.py
index faa2d04..bf3480c 100644
--- a/examples/evol/5_branchsite_cladetest.py
+++ b/examples/evol/5_branchsite_cladetest.py
@@ -26,12 +26,12 @@ tree.link_to_alignment ('data/S_example/alignment_S_measuring_evol.fasta')
print (tree)
print ('Tree and alignment loaded.')
-input ('Tree will be mark in order to contrast Gorilla and Chimpanzee as foreground \nspecies.')
+eval(input ('Tree will be mark in order to contrast Gorilla and Chimpanzee as foreground \nspecies.'))
marks = ['1', 3, '7']
tree.mark_tree (marks, ['#1'] * 3)
-print (tree.write ())
+print(tree.write ())
# display marked branches in orange
for node in tree.traverse ():
@@ -65,7 +65,7 @@ tree.run_model ('M1')
print ('''p-value that, in marked clade, we have one class of site
specifically evolving at a different rate:''')
-print (tree.get_most_likely ('bsC.137', 'M1'))
+print(tree.get_most_likely ('bsC.137', 'M1'))
#print ('p-value representing significance that omega is different of 1:')
#print (tree.get_most_likely ('bsD.137', 'M3'))