-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathquebrataxontree.m
282 lines (278 loc) · 15 KB
/
quebrataxontree.m
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
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
tree.all.NJ(1)=Phyqualyread('all100NJe.tree');
tree.all.NJ(2)=Phyqualyread('all200NJe.tree');
tree.all.NJ(3)=Phyqualyread('all300NJe.tree');
tree.all.NJ(4)=Phyqualyread('all400NJe.tree');
tree.all.NJ(5)=Phyqualyread('all500NJe.tree');
tree.all.NJ(6)=Phyqualyread('all600NJe.tree');
tree.all.NJ(7)=Phyqualyread('all700NJe.tree');
tree.all.NJ(8)=Phyqualyread('all800NJe.tree');
tree.all.NJ(9)=Phyqualyread('all900NJe.tree');
tree.all.NJ(10)=Phyqualyread('all1000NJe.tree');
tree.all.NJ(11)=Phyqualyread('all1100NJe.tree');
tree.all.NJ(12)=Phyqualyread('all1200NJe.tree');
tree.all.NJ(13)=Phyqualyread('all1300NJe.tree');
tree.all.NJ(14)=Phyqualyread('all1400NJe.tree');
tree.all.NJ(15)=Phyqualyread('all1500NJe.tree');
tree.all.NJ(16)=Phyqualyread('all1600NJe.tree');
tree.all.NJ(17)=Phyqualyread('all1700NJe.tree');
tree.all.NJ(18)=Phyqualyread('all1800NJe.tree');
tree.all.NJ(19)=Phyqualyread('all1900NJe.tree');
tree.all.NJ(20)=Phyqualyread('all2000NJe.tree');
tree.all.NJ(21)=Phyqualyread('all2100NJe.tree');
tree.all.NJ(22)=Phyqualyread('all2200NJe.tree');
tree.all.NJ(23)=Phyqualyread('all2300NJe.tree');
tree.all.NJ(24)=Phyqualyread('all2400NJe.tree');
tree.all.NJ(25)=Phyqualyread('all2500NJe.tree');
tree.all.NJ(26)=Phyqualyread('all2600NJe.tree');
tree.all.NJ(27)=Phyqualyread('all2700NJe.tree');
tree.all.NJ(28)=Phyqualyread('all2800NJe.tree');
tree.all.NJ(29)=Phyqualyread('all2900NJe.tree');
tree.all.NJ(30)=Phyqualyread('all3000NJe.tree');
tree.all.Ward(1)=Phyqualyread('all100Ward.tree');
tree.all.Ward(2)=Phyqualyread('all200Ward.tree');
tree.all.Ward(3)=Phyqualyread('all300Ward.tree');
tree.all.Ward(4)=Phyqualyread('all400Ward.tree');
tree.all.Ward(5)=Phyqualyread('all500Ward.tree');
tree.all.Ward(6)=Phyqualyread('all600Ward.tree');
tree.all.Ward(7)=Phyqualyread('all700Ward.tree');
tree.all.Ward(8)=Phyqualyread('all800Ward.tree');
tree.all.Ward(9)=Phyqualyread('all900Ward.tree');
tree.all.Ward(10)=Phyqualyread('all1000Ward.tree');
tree.all.Ward(11)=Phyqualyread('all1100Ward.tree');
tree.all.Ward(12)=Phyqualyread('all1200Ward.tree');
tree.all.Ward(13)=Phyqualyread('all1300Ward.tree');
tree.all.Ward(14)=Phyqualyread('all1400Ward.tree');
tree.all.Ward(15)=Phyqualyread('all1500Ward.tree');
tree.all.Ward(16)=Phyqualyread('all1600Ward.tree');
tree.all.Ward(17)=Phyqualyread('all1700Ward.tree');
tree.all.Ward(18)=Phyqualyread('all1800Ward.tree');
tree.all.Ward(19)=Phyqualyread('all1900Ward.tree');
tree.all.Ward(20)=Phyqualyread('all2000Ward.tree');
tree.all.Ward(21)=Phyqualyread('all2100Ward.tree');
tree.all.Ward(22)=Phyqualyread('all2200Ward.tree');
tree.all.Ward(23)=Phyqualyread('all2300Ward.tree');
tree.all.Ward(24)=Phyqualyread('all2400Ward.tree');
tree.all.Ward(25)=Phyqualyread('all2500Ward.tree');
tree.all.Ward(26)=Phyqualyread('all2600Ward.tree');
tree.all.Ward(27)=Phyqualyread('all2700Ward.tree');
tree.all.Ward(28)=Phyqualyread('all2800Ward.tree');
tree.all.Ward(29)=Phyqualyread('all2900Ward.tree');
tree.all.Ward(30)=Phyqualyread('all3000Ward.tree');
tree.all.UPGMA(1)=Phyqualyread('all100UPGMA.tree');
tree.all.UPGMA(2)=Phyqualyread('all200UPGMA.tree');
tree.all.UPGMA(3)=Phyqualyread('all300UPGMA.tree');
tree.all.UPGMA(4)=Phyqualyread('all400UPGMA.tree');
tree.all.UPGMA(5)=Phyqualyread('all500UPGMA.tree');
tree.all.UPGMA(6)=Phyqualyread('all600UPGMA.tree');
tree.all.UPGMA(7)=Phyqualyread('all700UPGMA.tree');
tree.all.UPGMA(8)=Phyqualyread('all800UPGMA.tree');
tree.all.UPGMA(9)=Phyqualyread('all900UPGMA.tree');
tree.all.UPGMA(10)=Phyqualyread('all1000UPGMA.tree');
tree.all.UPGMA(11)=Phyqualyread('all1100UPGMA.tree');
tree.all.UPGMA(12)=Phyqualyread('all1200UPGMA.tree');
tree.all.UPGMA(13)=Phyqualyread('all1300UPGMA.tree');
tree.all.UPGMA(14)=Phyqualyread('all1400UPGMA.tree');
tree.all.UPGMA(15)=Phyqualyread('all1500UPGMA.tree');
tree.all.UPGMA(16)=Phyqualyread('all1600UPGMA.tree');
tree.all.UPGMA(17)=Phyqualyread('all1700UPGMA.tree');
tree.all.UPGMA(18)=Phyqualyread('all1800UPGMA.tree');
tree.all.UPGMA(19)=Phyqualyread('all1900UPGMA.tree');
tree.all.UPGMA(20)=Phyqualyread('all2000UPGMA.tree');
tree.all.UPGMA(21)=Phyqualyread('all2100UPGMA.tree');
tree.all.UPGMA(22)=Phyqualyread('all2200UPGMA.tree');
tree.all.UPGMA(23)=Phyqualyread('all2300UPGMA.tree');
tree.all.UPGMA(24)=Phyqualyread('all2400UPGMA.tree');
tree.all.UPGMA(25)=Phyqualyread('all2500UPGMA.tree');
tree.all.UPGMA(26)=Phyqualyread('all2600UPGMA.tree');
tree.all.UPGMA(27)=Phyqualyread('all2700UPGMA.tree');
tree.all.UPGMA(28)=Phyqualyread('all2800UPGMA.tree');
tree.all.UPGMA(29)=Phyqualyread('all2900UPGMA.tree');
tree.all.UPGMA(30)=Phyqualyread('all3000UPGMA.tree');
tree.cloro.NJ(1)=Phyqualyread('cloro100NJe.tree');
tree.cloro.NJ(2)=Phyqualyread('cloro200NJe.tree');
tree.cloro.NJ(3)=Phyqualyread('cloro300NJe.tree');
tree.cloro.NJ(4)=Phyqualyread('cloro400NJe.tree');
tree.cloro.NJ(5)=Phyqualyread('cloro500NJe.tree');
tree.cloro.NJ(6)=Phyqualyread('cloro600NJe.tree');
tree.cloro.NJ(7)=Phyqualyread('cloro700NJe.tree');
tree.cloro.NJ(8)=Phyqualyread('cloro800NJe.tree');
tree.cloro.NJ(9)=Phyqualyread('cloro900NJe.tree');
tree.cloro.NJ(10)=Phyqualyread('cloro1000NJe.tree');
tree.cloro.NJ(11)=Phyqualyread('cloro1100NJe.tree');
tree.cloro.NJ(12)=Phyqualyread('cloro1200NJe.tree');
tree.cloro.NJ(13)=Phyqualyread('cloro1300NJe.tree');
tree.cloro.NJ(14)=Phyqualyread('cloro1400NJe.tree');
tree.cloro.NJ(15)=Phyqualyread('cloro1500NJe.tree');
tree.cloro.NJ(16)=Phyqualyread('cloro1600NJe.tree');
tree.cloro.NJ(17)=Phyqualyread('cloro1700NJe.tree');
tree.cloro.NJ(18)=Phyqualyread('cloro1800NJe.tree');
tree.cloro.NJ(19)=Phyqualyread('cloro1900NJe.tree');
tree.cloro.NJ(20)=Phyqualyread('cloro2000NJe.tree');
tree.cloro.NJ(21)=Phyqualyread('cloro2100NJe.tree');
tree.cloro.NJ(22)=Phyqualyread('cloro2200NJe.tree');
tree.cloro.NJ(23)=Phyqualyread('cloro2300NJe.tree');
tree.cloro.NJ(24)=Phyqualyread('cloro2400NJe.tree');
tree.cloro.NJ(25)=Phyqualyread('cloro2500NJe.tree');
tree.cloro.NJ(26)=Phyqualyread('cloro2600NJe.tree');
tree.cloro.NJ(27)=Phyqualyread('cloro2700NJe.tree');
tree.cloro.NJ(28)=Phyqualyread('cloro2800NJe.tree');
tree.cloro.NJ(29)=Phyqualyread('cloro2900NJe.tree');
tree.cloro.NJ(30)=Phyqualyread('cloro3000NJe.tree');
tree.cloro.Ward(1)=Phyqualyread('cloro100Ward.tree');
tree.cloro.Ward(2)=Phyqualyread('cloro200Ward.tree');
tree.cloro.Ward(3)=Phyqualyread('cloro300Ward.tree');
tree.cloro.Ward(4)=Phyqualyread('cloro400Ward.tree');
tree.cloro.Ward(5)=Phyqualyread('cloro500Ward.tree');
tree.cloro.Ward(6)=Phyqualyread('cloro600Ward.tree');
tree.cloro.Ward(7)=Phyqualyread('cloro700Ward.tree');
tree.cloro.Ward(8)=Phyqualyread('cloro800Ward.tree');
tree.cloro.Ward(9)=Phyqualyread('cloro900Ward.tree');
tree.cloro.Ward(10)=Phyqualyread('cloro1000Ward.tree');
tree.cloro.Ward(11)=Phyqualyread('cloro1100Ward.tree');
tree.cloro.Ward(12)=Phyqualyread('cloro1200Ward.tree');
tree.cloro.Ward(13)=Phyqualyread('cloro1300Ward.tree');
tree.cloro.Ward(14)=Phyqualyread('cloro1400Ward.tree');
tree.cloro.Ward(15)=Phyqualyread('cloro1500Ward.tree');
tree.cloro.Ward(16)=Phyqualyread('cloro1600Ward.tree');
tree.cloro.Ward(17)=Phyqualyread('cloro1700Ward.tree');
tree.cloro.Ward(18)=Phyqualyread('cloro1800Ward.tree');
tree.cloro.Ward(19)=Phyqualyread('cloro1900Ward.tree');
tree.cloro.Ward(20)=Phyqualyread('cloro2000Ward.tree');
tree.cloro.Ward(21)=Phyqualyread('cloro2100Ward.tree');
tree.cloro.Ward(22)=Phyqualyread('cloro2200Ward.tree');
tree.cloro.Ward(23)=Phyqualyread('cloro2300Ward.tree');
tree.cloro.Ward(24)=Phyqualyread('cloro2400Ward.tree');
tree.cloro.Ward(25)=Phyqualyread('cloro2500Ward.tree');
tree.cloro.Ward(26)=Phyqualyread('cloro2600Ward.tree');
tree.cloro.Ward(27)=Phyqualyread('cloro2700Ward.tree');
tree.cloro.Ward(28)=Phyqualyread('cloro2800Ward.tree');
tree.cloro.Ward(29)=Phyqualyread('cloro2900Ward.tree');
tree.cloro.Ward(30)=Phyqualyread('cloro3000Ward.tree');
tree.cloro.UPGMA(1)=Phyqualyread('cloro100UPGMA.tree');
tree.cloro.UPGMA(2)=Phyqualyread('cloro200UPGMA.tree');
tree.cloro.UPGMA(3)=Phyqualyread('cloro300UPGMA.tree');
tree.cloro.UPGMA(4)=Phyqualyread('cloro400UPGMA.tree');
tree.cloro.UPGMA(5)=Phyqualyread('cloro500UPGMA.tree');
tree.cloro.UPGMA(6)=Phyqualyread('cloro600UPGMA.tree');
tree.cloro.UPGMA(7)=Phyqualyread('cloro700UPGMA.tree');
tree.cloro.UPGMA(8)=Phyqualyread('cloro800UPGMA.tree');
tree.cloro.UPGMA(9)=Phyqualyread('cloro900UPGMA.tree');
tree.cloro.UPGMA(10)=Phyqualyread('cloro1000UPGMA.tree');
tree.cloro.UPGMA(11)=Phyqualyread('cloro1100UPGMA.tree');
tree.cloro.UPGMA(12)=Phyqualyread('cloro1200UPGMA.tree');
tree.cloro.UPGMA(13)=Phyqualyread('cloro1300UPGMA.tree');
tree.cloro.UPGMA(14)=Phyqualyread('cloro1400UPGMA.tree');
tree.cloro.UPGMA(15)=Phyqualyread('cloro1500UPGMA.tree');
tree.cloro.UPGMA(16)=Phyqualyread('cloro1600UPGMA.tree');
tree.cloro.UPGMA(17)=Phyqualyread('cloro1700UPGMA.tree');
tree.cloro.UPGMA(18)=Phyqualyread('cloro1800UPGMA.tree');
tree.cloro.UPGMA(19)=Phyqualyread('cloro1900UPGMA.tree');
tree.cloro.UPGMA(20)=Phyqualyread('cloro2000UPGMA.tree');
tree.cloro.UPGMA(21)=Phyqualyread('cloro2100UPGMA.tree');
tree.cloro.UPGMA(22)=Phyqualyread('cloro2200UPGMA.tree');
tree.cloro.UPGMA(23)=Phyqualyread('cloro2300UPGMA.tree');
tree.cloro.UPGMA(24)=Phyqualyread('cloro2400UPGMA.tree');
tree.cloro.UPGMA(25)=Phyqualyread('cloro2500UPGMA.tree');
tree.cloro.UPGMA(26)=Phyqualyread('cloro2600UPGMA.tree');
tree.cloro.UPGMA(27)=Phyqualyread('cloro2700UPGMA.tree');
tree.cloro.UPGMA(28)=Phyqualyread('cloro2800UPGMA.tree');
tree.cloro.UPGMA(29)=Phyqualyread('cloro2900UPGMA.tree');
tree.cloro.UPGMA(30)=Phyqualyread('cloro3000UPGMA.tree');
tree.plastid.NJ(1)=Phyqualyread('plastid100NJe.tree');
tree.plastid.NJ(2)=Phyqualyread('plastid200NJe.tree');
tree.plastid.NJ(3)=Phyqualyread('plastid300NJe.tree');
tree.plastid.NJ(4)=Phyqualyread('plastid400NJe.tree');
tree.plastid.NJ(5)=Phyqualyread('plastid500NJe.tree');
tree.plastid.NJ(6)=Phyqualyread('plastid600NJe.tree');
tree.plastid.NJ(7)=Phyqualyread('plastid700NJe.tree');
tree.plastid.NJ(8)=Phyqualyread('plastid800NJe.tree');
tree.plastid.NJ(9)=Phyqualyread('plastid900NJe.tree');
tree.plastid.NJ(10)=Phyqualyread('plastid1000NJe.tree');
tree.plastid.NJ(11)=Phyqualyread('plastid1100NJe.tree');
tree.plastid.NJ(12)=Phyqualyread('plastid1200NJe.tree');
tree.plastid.NJ(13)=Phyqualyread('plastid1300NJe.tree');
tree.plastid.NJ(14)=Phyqualyread('plastid1400NJe.tree');
tree.plastid.NJ(15)=Phyqualyread('plastid1500NJe.tree');
tree.plastid.NJ(16)=Phyqualyread('plastid1600NJe.tree');
tree.plastid.NJ(17)=Phyqualyread('plastid1700NJe.tree');
tree.plastid.NJ(18)=Phyqualyread('plastid1800NJe.tree');
tree.plastid.NJ(19)=Phyqualyread('plastid1900NJe.tree');
tree.plastid.NJ(20)=Phyqualyread('plastid2000NJe.tree');
tree.plastid.NJ(21)=Phyqualyread('plastid2100NJe.tree');
tree.plastid.NJ(22)=Phyqualyread('plastid2200NJe.tree');
tree.plastid.NJ(23)=Phyqualyread('plastid2300NJe.tree');
tree.plastid.NJ(24)=Phyqualyread('plastid2400NJe.tree');
tree.plastid.NJ(25)=Phyqualyread('plastid2500NJe.tree');
tree.plastid.NJ(26)=Phyqualyread('plastid2600NJe.tree');
tree.plastid.NJ(27)=Phyqualyread('plastid2700NJe.tree');
tree.plastid.NJ(28)=Phyqualyread('plastid2800NJe.tree');
tree.plastid.NJ(29)=Phyqualyread('plastid2900NJe.tree');
tree.plastid.NJ(30)=Phyqualyread('plastid3000NJe.tree');
tree.plastid.Ward(1)=Phyqualyread('plastid100Ward.tree');
tree.plastid.Ward(2)=Phyqualyread('plastid200Ward.tree');
tree.plastid.Ward(3)=Phyqualyread('plastid300Ward.tree');
tree.plastid.Ward(4)=Phyqualyread('plastid400Ward.tree');
tree.plastid.Ward(5)=Phyqualyread('plastid500Ward.tree');
tree.plastid.Ward(6)=Phyqualyread('plastid600Ward.tree');
tree.plastid.Ward(7)=Phyqualyread('plastid700Ward.tree');
tree.plastid.Ward(8)=Phyqualyread('plastid800Ward.tree');
tree.plastid.Ward(9)=Phyqualyread('plastid900Ward.tree');
tree.plastid.Ward(10)=Phyqualyread('plastid1000Ward.tree');
tree.plastid.Ward(11)=Phyqualyread('plastid1100Ward.tree');
tree.plastid.Ward(12)=Phyqualyread('plastid1200Ward.tree');
tree.plastid.Ward(13)=Phyqualyread('plastid1300Ward.tree');
tree.plastid.Ward(14)=Phyqualyread('plastid1400Ward.tree');
tree.plastid.Ward(15)=Phyqualyread('plastid1500Ward.tree');
tree.plastid.Ward(16)=Phyqualyread('plastid1600Ward.tree');
tree.plastid.Ward(17)=Phyqualyread('plastid1700Ward.tree');
tree.plastid.Ward(18)=Phyqualyread('plastid1800Ward.tree');
tree.plastid.Ward(19)=Phyqualyread('plastid1900Ward.tree');
tree.plastid.Ward(20)=Phyqualyread('plastid2000Ward.tree');
tree.plastid.Ward(21)=Phyqualyread('plastid2100Ward.tree');
tree.plastid.Ward(22)=Phyqualyread('plastid2200Ward.tree');
tree.plastid.Ward(23)=Phyqualyread('plastid2300Ward.tree');
tree.plastid.Ward(24)=Phyqualyread('plastid2400Ward.tree');
tree.plastid.Ward(25)=Phyqualyread('plastid2500Ward.tree');
tree.plastid.Ward(26)=Phyqualyread('plastid2600Ward.tree');
tree.plastid.Ward(27)=Phyqualyread('plastid2700Ward.tree');
tree.plastid.Ward(28)=Phyqualyread('plastid2800Ward.tree');
tree.plastid.Ward(29)=Phyqualyread('plastid2900Ward.tree');
tree.plastid.Ward(30)=Phyqualyread('plastid3000Ward.tree');
tree.plastid.UPGMA(1)=Phyqualyread('plastid100UPGMA.tree');
tree.plastid.UPGMA(2)=Phyqualyread('plastid200UPGMA.tree');
tree.plastid.UPGMA(3)=Phyqualyread('plastid300UPGMA.tree');
tree.plastid.UPGMA(4)=Phyqualyread('plastid400UPGMA.tree');
tree.plastid.UPGMA(5)=Phyqualyread('plastid500UPGMA.tree');
tree.plastid.UPGMA(6)=Phyqualyread('plastid600UPGMA.tree');
tree.plastid.UPGMA(7)=Phyqualyread('plastid700UPGMA.tree');
tree.plastid.UPGMA(8)=Phyqualyread('plastid800UPGMA.tree');
tree.plastid.UPGMA(9)=Phyqualyread('plastid900UPGMA.tree');
tree.plastid.UPGMA(10)=Phyqualyread('plastid1000UPGMA.tree');
tree.plastid.UPGMA(11)=Phyqualyread('plastid1100UPGMA.tree');
tree.plastid.UPGMA(12)=Phyqualyread('plastid1200UPGMA.tree');
tree.plastid.UPGMA(13)=Phyqualyread('plastid1300UPGMA.tree');
tree.plastid.UPGMA(14)=Phyqualyread('plastid1400UPGMA.tree');
tree.plastid.UPGMA(15)=Phyqualyread('plastid1500UPGMA.tree');
tree.plastid.UPGMA(16)=Phyqualyread('plastid1600UPGMA.tree');
tree.plastid.UPGMA(17)=Phyqualyread('plastid1700UPGMA.tree');
tree.plastid.UPGMA(18)=Phyqualyread('plastid1800UPGMA.tree');
tree.plastid.UPGMA(19)=Phyqualyread('plastid1900UPGMA.tree');
tree.plastid.UPGMA(20)=Phyqualyread('plastid2000UPGMA.tree');
tree.plastid.UPGMA(21)=Phyqualyread('plastid2100UPGMA.tree');
tree.plastid.UPGMA(22)=Phyqualyread('plastid2200UPGMA.tree');
tree.plastid.UPGMA(23)=Phyqualyread('plastid2300UPGMA.tree');
tree.plastid.UPGMA(24)=Phyqualyread('plastid2400UPGMA.tree');
tree.plastid.UPGMA(25)=Phyqualyread('plastid2500UPGMA.tree');
tree.plastid.UPGMA(26)=Phyqualyread('plastid2600UPGMA.tree');
tree.plastid.UPGMA(27)=Phyqualyread('plastid2700UPGMA.tree');
tree.plastid.UPGMA(28)=Phyqualyread('plastid2800UPGMA.tree');
tree.plastid.UPGMA(29)=Phyqualyread('plastid2900UPGMA.tree');
tree.plastid.UPGMA(30)=Phyqualyread('plastid3000UPGMA.tree');
cont=1;
for it=1:length(xfas)
if strfind(xfas(it).Header, 'Corydalis')
Corydalis(cont)=xfas(it);
cont=cont+1;
end
end
tree.plastid.NJ(31)=Phyqualyread('plastidNJeConsensus.tree');