Jess
unread,Nov 11, 2011, 11:12:06 AM11/11/11Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to raxml
Hello,
I've just been trying to create bootstrap consensus trees in RAxML.
7.2.8 and have come across some problems.
Firstly, when trying to create the three consensus trees, (in the -J
option - MR, MRE or STRICT), if I go to open them in figtree (or
treeedit) I get the error message: "Import Error: Error reading tree
file: Missing closing ')' in tree"
However, clearly this isn't the case, (there are the same number of
'(' and ')' brackets - and adding an extra ')' gives you the same
error message). Is this something to do with the way the consensus
tree allows polytomies at the root or something? The best trees all
open fine in figtree. It would be really great to be able to view
these consensus trees in figtree - but I can't work out how to correct
the consensus tree as it is produced from RAxML into something that
figtree will import. (I have pasted the MR tree below)
The second problem I have is that when drawing the bootstraps onto the
best tree (-f b) to create the bipartitions tree file, v7.2.8 produces
a bipartitions tree that is different to the best tree. In the
particular dataset I was using, it created a massive root branchlength
of 0.45. Calling exactly the same command on the old v7.0.4 version of
RAxML produced the correct best-tree with bootstraps.
I am happy to send you any of the files if that would help, otherwise
I have pasted the tree files below.
I would be very grateful if anyone is able to help with this. Many
thanks.
MR tree:
(D_torquatus_E_Kolyma2(79_layer12_M4D_6,49_layer12_M3_1,150_UK_Yngr_dryas,
(D_hudsonius_b,
((E023,144_layer16_C7_19,E015,E012,E050,E051,E029,E049,E018,E019,
(E052,E024,(E027,E020):1.0[83]):1.0[64]):1.0[75],(E026,(E036,E014):
1.0[68]):1.0[52]):1.0[62],
(((129_layer12_M4_6,134_Layer14_K6_9,146_layer16_M7_16,135_layer14_D7_16,133_layer14_K6_9,
(132_layer14_K7_9,145_layer16_M7_16):1.0[68]):1.0[85],
(130_12_M4_6_Haplo,136_layer14_D7_16):1.0[58]):1.0[57],(E047,
(131_layer12_M4_6,48_layer12_M3_9,47_layer12_M3_9):1.0[65]):1.0[64]):
1.0[78],(56_12_M4_5_Haplo,
46_layer12_M3_9,80_layer12_M4D_6,54_12_M4_5_Haplo,
86_layer12_M4D_6,73_layer12_M4D_6,62_layer12_M4_5,57_12_M4_5_Haplo,
66_layer12_M4_5,(65_layer12_M4_5,71_layer12_M4D_6):1.0[52]):1.0[76]):
1.0[66],(126_layer12_M3_2,(51_layer12_M3_1,128_12_M3_2_Haplo):1.0[56]):
1.0[65],(D_torquatus_Yamal,(D_torquatus_Pechora,D_torquatus_Kanin):
1.0[67]):1.0[61],(D_torquatus_Taymyr,D_torquatus_Taymyr2):1.0[99],
(153_UK_Yngr_dryas,149_UK_dryas_Haplo):1.0[73],
(D_torquatus_W_Kolyma,D_torquatus_Olenski):1.0[73]):1.0[61],
(D_torquatus_E_Kolyma,D_torquatus_E_Kolyma3):1.0[100]);
bipartition-bootstrap analysis called: ./raxmlHPC -m GTRGAMMA -f b -t
RAxML_bestTree.dicro.best.tre -z RAxML_bootstrap.dicro.boots -n
dicro.bip.2 -o D_hudsonius_b
besttree: ((((((48_layer12_M3_9:0.00000363619534575013,
(47_layer12_M3_9:0.00129942266970888877,131_layer12_M4_6:0.00131404211215039321):
0.00000363619534575013):
0.00132730118271114994,E047:0.00130513155135699973):
0.00134846575657471333,((((133_layer14_K6_9:0.00134145769075794011,
(135_layer14_D7_16:0.00271945986349485073,
(129_layer12_M4_6:0.00546163750744719286,146_layer16_M7_16:0.00000363619534575013):
0.00133568851414785777):0.00000363619534575013):
0.00000363619534575013,134_Layer14_K6_9:0.00000363619534575013):
0.00000363619534575013,
(145_layer16_M7_16:0.00132774593050852567,132_layer14_K7_9:0.00000363619534575013):
0.00133713032142029464):0.00278647994986401396,
(136_layer14_D7_16:0.00000363619534575013,130_12_M4_6_Haplo:
0.00131994653001272968):0.00141651841126976832):
0.00125910185854779654):0.00403005260407361416,
((62_layer12_M4_5:0.00126937068342301146,
((86_layer12_M4D_6:0.00128584289693581526,
(73_layer12_M4D_6:0.00128363904623272625,57_12_M4_5_Haplo:
0.00000363619534575013):0.00000363619534575013):0.00000363619534575013,
(((66_layer12_M4_5:0.00128838377089412386,54_12_M4_5_Haplo:
0.00128853886813459946):0.00000363619534575013,
(((65_layer12_M4_5:0.00128884461869839973,71_layer12_M4D_6:0.00000363619534575013):
0.00000363619534575013,46_layer12_M3_9:0.00000363619534575013):
0.00000363619534575013,56_12_M4_5_Haplo:0.00000363619534575013):
0.00000363619534575013):
0.00000363619534575013,80_layer12_M4D_6:0.00129274952360672085):
0.00128747536906372701):0.00130570538219615829):0.00817091584735905874,
(49_layer12_M3_1:0.00256999670460223003,(((((D_torquatus_E_Kolyma:
0.00120531358297111383,D_torquatus_E_Kolyma3:0.00000363619534575013):
0.00771015462378700587,D_torquatus_E_Kolyma2:0.00497959377899245825):
0.00395369005417124230,79_layer12_M4D_6:0.00124108783173599549):
0.00000363619534575013,(126_layer12_M3_2:0.00125541203451745163,
(51_layer12_M3_1:0.00000363619534575013,128_12_M3_2_Haplo:
0.00000363619534575013):0.00000363619534575013):
0.00124131096933333888):0.00000363619534575013,(((D_torquatus_Taymyr:
0.00000363619534575013,D_torquatus_Taymyr2:0.00000363619534575013):
0.00517034449906722277,((D_torquatus_Yamal:0.00265788371535682563,
(D_torquatus_Pechora:0.00260245233305352475,D_torquatus_Kanin:
0.00522030069062469609):0.00128598315030102996):0.00128459424545687923,
(150_UK_Yngr_dryas:0.00520455735781596464,(149_UK_dryas_Haplo:
0.00125427769367385147,153_UK_Yngr_dryas:0.00000363619534575013):
0.00000363619534575013):0.00000363619534575013):
0.00000363619534575013):0.00127745439282606044,(D_torquatus_W_Kolyma:
0.00259264052177891153,D_torquatus_Olenski:0.00122822421295203338):
0.00261310670036057598):0.00124373314667413411):
0.00000363619534575013):0.00403594532436997271):
0.00106724494576406663):0.00266572786165118845,
((E026:0.00206182631963185582,
(E036:0.00000363619534575013,E014:0.00262369135172486109):
0.00185264307150822920):0.00258387553564366377,
(((E019:0.00131795249856911140,((E049:0.00267010389584985573,
(E029:0.00000363619534575013,(E051:0.00844430881099344881,
(E050:0.00142855263970429001,
((E020:0.00144999712504058635,E027:0.00983731550324775349):
0.00539467881618129093,
(E052:0.00264611972331375844,E024:0.00000363619534575013):
0.00000363619534575013):0.00411087587953217360):
0.00278091458703202867):0.00130152933468939518):
0.00135011154385100045):
0.00131564916088999112,E018:0.00000363619534575013):
0.00000363619534575013):
0.00130834356062447198,E023:0.00000363619534575013):
0.00000363619534575013,
((144_layer16_C7_19:0.01276570132063204431,E015:0.00866622710162602734):
0.00567997135407338500,E012:0.00234361438279412837):
0.00165979415839065956):0.00761133019019519068):
0.00198464714858766602):0.02682618930035064048,D_hudsonius_b:
0.02682618930035064048);
bipartition tree:
((((((48_layer12_M3_9:0.00000363619534575013,
(47_layer12_M3_9:0.00129942266970888877,131_layer12_M4_6:0.00131404211215039321)16:0.00000363619534575013)65:0.00132730118271114994,E047:0.00130513155135699973)64:0.00134846575657471333,
((((133_layer14_K6_9:0.00134145769075794011,
(135_layer14_D7_16:0.00271945986349485073,
(129_layer12_M4_6:0.00546163750744719286,146_layer16_M7_16:0.00000363619534575013)47:0.00133568851414785777)17:0.00000363619534575013)4:0.00000363619534575013,134_Layer14_K6_9:0.00000363619534575013)12:0.00000363619534575013,
(145_layer16_M7_16:0.00132774593050852567,132_layer14_K7_9:0.00000363619534575013)68:0.00133713032142029464)85:0.00278647994986401396,
(136_layer14_D7_16:0.00000363619534575013,130_12_M4_6_Haplo:
0.00131994653001272968)58:0.00141651841126976832)57:0.00125910185854779654)78:0.00403005260407361416,
((62_layer12_M4_5:0.00126937068342301146,
((86_layer12_M4D_6:0.00128584289693581526,
(73_layer12_M4D_6:0.00128363904623272625,57_12_M4_5_Haplo:
0.00000363619534575013)19:0.00000363619534575013)26:0.00000363619534575013,
(((66_layer12_M4_5:0.00128838377089412386,54_12_M4_5_Haplo:
0.00128853886813459946)7:0.00000363619534575013,
(((65_layer12_M4_5:0.00128884461869839973,71_layer12_M4D_6:0.00000363619534575013)52:0.00000363619534575013,46_layer12_M3_9:0.00000363619534575013)4:0.00000363619534575013,56_12_M4_5_Haplo:
0.00000363619534575013)8:0.00000363619534575013)12:0.00000363619534575013,80_layer12_M4D_6:0.00129274952360672085)26:0.00128747536906372701)38:0.00130570538219615829)76:0.00817091584735905874,
(49_layer12_M3_1:0.00256999670460223003,(((((D_torquatus_E_Kolyma:
0.00120531358297111383,D_torquatus_E_Kolyma3:0.00000363619534575013)100:0.00771015462378700587,D_torquatus_E_Kolyma2:0.00497959377899245825)61:0.00395369005417124230,79_layer12_M4D_6:0.00124108783173599549)11:0.00000363619534575013,
(126_layer12_M3_2:0.00125541203451745163,
(51_layer12_M3_1:0.00000363619534575013,128_12_M3_2_Haplo:
0.00000363619534575013)56:0.00000363619534575013)65:0.00124131096933333888)5:0.00000363619534575013,
(((D_torquatus_Taymyr:
0.00000363619534575013,D_torquatus_Taymyr2:0.00000363619534575013)99:0.00517034449906722277,
((D_torquatus_Yamal:0.00265788371535682563,(D_torquatus_Pechora:
0.00260245233305352475,D_torquatus_Kanin:
0.00522030069062469609)67:0.00128598315030102996)61:0.00128459424545687923,
(150_UK_Yngr_dryas:0.00520455735781596464,(149_UK_dryas_Haplo:
0.00125427769367385147,153_UK_Yngr_dryas:
0.00000363619534575013)73:0.00000363619534575013)26:0.00000363619534575013)8:0.00000363619534575013)31:0.00127745439282606044,
(D_torquatus_W_Kolyma:0.00259264052177891153,D_torquatus_Olenski:
0.00122822421295203338)73:0.00261310670036057598)35:0.00124373314667413411)11:0.00000363619534575013)66:0.00403594532436997271)24:0.00106724494576406663)33:0.00266572786165118845,
((E026:0.00206182631963185582,
(E036:0.00000363619534575013,E014:0.00262369135172486109)68:0.00185264307150822920)52:0.00258387553564366377,
(((E019:0.00131795249856911140,((E049:0.00267010389584985573,
(E029:0.00000363619534575013,(E051:0.00844430881099344881,
(E050:0.00142855263970429001,
((E020:0.00144999712504058635,E027:0.00983731550324775349)83:0.00539467881618129093,
(E052:0.00264611972331375844,E024:0.00000363619534575013)34:0.00000363619534575013)64:0.00411087587953217360)49:0.00278091458703202867)23:0.00130152933468939518)27:0.00135011154385100045)30:0.00131564916088999112,E018:0.00000363619534575013)4:0.00000363619534575013)17:0.00130834356062447198,E023:0.00000363619534575013)12:0.00000363619534575013,
((144_layer16_C7_19:0.01276570132063204431,E015:0.00866622710162602734)48:0.00567997135407338500,E012:0.00234361438279412837)22:0.00165979415839065956)75:0.00761133019019519068)62:0.00198464714858766602):
0.45000000000000001110,D_hudsonius_b:0.45000000000000001110);