-
Notifications
You must be signed in to change notification settings - Fork 2
/
draw.sh
75 lines (63 loc) · 1.32 KB
/
draw.sh
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
#!/bin/bash
TREE=Mbasal_mod1.bt.rescaled.nex
ASSOC=A,B,G,M
LOG=run2/HostFungusAssociations.txt.unconstrained.log
STATES=states.tsv
DATA=HostFungusAssociations.txt.tsv
PHYLA=phyla.tsv
WIDTH=3000
HEIGHT=4000
BURNIN=500000
# HERE WE MAKE THE FIGURE FOR IN THE MANUSCRIPT
# to make a phylogram, set this to 'phylo'
#MODE=clado
MODE=phylo
# to change node radius, e.g. to hide pies, set to 1
#RADIUS=12
RADIUS=0
# to hide pies, set to 'no'
#PIES=yes
PIES=no
# to hide tips, set to 'no'
#TIPS=yes
TIPS=no
perl ../../script/visualize.pl \
-verbose \
-assoc $ASSOC \
-tree $TREE \
-states $STATES \
-log $LOG \
-data $DATA \
-width $WIDTH \
-height $HEIGHT \
-mode $MODE \
-taxa $PHYLA \
-config radius=$RADIUS \
-config pies=$PIES \
-config tips=$TIPS \
-config width=6 \
-burnin $BURNIN > $TREE.ms.svg
# HERE WE MAKE THE SUPPLEMENTARY FIGURE
# to make a phylogram, set this to 'phylo'
MODE=clado
# to change node radius, e.g. to hide pies, set to 1
RADIUS=12
# to hide pies, set to 'no'
PIES=yes
# to hide tips, set to 'no'
TIPS=yes
perl ../../script/visualize.pl \
-verbose \
-assoc $ASSOC \
-tree $TREE \
-states $STATES \
-log $LOG \
-data $DATA \
-width $WIDTH \
-height $HEIGHT \
-mode $MODE \
-config radius=$RADIUS \
-config pies=$PIES \
-config tips=$TIPS \
-config width=2 \
-burnin $BURNIN > $TREE.supp.svg