'Clockifying' maximum likelihood trees

Load the libraries

library(ape)
library(stringr)
library(phangorn)
library(magrittr)
source("chronos.R")

Preprocessing the data

We need is a tree, plus how long the sequences were used to construct the tree. I will use a tree generated using ExaML, but in principle, you could use any tree, including one generated using optim.pml in phangorn.

stub <- "village_examl"
seq.len <- 6987
tr.fn <- paste(stub,".tre",sep="")
td.fn <- paste(stub,".td",sep="")
tr <- read.tree(tr.fn)

The sample times are in the sequence names at the end.

tipnames <- tr$tip.label
tipdates <- tipnames %>% strsplit(.,"_",fixed=TRUE) %>% lapply(.,tail,1) %>% unlist %>% as.double

Now save the sequence names and tip dates to a file.

write.table(rbind(c(length(tipnames),""),cbind(tipnames,tipdates)),td.fn,col.names=FALSE,row.names=FALSE,quote=FALSE)

Root to tip regression

Using the tree and the tip dates (in the same order as the tip labels in the tree), we can root the tree using rtt.

tr.rtt <- rtt(tr,tipdates)
tr.rtt$edge.length[tr.rtt$edge.length<0.00000001] <- 0.00000001
tr.rtt.tipnames <- tr.rtt$tip.label
tr.rtt.tipdates <- tr.rtt.tipnames %>% strsplit(.,"_",fixed=TRUE) %>% lapply(.,tail,1) %>% unlist %>% as.double

Now the tree is rooted, we can calculate the evolutionary rate, as the slope of the regression of root-to-tip distances against sampling times.

rootdistance <- distRoot(tr.rtt)
pathlm <- lm(rootdistance~tr.rtt.tipdates)
rate <- coef(pathlm)[2]
rate
## tr.rtt.tipdates 
##     0.002587992

The time to the most recent common ancestor can then be obtained using the intercept term and the slope, as it is the point at which the regression line crosses zero.

tmrca <- unname(-coef(pathlm)[1]/coef(pathlm)[2])
tmrca
## [1] 0.732083

Chronos

For us to be able to run chronos, we have to set up a quite complex calibration structure.

max.time <- max(tipdates)
ncat <- 1
strict.clock.ctrl <- chronos.control(nb.rate.cat=as.integer(ncat))
calibrating.values <- makeChronosCalib(tr)
calibrating.values$age.min <- max.time - tmrca
calibrating.values$age.max <- max.time - tmrca
# pins the tips to sampling years
calibrating.values <- rbind(calibrating.values,
                            data.frame(node=seq(1,length(td)),
                                       age.min=max.time - td,
                                       age.max=max.time - td,
                                       soft.bounds=FALSE))

We now set options for the model.

ncat <- 1
sub.rate.model <- "discrete"
strict.clock.ctrl <- chronos.control(nb.rate.cat=as.integer(ncat))

We can now fit a dated tree.

chronos.tree <- RLchronos(tr, 
                     lambda=1, 
                     model=sub.rate.model, 
                     calibration=calibrating.values,
                     control=strict.clock.ctrl,
                     quiet=TRUE)
## Initial dates set.
plot(chronos.tree,show.tip.label=FALSE)
add.scale.bar()

plot of chunk unnamed-chunk-12
Figure: plot of chunk unnamed-chunk-12

We can write the tree to a file as a Newick string with write.tree.

write.tree(dated.tree,file="village_chronos.tre")

Least squares dating

lsd is an external application, which can be run using a system call. It requires:

  • The filename of the tree
  • The filename of the sampling times
  • A variance parameter (default 10)
  • The sequence length
lsd.cmd <- sprintf("lsd -i %s -d %s -c -n 1 -r -b %s -s %s -v",tr.fn,td.fn,paste(10),seq.len)
lsd.cmd
## [1] "lsd -i village_examl.tre -d village_examl.td -c -n 1 -r -b 10 -s 6987 -v"
lsd <- system(lsd.cmd,intern=TRUE)

There is lots of information in the result.

lsd
##    [1] "Reading the tree 1..."                                                         
##    [2] "Pre-estimating the root without constraints ..."                               
##    [3] "Optimizing the root position on the branch 1 ... 0.0100313653"                 
##    [4] "Optimizing the root position on the branch 2 ... 0.0090787440"                 
##    [5] "Optimizing the root position on the branch 3 ... 0.0090534471"                 
##    [6] "Optimizing the root position on the branch 4 ... 0.0147294111"                 
##    [7] "Optimizing the root position on the branch 5 ... 0.0150153208"                 
##    [8] "Optimizing the root position on the branch 6 ... 0.0180342072"                 
##    [9] "Optimizing the root position on the branch 7 ... 0.0090998256"                 
##   [10] "Optimizing the root position on the branch 8 ... 0.0090896480"                 
##   [11] "Optimizing the root position on the branch 9 ... 0.0090891878"                 
##   [12] "Optimizing the root position on the branch 10 ... 0.0099988852"                
##   [13] "Optimizing the root position on the branch 11 ... 0.0102494605"                
##   [14] "Optimizing the root position on the branch 12 ... 0.0102921770"                
##   [15] "Optimizing the root position on the branch 13 ... 0.0103728229"                
##   [16] "Optimizing the root position on the branch 14 ... 0.0122575648"                
##   [17] "Optimizing the root position on the branch 15 ... 0.0124530444"                
##   [18] "Optimizing the root position on the branch 16 ... 0.0126307471"                
##   [19] "Optimizing the root position on the branch 17 ... 0.0126357945"                
##   [20] "Optimizing the root position on the branch 18 ... 0.0103655914"                
##   [21] "Optimizing the root position on the branch 19 ... 0.0101469853"                
##   [22] "Optimizing the root position on the branch 20 ... 0.0116820327"                
##   [23] "Optimizing the root position on the branch 21 ... 0.0118149356"                
##   [24] "Optimizing the root position on the branch 22 ... 0.0116291749"                
##   [25] "Optimizing the root position on the branch 23 ... 0.0118342906"                
##   [26] "Optimizing the root position on the branch 24 ... 0.0119317621"                
##   [27] "Optimizing the root position on the branch 25 ... 0.0131969372"                
##   [28] "Optimizing the root position on the branch 26 ... 0.0133652061"                
##   [29] "Optimizing the root position on the branch 27 ... 0.0134559361"                
##   [30] "Optimizing the root position on the branch 28 ... 0.0134621686"                
##   [31] "Optimizing the root position on the branch 29 ... 0.0090664546"                
##   [32] "Optimizing the root position on the branch 30 ... 0.0089555290"                
##   [33] "Optimizing the root position on the branch 31 ... 0.0089542341"                
##   [34] "Optimizing the root position on the branch 32 ... 0.0089881570"                
##   [35] "Optimizing the root position on the branch 33 ... 0.0089909891"                
##   [36] "Optimizing the root position on the branch 34 ... 0.0089814172"                
##   [37] "Optimizing the root position on the branch 35 ... 0.0090408912"                
##   [38] "Optimizing the root position on the branch 36 ... 0.0091707028"                
##   [39] "Optimizing the root position on the branch 37 ... 0.0097532962"                
##   [40] "Optimizing the root position on the branch 38 ... 0.0136307627"                
##   [41] "Optimizing the root position on the branch 39 ... 0.0136325525"                
##   [42] "Optimizing the root position on the branch 40 ... 0.0141591168"                
##   [43] "Optimizing the root position on the branch 41 ... 0.0143344263"                
##   [44] "Optimizing the root position on the branch 42 ... 0.0147172654"                
##   [45] "Optimizing the root position on the branch 43 ... 0.0147644562"                
##   [46] "Optimizing the root position on the branch 44 ... 0.0147629382"                
##   [47] "Optimizing the root position on the branch 45 ... 0.0097481810"                
##   [48] "Optimizing the root position on the branch 46 ... 0.0090315389"                
##   [49] "Optimizing the root position on the branch 47 ... 0.0091102437"                
##   [50] "Optimizing the root position on the branch 48 ... 0.0091402305"                
##   [51] "Optimizing the root position on the branch 49 ... 0.0091468662"                
##   [52] "Optimizing the root position on the branch 50 ... 0.0092077595"                
##   [53] "Optimizing the root position on the branch 51 ... 0.0100218049"                
##   [54] "Optimizing the root position on the branch 52 ... 0.0103080343"                
##   [55] "Optimizing the root position on the branch 53 ... 0.0106215815"                
##   [56] "Optimizing the root position on the branch 54 ... 0.0106540804"                
##   [57] "Optimizing the root position on the branch 55 ... 0.0106703013"                
##   [58] "Optimizing the root position on the branch 56 ... 0.0106699537"                
##   [59] "Optimizing the root position on the branch 57 ... 0.0106689158"                
##   [60] "Optimizing the root position on the branch 58 ... 0.0106643950"                
##   [61] "Optimizing the root position on the branch 59 ... 0.0106561076"                
##   [62] "Optimizing the root position on the branch 60 ... 0.0107286516"                
##   [63] "Optimizing the root position on the branch 61 ... 0.0112546545"                
##   [64] "Optimizing the root position on the branch 62 ... 0.0112721132"                
##   [65] "Optimizing the root position on the branch 63 ... 0.0107547676"                
##   [66] "Optimizing the root position on the branch 64 ... 0.0102841847"                
##   [67] "Optimizing the root position on the branch 65 ... 0.0104204325"                
##   [68] "Optimizing the root position on the branch 66 ... 0.0115541718"                
##   [69] "Optimizing the root position on the branch 67 ... 0.0121582956"                
##   [70] "Optimizing the root position on the branch 68 ... 0.0100022054"                
##   [71] "Optimizing the root position on the branch 69 ... 0.0108469995"                
##   [72] "Optimizing the root position on the branch 70 ... 0.0118839960"                
##   [73] "Optimizing the root position on the branch 71 ... 0.0120335805"                
##   [74] "Optimizing the root position on the branch 72 ... 0.0120213680"                
##   [75] "Optimizing the root position on the branch 73 ... 0.0121075965"                
##   [76] "Optimizing the root position on the branch 74 ... 0.0121022654"                
##   [77] "Optimizing the root position on the branch 75 ... 0.0121425227"                
##   [78] "Optimizing the root position on the branch 76 ... 0.0121441439"                
##   [79] "Optimizing the root position on the branch 77 ... 0.0118517082"                
##   [80] "Optimizing the root position on the branch 78 ... 0.0092370043"                
##   [81] "Optimizing the root position on the branch 79 ... 0.0143984859"                
##   [82] "Optimizing the root position on the branch 80 ... 0.0151305305"                
##   [83] "Optimizing the root position on the branch 81 ... 0.0154358188"                
##   [84] "Optimizing the root position on the branch 82 ... 0.0154796980"                
##   [85] "Optimizing the root position on the branch 83 ... 0.0154714234"                
##   [86] "Optimizing the root position on the branch 84 ... 0.0150541348"                
##   [87] "Optimizing the root position on the branch 85 ... 0.0091428930"                
##   [88] "Optimizing the root position on the branch 86 ... 0.0091820010"                
##   [89] "Optimizing the root position on the branch 87 ... 0.0095917579"                
##   [90] "Optimizing the root position on the branch 88 ... 0.0111226057"                
##   [91] "Optimizing the root position on the branch 89 ... 0.0096561958"                
##   [92] "Optimizing the root position on the branch 90 ... 0.0125992199"                
##   [93] "Optimizing the root position on the branch 91 ... 0.0126108172"                
##   [94] "Optimizing the root position on the branch 92 ... 0.0092102482"                
##   [95] "Optimizing the root position on the branch 93 ... 0.0143158341"                
##   [96] "Optimizing the root position on the branch 94 ... 0.0150595795"                
##   [97] "Optimizing the root position on the branch 95 ... 0.0089884614"                
##   [98] "Optimizing the root position on the branch 96 ... 0.0089983392"                
##   [99] "Optimizing the root position on the branch 97 ... 0.0090000299"                
##  [100] "Optimizing the root position on the branch 98 ... 0.0090078647"                
##  [101] "Optimizing the root position on the branch 99 ... 0.0098816650"                
##  [102] "Optimizing the root position on the branch 100 ... 0.0099839750"               
##  [103] "Optimizing the root position on the branch 101 ... 0.0102946418"               
##  [104] "Optimizing the root position on the branch 102 ... 0.0102878193"               
##  [105] "Optimizing the root position on the branch 103 ... 0.0100023907"               
##  [106] "Optimizing the root position on the branch 104 ... 0.0101687462"               
##  [107] "Optimizing the root position on the branch 105 ... 0.0101668688"               
##  [108] "Optimizing the root position on the branch 106 ... 0.0105930607"               
##  [109] "Optimizing the root position on the branch 107 ... 0.0108108207"               
##  [110] "Optimizing the root position on the branch 108 ... 0.0107983898"               
##  [111] "Optimizing the root position on the branch 109 ... 0.0109693939"               
##  [112] "Optimizing the root position on the branch 110 ... 0.0110166373"               
##  [113] "Optimizing the root position on the branch 111 ... 0.0110847196"               
##  [114] "Optimizing the root position on the branch 112 ... 0.0110945603"               
##  [115] "Optimizing the root position on the branch 113 ... 0.0110952845"               
##  [116] "Optimizing the root position on the branch 114 ... 0.0110309748"               
##  [117] "Optimizing the root position on the branch 115 ... 0.0112825352"               
##  [118] "Optimizing the root position on the branch 116 ... 0.0105998122"               
##  [119] "Optimizing the root position on the branch 117 ... 0.0107958636"               
##  [120] "Optimizing the root position on the branch 118 ... 0.0107786999"               
##  [121] "Optimizing the root position on the branch 119 ... 0.0101115931"               
##  [122] "Optimizing the root position on the branch 120 ... 0.0113140901"               
##  [123] "Optimizing the root position on the branch 121 ... 0.0129475561"               
##  [124] "Optimizing the root position on the branch 122 ... 0.0129912503"               
##  [125] "Optimizing the root position on the branch 123 ... 0.0130149302"               
##  [126] "Optimizing the root position on the branch 124 ... 0.0130063751"               
##  [127] "Optimizing the root position on the branch 125 ... 0.0113446206"               
##  [128] "Optimizing the root position on the branch 126 ... 0.0115703411"               
##  [129] "Optimizing the root position on the branch 127 ... 0.0098138148"               
##  [130] "Optimizing the root position on the branch 128 ... 0.0090012002"               
##  [131] "Optimizing the root position on the branch 129 ... 0.0099239516"               
##  [132] "Optimizing the root position on the branch 130 ... 0.0162807180"               
##  [133] "Optimizing the root position on the branch 131 ... 0.0090028153"               
##  [134] "Optimizing the root position on the branch 132 ... 0.0090036913"               
##  [135] "Optimizing the root position on the branch 133 ... 0.0091867603"               
##  [136] "Optimizing the root position on the branch 134 ... 0.0095484667"               
##  [137] "Optimizing the root position on the branch 135 ... 0.0095559312"               
##  [138] "Optimizing the root position on the branch 136 ... 0.0095390563"               
##  [139] "Optimizing the root position on the branch 137 ... 0.0095261688"               
##  [140] "Optimizing the root position on the branch 138 ... 0.0095302096"               
##  [141] "Optimizing the root position on the branch 139 ... 0.0096476150"               
##  [142] "Optimizing the root position on the branch 140 ... 0.0096563151"               
##  [143] "Optimizing the root position on the branch 141 ... 0.0121262975"               
##  [144] "Optimizing the root position on the branch 142 ... 0.0121845384"               
##  [145] "Optimizing the root position on the branch 143 ... 0.0121563055"               
##  [146] "Optimizing the root position on the branch 144 ... 0.0132819908"               
##  [147] "Optimizing the root position on the branch 145 ... 0.0096597558"               
##  [148] "Optimizing the root position on the branch 146 ... 0.0102963093"               
##  [149] "Optimizing the root position on the branch 147 ... 0.0104228704"               
##  [150] "Optimizing the root position on the branch 148 ... 0.0105026739"               
##  [151] "Optimizing the root position on the branch 149 ... 0.0125016154"               
##  [152] "Optimizing the root position on the branch 150 ... 0.0104225789"               
##  [153] "Optimizing the root position on the branch 151 ... 0.0116570929"               
##  [154] "Optimizing the root position on the branch 152 ... 0.0116971440"               
##  [155] "Optimizing the root position on the branch 153 ... 0.0117713472"               
##  [156] "Optimizing the root position on the branch 154 ... 0.0119060868"               
##  [157] "Optimizing the root position on the branch 155 ... 0.0117675200"               
##  [158] "Optimizing the root position on the branch 156 ... 0.0119575240"               
##  [159] "Optimizing the root position on the branch 157 ... 0.0119817600"               
##  [160] "Optimizing the root position on the branch 158 ... 0.0103000389"               
##  [161] "Optimizing the root position on the branch 159 ... 0.0096545526"               
##  [162] "Optimizing the root position on the branch 160 ... 0.0097829066"               
##  [163] "Optimizing the root position on the branch 161 ... 0.0099170141"               
##  [164] "Optimizing the root position on the branch 162 ... 0.0107445471"               
##  [165] "Optimizing the root position on the branch 163 ... 0.0110217897"               
##  [166] "Optimizing the root position on the branch 164 ... 0.0113050370"               
##  [167] "Optimizing the root position on the branch 165 ... 0.0114966166"               
##  [168] "Optimizing the root position on the branch 166 ... 0.0114882839"               
##  [169] "Optimizing the root position on the branch 167 ... 0.0115456052"               
##  [170] "Optimizing the root position on the branch 168 ... 0.0115946773"               
##  [171] "Optimizing the root position on the branch 169 ... 0.0099104755"               
##  [172] "Optimizing the root position on the branch 170 ... 0.0101124436"               
##  [173] "Optimizing the root position on the branch 171 ... 0.0104906374"               
##  [174] "Optimizing the root position on the branch 172 ... 0.0105416812"               
##  [175] "Optimizing the root position on the branch 173 ... 0.0106173035"               
##  [176] "Optimizing the root position on the branch 174 ... 0.0106976331"               
##  [177] "Optimizing the root position on the branch 175 ... 0.0104907621"               
##  [178] "Optimizing the root position on the branch 176 ... 0.0097652716"               
##  [179] "Optimizing the root position on the branch 177 ... 0.0103103657"               
##  [180] "Optimizing the root position on the branch 178 ... 0.0094932470"               
##  [181] "Optimizing the root position on the branch 179 ... 0.0095528955"               
##  [182] "Optimizing the root position on the branch 180 ... 0.0095642079"               
##  [183] "Optimizing the root position on the branch 181 ... 0.0092228832"               
##  [184] "Optimizing the root position on the branch 182 ... 0.0117182460"               
##  [185] "Optimizing the root position on the branch 183 ... 0.0115296489"               
##  [186] "Optimizing the root position on the branch 184 ... 0.0152075174"               
##  [187] "Optimizing the root position on the branch 185 ... 0.0155696538"               
##  [188] "Optimizing the root position on the branch 186 ... 0.0090044626"               
##  [189] "Optimizing the root position on the branch 187 ... 0.0090075279"               
##  [190] "Optimizing the root position on the branch 188 ... 0.0097358324"               
##  [191] "Optimizing the root position on the branch 189 ... 0.0097757366"               
##  [192] "Optimizing the root position on the branch 190 ... 0.0121410001"               
##  [193] "Optimizing the root position on the branch 191 ... 0.0097494935"               
##  [194] "Optimizing the root position on the branch 192 ... 0.0099405533"               
##  [195] "Optimizing the root position on the branch 193 ... 0.0101947982"               
##  [196] "Optimizing the root position on the branch 194 ... 0.0125783665"               
##  [197] "Optimizing the root position on the branch 195 ... 0.0126816630"               
##  [198] "Optimizing the root position on the branch 196 ... 0.0129115272"               
##  [199] "Optimizing the root position on the branch 197 ... 0.0129652687"               
##  [200] "Optimizing the root position on the branch 198 ... 0.0101767885"               
##  [201] "Optimizing the root position on the branch 199 ... 0.0102838409"               
##  [202] "Optimizing the root position on the branch 200 ... 0.0135238139"               
##  [203] "Optimizing the root position on the branch 201 ... 0.0102920709"               
##  [204] "Optimizing the root position on the branch 202 ... 0.0110416395"               
##  [205] "Optimizing the root position on the branch 203 ... 0.0117183849"               
##  [206] "Optimizing the root position on the branch 204 ... 0.0117787935"               
##  [207] "Optimizing the root position on the branch 205 ... 0.0117890125"               
##  [208] "Optimizing the root position on the branch 206 ... 0.0120117291"               
##  [209] "Optimizing the root position on the branch 207 ... 0.0117049804"               
##  [210] "Optimizing the root position on the branch 208 ... 0.0118930339"               
##  [211] "Optimizing the root position on the branch 209 ... 0.0119703283"               
##  [212] "Optimizing the root position on the branch 210 ... 0.0119690311"               
##  [213] "Optimizing the root position on the branch 211 ... 0.0099109962"               
##  [214] "Optimizing the root position on the branch 212 ... 0.0102569009"               
##  [215] "Optimizing the root position on the branch 213 ... 0.0102768737"               
##  [216] "Optimizing the root position on the branch 214 ... 0.0109697119"               
##  [217] "Optimizing the root position on the branch 215 ... 0.0113326081"               
##  [218] "Optimizing the root position on the branch 216 ... 0.0113568564"               
##  [219] "Optimizing the root position on the branch 217 ... 0.0114979939"               
##  [220] "Optimizing the root position on the branch 218 ... 0.0109565935"               
##  [221] "Optimizing the root position on the branch 219 ... 0.0122759012"               
##  [222] "Optimizing the root position on the branch 220 ... 0.0096749090"               
##  [223] "Optimizing the root position on the branch 221 ... 0.0098731743"               
##  [224] "Optimizing the root position on the branch 222 ... 0.0099075101"               
##  [225] "Optimizing the root position on the branch 223 ... 0.0128107455"               
##  [226] "Optimizing the root position on the branch 224 ... 0.0098711367"               
##  [227] "Optimizing the root position on the branch 225 ... 0.0101420266"               
##  [228] "Optimizing the root position on the branch 226 ... 0.0103926174"               
##  [229] "Optimizing the root position on the branch 227 ... 0.0111457939"               
##  [230] "Optimizing the root position on the branch 228 ... 0.0112593070"               
##  [231] "Optimizing the root position on the branch 229 ... 0.0112846521"               
##  [232] "Optimizing the root position on the branch 230 ... 0.0113201903"               
##  [233] "Optimizing the root position on the branch 231 ... 0.0113487075"               
##  [234] "Optimizing the root position on the branch 232 ... 0.0113502007"               
##  [235] "Optimizing the root position on the branch 233 ... 0.0113194970"               
##  [236] "Optimizing the root position on the branch 234 ... 0.0120860280"               
##  [237] "Optimizing the root position on the branch 235 ... 0.0121143750"               
##  [238] "Optimizing the root position on the branch 236 ... 0.0121526095"               
##  [239] "Optimizing the root position on the branch 237 ... 0.0121549406"               
##  [240] "Optimizing the root position on the branch 238 ... 0.0121919209"               
##  [241] "Optimizing the root position on the branch 239 ... 0.0112590113"               
##  [242] "Optimizing the root position on the branch 240 ... 0.0111294054"               
##  [243] "Optimizing the root position on the branch 241 ... 0.0104082794"               
##  [244] "Optimizing the root position on the branch 242 ... 0.0114459017"               
##  [245] "Optimizing the root position on the branch 243 ... 0.0090080087"               
##  [246] "Optimizing the root position on the branch 244 ... 0.0097486102"               
##  [247] "Optimizing the root position on the branch 245 ... 0.0097445097"               
##  [248] "Optimizing the root position on the branch 246 ... 0.0090017371"               
##  [249] "Optimizing the root position on the branch 247 ... 0.0090073054"               
##  [250] "Optimizing the root position on the branch 248 ... 0.0089885047"               
##  [251] "Optimizing the root position on the branch 249 ... 0.0090274724"               
##  [252] "Optimizing the root position on the branch 250 ... 0.0090278803"               
##  [253] "Optimizing the root position on the branch 251 ... 0.0096326532"               
##  [254] "Optimizing the root position on the branch 252 ... 0.0096549325"               
##  [255] "Optimizing the root position on the branch 253 ... 0.0096651448"               
##  [256] "Optimizing the root position on the branch 254 ... 0.0096959088"               
##  [257] "Optimizing the root position on the branch 255 ... 0.0096950401"               
##  [258] "Optimizing the root position on the branch 256 ... 0.0097297304"               
##  [259] "Optimizing the root position on the branch 257 ... 0.0097812456"               
##  [260] "Optimizing the root position on the branch 258 ... 0.0121524904"               
##  [261] "Optimizing the root position on the branch 259 ... 0.0097827129"               
##  [262] "Optimizing the root position on the branch 260 ... 0.0100748480"               
##  [263] "Optimizing the root position on the branch 261 ... 0.0111620591"               
##  [264] "Optimizing the root position on the branch 262 ... 0.0115477465"               
##  [265] "Optimizing the root position on the branch 263 ... 0.0116776012"               
##  [266] "Optimizing the root position on the branch 264 ... 0.0116615748"               
##  [267] "Optimizing the root position on the branch 265 ... 0.0117884120"               
##  [268] "Optimizing the root position on the branch 266 ... 0.0118012797"               
##  [269] "Optimizing the root position on the branch 267 ... 0.0096901911"               
##  [270] "Optimizing the root position on the branch 268 ... 0.0097099032"               
##  [271] "Optimizing the root position on the branch 269 ... 0.0097173435"               
##  [272] "Optimizing the root position on the branch 270 ... 0.0123732520"               
##  [273] "Optimizing the root position on the branch 271 ... 0.0131277869"               
##  [274] "Optimizing the root position on the branch 272 ... 0.0097146428"               
##  [275] "Optimizing the root position on the branch 273 ... 0.0099037942"               
##  [276] "Optimizing the root position on the branch 274 ... 0.0097033878"               
##  [277] "Optimizing the root position on the branch 275 ... 0.0097092876"               
##  [278] "Optimizing the root position on the branch 276 ... 0.0100809393"               
##  [279] "Optimizing the root position on the branch 277 ... 0.0097075392"               
##  [280] "Optimizing the root position on the branch 278 ... 0.0097351367"               
##  [281] "Optimizing the root position on the branch 279 ... 0.0097334133"               
##  [282] "Optimizing the root position on the branch 280 ... 0.0097060839"               
##  [283] "Optimizing the root position on the branch 281 ... 0.0101738206"               
##  [284] "Optimizing the root position on the branch 282 ... 0.0097169702"               
##  [285] "Optimizing the root position on the branch 283 ... 0.0139126141"               
##  [286] "Optimizing the root position on the branch 284 ... 0.0095644324"               
##  [287] "Optimizing the root position on the branch 285 ... 0.0090383350"               
##  [288] "Optimizing the root position on the branch 286 ... 0.0094153413"               
##  [289] "Optimizing the root position on the branch 287 ... 0.0108249279"               
##  [290] "Optimizing the root position on the branch 288 ... 0.0113750314"               
##  [291] "Optimizing the root position on the branch 289 ... 0.0138778828"               
##  [292] "Optimizing the root position on the branch 290 ... 0.0142375410"               
##  [293] "Optimizing the root position on the branch 291 ... 0.0113867421"               
##  [294] "Optimizing the root position on the branch 292 ... 0.0115069795"               
##  [295] "Optimizing the root position on the branch 293 ... 0.0123699907"               
##  [296] "Optimizing the root position on the branch 294 ... 0.0115237247"               
##  [297] "Optimizing the root position on the branch 295 ... 0.0122977559"               
##  [298] "Optimizing the root position on the branch 296 ... 0.0123363253"               
##  [299] "Optimizing the root position on the branch 297 ... 0.0123371899"               
##  [300] "Optimizing the root position on the branch 298 ... 0.0094147706"               
##  [301] "Optimizing the root position on the branch 299 ... 0.0094648331"               
##  [302] "Optimizing the root position on the branch 300 ... 0.0094827183"               
##  [303] "Optimizing the root position on the branch 301 ... 0.0095234635"               
##  [304] "Optimizing the root position on the branch 302 ... 0.0103778730"               
##  [305] "Optimizing the root position on the branch 303 ... 0.0095261732"               
##  [306] "Optimizing the root position on the branch 304 ... 0.0109030576"               
##  [307] "Optimizing the root position on the branch 305 ... 0.0134048541"               
##  [308] "Optimizing the root position on the branch 306 ... 0.0135799392"               
##  [309] "Optimizing the root position on the branch 307 ... 0.0135601256"               
##  [310] "Optimizing the root position on the branch 308 ... 0.0109868659"               
##  [311] "Optimizing the root position on the branch 309 ... 0.0114281738"               
##  [312] "Optimizing the root position on the branch 310 ... 0.0114089508"               
##  [313] "Optimizing the root position on the branch 311 ... 0.0094723078"               
##  [314] "Optimizing the root position on the branch 312 ... 0.0095504635"               
##  [315] "Optimizing the root position on the branch 313 ... 0.0110563398"               
##  [316] "Optimizing the root position on the branch 314 ... 0.0119149097"               
##  [317] "Optimizing the root position on the branch 315 ... 0.0120244575"               
##  [318] "Optimizing the root position on the branch 316 ... 0.0125879890"               
##  [319] "Optimizing the root position on the branch 317 ... 0.0126127014"               
##  [320] "Optimizing the root position on the branch 318 ... 0.0126445427"               
##  [321] "Optimizing the root position on the branch 319 ... 0.0126699620"               
##  [322] "Optimizing the root position on the branch 320 ... 0.0126725874"               
##  [323] "Optimizing the root position on the branch 321 ... 0.0126471584"               
##  [324] "Optimizing the root position on the branch 322 ... 0.0120170471"               
##  [325] "Optimizing the root position on the branch 323 ... 0.0121732214"               
##  [326] "Optimizing the root position on the branch 324 ... 0.0110374032"               
##  [327] "Optimizing the root position on the branch 325 ... 0.0131219846"               
##  [328] "Optimizing the root position on the branch 326 ... 0.0095688526"               
##  [329] "Optimizing the root position on the branch 327 ... 0.0127051617"               
##  [330] "Optimizing the root position on the branch 328 ... 0.0127607645"               
##  [331] "Optimizing the root position on the branch 329 ... 0.0130509786"               
##  [332] "Optimizing the root position on the branch 330 ... 0.0133716713"               
##  [333] "Optimizing the root position on the branch 331 ... 0.0134631663"               
##  [334] "Optimizing the root position on the branch 332 ... 0.0135151718"               
##  [335] "Optimizing the root position on the branch 333 ... 0.0135119763"               
##  [336] "Optimizing the root position on the branch 334 ... 0.0130943497"               
##  [337] "Optimizing the root position on the branch 335 ... 0.0090320982"               
##  [338] "Optimizing the root position on the branch 336 ... 0.0090441579"               
##  [339] "Optimizing the root position on the branch 337 ... 0.0091292876"               
##  [340] "Optimizing the root position on the branch 338 ... 0.0091288697"               
##  [341] "Optimizing the root position on the branch 339 ... 0.0092260395"               
##  [342] "Optimizing the root position on the branch 340 ... 0.0098760557"               
##  [343] "Optimizing the root position on the branch 341 ... 0.0099749405"               
##  [344] "Optimizing the root position on the branch 342 ... 0.0108388580"               
##  [345] "Optimizing the root position on the branch 343 ... 0.0112531888"               
##  [346] "Optimizing the root position on the branch 344 ... 0.0114903670"               
##  [347] "Optimizing the root position on the branch 345 ... 0.0115976468"               
##  [348] "Optimizing the root position on the branch 346 ... 0.0116078552"               
##  [349] "Optimizing the root position on the branch 347 ... 0.0116124257"               
##  [350] "Optimizing the root position on the branch 348 ... 0.0116182108"               
##  [351] "Optimizing the root position on the branch 349 ... 0.0115893257"               
##  [352] "Optimizing the root position on the branch 350 ... 0.0112650676"               
##  [353] "Optimizing the root position on the branch 351 ... 0.0092193456"               
##  [354] "Optimizing the root position on the branch 352 ... 0.0093530930"               
##  [355] "Optimizing the root position on the branch 353 ... 0.0094642504"               
##  [356] "Optimizing the root position on the branch 354 ... 0.0116933473"               
##  [357] "Optimizing the root position on the branch 355 ... 0.0119296733"               
##  [358] "Optimizing the root position on the branch 356 ... 0.0119978940"               
##  [359] "Optimizing the root position on the branch 357 ... 0.0093448721"               
##  [360] "Optimizing the root position on the branch 358 ... 0.0091337192"               
##  [361] "Optimizing the root position on the branch 359 ... 0.0146407581"               
##  [362] "Optimizing the root position on the branch 360 ... 0.0091182878"               
##  [363] "Optimizing the root position on the branch 361 ... 0.0114046821"               
##  [364] "Optimizing the root position on the branch 362 ... 0.0090381046"               
##  [365] "Optimizing the root position on the branch 363 ... 0.0092816694"               
##  [366] "Optimizing the root position on the branch 364 ... 0.0092960081"               
##  [367] "Optimizing the root position on the branch 365 ... 0.0093499041"               
##  [368] "Optimizing the root position on the branch 366 ... 0.0096849559"               
##  [369] "Optimizing the root position on the branch 367 ... 0.0111015044"               
##  [370] "Optimizing the root position on the branch 368 ... 0.0110515435"               
##  [371] "Optimizing the root position on the branch 369 ... 0.0092962856"               
##  [372] "Optimizing the root position on the branch 370 ... 0.0093248217"               
##  [373] "Optimizing the root position on the branch 371 ... 0.0102954333"               
##  [374] "Optimizing the root position on the branch 372 ... 0.0123930641"               
##  [375] "Optimizing the root position on the branch 373 ... 0.0155482337"               
##  [376] "Optimizing the root position on the branch 374 ... 0.0103100538"               
##  [377] "Optimizing the root position on the branch 375 ... 0.0089851495"               
##  [378] "Optimizing the root position on the branch 376 ... 0.0090939366"               
##  [379] "Optimizing the root position on the branch 377 ... 0.0091496924"               
##  [380] "Optimizing the root position on the branch 378 ... 0.0091616870"               
##  [381] "Optimizing the root position on the branch 379 ... 0.0091892588"               
##  [382] "Optimizing the root position on the branch 380 ... 0.0095165019"               
##  [383] "Optimizing the root position on the branch 381 ... 0.0118953434"               
##  [384] "Optimizing the root position on the branch 382 ... 0.0124647221"               
##  [385] "Optimizing the root position on the branch 383 ... 0.0119005459"               
##  [386] "Optimizing the root position on the branch 384 ... 0.0091984199"               
##  [387] "Optimizing the root position on the branch 385 ... 0.0092242761"               
##  [388] "Optimizing the root position on the branch 386 ... 0.0097027448"               
##  [389] "Optimizing the root position on the branch 387 ... 0.0124199650"               
##  [390] "Optimizing the root position on the branch 388 ... 0.0124530548"               
##  [391] "Optimizing the root position on the branch 389 ... 0.0126976003"               
##  [392] "Optimizing the root position on the branch 390 ... 0.0130481202"               
##  [393] "Optimizing the root position on the branch 391 ... 0.0097388126"               
##  [394] "Optimizing the root position on the branch 392 ... 0.0122977145"               
##  [395] "Optimizing the root position on the branch 393 ... 0.0127725302"               
##  [396] "Optimizing the root position on the branch 394 ... 0.0132814242"               
##  [397] "Optimizing the root position on the branch 395 ... 0.0123100580"               
##  [398] "Optimizing the root position on the branch 396 ... 0.0092263214"               
##  [399] "Optimizing the root position on the branch 397 ... 0.0095211124"               
##  [400] "Optimizing the root position on the branch 398 ... 0.0097318421"               
##  [401] "Optimizing the root position on the branch 399 ... 0.0097483474"               
##  [402] "Optimizing the root position on the branch 400 ... 0.0113848289"               
##  [403] "Optimizing the root position on the branch 401 ... 0.0113658745"               
##  [404] "Optimizing the root position on the branch 402 ... 0.0115377037"               
##  [405] "Optimizing the root position on the branch 403 ... 0.0115809965"               
##  [406] "Optimizing the root position on the branch 404 ... 0.0117734135"               
##  [407] "Optimizing the root position on the branch 405 ... 0.0118942092"               
##  [408] "Optimizing the root position on the branch 406 ... 0.0118988317"               
##  [409] "Optimizing the root position on the branch 407 ... 0.0097641196"               
##  [410] "Optimizing the root position on the branch 408 ... 0.0112158747"               
##  [411] "Optimizing the root position on the branch 409 ... 0.0112260742"               
##  [412] "Optimizing the root position on the branch 410 ... 0.0117125855"               
##  [413] "Optimizing the root position on the branch 411 ... 0.0117664135"               
##  [414] "Optimizing the root position on the branch 412 ... 0.0117787518"               
##  [415] "Optimizing the root position on the branch 413 ... 0.0123447810"               
##  [416] "Optimizing the root position on the branch 414 ... 0.0117802539"               
##  [417] "Optimizing the root position on the branch 415 ... 0.0125131270"               
##  [418] "Optimizing the root position on the branch 416 ... 0.0097027939"               
##  [419] "Optimizing the root position on the branch 417 ... 0.0102186360"               
##  [420] "Optimizing the root position on the branch 418 ... 0.0095176605"               
##  [421] "Optimizing the root position on the branch 419 ... 0.0097017689"               
##  [422] "Optimizing the root position on the branch 420 ... 0.0128692038"               
##  [423] "Optimizing the root position on the branch 421 ... 0.0091602688"               
##  [424] "Optimizing the root position on the branch 422 ... 0.0091784251"               
##  [425] "Optimizing the root position on the branch 423 ... 0.0091918349"               
##  [426] "Optimizing the root position on the branch 424 ... 0.0091905329"               
##  [427] "Optimizing the root position on the branch 425 ... 0.0103651535"               
##  [428] "Optimizing the root position on the branch 426 ... 0.0104491982"               
##  [429] "Optimizing the root position on the branch 427 ... 0.0104798775"               
##  [430] "Optimizing the root position on the branch 428 ... 0.0104780414"               
##  [431] "Optimizing the root position on the branch 429 ... 0.0106148397"               
##  [432] "Optimizing the root position on the branch 430 ... 0.0119368917"               
##  [433] "Optimizing the root position on the branch 431 ... 0.0120776606"               
##  [434] "Optimizing the root position on the branch 432 ... 0.0119188033"               
##  [435] "Optimizing the root position on the branch 433 ... 0.0120883289"               
##  [436] "Optimizing the root position on the branch 434 ... 0.0104798799"               
##  [437] "Optimizing the root position on the branch 435 ... 0.0091839072"               
##  [438] "Optimizing the root position on the branch 436 ... 0.0152281354"               
##  [439] "Optimizing the root position on the branch 437 ... 0.0156276494"               
##  [440] "Optimizing the root position on the branch 438 ... 0.0091770345"               
##  [441] "Optimizing the root position on the branch 439 ... 0.0091813085"               
##  [442] "Optimizing the root position on the branch 440 ... 0.0091816221"               
##  [443] "Optimizing the root position on the branch 441 ... 0.0091595114"               
##  [444] "Optimizing the root position on the branch 442 ... 0.0096135462"               
##  [445] "Optimizing the root position on the branch 443 ... 0.0134485069"               
##  [446] "Optimizing the root position on the branch 444 ... 0.0139027810"               
##  [447] "Optimizing the root position on the branch 445 ... 0.0134262104"               
##  [448] "Optimizing the root position on the branch 446 ... 0.0095884208"               
##  [449] "Optimizing the root position on the branch 447 ... 0.0091749542"               
##  [450] "Optimizing the root position on the branch 448 ... 0.0091746489"               
##  [451] "Optimizing the root position on the branch 449 ... 0.0091617850"               
##  [452] "Optimizing the root position on the branch 450 ... 0.0091577675"               
##  [453] "Optimizing the root position on the branch 451 ... 0.0093711711"               
##  [454] "Optimizing the root position on the branch 452 ... 0.0165103400"               
##  [455] "Optimizing the root position on the branch 453 ... 0.0168844212"               
##  [456] "Optimizing the root position on the branch 454 ... 0.0172958338"               
##  [457] "Optimizing the root position on the branch 455 ... 0.0093548876"               
##  [458] "Optimizing the root position on the branch 456 ... 0.0091578430"               
##  [459] "Optimizing the root position on the branch 457 ... 0.0091568303"               
##  [460] "Optimizing the root position on the branch 458 ... 0.0091897731"               
##  [461] "Optimizing the root position on the branch 459 ... 0.0092801273"               
##  [462] "Optimizing the root position on the branch 460 ... 0.0116070522"               
##  [463] "Optimizing the root position on the branch 461 ... 0.0127570123"               
##  [464] "Optimizing the root position on the branch 462 ... 0.0116611836"               
##  [465] "Optimizing the root position on the branch 463 ... 0.0136874238"               
##  [466] "Optimizing the root position on the branch 464 ... 0.0136270675"               
##  [467] "Optimizing the root position on the branch 465 ... 0.0091497069"               
##  [468] "Optimizing the root position on the branch 466 ... 0.0091542844"               
##  [469] "Optimizing the root position on the branch 467 ... 0.0093150491"               
##  [470] "Optimizing the root position on the branch 468 ... 0.0096155692"               
##  [471] "Optimizing the root position on the branch 469 ... 0.0096242622"               
##  [472] "Optimizing the root position on the branch 470 ... 0.0096130171"               
##  [473] "Optimizing the root position on the branch 471 ... 0.0097629848"               
##  [474] "Optimizing the root position on the branch 472 ... 0.0095487108"               
##  [475] "Optimizing the root position on the branch 473 ... 0.0095969754"               
##  [476] "Optimizing the root position on the branch 474 ... 0.0097127205"               
##  [477] "Optimizing the root position on the branch 475 ... 0.0101571271"               
##  [478] "Optimizing the root position on the branch 476 ... 0.0101803951"               
##  [479] "Optimizing the root position on the branch 477 ... 0.0102664946"               
##  [480] "Optimizing the root position on the branch 478 ... 0.0101733774"               
##  [481] "Optimizing the root position on the branch 479 ... 0.0104492348"               
##  [482] "Optimizing the root position on the branch 480 ... 0.0104556553"               
##  [483] "Optimizing the root position on the branch 481 ... 0.0104572121"               
##  [484] "Optimizing the root position on the branch 482 ... 0.0093256347"               
##  [485] "Optimizing the root position on the branch 483 ... 0.0102386421"               
##  [486] "Optimizing the root position on the branch 484 ... 0.0107528185"               
##  [487] "Optimizing the root position on the branch 485 ... 0.0107851689"               
##  [488] "Optimizing the root position on the branch 486 ... 0.0145970412"               
##  [489] "Optimizing the root position on the branch 487 ... 0.0091555828"               
##  [490] "Optimizing the root position on the branch 488 ... 0.0091470614"               
##  [491] "Optimizing the root position on the branch 489 ... 0.0106674443"               
##  [492] "Optimizing the root position on the branch 490 ... 0.0125264771"               
##  [493] "Optimizing the root position on the branch 491 ... 0.0125546704"               
##  [494] "Optimizing the root position on the branch 492 ... 0.0126312211"               
##  [495] "Optimizing the root position on the branch 493 ... 0.0129470138"               
##  [496] "Optimizing the root position on the branch 494 ... 0.0132737736"               
##  [497] "Optimizing the root position on the branch 495 ... 0.0132942999"               
##  [498] "Optimizing the root position on the branch 496 ... 0.0132991833"               
##  [499] "Optimizing the root position on the branch 497 ... 0.0129659981"               
##  [500] "Optimizing the root position on the branch 498 ... 0.0090898042"               
##  [501] "Optimizing the root position on the branch 499 ... 0.0095704184"               
##  [502] "Optimizing the root position on the branch 500 ... 0.0101759996"               
##  [503] "Optimizing the root position on the branch 501 ... 0.0106541047"               
##  [504] "Optimizing the root position on the branch 502 ... 0.0116197409"               
##  [505] "Optimizing the root position on the branch 503 ... 0.0117878574"               
##  [506] "Optimizing the root position on the branch 504 ... 0.0119642496"               
##  [507] "Optimizing the root position on the branch 505 ... 0.0120035860"               
##  [508] "Optimizing the root position on the branch 506 ... 0.0120401157"               
##  [509] "Optimizing the root position on the branch 507 ... 0.0119854257"               
##  [510] "Optimizing the root position on the branch 508 ... 0.0116050472"               
##  [511] "Optimizing the root position on the branch 509 ... 0.0128562359"               
##  [512] "Optimizing the root position on the branch 510 ... 0.0106737811"               
##  [513] "Optimizing the root position on the branch 511 ... 0.0145444733"               
##  [514] "Optimizing the root position on the branch 512 ... 0.0101551405"               
##  [515] "Optimizing the root position on the branch 513 ... 0.0113714436"               
##  [516] "Optimizing the root position on the branch 514 ... 0.0119463910"               
##  [517] "Optimizing the root position on the branch 515 ... 0.0124278890"               
##  [518] "Optimizing the root position on the branch 516 ... 0.0119646729"               
##  [519] "Optimizing the root position on the branch 517 ... 0.0113127310"               
##  [520] "Optimizing the root position on the branch 518 ... 0.0095891302"               
##  [521] "Optimizing the root position on the branch 519 ... 0.0089511653"               
##  [522] "Optimizing the root position on the branch 520 ... 0.0089563697"               
##  [523] "Optimizing the root position on the branch 521 ... 0.0089502245"               
##  [524] "Optimizing the root position on the branch 522 ... 0.0089510050"               
##  [525] "Optimizing the root position on the branch 523 ... 0.0089763406"               
##  [526] "Optimizing the root position on the branch 524 ... 0.0092084893"               
##  [527] "Optimizing the root position on the branch 525 ... 0.0092057147"               
##  [528] "Optimizing the root position on the branch 526 ... 0.0098398089"               
##  [529] "Optimizing the root position on the branch 527 ... 0.0099395500"               
##  [530] "Optimizing the root position on the branch 528 ... 0.0102051530"               
##  [531] "Optimizing the root position on the branch 529 ... 0.0134093475"               
##  [532] "Optimizing the root position on the branch 530 ... 0.0140669383"               
##  [533] "Optimizing the root position on the branch 531 ... 0.0101940539"               
##  [534] "Optimizing the root position on the branch 532 ... 0.0102693302"               
##  [535] "Optimizing the root position on the branch 533 ... 0.0112261683"               
##  [536] "Optimizing the root position on the branch 534 ... 0.0112650642"               
##  [537] "Optimizing the root position on the branch 535 ... 0.0120261132"               
##  [538] "Optimizing the root position on the branch 536 ... 0.0121542701"               
##  [539] "Optimizing the root position on the branch 537 ... 0.0121659244"               
##  [540] "Optimizing the root position on the branch 538 ... 0.0121670145"               
##  [541] "Optimizing the root position on the branch 539 ... 0.0112605947"               
##  [542] "Optimizing the root position on the branch 540 ... 0.0113140849"               
##  [543] "Optimizing the root position on the branch 541 ... 0.0113458416"               
##  [544] "Optimizing the root position on the branch 542 ... 0.0113717216"               
##  [545] "Optimizing the root position on the branch 543 ... 0.0113788332"               
##  [546] "Optimizing the root position on the branch 544 ... 0.0118266512"               
##  [547] "Optimizing the root position on the branch 545 ... 0.0102628659"               
##  [548] "Optimizing the root position on the branch 546 ... 0.0102760473"               
##  [549] "Optimizing the root position on the branch 547 ... 0.0097866124"               
##  [550] "Optimizing the root position on the branch 548 ... 0.0119680561"               
##  [551] "Optimizing the root position on the branch 549 ... 0.0120785161"               
##  [552] "Optimizing the root position on the branch 550 ... 0.0124026303"               
##  [553] "Optimizing the root position on the branch 551 ... 0.0092105402"               
##  [554] "Optimizing the root position on the branch 552 ... 0.0092027279"               
##  [555] "Optimizing the root position on the branch 553 ... 0.0092001389"               
##  [556] "Optimizing the root position on the branch 554 ... 0.0092427154"               
##  [557] "Optimizing the root position on the branch 555 ... 0.0095746002"               
##  [558] "Optimizing the root position on the branch 556 ... 0.0097628925"               
##  [559] "Optimizing the root position on the branch 557 ... 0.0109164999"               
##  [560] "Optimizing the root position on the branch 558 ... 0.0109556510"               
##  [561] "Optimizing the root position on the branch 559 ... 0.0096070171"               
##  [562] "Optimizing the root position on the branch 560 ... 0.0100807349"               
##  [563] "Optimizing the root position on the branch 561 ... 0.0091921231"               
##  [564] "Optimizing the root position on the branch 562 ... 0.0092074941"               
##  [565] "Optimizing the root position on the branch 563 ... 0.0095450772"               
##  [566] "Optimizing the root position on the branch 564 ... 0.0091737561"               
##  [567] "Optimizing the root position on the branch 565 ... 0.0119341278"               
##  [568] "Optimizing the root position on the branch 566 ... 0.0126664581"               
##  [569] "Optimizing the root position on the branch 567 ... 0.0126634575"               
##  [570] "Optimizing the root position on the branch 568 ... 0.0089803708"               
##  [571] "Optimizing the root position on the branch 569 ... 0.0089955596"               
##  [572] "Optimizing the root position on the branch 570 ... 0.0094483436"               
##  [573] "Optimizing the root position on the branch 571 ... 0.0096096823"               
##  [574] "Optimizing the root position on the branch 572 ... 0.0096201123"               
##  [575] "Optimizing the root position on the branch 573 ... 0.0095697032"               
##  [576] "Optimizing the root position on the branch 574 ... 0.0095765386"               
##  [577] "Optimizing the root position on the branch 575 ... 0.0164897129"               
##  [578] "Optimizing the root position on the branch 576 ... 0.0174402462"               
##  [579] "Optimizing the root position on the branch 577 ... 0.0165098981"               
##  [580] "Optimizing the root position on the branch 578 ... 0.0173148396"               
##  [581] "Optimizing the root position on the branch 579 ... 0.0173926985"               
##  [582] "Optimizing the root position on the branch 580 ... 0.0174222085"               
##  [583] "Optimizing the root position on the branch 581 ... 0.0096152492"               
##  [584] "Optimizing the root position on the branch 582 ... 0.0097254789"               
##  [585] "Optimizing the root position on the branch 583 ... 0.0099180377"               
##  [586] "Optimizing the root position on the branch 584 ... 0.0099950201"               
##  [587] "Optimizing the root position on the branch 585 ... 0.0106982848"               
##  [588] "Optimizing the root position on the branch 586 ... 0.0110046385"               
##  [589] "Optimizing the root position on the branch 587 ... 0.0110440478"               
##  [590] "Optimizing the root position on the branch 588 ... 0.0112706519"               
##  [591] "Optimizing the root position on the branch 589 ... 0.0113611176"               
##  [592] "Optimizing the root position on the branch 590 ... 0.0113601245"               
##  [593] "Optimizing the root position on the branch 591 ... 0.0118992957"               
##  [594] "Optimizing the root position on the branch 592 ... 0.0110035655"               
##  [595] "Optimizing the root position on the branch 593 ... 0.0096035650"               
##  [596] "Optimizing the root position on the branch 594 ... 0.0096440100"               
##  [597] "Optimizing the root position on the branch 595 ... 0.0096472416"               
##  [598] "Optimizing the root position on the branch 596 ... 0.0101405342"               
##  [599] "Optimizing the root position on the branch 597 ... 0.0107050095"               
##  [600] "Optimizing the root position on the branch 598 ... 0.0107956532"               
##  [601] "Optimizing the root position on the branch 599 ... 0.0115970185"               
##  [602] "Optimizing the root position on the branch 600 ... 0.0115920865"               
##  [603] "Optimizing the root position on the branch 601 ... 0.0117101936"               
##  [604] "Optimizing the root position on the branch 602 ... 0.0107113388"               
##  [605] "Optimizing the root position on the branch 603 ... 0.0120734862"               
##  [606] "Optimizing the root position on the branch 604 ... 0.0121751753"               
##  [607] "Optimizing the root position on the branch 605 ... 0.0122046284"               
##  [608] "Optimizing the root position on the branch 606 ... 0.0100637656"               
##  [609] "Optimizing the root position on the branch 607 ... 0.0127348032"               
##  [610] "Optimizing the root position on the branch 608 ... 0.0096470638"               
##  [611] "Optimizing the root position on the branch 609 ... 0.0096727548"               
##  [612] "Optimizing the root position on the branch 610 ... 0.0099450124"               
##  [613] "Optimizing the root position on the branch 611 ... 0.0100705915"               
##  [614] "Optimizing the root position on the branch 612 ... 0.0099199236"               
##  [615] "Optimizing the root position on the branch 613 ... 0.0107477702"               
##  [616] "Optimizing the root position on the branch 614 ... 0.0109459004"               
##  [617] "Optimizing the root position on the branch 615 ... 0.0109264596"               
##  [618] "Optimizing the root position on the branch 616 ... 0.0109574619"               
##  [619] "Optimizing the root position on the branch 617 ... 0.0110167386"               
##  [620] "Optimizing the root position on the branch 618 ... 0.0110264447"               
##  [621] "Optimizing the root position on the branch 619 ... 0.0110627070"               
##  [622] "Optimizing the root position on the branch 620 ... 0.0110643270"               
##  [623] "Optimizing the root position on the branch 621 ... 0.0110158646"               
##  [624] "Optimizing the root position on the branch 622 ... 0.0120705773"               
##  [625] "Optimizing the root position on the branch 623 ... 0.0110201834"               
##  [626] "Optimizing the root position on the branch 624 ... 0.0115856642"               
##  [627] "Optimizing the root position on the branch 625 ... 0.0109495966"               
##  [628] "Optimizing the root position on the branch 626 ... 0.0111481241"               
##  [629] "Optimizing the root position on the branch 627 ... 0.0111959963"               
##  [630] "Optimizing the root position on the branch 628 ... 0.0111988185"               
##  [631] "Optimizing the root position on the branch 629 ... 0.0111979269"               
##  [632] "Optimizing the root position on the branch 630 ... 0.0112276424"               
##  [633] "Optimizing the root position on the branch 631 ... 0.0112769189"               
##  [634] "Optimizing the root position on the branch 632 ... 0.0096834682"               
##  [635] "Optimizing the root position on the branch 633 ... 0.0127924165"               
##  [636] "Optimizing the root position on the branch 634 ... 0.0132436036"               
##  [637] "Optimizing the root position on the branch 635 ... 0.0135232077"               
##  [638] "Optimizing the root position on the branch 636 ... 0.0135562483"               
##  [639] "Optimizing the root position on the branch 637 ... 0.0096066558"               
##  [640] "Optimizing the root position on the branch 638 ... 0.0095827600"               
##  [641] "Optimizing the root position on the branch 639 ... 0.0096540188"               
##  [642] "Optimizing the root position on the branch 640 ... 0.0096632959"               
##  [643] "Optimizing the root position on the branch 641 ... 0.0096689382"               
##  [644] "Optimizing the root position on the branch 642 ... 0.0111901994"               
##  [645] "Optimizing the root position on the branch 643 ... 0.0132403871"               
##  [646] "Optimizing the root position on the branch 644 ... 0.0136298194"               
##  [647] "Optimizing the root position on the branch 645 ... 0.0096679127"               
##  [648] "Optimizing the root position on the branch 646 ... 0.0101313215"               
##  [649] "Optimizing the root position on the branch 647 ... 0.0102196310"               
##  [650] "Optimizing the root position on the branch 648 ... 0.0132330152"               
##  [651] "Optimizing the root position on the branch 649 ... 0.0133388693"               
##  [652] "Optimizing the root position on the branch 650 ... 0.0133574675"               
##  [653] "Optimizing the root position on the branch 651 ... 0.0133609305"               
##  [654] "Optimizing the root position on the branch 652 ... 0.0101858701"               
##  [655] "Optimizing the root position on the branch 653 ... 0.0103842211"               
##  [656] "Optimizing the root position on the branch 654 ... 0.0104205363"               
##  [657] "Optimizing the root position on the branch 655 ... 0.0104268134"               
##  [658] "Optimizing the root position on the branch 656 ... 0.0118987043"               
##  [659] "Optimizing the root position on the branch 657 ... 0.0119781255"               
##  [660] "Optimizing the root position on the branch 658 ... 0.0118238856"               
##  [661] "Optimizing the root position on the branch 659 ... 0.0096622464"               
##  [662] "Optimizing the root position on the branch 660 ... 0.0096731781"               
##  [663] "Optimizing the root position on the branch 661 ... 0.0096843231"               
##  [664] "Optimizing the root position on the branch 662 ... 0.0102941119"               
##  [665] "Optimizing the root position on the branch 663 ... 0.0106914027"               
##  [666] "Optimizing the root position on the branch 664 ... 0.0122320208"               
##  [667] "Optimizing the root position on the branch 665 ... 0.0126107329"               
##  [668] "Optimizing the root position on the branch 666 ... 0.0126503682"               
##  [669] "Optimizing the root position on the branch 667 ... 0.0127074368"               
##  [670] "Optimizing the root position on the branch 668 ... 0.0107174042"               
##  [671] "Optimizing the root position on the branch 669 ... 0.0110153070"               
##  [672] "Optimizing the root position on the branch 670 ... 0.0116185285"               
##  [673] "Optimizing the root position on the branch 671 ... 0.0117845043"               
##  [674] "Optimizing the root position on the branch 672 ... 0.0118098437"               
##  [675] "Optimizing the root position on the branch 673 ... 0.0118158254"               
##  [676] "Optimizing the root position on the branch 674 ... 0.0118186043"               
##  [677] "Optimizing the root position on the branch 675 ... 0.0118593894"               
##  [678] "Optimizing the root position on the branch 676 ... 0.0119039565"               
##  [679] "Optimizing the root position on the branch 677 ... 0.0119485459"               
##  [680] "Optimizing the root position on the branch 678 ... 0.0119723539"               
##  [681] "Optimizing the root position on the branch 679 ... 0.0118592243"               
##  [682] "Optimizing the root position on the branch 680 ... 0.0118206881"               
##  [683] "Optimizing the root position on the branch 681 ... 0.0118775973"               
##  [684] "Optimizing the root position on the branch 682 ... 0.0117958985"               
##  [685] "Optimizing the root position on the branch 683 ... 0.0110070340"               
##  [686] "Optimizing the root position on the branch 684 ... 0.0110984603"               
##  [687] "Optimizing the root position on the branch 685 ... 0.0096726168"               
##  [688] "Optimizing the root position on the branch 686 ... 0.0105941726"               
##  [689] "Optimizing the root position on the branch 687 ... 0.0130518665"               
##  [690] "Optimizing the root position on the branch 688 ... 0.0105685439"               
##  [691] "Optimizing the root position on the branch 689 ... 0.0125603916"               
##  [692] "Optimizing the root position on the branch 690 ... 0.0127232171"               
##  [693] "Optimizing the root position on the branch 691 ... 0.0125420835"               
##  [694] "Optimizing the root position on the branch 692 ... 0.0135815585"               
##  [695] "Optimizing the root position on the branch 693 ... 0.0135583359"               
##  [696] "Optimizing the root position on the branch 694 ... 0.0090016514"               
##  [697] "Optimizing the root position on the branch 695 ... 0.0093585766"               
##  [698] "Optimizing the root position on the branch 696 ... 0.0154703267"               
##  [699] "Optimizing the root position on the branch 697 ... 0.0160190891"               
##  [700] "Optimizing the root position on the branch 698 ... 0.0172242083"               
##  [701] "Optimizing the root position on the branch 699 ... 0.0173283037"               
##  [702] "Optimizing the root position on the branch 700 ... 0.0173234109"               
##  [703] "Optimizing the root position on the branch 701 ... 0.0160162869"               
##  [704] "Optimizing the root position on the branch 702 ... 0.0094126227"               
##  [705] "Optimizing the root position on the branch 703 ... 0.0095863620"               
##  [706] "Optimizing the root position on the branch 704 ... 0.0096657694"               
##  [707] "Optimizing the root position on the branch 705 ... 0.0096834535"               
##  [708] "Optimizing the root position on the branch 706 ... 0.0110594120"               
##  [709] "Optimizing the root position on the branch 707 ... 0.0130681413"               
##  [710] "Optimizing the root position on the branch 708 ... 0.0134124816"               
##  [711] "Optimizing the root position on the branch 709 ... 0.0135360478"               
##  [712] "Optimizing the root position on the branch 710 ... 0.0135437985"               
##  [713] "Optimizing the root position on the branch 711 ... 0.0131069701"               
##  [714] "Optimizing the root position on the branch 712 ... 0.0136240725"               
##  [715] "Optimizing the root position on the branch 713 ... 0.0136672676"               
##  [716] "Optimizing the root position on the branch 714 ... 0.0136979787"               
##  [717] "Optimizing the root position on the branch 715 ... 0.0096593326"               
##  [718] "Optimizing the root position on the branch 716 ... 0.0144875466"               
##  [719] "Optimizing the root position on the branch 717 ... 0.0147307459"               
##  [720] "Optimizing the root position on the branch 718 ... 0.0089499809"               
##  [721] "Optimizing the root position on the branch 719 ... 0.0089508606"               
##  [722] "Optimizing the root position on the branch 720 ... 0.0089549782"               
##  [723] "Optimizing the root position on the branch 721 ... 0.0092129218"               
##  [724] "Optimizing the root position on the branch 722 ... 0.0095874230"               
##  [725] "Optimizing the root position on the branch 723 ... 0.0105791461"               
##  [726] "Optimizing the root position on the branch 724 ... 0.0125358974"               
##  [727] "Optimizing the root position on the branch 725 ... 0.0127574622"               
##  [728] "Optimizing the root position on the branch 726 ... 0.0128139385"               
##  [729] "Optimizing the root position on the branch 727 ... 0.0129732589"               
##  [730] "Optimizing the root position on the branch 728 ... 0.0129890494"               
##  [731] "Optimizing the root position on the branch 729 ... 0.0129964065"               
##  [732] "Optimizing the root position on the branch 730 ... 0.0130049832"               
##  [733] "Optimizing the root position on the branch 731 ... 0.0130033119"               
##  [734] "Optimizing the root position on the branch 732 ... 0.0128085018"               
##  [735] "Optimizing the root position on the branch 733 ... 0.0129397105"               
##  [736] "Optimizing the root position on the branch 734 ... 0.0130840233"               
##  [737] "Optimizing the root position on the branch 735 ... 0.0131390067"               
##  [738] "Optimizing the root position on the branch 736 ... 0.0131574814"               
##  [739] "Optimizing the root position on the branch 737 ... 0.0106445140"               
##  [740] "Optimizing the root position on the branch 738 ... 0.0115088944"               
##  [741] "Optimizing the root position on the branch 739 ... 0.0114670568"               
##  [742] "Optimizing the root position on the branch 740 ... 0.0095520957"               
##  [743] "Optimizing the root position on the branch 741 ... 0.0097434061"               
##  [744] "Optimizing the root position on the branch 742 ... 0.0092479048"               
##  [745] "Optimizing the root position on the branch 743 ... 0.0089568969"               
##  [746] "Optimizing the root position on the branch 744 ... 0.0140786926"               
##  [747] "Optimizing the root position on the branch 745 ... 0.0146907496"               
##  [748] "Optimizing the root position on the branch 746 ... 0.0157102429"               
##  [749] "Optimizing the root position on the branch 747 ... 0.0156815208"               
##  [750] "Optimizing the root position on the branch 748 ... 0.0146346020"               
##  [751] "Optimizing the root position on the branch 749 ... 0.0141418231"               
##  [752] "Optimizing the root position on the branch 750 ... 0.0179697986"               
##  [753] "Optimizing the root position on the branch 751 ... 0.0182763369"               
##  [754] "Optimizing the root position on the branch 752 ... 0.0179729462"               
##  [755] "Optimizing the root position on the branch 753 ... 0.0184976795"               
##  [756] "Optimizing the root position on the branch 754 ... 0.0194413189"               
##  [757] "Optimizing the root position on the branch 755 ... 0.0195166210"               
##  [758] "Optimizing the root position on the branch 756 ... 0.0089504556"               
##  [759] "Optimizing the root position on the branch 757 ... 0.0196572126"               
##  [760] "Optimizing the root position on the branch 758 ... 0.0204068782"               
##  [761] "Optimizing the root position on the branch 759 ... 0.0089433634"               
##  [762] "Optimizing the root position on the branch 760 ... 0.0093068477"               
##  [763] "Optimizing the root position on the branch 761 ... 0.0092987398"               
##  [764] "Optimizing the root position on the branch 762 ... 0.0093563773"               
##  [765] "Optimizing the root position on the branch 763 ... 0.0097576328"               
##  [766] "Optimizing the root position on the branch 764 ... 0.0097480774"               
##  [767] "Optimizing the root position on the branch 765 ... 0.0093625679"               
##  [768] "Optimizing the root position on the branch 766 ... 0.0102725737"               
##  [769] "Optimizing the root position on the branch 767 ... 0.0092833408"               
##  [770] "Optimizing the root position on the branch 768 ... 0.0092902478"               
##  [771] "Optimizing the root position on the branch 769 ... 0.0094163673"               
##  [772] "Optimizing the root position on the branch 770 ... 0.0100382213"               
##  [773] "Optimizing the root position on the branch 771 ... 0.0113352271"               
##  [774] "Optimizing the root position on the branch 772 ... 0.0114272324"               
##  [775] "Optimizing the root position on the branch 773 ... 0.0117077776"               
##  [776] "Optimizing the root position on the branch 774 ... 0.0117159750"               
##  [777] "Optimizing the root position on the branch 775 ... 0.0117379343"               
##  [778] "Optimizing the root position on the branch 776 ... 0.0117989010"               
##  [779] "Optimizing the root position on the branch 777 ... 0.0117336366"               
##  [780] "Optimizing the root position on the branch 778 ... 0.0119036194"               
##  [781] "Optimizing the root position on the branch 779 ... 0.0119779901"               
##  [782] "Optimizing the root position on the branch 780 ... 0.0119946969"               
##  [783] "Optimizing the root position on the branch 781 ... 0.0121592432"               
##  [784] "Optimizing the root position on the branch 782 ... 0.0121642332"               
##  [785] "Optimizing the root position on the branch 783 ... 0.0122852274"               
##  [786] "Optimizing the root position on the branch 784 ... 0.0117170686"               
##  [787] "Optimizing the root position on the branch 785 ... 0.0114714955"               
##  [788] "Optimizing the root position on the branch 786 ... 0.0116145349"               
##  [789] "Optimizing the root position on the branch 787 ... 0.0120561662"               
##  [790] "Optimizing the root position on the branch 788 ... 0.0121091503"               
##  [791] "Optimizing the root position on the branch 789 ... 0.0122811810"               
##  [792] "Optimizing the root position on the branch 790 ... 0.0123547000"               
##  [793] "Optimizing the root position on the branch 791 ... 0.0121240331"               
##  [794] "Optimizing the root position on the branch 792 ... 0.0116371119"               
##  [795] "Optimizing the root position on the branch 793 ... 0.0121219026"               
##  [796] "Optimizing the root position on the branch 794 ... 0.0122610198"               
##  [797] "Optimizing the root position on the branch 795 ... 0.0123454704"               
##  [798] "Optimizing the root position on the branch 796 ... 0.0126417786"               
##  [799] "Optimizing the root position on the branch 797 ... 0.0126539375"               
##  [800] "Optimizing the root position on the branch 798 ... 0.0121133095"               
##  [801] "Optimizing the root position on the branch 799 ... 0.0123406254"               
##  [802] "Optimizing the root position on the branch 800 ... 0.0113003888"               
##  [803] "Optimizing the root position on the branch 801 ... 0.0120772404"               
##  [804] "Optimizing the root position on the branch 802 ... 0.0129104483"               
##  [805] "Optimizing the root position on the branch 803 ... 0.0120460755"               
##  [806] "Optimizing the root position on the branch 804 ... 0.0094320597"               
##  [807] "Optimizing the root position on the branch 805 ... 0.0099820262"               
##  [808] "Optimizing the root position on the branch 806 ... 0.0101768845"               
##  [809] "Optimizing the root position on the branch 807 ... 0.0111795178"               
##  [810] "Optimizing the root position on the branch 808 ... 0.0098668638"               
##  [811] "Optimizing the root position on the branch 809 ... 0.0092867681"               
##  [812] "Optimizing the root position on the branch 810 ... 0.0095175234"               
##  [813] "Optimizing the root position on the branch 811 ... 0.0094936520"               
##  [814] "Optimizing the root position on the branch 812 ... 0.0094954106"               
##  [815] "Optimizing the root position on the branch 813 ... 0.0100404977"               
##  [816] "Optimizing the root position on the branch 814 ... 0.0102547226"               
##  [817] "Optimizing the root position on the branch 815 ... 0.0100177190"               
##  [818] "Optimizing the root position on the branch 816 ... 0.0095247678"               
##  [819] "Optimizing the root position on the branch 817 ... 0.0096178518"               
##  [820] "Optimizing the root position on the branch 818 ... 0.0089569632"               
##  [821] "Optimizing the root position on the branch 819 ... 0.0099492269"               
##  [822] "Optimizing the root position on the branch 820 ... 0.0100716095"               
##  [823] "Optimizing the root position on the branch 821 ... 0.0101394985"               
##  [824] "Optimizing the root position on the branch 822 ... 0.0121086264"               
##  [825] "Optimizing the root position on the branch 823 ... 0.0100775627"               
##  [826] "Optimizing the root position on the branch 824 ... 0.0129791703"               
##  [827] "Optimizing the root position on the branch 825 ... 0.0099949115"               
##  [828] "Optimizing the root position on the branch 826 ... 0.0090585938"               
##  [829] "Optimizing the root position on the branch 827 ... 0.0091687148"               
##  [830] "Optimizing the root position on the branch 828 ... 0.0095096938"               
##  [831] "Optimizing the root position on the branch 829 ... 0.0095306824"               
##  [832] "Optimizing the root position on the branch 830 ... 0.0095636141"               
##  [833] "Optimizing the root position on the branch 831 ... 0.0096861674"               
##  [834] "Optimizing the root position on the branch 832 ... 0.0105866177"               
##  [835] "Optimizing the root position on the branch 833 ... 0.0109264943"               
##  [836] "Optimizing the root position on the branch 834 ... 0.0109107672"               
##  [837] "Optimizing the root position on the branch 835 ... 0.0111501901"               
##  [838] "Optimizing the root position on the branch 836 ... 0.0119666268"               
##  [839] "Optimizing the root position on the branch 837 ... 0.0120469706"               
##  [840] "Optimizing the root position on the branch 838 ... 0.0105170287"               
##  [841] "Optimizing the root position on the branch 839 ... 0.0109974679"               
##  [842] "Optimizing the root position on the branch 840 ... 0.0118646993"               
##  [843] "Optimizing the root position on the branch 841 ... 0.0096821594"               
##  [844] "Optimizing the root position on the branch 842 ... 0.0098977971"               
##  [845] "Optimizing the root position on the branch 843 ... 0.0099154625"               
##  [846] "Optimizing the root position on the branch 844 ... 0.0106568877"               
##  [847] "Optimizing the root position on the branch 845 ... 0.0098948163"               
##  [848] "Optimizing the root position on the branch 846 ... 0.0095742927"               
##  [849] "Optimizing the root position on the branch 847 ... 0.0181121677"               
##  [850] "Optimizing the root position on the branch 848 ... 0.0095266093"               
##  [851] "Optimizing the root position on the branch 849 ... 0.0095379882"               
##  [852] "Optimizing the root position on the branch 850 ... 0.0095526899"               
##  [853] "Optimizing the root position on the branch 851 ... 0.0153781562"               
##  [854] "Optimizing the root position on the branch 852 ... 0.0162644245"               
##  [855] "Optimizing the root position on the branch 853 ... 0.0095428804"               
##  [856] "Optimizing the root position on the branch 854 ... 0.0130318085"               
##  [857] "Optimizing the root position on the branch 855 ... 0.0095426435"               
##  [858] "Optimizing the root position on the branch 856 ... 0.0097611942"               
##  [859] "Optimizing the root position on the branch 857 ... 0.0135889290"               
##  [860] "Optimizing the root position on the branch 858 ... 0.0145333233"               
##  [861] "Optimizing the root position on the branch 859 ... 0.0145592042"               
##  [862] "Optimizing the root position on the branch 860 ... 0.0097790544"               
##  [863] "Optimizing the root position on the branch 861 ... 0.0094853965"               
##  [864] "Optimizing the root position on the branch 862 ... 0.0091452588"               
##  [865] "Optimizing the root position on the branch 863 ... 0.0108135923"               
##  [866] "Optimizing the root position on the branch 864 ... 0.0108252352"               
##  [867] "Optimizing the root position on the branch 865 ... 0.0111749962"               
##  [868] "Optimizing the root position on the branch 866 ... 0.0117672807"               
##  [869] "Optimizing the root position on the branch 867 ... 0.0124733503"               
##  [870] "Optimizing the root position on the branch 868 ... 0.0131144994"               
##  [871] "Optimizing the root position on the branch 869 ... 0.0090985168"               
##  [872] "Optimizing the root position on the branch 870 ... 0.0091489955"               
##  [873] "Optimizing the root position on the branch 871 ... 0.0092245982"               
##  [874] "Optimizing the root position on the branch 872 ... 0.0092451712"               
##  [875] "Optimizing the root position on the branch 873 ... 0.0092719044"               
##  [876] "Optimizing the root position on the branch 874 ... 0.0093462695"               
##  [877] "Optimizing the root position on the branch 875 ... 0.0101220137"               
##  [878] "Optimizing the root position on the branch 876 ... 0.0103951654"               
##  [879] "Optimizing the root position on the branch 877 ... 0.0104152674"               
##  [880] "Optimizing the root position on the branch 878 ... 0.0104368712"               
##  [881] "Optimizing the root position on the branch 879 ... 0.0107123630"               
##  [882] "Optimizing the root position on the branch 880 ... 0.0107181687"               
##  [883] "Optimizing the root position on the branch 881 ... 0.0112062311"               
##  [884] "Optimizing the root position on the branch 882 ... 0.0116780381"               
##  [885] "Optimizing the root position on the branch 883 ... 0.0107114752"               
##  [886] "Optimizing the root position on the branch 884 ... 0.0107991331"               
##  [887] "Optimizing the root position on the branch 885 ... 0.0108075351"               
##  [888] "Optimizing the root position on the branch 886 ... 0.0108183615"               
##  [889] "Optimizing the root position on the branch 887 ... 0.0111712965"               
##  [890] "Optimizing the root position on the branch 888 ... 0.0108161906"               
##  [891] "Optimizing the root position on the branch 889 ... 0.0108085181"               
##  [892] "Optimizing the root position on the branch 890 ... 0.0109219643"               
##  [893] "Optimizing the root position on the branch 891 ... 0.0109272080"               
##  [894] "Optimizing the root position on the branch 892 ... 0.0103508136"               
##  [895] "Optimizing the root position on the branch 893 ... 0.0129996175"               
##  [896] "Optimizing the root position on the branch 894 ... 0.0101531523"               
##  [897] "Optimizing the root position on the branch 895 ... 0.0092849448"               
##  [898] "Optimizing the root position on the branch 896 ... 0.0092946674"               
##  [899] "Optimizing the root position on the branch 897 ... 0.0094143559"               
##  [900] "Optimizing the root position on the branch 898 ... 0.0094527509"               
##  [901] "Optimizing the root position on the branch 899 ... 0.0094533525"               
##  [902] "Optimizing the root position on the branch 900 ... 0.0110691821"               
##  [903] "Optimizing the root position on the branch 901 ... 0.0117023513"               
##  [904] "Optimizing the root position on the branch 902 ... 0.0118227929"               
##  [905] "Optimizing the root position on the branch 903 ... 0.0118501675"               
##  [906] "Optimizing the root position on the branch 904 ... 0.0118816208"               
##  [907] "Optimizing the root position on the branch 905 ... 0.0119144362"               
##  [908] "Optimizing the root position on the branch 906 ... 0.0119202176"               
##  [909] "Optimizing the root position on the branch 907 ... 0.0118794543"               
##  [910] "Optimizing the root position on the branch 908 ... 0.0118517398"               
##  [911] "Optimizing the root position on the branch 909 ... 0.0118384320"               
##  [912] "Optimizing the root position on the branch 910 ... 0.0116760550"               
##  [913] "Optimizing the root position on the branch 911 ... 0.0126074479"               
##  [914] "Optimizing the root position on the branch 912 ... 0.0094526925"               
##  [915] "Optimizing the root position on the branch 913 ... 0.0095641528"               
##  [916] "Optimizing the root position on the branch 914 ... 0.0096823426"               
##  [917] "Optimizing the root position on the branch 915 ... 0.0110242974"               
##  [918] "Optimizing the root position on the branch 916 ... 0.0111411925"               
##  [919] "Optimizing the root position on the branch 917 ... 0.0115332703"               
##  [920] "Optimizing the root position on the branch 918 ... 0.0115797532"               
##  [921] "Optimizing the root position on the branch 919 ... 0.0115794709"               
##  [922] "Optimizing the root position on the branch 920 ... 0.0117080009"               
##  [923] "Optimizing the root position on the branch 921 ... 0.0096449675"               
##  [924] "Optimizing the root position on the branch 922 ... 0.0140697856"               
##  [925] "Optimizing the root position on the branch 923 ... 0.0095582370"               
##  [926] "Optimizing the root position on the branch 924 ... 0.0096059207"               
##  [927] "Optimizing the root position on the branch 925 ... 0.0138095378"               
##  [928] "Optimizing the root position on the branch 926 ... 0.0140717466"               
##  [929] "Optimizing the root position on the branch 927 ... 0.0096004608"               
##  [930] "Optimizing the root position on the branch 928 ... 0.0108714009"               
##  [931] "Optimizing the root position on the branch 929 ... 0.0094301945"               
##  [932] "Optimizing the root position on the branch 930 ... 0.0112976064"               
##  [933] "Optimizing the root position on the branch 931 ... 0.0114139429"               
##  [934] "Optimizing the root position on the branch 932 ... 0.0114688626"               
##  [935] "Optimizing the root position on the branch 933 ... 0.0116282110"               
##  [936] "Optimizing the root position on the branch 934 ... 0.0117195291"               
##  [937] "Optimizing the root position on the branch 935 ... 0.0114557729"               
##  [938] "Optimizing the root position on the branch 936 ... 0.0112873714"               
##  [939] "Optimizing the root position on the branch 937 ... 0.0118394841"               
##  [940] "Optimizing the root position on the branch 938 ... 0.0094158324"               
##  [941] "Optimizing the root position on the branch 939 ... 0.0092979954"               
##  [942] "Optimizing the root position on the branch 940 ... 0.0092878075"               
##  [943] "Optimizing the root position on the branch 941 ... 0.0092736556"               
##  [944] "Optimizing the root position on the branch 942 ... 0.0092909890"               
##  [945] "Optimizing the root position on the branch 943 ... 0.0094423865"               
##  [946] "Optimizing the root position on the branch 944 ... 0.0096351005"               
##  [947] "Optimizing the root position on the branch 945 ... 0.0094184718"               
##  [948] "Optimizing the root position on the branch 946 ... 0.0092564143"               
##  [949] "Optimizing the root position on the branch 947 ... 0.0111941601"               
##  [950] "Optimizing the root position on the branch 948 ... 0.0117399824"               
##  [951] "Optimizing the root position on the branch 949 ... 0.0119077971"               
##  [952] "Optimizing the root position on the branch 950 ... 0.0124122369"               
##  [953] "Optimizing the root position on the branch 951 ... 0.0124630827"               
##  [954] "Optimizing the root position on the branch 952 ... 0.0130438308"               
##  [955] "Optimizing the root position on the branch 953 ... 0.0119218124"               
##  [956] "Optimizing the root position on the branch 954 ... 0.0119793103"               
##  [957] "Optimizing the root position on the branch 955 ... 0.0120508879"               
##  [958] "Optimizing the root position on the branch 956 ... 0.0119671449"               
##  [959] "Optimizing the root position on the branch 957 ... 0.0121375746"               
##  [960] "Optimizing the root position on the branch 958 ... 0.0117164879"               
##  [961] "Optimizing the root position on the branch 959 ... 0.0127279316"               
##  [962] "Optimizing the root position on the branch 960 ... 0.0112376619"               
##  [963] "Optimizing the root position on the branch 961 ... 0.0092160081"               
##  [964] "Optimizing the root position on the branch 962 ... 0.0110774689"               
##  [965] "Optimizing the root position on the branch 963 ... 0.0119082229"               
##  [966] "Optimizing the root position on the branch 964 ... 0.0091480617"               
##  [967] "Optimizing the root position on the branch 965 ... 0.0092352783"               
##  [968] "Optimizing the root position on the branch 966 ... 0.0114802755"               
##  [969] "Optimizing the root position on the branch 967 ... 0.0127714380"               
##  [970] "Optimizing the root position on the branch 968 ... 0.0131109261"               
##  [971] "Optimizing the root position on the branch 969 ... 0.0133235850"               
##  [972] "Optimizing the root position on the branch 970 ... 0.0133360096"               
##  [973] "Optimizing the root position on the branch 971 ... 0.0134173441"               
##  [974] "Optimizing the root position on the branch 972 ... 0.0092247899"               
##  [975] "Optimizing the root position on the branch 973 ... 0.0095518543"               
##  [976] "Optimizing the root position on the branch 974 ... 0.0103665163"               
##  [977] "Optimizing the root position on the branch 975 ... 0.0106064870"               
##  [978] "Optimizing the root position on the branch 976 ... 0.0107310071"               
##  [979] "Optimizing the root position on the branch 977 ... 0.0105668711"               
##  [980] "Optimizing the root position on the branch 978 ... 0.0095527307"               
##  [981] "Optimizing the root position on the branch 979 ... 0.0108250552"               
##  [982] "Optimizing the root position on the branch 980 ... 0.0108472424"               
##  [983] "Optimizing the root position on the branch 981 ... 0.0126962731"               
##  [984] "Optimizing the root position on the branch 982 ... 0.0199091959"               
##  [985] "Optimizing the root position on the branch 983 ... 0.0129890471"               
##  [986] "Optimizing the root position on the branch 984 ... 0.0129571021"               
##  [987] "Optimizing the root position on the branch 985 ... 0.0125145273"               
##  [988] "Optimizing the root position on the branch 986 ... 0.0118765946"               
##  [989] "Optimizing the root position on the branch 987 ... 0.0119208932"               
##  [990] "Optimizing the root position on the branch 988 ... 0.0103252389"               
##  [991] "Optimizing the root position on the branch 989 ... 0.0143237843"               
##  [992] "Optimizing the root position on the branch 990 ... 0.0142614703"               
##  [993] "Optimizing the root position on the branch 991 ... 0.0129577183"               
##  [994] "Optimizing the root position on the branch 992 ... 0.0129533386"               
##  [995] "Optimizing the root position on the branch 993 ... 0.0106899129"               
##  [996] "Optimizing the root position on the branch 994 ... 0.0114776537"               
##  [997] "Optimizing the root position on the branch 995 ... 0.0126884609"               
##  [998] "Optimizing the root position on the branch 996 ... 0.0130431647"               
##  [999] "Optimizing the root position on the branch 997 ... 0.0134243251"               
## [1000] "Optimizing the root position on the branch 998 ... 0.0141050977"               
## [1001] "Optimizing the root position on the branch 999 ... 0.0141105470"               
## [1002] "Optimizing the root position on the branch 1000 ... 0.0138195140"              
## [1003] "Optimizing the root position on the branch 1001 ... 0.0138237633"              
## [1004] "Optimizing the root position on the branch 1002 ... 0.0109121895"              
## [1005] "Optimizing the root position on the branch 1003 ... 0.0119753380"              
## [1006] "Optimizing the root position on the branch 1004 ... 0.0122297921"              
## [1007] "Optimizing the root position on the branch 1005 ... 0.0121944296"              
## [1008] "Optimizing the root position on the branch 1006 ... 0.0118713033"              
## [1009] "Optimizing the root position on the branch 1007 ... 0.0118862462"              
## [1010] "Optimizing the root position on the branch 1008 ... 0.0142830811"              
## [1011] "Optimizing the root position on the branch 1009 ... 0.0142879118"              
## [1012] "Optimizing the root position on the branch 1010 ... 0.0126656385"              
## [1013] "Optimizing the root position on the branch 1011 ... 0.0121812377"              
## [1014] "Optimizing the root position on the branch 1012 ... 0.0121721970"              
## [1015] "Optimizing the root position on the branch 1013 ... 0.0121128137"              
## [1016] "Optimizing the root position on the branch 1014 ... 0.0121642008"              
## [1017] "Optimizing the root position on the branch 1015 ... 0.0122036811"              
## [1018] "Optimizing the root position on the branch 1016 ... 0.0120529428"              
## [1019] "Optimizing the root position on the branch 1017 ... 0.0123679623"              
## [1020] "Optimizing the root position on the branch 1018 ... 0.0124487259"              
## [1021] "Optimizing the root position on the branch 1019 ... 0.0130248961"              
## [1022] "Optimizing the root position on the branch 1020 ... 0.0131885697"              
## [1023] "Optimizing the root position on the branch 1021 ... 0.0131766061"              
## [1024] "Optimizing the root position on the branch 1022 ... 0.0128656530"              
## [1025] "Optimizing the root position on the branch 1023 ... 0.0127335500"              
## [1026] "Optimizing the root position on the branch 1024 ... 0.0096125950"              
## [1027] "Optimizing the root position on the branch 1025 ... 0.0097444559"              
## [1028] "Optimizing the root position on the branch 1026 ... 0.0097381362"              
## [1029] "Optimizing the root position on the branch 1027 ... 0.0092580202"              
## [1030] "Optimizing the root position on the branch 1028 ... 0.0174948256"              
## [1031] "Optimizing the root position on the branch 1029 ... 0.0175004167"              
## [1032] "Optimizing the root position on the branch 1030 ... 0.0107643293"              
## [1033] "Optimizing the root position on the branch 1031 ... 0.0106839815"              
## [1034] "Optimizing the root position on the branch 1032 ... 0.0117603039"              
## [1035] "Optimizing the root position on the branch 1033 ... 0.0119529524"              
## [1036] "Optimizing the root position on the branch 1034 ... 0.0119725298"              
## [1037] "Optimizing the root position on the branch 1035 ... 0.0114101053"              
## [1038] "Optimizing the root position on the branch 1036 ... 0.0116944330"              
## [1039] "Optimizing the root position on the branch 1037 ... 0.0116796462"              
## [1040] "Optimizing the root position on the branch 1038 ... 0.0115951925"              
## [1041] "Optimizing the root position on the branch 1039 ... 0.0117069093"              
## [1042] "Optimizing the root position on the branch 1040 ... 0.0128850817"              
## [1043] "Optimizing the root position on the branch 1041 ... 0.0129019616"              
## [1044] "Optimizing the root position on the branch 1042 ... 0.0107110952"              
## [1045] "Optimizing the root position on the branch 1043 ... 0.0110007922"              
## [1046] "Optimizing the root position on the branch 1044 ... 0.0110347087"              
## [1047] "Optimizing the root position on the branch 1045 ... 0.0140158378"              
## [1048] "Optimizing the root position on the branch 1046 ... 0.0143087830"              
## [1049] "Optimizing the root position on the branch 1047 ... 0.0143175142"              
## [1050] "Optimizing the root position on the branch 1048 ... 0.0135848359"              
## [1051] "Optimizing the root position on the branch 1049 ... 0.0139004385"              
## [1052] "Optimizing the root position on the branch 1050 ... 0.0143623450"              
## [1053] "Optimizing the root position on the branch 1051 ... 0.0143364737"              
## [1054] "Optimizing the root position on the branch 1052 ... 0.0109750077"              
## [1055] "Optimizing the root position on the branch 1053 ... 0.0111172364"              
## [1056] "Optimizing the root position on the branch 1054 ... 0.0114902750"              
## [1057] "Optimizing the root position on the branch 1055 ... 0.0117409036"              
## [1058] "Optimizing the root position on the branch 1056 ... 0.0117391231"              
## [1059] "Optimizing the root position on the branch 1057 ... 0.0116985888"              
## [1060] "Optimizing the root position on the branch 1058 ... 0.0117539737"              
## [1061] "Optimizing the root position on the branch 1059 ... 0.0117444592"              
## [1062] "Optimizing the root position on the branch 1060 ... 0.0125952156"              
## [1063] "Optimizing the root position on the branch 1061 ... 0.0127670514"              
## [1064] "Optimizing the root position on the branch 1062 ... 0.0127609337"              
## [1065] "Optimizing the root position on the branch 1063 ... 0.0118904892"              
## [1066] "Optimizing the root position on the branch 1064 ... 0.0118832352"              
## [1067] "Optimizing the root position on the branch 1065 ... 0.0118780894"              
## [1068] "Optimizing the root position on the branch 1066 ... 0.0118643294"              
## [1069] "Optimizing the root position on the branch 1067 ... 0.0118886281"              
## [1070] "Optimizing the root position on the branch 1068 ... 0.0118894861"              
## [1071] "Optimizing the root position on the branch 1069 ... 0.0119110591"              
## [1072] "Optimizing the root position on the branch 1070 ... 0.0119222093"              
## [1073] "Optimizing the root position on the branch 1071 ... 0.0120612279"              
## [1074] "Optimizing the root position on the branch 1072 ... 0.0120632579"              
## [1075] "Optimizing the root position on the branch 1073 ... 0.0109901688"              
## [1076] "Optimizing the root position on the branch 1074 ... 0.0128843617"              
## [1077] "Optimizing the root position on the branch 1075 ... 0.0128911004"              
## [1078] "Optimizing the root position on the branch 1076 ... 0.0132576442"              
## [1079] "Optimizing the root position on the branch 1077 ... 0.0132477933"              
## [1080] "Optimizing the root position on the branch 1078 ... 0.0127486284"              
## [1081] "Optimizing the root position on the branch 1079 ... 0.0106506387"              
## [1082] "Optimizing the root position on the branch 1080 ... 0.0107671698"              
## [1083] "Optimizing the root position on the branch 1081 ... 0.0109372201"              
## [1084] "Optimizing the root position on the branch 1082 ... 0.0109310860"              
## [1085] "Optimizing the root position on the branch 1083 ... 0.0109474088"              
## [1086] "Optimizing the root position on the branch 1084 ... 0.0109534509"              
## [1087] "Optimizing the root position on the branch 1085 ... 0.0110652705"              
## [1088] "Optimizing the root position on the branch 1086 ... 0.0110508103"              
## [1089] "Optimizing the root position on the branch 1087 ... 0.0117629591"              
## [1090] "Optimizing the root position on the branch 1088 ... 0.0117492198"              
## [1091] "Optimizing the root position on the branch 1089 ... 0.0111379025"              
## [1092] "Optimizing the root position on the branch 1090 ... 0.0111655371"              
## [1093] "Optimizing the root position on the branch 1091 ... 0.0116587979"              
## [1094] "Optimizing the root position on the branch 1092 ... 0.0119136890"              
## [1095] "Optimizing the root position on the branch 1093 ... 0.0119065429"              
## [1096] "Optimizing the root position on the branch 1094 ... 0.0104611841"              
## [1097] "Optimizing the root position on the branch 1095 ... 0.0104138882"              
## [1098] "Optimizing the root position on the branch 1096 ... 0.0093080628"              
## [1099] "Optimizing the root position on the branch 1097 ... 0.0110233241"              
## [1100] "Optimizing the root position on the branch 1098 ... 0.0117744462"              
## [1101] "Optimizing the root position on the branch 1099 ... 0.0135551251"              
## [1102] "Optimizing the root position on the branch 1100 ... 0.0135516487"              
## [1103] "Optimizing the root position on the branch 1101 ... 0.0130639700"              
## [1104] "Optimizing the root position on the branch 1102 ... 0.0124738428"              
## [1105] "Optimizing the root position on the branch 1103 ... 0.0117873348"              
## [1106] "Optimizing the root position on the branch 1104 ... 0.0117818068"              
## [1107] "Optimizing the root position on the branch 1105 ... 0.0136138801"              
## [1108] "Optimizing the root position on the branch 1106 ... 0.0135503523"              
## [1109] "Optimizing the root position on the branch 1107 ... 0.0102393255"              
## [1110] "Optimizing the root position on the branch 1108 ... 0.0103441549"              
## [1111] "Optimizing the root position on the branch 1109 ... 0.0132233038"              
## [1112] "Optimizing the root position on the branch 1110 ... 0.0148237845"              
## [1113] "Optimizing the root position on the branch 1111 ... 0.0148309806"              
## [1114] "Optimizing the root position on the branch 1112 ... 0.0146489615"              
## [1115] "Optimizing the root position on the branch 1113 ... 0.0146519870"              
## [1116] "Optimizing the root position on the branch 1114 ... 0.0128897966"              
## [1117] "Optimizing the root position on the branch 1115 ... 0.0132247357"              
## [1118] "Optimizing the root position on the branch 1116 ... 0.0132390863"              
## [1119] "Optimizing the root position on the branch 1117 ... 0.0149825462"              
## [1120] "Optimizing the root position on the branch 1118 ... 0.0161940217"              
## [1121] "Optimizing the root position on the branch 1119 ... 0.0163520503"              
## [1122] "Optimizing the root position on the branch 1120 ... 0.0163383281"              
## [1123] "Optimizing the root position on the branch 1121 ... 0.0183695252"              
## [1124] "Optimizing the root position on the branch 1122 ... 0.0184033742"              
## [1125] "Optimizing the root position on the branch 1123 ... 0.0176202217"              
## [1126] "Optimizing the root position on the branch 1124 ... 0.0099260331"              
## [1127] "Optimizing the root position on the branch 1125 ... 0.0099121957"              
## [1128] "Optimizing the root position on the branch 1126 ... 0.0099178785"              
## [1129] "Optimizing the root position on the branch 1127 ... 0.0106476479"              
## [1130] "Optimizing the root position on the branch 1128 ... 0.0117215648"              
## [1131] "Optimizing the root position on the branch 1129 ... 0.0117224905"              
## [1132] "Optimizing the root position on the branch 1130 ... 0.0109439175"              
## [1133] "Optimizing the root position on the branch 1131 ... 0.0118554161"              
## [1134] "Optimizing the root position on the branch 1132 ... 0.0120787952"              
## [1135] "Optimizing the root position on the branch 1133 ... 0.0120492291"              
## [1136] "Optimizing the root position on the branch 1134 ... 0.0111562220"              
## [1137] "Optimizing the root position on the branch 1135 ... 0.0126376832"              
## [1138] "Optimizing the root position on the branch 1136 ... 0.0126278623"              
## [1139] "Optimizing the root position on the branch 1137 ... 0.0120532972"              
## [1140] "Optimizing the root position on the branch 1138 ... 0.0118781212"              
## [1141] "Optimizing the root position on the branch 1139 ... 0.0115114036"              
## [1142] "Optimizing the root position on the branch 1140 ... 0.0114530339"              
## [1143] "Optimizing the root position on the branch 1141 ... 0.0100793864"              
## [1144] "Optimizing the root position on the branch 1142 ... 0.0100907535"              
## [1145] "Optimizing the root position on the branch 1143 ... 0.0131382613"              
## [1146] "Optimizing the root position on the branch 1144 ... 0.0130986430"              
## [1147] "Optimizing the root position on the branch 1145 ... 0.0126403026"              
## [1148] "Optimizing the root position on the branch 1146 ... 0.0120097449"              
## [1149] "Optimizing the root position on the branch 1147 ... 0.0122878811"              
## [1150] "Optimizing the root position on the branch 1148 ... 0.0122805059"              
## [1151] "Optimizing the root position on the branch 1149 ... 0.0101019219"              
## [1152] "Optimizing the root position on the branch 1150 ... 0.0095747491"              
## [1153] "Optimizing the root position on the branch 1151 ... 0.0162203196"              
## [1154] "Optimizing the root position on the branch 1152 ... 0.0163570164"              
## [1155] "Optimizing the root position on the branch 1153 ... 0.0125106754"              
## [1156] "Optimizing the root position on the branch 1154 ... 0.0126434179"              
## [1157] "Optimizing the root position on the branch 1155 ... 0.0102717622"              
## [1158] "Optimizing the root position on the branch 1156 ... 0.0102501579"              
## [1159] "Optimizing the root position on the branch 1157 ... 0.0102339900"              
## [1160] "Optimizing the root position on the branch 1158 ... 0.0094930319"              
## [1161] "Optimizing the root position on the branch 1159 ... 0.0095067389"              
## [1162] "Optimizing the root position on the branch 1160 ... 0.0137331875"              
## [1163] "Optimizing the root position on the branch 1161 ... 0.0134347774"              
## [1164] "Optimizing the root position on the branch 1162 ... 0.0103028889"              
## [1165] "Optimizing the root position on the branch 1163 ... 0.0111734125"              
## [1166] "Optimizing the root position on the branch 1164 ... 0.0112331774"              
## [1167] "Optimizing the root position on the branch 1165 ... 0.0112974254"              
## [1168] "Optimizing the root position on the branch 1166 ... 0.0100181391"              
## [1169] "Optimizing the root position on the branch 1167 ... 0.0141339750"              
## [1170] "Optimizing the root position on the branch 1168 ... 0.0141628423"              
## [1171] "Optimizing the root position on the branch 1169 ... 0.0128835090"              
## [1172] "Optimizing the root position on the branch 1170 ... 0.0139739243"              
## [1173] "Optimizing the root position on the branch 1171 ... 0.0139878602"              
## [1174] "Optimizing the root position on the branch 1172 ... 0.0123074427"              
## [1175] "Optimizing the root position on the branch 1173 ... 0.0131115312"              
## [1176] "Optimizing the root position on the branch 1174 ... 0.0130974005"              
## [1177] "Optimizing the root position on the branch 1175 ... 0.0122646194"              
## [1178] "Optimizing the root position on the branch 1176 ... 0.0123524822"              
## [1179] "Optimizing the root position on the branch 1177 ... 0.0126070541"              
## [1180] "Optimizing the root position on the branch 1178 ... 0.0126625588"              
## [1181] "Optimizing the root position on the branch 1179 ... 0.0126626004"              
## [1182] "Optimizing the root position on the branch 1180 ... 0.0126526330"              
## [1183] "Optimizing the root position on the branch 1181 ... 0.0120497171"              
## [1184] "Optimizing the root position on the branch 1182 ... 0.0123000626"              
## [1185] "Optimizing the root position on the branch 1183 ... 0.0123228405"              
## [1186] "Optimizing the root position on the branch 1184 ... 0.0123362916"              
## [1187] "Optimizing the root position on the branch 1185 ... 0.0124262971"              
## [1188] "Optimizing the root position on the branch 1186 ... 0.0124196377"              
## [1189] "Optimizing the root position on the branch 1187 ... 0.0122690508"              
## [1190] "Optimizing the root position on the branch 1188 ... 0.0116130782"              
## [1191] "Optimizing the root position on the branch 1189 ... 0.0117343300"              
## [1192] "Optimizing the root position on the branch 1190 ... 0.0117400620"              
## [1193] "Optimizing the root position on the branch 1191 ... 0.0118847935"              
## [1194] "Optimizing the root position on the branch 1192 ... 0.0119824574"              
## [1195] "Optimizing the root position on the branch 1193 ... 0.0119974577"              
## [1196] "Optimizing the root position on the branch 1194 ... 0.0122899361"              
## [1197] "Optimizing the root position on the branch 1195 ... 0.0123242360"              
## [1198] "Optimizing the root position on the branch 1196 ... 0.0123230387"              
## [1199] "Optimizing the root position on the branch 1197 ... 0.0121642470"              
## [1200] "Optimizing the root position on the branch 1198 ... 0.0121383424"              
## [1201] "Optimizing the root position on the branch 1199 ... 0.0129540669"              
## [1202] "Optimizing the root position on the branch 1200 ... 0.0129332805"              
## [1203] "Optimizing the root position on the branch 1201 ... 0.0118114322"              
## [1204] "Optimizing the root position on the branch 1202 ... 0.0101751784"              
## [1205] "Optimizing the root position on the branch 1203 ... 0.0123243214"              
## [1206] "Optimizing the root position on the branch 1204 ... 0.0122043139"              
## [1207] "Optimizing the root position on the branch 1205 ... 0.0152432225"              
## [1208] "Optimizing the root position on the branch 1206 ... 0.0153533706"              
## [1209] "Optimizing the root position on the branch 1207 ... 0.0098133257"              
## [1210] "Optimizing the root position on the branch 1208 ... 0.0098122942"              
## [1211] "Optimizing the root position on the branch 1209 ... 0.0092033997"              
## [1212] "Optimizing the root position on the branch 1210 ... 0.0191530782"              
## [1213] "Optimizing the root position on the branch 1211 ... 0.0203807728"              
## [1214] "Optimizing the root position on the branch 1212 ... 0.0212776576"              
## [1215] "Optimizing the root position on the branch 1213 ... 0.0212852160"              
## [1216] "Optimizing the root position on the branch 1214 ... 0.0193738205"              
## [1217] "Optimizing the root position on the branch 1215 ... 0.0195093090"              
## [1218] "Optimizing the root position on the branch 1216 ... 0.0195519744"              
## [1219] "Optimizing the root position on the branch 1217 ... 0.0195533619"              
## [1220] "Optimizing the root position on the branch 1218 ... 0.0183906109"              
## [1221] "Optimizing the root position on the branch 1219 ... 0.0183416563"              
## [1222] "Optimizing the root position on the branch 1220 ... 0.0183258010"              
## [1223] "Optimizing the root position on the branch 1221 ... 0.0182239240"              
## [1224] "Optimizing the root position on the branch 1222 ... 0.0148376504"              
## [1225] "Optimizing the root position on the branch 1223 ... 0.0147378439"              
## [1226] "Optimizing the root position on the branch 1224 ... 0.0178088185"              
## [1227] "Optimizing the root position on the branch 1225 ... 0.0178132797"              
## [1228] "Optimizing the root position on the branch 1226 ... 0.0158091007"              
## [1229] "Optimizing the root position on the branch 1227 ... 0.0157861741"              
## [1230] "Optimizing the root position on the branch 1228 ... 0.0158602253"              
## [1231] "Optimizing the root position on the branch 1229 ... 0.0158069632"              
## [1232] "Optimizing the root position on the branch 1230 ... 0.0096528222"              
## [1233] "Optimizing the root position on the branch 1231 ... 0.0136676103"              
## [1234] "Optimizing the root position on the branch 1232 ... 0.0139891491"              
## [1235] "Optimizing the root position on the branch 1233 ... 0.0137649674"              
## [1236] "Optimizing the root position on the branch 1234 ... 0.0138134535"              
## [1237] "Optimizing the root position on the branch 1235 ... 0.0124369064"              
## [1238] "Optimizing the root position on the branch 1236 ... 0.0124410667"              
## [1239] "Optimizing the root position on the branch 1237 ... 0.0124081322"              
## [1240] "Optimizing the root position on the branch 1238 ... 0.0127434974"              
## [1241] "Optimizing the root position on the branch 1239 ... 0.0129134567"              
## [1242] "Optimizing the root position on the branch 1240 ... 0.0130778830"              
## [1243] "Optimizing the root position on the branch 1241 ... 0.0131340776"              
## [1244] "Optimizing the root position on the branch 1242 ... 0.0131503603"              
## [1245] "Optimizing the root position on the branch 1243 ... 0.0131655885"              
## [1246] "Optimizing the root position on the branch 1244 ... 0.0131680724"              
## [1247] "Optimizing the root position on the branch 1245 ... 0.0129942419"              
## [1248] "Optimizing the root position on the branch 1246 ... 0.0130115444"              
## [1249] "Optimizing the root position on the branch 1247 ... 0.0130092326"              
## [1250] "Optimizing the root position on the branch 1248 ... 0.0130838378"              
## [1251] "Optimizing the root position on the branch 1249 ... 0.0130944811"              
## [1252] "Optimizing the root position on the branch 1250 ... 0.0136422091"              
## [1253] "Optimizing the root position on the branch 1251 ... 0.0136512042"              
## [1254] "Optimizing the root position on the branch 1252 ... 0.0143242979"              
## [1255] "Optimizing the root position on the branch 1253 ... 0.0146592947"              
## [1256] "Optimizing the root position on the branch 1254 ... 0.0147873898"              
## [1257] "Optimizing the root position on the branch 1255 ... 0.0147702326"              
## [1258] "Optimizing the root position on the branch 1256 ... 0.0110210147"              
## [1259] "Optimizing the root position on the branch 1257 ... 0.0135812986"              
## [1260] "Optimizing the root position on the branch 1258 ... 0.0136417399"              
## [1261] "Optimizing the root position on the branch 1259 ... 0.0136926911"              
## [1262] "Optimizing the root position on the branch 1260 ... 0.0137234664"              
## [1263] "Optimizing the root position on the branch 1261 ... 0.0137226357"              
## [1264] "Optimizing the root position on the branch 1262 ... 0.0133391531"              
## [1265] "Optimizing the root position on the branch 1263 ... 0.0135135549"              
## [1266] "Optimizing the root position on the branch 1264 ... 0.0137930593"              
## [1267] "Optimizing the root position on the branch 1265 ... 0.0137747306"              
## [1268] "Optimizing the root position on the branch 1266 ... 0.0135422384"              
## [1269] "Optimizing the root position on the branch 1267 ... 0.0097287075"              
## [1270] "Optimizing the root position on the branch 1268 ... 0.0097283724"              
## [1271] "Optimizing the root position on the branch 1269 ... 0.0151106024"              
## [1272] "Optimizing the root position on the branch 1270 ... 0.0162339577"              
## [1273] "Optimizing the root position on the branch 1271 ... 0.0162685632"              
## [1274] "Optimizing the root position on the branch 1272 ... 0.0175377461"              
## [1275] "Optimizing the root position on the branch 1273 ... 0.0175423690"              
## [1276] "Optimizing the root position on the branch 1274 ... 0.0174135415"              
## [1277] "Optimizing the root position on the branch 1275 ... 0.0174071597"              
## [1278] "Optimizing the root position on the branch 1276 ... 0.0171186999"              
## [1279] "Optimizing the root position on the branch 1277 ... 0.0096069325"              
## [1280] "Optimizing the root position on the branch 1278 ... 0.0139029902"              
## [1281] "Optimizing the root position on the branch 1279 ... 0.0138927217"              
## [1282] "Optimizing the root position on the branch 1280 ... 0.0137810379"              
## [1283] "Optimizing the root position on the branch 1281 ... 0.0138206130"              
## [1284] "Optimizing the root position on the branch 1282 ... 0.0129174144"              
## [1285] "Optimizing the root position on the branch 1283 ... 0.0129477826"              
## [1286] "Optimizing the root position on the branch 1284 ... 0.0126989835"              
## [1287] "Optimizing the root position on the branch 1285 ... 0.0129740576"              
## [1288] "Optimizing the root position on the branch 1286 ... 0.0142835866"              
## [1289] "Optimizing the root position on the branch 1287 ... 0.0143727587"              
## [1290] "Optimizing the root position on the branch 1288 ... 0.0096895499"              
## [1291] "Optimizing the root position on the branch 1289 ... 0.0102177971"              
## [1292] "Optimizing the root position on the branch 1290 ... 0.0110716729"              
## [1293] "Optimizing the root position on the branch 1291 ... 0.0116951774"              
## [1294] "Optimizing the root position on the branch 1292 ... 0.0117801646"              
## [1295] "Optimizing the root position on the branch 1293 ... 0.0115379064"              
## [1296] "Optimizing the root position on the branch 1294 ... 0.0119979446"              
## [1297] "Optimizing the root position on the branch 1295 ... 0.0119945961"              
## [1298] "Optimizing the root position on the branch 1296 ... 0.0119273275"              
## [1299] "Optimizing the root position on the branch 1297 ... 0.0119187309"              
## [1300] "Optimizing the root position on the branch 1298 ... 0.0118823847"              
## [1301] "Optimizing the root position on the branch 1299 ... 0.0119703913"              
## [1302] "Optimizing the root position on the branch 1300 ... 0.0119725051"              
## [1303] "Optimizing the root position on the branch 1301 ... 0.0119000816"              
## [1304] "Optimizing the root position on the branch 1302 ... 0.0119428179"              
## [1305] "Optimizing the root position on the branch 1303 ... 0.0119688582"              
## [1306] "Optimizing the root position on the branch 1304 ... 0.0120917283"              
## [1307] "Optimizing the root position on the branch 1305 ... 0.0120924628"              
## [1308] "Optimizing the root position on the branch 1306 ... 0.0121249052"              
## [1309] "Optimizing the root position on the branch 1307 ... 0.0121243598"              
## [1310] "Optimizing the root position on the branch 1308 ... 0.0120757943"              
## [1311] "Optimizing the root position on the branch 1309 ... 0.0125926761"              
## [1312] "Optimizing the root position on the branch 1310 ... 0.0126533771"              
## [1313] "Optimizing the root position on the branch 1311 ... 0.0127085321"              
## [1314] "Optimizing the root position on the branch 1312 ... 0.0127198461"              
## [1315] "Optimizing the root position on the branch 1313 ... 0.0127159140"              
## [1316] "Optimizing the root position on the branch 1314 ... 0.0104375815"              
## [1317] "Optimizing the root position on the branch 1315 ... 0.0121602397"              
## [1318] "Optimizing the root position on the branch 1316 ... 0.0121674396"              
## [1319] "Optimizing the root position on the branch 1317 ... 0.0119583808"              
## [1320] "Optimizing the root position on the branch 1318 ... 0.0119882339"              
## [1321] "Optimizing the root position on the branch 1319 ... 0.0119756756"              
## [1322] "Optimizing the root position on the branch 1320 ... 0.0117489199"              
## [1323] "Optimizing the root position on the branch 1321 ... 0.0118739714"              
## [1324] "Optimizing the root position on the branch 1322 ... 0.0128982412"              
## [1325] "Optimizing the root position on the branch 1323 ... 0.0133186358"              
## [1326] "Optimizing the root position on the branch 1324 ... 0.0140110758"              
## [1327] "Optimizing the root position on the branch 1325 ... 0.0140357323"              
## [1328] "Optimizing the root position on the branch 1326 ... 0.0133831593"              
## [1329] "Optimizing the root position on the branch 1327 ... 0.0133901880"              
## [1330] "Optimizing the root position on the branch 1328 ... 0.0101970875"              
## [1331] "Optimizing the root position on the branch 1329 ... 0.0111966513"              
## [1332] "Optimizing the root position on the branch 1330 ... 0.0131435142"              
## [1333] "Optimizing the root position on the branch 1331 ... 0.0135528593"              
## [1334] "Optimizing the root position on the branch 1332 ... 0.0136556014"              
## [1335] "Optimizing the root position on the branch 1333 ... 0.0136588710"              
## [1336] "Optimizing the root position on the branch 1334 ... 0.0156653397"              
## [1337] "Optimizing the root position on the branch 1335 ... 0.0157291726"              
## [1338] "Optimizing the root position on the branch 1336 ... 0.0127274849"              
## [1339] "Optimizing the root position on the branch 1337 ... 0.0131908788"              
## [1340] "Optimizing the root position on the branch 1338 ... 0.0136595425"              
## [1341] "Optimizing the root position on the branch 1339 ... 0.0136594968"              
## [1342] "Optimizing the root position on the branch 1340 ... 0.0135580114"              
## [1343] "Optimizing the root position on the branch 1341 ... 0.0134711708"              
## [1344] "Optimizing the root position on the branch 1342 ... 0.0111635718"              
## [1345] "Optimizing the root position on the branch 1343 ... 0.0111567639"              
## [1346] "Optimizing the root position on the branch 1344 ... 0.0114342443"              
## [1347] "Optimizing the root position on the branch 1345 ... 0.0114272880"              
## [1348] "Optimizing the root position on the branch 1346 ... 0.0112775679"              
## [1349] "Optimizing the root position on the branch 1347 ... 0.0112237989"              
## [1350] "Optimizing the root position on the branch 1348 ... 0.0112055621"              
## [1351] "Optimizing the root position on the branch 1349 ... 0.0111997510"              
## [1352] "Optimizing the root position on the branch 1350 ... 0.0120213505"              
## [1353] "Optimizing the root position on the branch 1351 ... 0.0120201294"              
## [1354] "Optimizing the root position on the branch 1352 ... 0.0115115561"              
## [1355] "Optimizing the root position on the branch 1353 ... 0.0123845755"              
## [1356] "Optimizing the root position on the branch 1354 ... 0.0123941289"              
## [1357] "Optimizing the root position on the branch 1355 ... 0.0119973085"              
## [1358] "Optimizing the root position on the branch 1356 ... 0.0113939671"              
## [1359] "Optimizing the root position on the branch 1357 ... 0.0113886970"              
## [1360] "Optimizing the root position on the branch 1358 ... 0.0112256364"              
## [1361] "Optimizing the root position on the branch 1359 ... 0.0111886740"              
## [1362] "Optimizing the root position on the branch 1360 ... 0.0131441589"              
## [1363] "Optimizing the root position on the branch 1361 ... 0.0131635781"              
## [1364] "Optimizing the root position on the branch 1362 ... 0.0105682433"              
## [1365] "Optimizing the root position on the branch 1363 ... 0.0100762286"              
## [1366] "Optimizing the root position on the branch 1364 ... 0.0145543751"              
## [1367] "Optimizing the root position on the branch 1365 ... 0.0145851425"              
## [1368] "Optimizing the root position on the branch 1366 ... 0.0125790313"              
## [1369] "Optimizing the root position on the branch 1367 ... 0.0145480311"              
## [1370] "Optimizing the root position on the branch 1368 ... 0.0145108828"              
## [1371] "Optimizing the root position on the branch 1369 ... 0.0125806062"              
## [1372] "Optimizing the root position on the branch 1370 ... 0.0125790319"              
## [1373] "Optimizing the root position on the branch 1371 ... 0.0122538233"              
## [1374] "Optimizing the root position on the branch 1372 ... 0.0135215305"              
## [1375] "Optimizing the root position on the branch 1373 ... 0.0135399824"              
## [1376] "Optimizing the root position on the branch 1374 ... 0.0107861464"              
## [1377] "Optimizing the root position on the branch 1375 ... 0.0117005237"              
## [1378] "Optimizing the root position on the branch 1376 ... 0.0122767062"              
## [1379] "Optimizing the root position on the branch 1377 ... 0.0122559888"              
## [1380] "Optimizing the root position on the branch 1378 ... 0.0119058109"              
## [1381] "Optimizing the root position on the branch 1379 ... 0.0119817565"              
## [1382] "Optimizing the root position on the branch 1380 ... 0.0099113857"              
## [1383] "Optimizing the root position on the branch 1381 ... 0.0100178963"              
## [1384] "Optimizing the root position on the branch 1382 ... 0.0110385218"              
## [1385] "Optimizing the root position on the branch 1383 ... 0.0110455077"              
## [1386] "Optimizing the root position on the branch 1384 ... 0.0110501974"              
## [1387] "Optimizing the root position on the branch 1385 ... 0.0112363264"              
## [1388] "Optimizing the root position on the branch 1386 ... 0.0118456422"              
## [1389] "Optimizing the root position on the branch 1387 ... 0.0125789441"              
## [1390] "Optimizing the root position on the branch 1388 ... 0.0125631506"              
## [1391] "Optimizing the root position on the branch 1389 ... 0.0113679000"              
## [1392] "Optimizing the root position on the branch 1390 ... 0.0113779939"              
## [1393] "Optimizing the root position on the branch 1391 ... 0.0105906538"              
## [1394] "Optimizing the root position on the branch 1392 ... 0.0096877590"              
## [1395] "Optimizing the root position on the branch 1393 ... 0.0095467971"              
## [1396] "Optimizing the root position on the branch 1394 ... 0.0172008398"              
## [1397] "Optimizing the root position on the branch 1395 ... 0.0173887976"              
## [1398] "Optimizing the root position on the branch 1396 ... 0.0174734206"              
## [1399] "Optimizing the root position on the branch 1397 ... 0.0174671146"              
## [1400] "Optimizing the root position on the branch 1398 ... 0.0174082624"              
## [1401] "Optimizing the root position on the branch 1399 ... 0.0178434301"              
## [1402] "Optimizing the root position on the branch 1400 ... 0.0178420343"              
## [1403] "Optimizing the root position on the branch 1401 ... 0.0173622625"              
## [1404] "Optimizing the root position on the branch 1402 ... 0.0094210631"              
## [1405] "Optimizing the root position on the branch 1403 ... 0.0117976108"              
## [1406] "Optimizing the root position on the branch 1404 ... 0.0157874324"              
## [1407] "Optimizing the root position on the branch 1405 ... 0.0157393331"              
## [1408] "Optimizing the root position on the branch 1406 ... 0.0163078656"              
## [1409] "Optimizing the root position on the branch 1407 ... 0.0162724486"              
## [1410] "Optimizing the root position on the branch 1408 ... 0.0095450029"              
## [1411] "Optimizing the root position on the branch 1409 ... 0.0216123929"              
## [1412] "Optimizing the root position on the branch 1410 ... 0.0217835458"              
## [1413] "Optimizing the root position on the branch 1411 ... 0.0144586841"              
## [1414] "Optimizing the root position on the branch 1412 ... 0.0146732360"              
## [1415] "Optimizing the root position on the branch 1413 ... 0.0100222663"              
## [1416] "Optimizing the root position on the branch 1414 ... 0.0107170775"              
## [1417] "Optimizing the root position on the branch 1415 ... 0.0107039368"              
## [1418] "Optimizing the root position on the branch 1416 ... 0.0154043789"              
## [1419] "Optimizing the root position on the branch 1417 ... 0.0154364453"              
## [1420] "Optimizing the root position on the branch 1418 ... 0.0146126443"              
## [1421] "Optimizing the root position on the branch 1419 ... 0.0145506270"              
## [1422] "Optimizing the root position on the branch 1420 ... 0.0097060876"              
## [1423] "Optimizing the root position on the branch 1421 ... 0.0092622103"              
## [1424] "Optimizing the root position on the branch 1422 ... 0.0118765634"              
## [1425] "Optimizing the root position on the branch 1423 ... 0.0124770762"              
## [1426] "Optimizing the root position on the branch 1424 ... 0.0124253776"              
## [1427] "Optimizing the root position on the branch 1425 ... 0.0123711705"              
## [1428] "Optimizing the root position on the branch 1426 ... 0.0120675655"              
## [1429] "Optimizing the root position on the branch 1427 ... 0.0116168354"              
## [1430] "Optimizing the root position on the branch 1428 ... 0.0116780559"              
## [1431] "Optimizing the root position on the branch 1429 ... 0.0102756871"              
## [1432] "Optimizing the root position on the branch 1430 ... 0.0111262298"              
## [1433] "Optimizing the root position on the branch 1431 ... 0.0113001332"              
## [1434] "Optimizing the root position on the branch 1432 ... 0.0113469297"              
## [1435] "Optimizing the root position on the branch 1433 ... 0.0119112625"              
## [1436] "Optimizing the root position on the branch 1434 ... 0.0119300622"              
## [1437] "Optimizing the root position on the branch 1435 ... 0.0118102259"              
## [1438] "Optimizing the root position on the branch 1436 ... 0.0124841161"              
## [1439] "Optimizing the root position on the branch 1437 ... 0.0124995946"              
## [1440] "Optimizing the root position on the branch 1438 ... 0.0119713792"              
## [1441] "Optimizing the root position on the branch 1439 ... 0.0121350750"              
## [1442] "Optimizing the root position on the branch 1440 ... 0.0121682524"              
## [1443] "Optimizing the root position on the branch 1441 ... 0.0121673129"              
## [1444] "Optimizing the root position on the branch 1442 ... 0.0121650505"              
## [1445] "Optimizing the root position on the branch 1443 ... 0.0121587983"              
## [1446] "Optimizing the root position on the branch 1444 ... 0.0133038304"              
## [1447] "Optimizing the root position on the branch 1445 ... 0.0139760805"              
## [1448] "Optimizing the root position on the branch 1446 ... 0.0142890185"              
## [1449] "Optimizing the root position on the branch 1447 ... 0.0142816764"              
## [1450] "Optimizing the root position on the branch 1448 ... 0.0099749850"              
## [1451] "Optimizing the root position on the branch 1449 ... 0.0154577118"              
## [1452] "Optimizing the root position on the branch 1450 ... 0.0149754566"              
## [1453] "Optimizing the root position on the branch 1451 ... 0.0123906659"              
## [1454] "Optimizing the root position on the branch 1452 ... 0.0122001342"              
## [1455] "Optimizing the root position on the branch 1453 ... 0.0132345664"              
## [1456] "Optimizing the root position on the branch 1454 ... 0.0132536682"              
## [1457] "Optimizing the root position on the branch 1455 ... 0.0123776822"              
## [1458] "Optimizing the root position on the branch 1456 ... 0.0128483282"              
## [1459] "Optimizing the root position on the branch 1457 ... 0.0128524740"              
## [1460] "Optimizing the root position on the branch 1458 ... 0.0142650088"              
## [1461] "Optimizing the root position on the branch 1459 ... 0.0159224331"              
## [1462] "Optimizing the root position on the branch 1460 ... 0.0159165178"              
## [1463] "Optimizing the root position on the branch 1461 ... 0.0127557929"              
## [1464] "Optimizing the root position on the branch 1462 ... 0.0138344619"              
## [1465] "Optimizing the root position on the branch 1463 ... 0.0138331836"              
## [1466] "Optimizing the root position on the branch 1464 ... 0.0120699835"              
## [1467] "Optimizing the root position on the branch 1465 ... 0.0121056207"              
## [1468] "Optimizing the root position on the branch 1466 ... 0.0119831214"              
## [1469] "Optimizing the root position on the branch 1467 ... 0.0120368921"              
## [1470] "Optimizing the root position on the branch 1468 ... 0.0123936496"              
## [1471] "Optimizing the root position on the branch 1469 ... 0.0123884091"              
## [1472] "Optimizing the root position on the branch 1470 ... 0.0117754055"              
## [1473] "Optimizing the root position on the branch 1471 ... 0.0091420323"              
## [1474] "Optimizing the root position on the branch 1472 ... 0.0104887641"              
## [1475] "Optimizing the root position on the branch 1473 ... 0.0122852259"              
## [1476] "Optimizing the root position on the branch 1474 ... 0.0125703493"              
## [1477] "Optimizing the root position on the branch 1475 ... 0.0138068624"              
## [1478] "Optimizing the root position on the branch 1476 ... 0.0137921760"              
## [1479] "Optimizing the root position on the branch 1477 ... 0.0132335318"              
## [1480] "Optimizing the root position on the branch 1478 ... 0.0133201222"              
## [1481] "Optimizing the root position on the branch 1479 ... 0.0133267787"              
## [1482] "Optimizing the root position on the branch 1480 ... 0.0133101683"              
## [1483] "Optimizing the root position on the branch 1481 ... 0.0133138458"              
## [1484] "Optimizing the root position on the branch 1482 ... 0.0126111874"              
## [1485] "Optimizing the root position on the branch 1483 ... 0.0102906105"              
## [1486] "Optimizing the root position on the branch 1484 ... 0.0144044711"              
## [1487] "Optimizing the root position on the branch 1485 ... 0.0149933536"              
## [1488] "Optimizing the root position on the branch 1486 ... 0.0149795038"              
## [1489] "Optimizing the root position on the branch 1487 ... 0.0108128865"              
## [1490] "Optimizing the root position on the branch 1488 ... 0.0108406784"              
## [1491] "Optimizing the root position on the branch 1489 ... 0.0095983634"              
## [1492] "Optimizing the root position on the branch 1490 ... 0.0096783017"              
## [1493] "Optimizing the root position on the branch 1491 ... 0.0100882181"              
## [1494] "Optimizing the root position on the branch 1492 ... 0.0103877784"              
## [1495] "Optimizing the root position on the branch 1493 ... 0.0109761376"              
## [1496] "Optimizing the root position on the branch 1494 ... 0.0109681348"              
## [1497] "Optimizing the root position on the branch 1495 ... 0.0107752817"              
## [1498] "Optimizing the root position on the branch 1496 ... 0.0108376550"              
## [1499] "Optimizing the root position on the branch 1497 ... 0.0102352144"              
## [1500] "Optimizing the root position on the branch 1498 ... 0.0131325774"              
## [1501] "Optimizing the root position on the branch 1499 ... 0.0131104837"              
## [1502] "Optimizing the root position on the branch 1500 ... 0.0184745659"              
## [1503] "Optimizing the root position on the branch 1501 ... 0.0186663169"              
## [1504] "Optimizing the root position on the branch 1502 ... 0.0097705569"              
## [1505] "Optimizing the root position on the branch 1503 ... 0.0115354433"              
## [1506] "Optimizing the root position on the branch 1504 ... 0.0115900959"              
## [1507] "Optimizing the root position on the branch 1505 ... 0.0091970105"              
## [1508] "Optimizing the root position on the branch 1506 ... 0.0092952721"              
## [1509] "Optimizing the root position on the branch 1507 ... 0.0147368287"              
## [1510] "Optimizing the root position on the branch 1508 ... 0.0147444058"              
## [1511] "Optimizing the root position on the branch 1509 ... 0.0140311541"              
## [1512] "Optimizing the root position on the branch 1510 ... 0.0140360733"              
## [1513] "Optimizing the root position on the branch 1511 ... 0.0125212663"              
## [1514] "Optimizing the root position on the branch 1512 ... 0.0138874791"              
## [1515] "Optimizing the root position on the branch 1513 ... 0.0139276957"              
## [1516] "Optimizing the root position on the branch 1514 ... 0.0142559364"              
## [1517] "Optimizing the root position on the branch 1515 ... 0.0142554720"              
## [1518] "Optimizing the root position on the branch 1516 ... 0.0160634221"              
## [1519] "Optimizing the root position on the branch 1517 ... 0.0168651661"              
## [1520] "Optimizing the root position on the branch 1518 ... 0.0172724006"              
## [1521] "Optimizing the root position on the branch 1519 ... 0.0174086941"              
## [1522] "Optimizing the root position on the branch 1520 ... 0.0173906294"              
## [1523] "Optimizing the root position on the branch 1521 ... 0.0120821846"              
## [1524] "Optimizing the root position on the branch 1522 ... 0.0120160810"              
## [1525] "Optimizing the root position on the branch 1523 ... 0.0101310338"              
## [1526] "Optimizing the root position on the branch 1524 ... 0.0101805296"              
## [1527] "Optimizing the root position on the branch 1525 ... 0.0153393928"              
## [1528] "Optimizing the root position on the branch 1526 ... 0.0154769330"              
## [1529] "Optimizing the root position on the branch 1527 ... 0.0153072952"              
## [1530] "Optimizing the root position on the branch 1528 ... 0.0153679177"              
## [1531] "Optimizing the root position on the branch 1529 ... 0.0138470629"              
## [1532] "Optimizing the root position on the branch 1530 ... 0.0147435496"              
## [1533] "Optimizing the root position on the branch 1531 ... 0.0147684054"              
## [1534] "Optimizing the root position on the branch 1532 ... 0.0091713091"              
## [1535] "Optimizing the root position on the branch 1533 ... 0.0091958329"              
## [1536] "Optimizing the root position on the branch 1534 ... 0.0148624419"              
## [1537] "Optimizing the root position on the branch 1535 ... 0.0156212294"              
## [1538] "Optimizing the root position on the branch 1536 ... 0.0158838859"              
## [1539] "Optimizing the root position on the branch 1537 ... 0.0158846726"              
## [1540] "Optimizing the root position on the branch 1538 ... 0.0141116594"              
## [1541] "Optimizing the root position on the branch 1539 ... 0.0140179448"              
## [1542] "Optimizing the root position on the branch 1540 ... 0.0106128475"              
## [1543] "Optimizing the root position on the branch 1541 ... 0.0120583795"              
## [1544] "Optimizing the root position on the branch 1542 ... 0.0121458793"              
## [1545] "Optimizing the root position on the branch 1543 ... 0.0121325771"              
## [1546] "Optimizing the root position on the branch 1544 ... 0.0121641885"              
## [1547] "Optimizing the root position on the branch 1545 ... 0.0121436621"              
## [1548] "Optimizing the root position on the branch 1546 ... 0.0120560482"              
## [1549] "Optimizing the root position on the branch 1547 ... 0.0120587361"              
## [1550] "Optimizing the root position on the branch 1548 ... 0.0122093819"              
## [1551] "Optimizing the root position on the branch 1549 ... 0.0102473661"              
## [1552] "Optimizing the root position on the branch 1550 ... 0.0096544544"              
## [1553] "Optimizing the root position on the branch 1551 ... 0.0127964924"              
## [1554] "Optimizing the root position on the branch 1552 ... 0.0131560711"              
## [1555] "Optimizing the root position on the branch 1553 ... 0.0131819378"              
## [1556] "Optimizing the root position on the branch 1554 ... 0.0101405673"              
## [1557] "Optimizing the root position on the branch 1555 ... 0.0130749516"              
## [1558] "Optimizing the root position on the branch 1556 ... 0.0130908374"              
## [1559] "Optimizing the root position on the branch 1557 ... 0.0116331372"              
## [1560] "Optimizing the root position on the branch 1558 ... 0.0117476903"              
## [1561] "Optimizing the root position on the branch 1559 ... 0.0124710802"              
## [1562] "Optimizing the root position on the branch 1560 ... 0.0125977632"              
## [1563] "Optimizing the root position on the branch 1561 ... 0.0125935974"              
## [1564] "Optimizing the root position on the branch 1562 ... 0.0123541318"              
## [1565] "Optimizing the root position on the branch 1563 ... 0.0123647002"              
## [1566] "Optimizing the root position on the branch 1564 ... 0.0122988641"              
## [1567] "Optimizing the root position on the branch 1565 ... 0.0124919651"              
## [1568] "Optimizing the root position on the branch 1566 ... 0.0125418868"              
## [1569] "Optimizing the root position on the branch 1567 ... 0.0115137155"              
## [1570] "Optimizing the root position on the branch 1568 ... 0.0115764498"              
## [1571] "Optimizing the root position on the branch 1569 ... 0.0117602864"              
## [1572] "Optimizing the root position on the branch 1570 ... 0.0119757360"              
## [1573] "Optimizing the root position on the branch 1571 ... 0.0119967563"              
## [1574] "Optimizing the root position on the branch 1572 ... 0.0119548174"              
## [1575] "Optimizing the root position on the branch 1573 ... 0.0119423418"              
## [1576] "Optimizing the root position on the branch 1574 ... 0.0132088410"              
## [1577] "Optimizing the root position on the branch 1575 ... 0.0131556405"              
## [1578] "Optimizing the root position on the branch 1576 ... 0.0135075664"              
## [1579] "Optimizing the root position on the branch 1577 ... 0.0135857724"              
## [1580] "Optimizing the root position on the branch 1578 ... 0.0132018920"              
## [1581] "Optimizing the root position on the branch 1579 ... 0.0135157434"              
## [1582] "Optimizing the root position on the branch 1580 ... 0.0135039730"              
## [1583] "Optimizing the root position on the branch 1581 ... 0.0127430665"              
## [1584] "Optimizing the root position on the branch 1582 ... 0.0120746288"              
## [1585] "Optimizing the root position on the branch 1583 ... 0.0126250917"              
## [1586] "Optimizing the root position on the branch 1584 ... 0.0130240686"              
## [1587] "Optimizing the root position on the branch 1585 ... 0.0134120864"              
## [1588] "Optimizing the root position on the branch 1586 ... 0.0134001324"              
## [1589] "Optimizing the root position on the branch 1587 ... 0.0124902074"              
## [1590] "Optimizing the root position on the branch 1588 ... 0.0094585667"              
## [1591] "Optimizing the root position on the branch 1589 ... 0.0147329273"              
## [1592] "Optimizing the root position on the branch 1590 ... 0.0146604705"              
## [1593] "Optimizing the root position on the branch 1591 ... 0.0124198992"              
## [1594] "Optimizing the root position on the branch 1592 ... 0.0133098301"              
## [1595] "Optimizing the root position on the branch 1593 ... 0.0133405231"              
## [1596] "Optimizing the root position on the branch 1594 ... 0.0092506162"              
## [1597] "Optimizing the root position on the branch 1595 ... 0.0093283507"              
## [1598] "Optimizing the root position on the branch 1596 ... 0.0173792224"              
## [1599] "Optimizing the root position on the branch 1597 ... 0.0173413484"              
## [1600] "Optimizing the root position on the branch 1598 ... 0.0122548328"              
## [1601] "Optimizing the root position on the branch 1599 ... 0.0153624101"              
## [1602] "Optimizing the root position on the branch 1600 ... 0.0156919895"              
## [1603] "Optimizing the root position on the branch 1601 ... 0.0156943540"              
## [1604] "Optimizing the root position on the branch 1602 ... 0.0093434856"              
## [1605] "Optimizing the root position on the branch 1603 ... 0.0123433880"              
## [1606] "Optimizing the root position on the branch 1604 ... 0.0124138490"              
## [1607] "Optimizing the root position on the branch 1605 ... 0.0121875863"              
## [1608] "Optimizing the root position on the branch 1606 ... 0.0122468381"              
## [1609] "Optimizing the root position on the branch 1607 ... 0.0096929842"              
## [1610] "Optimizing the root position on the branch 1608 ... 0.0114318551"              
## [1611] "Optimizing the root position on the branch 1609 ... 0.0132924181"              
## [1612] "Optimizing the root position on the branch 1610 ... 0.0133537163"              
## [1613] "Optimizing the root position on the branch 1611 ... 0.0142626041"              
## [1614] "Optimizing the root position on the branch 1612 ... 0.0148857074"              
## [1615] "Optimizing the root position on the branch 1613 ... 0.0147880396"              
## [1616] "Optimizing the root position on the branch 1614 ... 0.0151834168"              
## [1617] "Optimizing the root position on the branch 1615 ... 0.0149685843"              
## [1618] "Optimizing the root position on the branch 1616 ... 0.0094327809"              
## [1619] "Optimizing the root position on the branch 1617 ... 0.0118514362"              
## [1620] "Optimizing the root position on the branch 1618 ... 0.0123003549"              
## [1621] "Optimizing the root position on the branch 1619 ... 0.0122773671"              
## [1622] "Optimizing the root position on the branch 1620 ... 0.0120364137"              
## [1623] "Optimizing the root position on the branch 1621 ... 0.0116194247"              
## [1624] "Optimizing the root position on the branch 1622 ... 0.0098698472"              
## [1625] "Optimizing the root position on the branch 1623 ... 0.0099573483"              
## [1626] "Optimizing the root position on the branch 1624 ... 0.0130406532"              
## [1627] "Optimizing the root position on the branch 1625 ... 0.0130047906"              
## [1628] "Optimizing the root position on the branch 1626 ... 0.0114721773"              
## [1629] "Optimizing the root position on the branch 1627 ... 0.0116152050"              
## [1630] "Optimizing the root position on the branch 1628 ... 0.0116066674"              
## [1631] "Optimizing the root position on the branch 1629 ... 0.0116128206"              
## [1632] "Optimizing the root position on the branch 1630 ... 0.0116198924"              
## [1633] "Optimizing the root position on the branch 1631 ... 0.0117672948"              
## [1634] "Optimizing the root position on the branch 1632 ... 0.0117672100"              
## [1635] "Optimizing the root position on the branch 1633 ... 0.0116060645"              
## [1636] "Optimizing the root position on the branch 1634 ... 0.0107371842"              
## [1637] "Optimizing the root position on the branch 1635 ... 0.0137986962"              
## [1638] "Optimizing the root position on the branch 1636 ... 0.0138602459"              
## [1639] "Optimizing the root position on the branch 1637 ... 0.0133200434"              
## [1640] "Optimizing the root position on the branch 1638 ... 0.0134461579"              
## [1641] "Optimizing the root position on the branch 1639 ... 0.0136373044"              
## [1642] "Optimizing the root position on the branch 1640 ... 0.0136274927"              
## [1643] "Optimizing the root position on the branch 1641 ... 0.0135740366"              
## [1644] "Optimizing the root position on the branch 1642 ... 0.0135768786"              
## [1645] "Optimizing the root position on the branch 1643 ... 0.0143211424"              
## [1646] "Optimizing the root position on the branch 1644 ... 0.0142793388"              
## [1647] "Optimizing the root position on the branch 1645 ... 0.0129929247"              
## [1648] "Optimizing the root position on the branch 1646 ... 0.0137741133"              
## [1649] "Optimizing the root position on the branch 1647 ... 0.0137552155"              
## [1650] "Optimizing the root position on the branch 1648 ... 0.0118248754"              
## [1651] "Optimizing the root position on the branch 1649 ... 0.0121695275"              
## [1652] "Optimizing the root position on the branch 1650 ... 0.0123467027"              
## [1653] "Optimizing the root position on the branch 1651 ... 0.0123280633"              
## [1654] "Optimizing the root position on the branch 1652 ... 0.0128554272"              
## [1655] "Optimizing the root position on the branch 1653 ... 0.0128434850"              
## [1656] "Optimizing the root position on the branch 1654 ... 0.0126642662"              
## [1657] "Optimizing the root position on the branch 1655 ... 0.0126716180"              
## [1658] "Optimizing the root position on the branch 1656 ... 0.0126931877"              
## [1659] "Optimizing the root position on the branch 1657 ... 0.0126971756"              
## [1660] "Optimizing the root position on the branch 1658 ... 0.0128332435"              
## [1661] "Optimizing the root position on the branch 1659 ... 0.0128297942"              
## [1662] "Optimizing the root position on the branch 1660 ... 0.0094664769"              
## [1663] "Optimizing the root position on the branch 1661 ... 0.0132529404"              
## [1664] "Optimizing the root position on the branch 1662 ... 0.0132450441"              
## [1665] "Optimizing the root position on the branch 1663 ... 0.0122976953"              
## [1666] "Optimizing the root position on the branch 1664 ... 0.0123321218"              
## [1667] "Optimizing the root position on the branch 1665 ... 0.0131943625"              
## [1668] "Optimizing the root position on the branch 1666 ... 0.0135873950"              
## [1669] "Optimizing the root position on the branch 1667 ... 0.0136134096"              
## [1670] "Optimizing the root position on the branch 1668 ... 0.0137440178"              
## [1671] "Optimizing the root position on the branch 1669 ... 0.0137552300"              
## [1672] "Optimizing the root position on the branch 1670 ... 0.0103870044"              
## [1673] "Optimizing the root position on the branch 1671 ... 0.0152085534"              
## [1674] "Optimizing the root position on the branch 1672 ... 0.0152017458"              
## [1675] "Optimizing the root position on the branch 1673 ... 0.0107791155"              
## [1676] "Optimizing the root position on the branch 1674 ... 0.0122361047"              
## [1677] "Optimizing the root position on the branch 1675 ... 0.0124330612"              
## [1678] "Optimizing the root position on the branch 1676 ... 0.0124287541"              
## [1679] "Optimizing the root position on the branch 1677 ... 0.0136312018"              
## [1680] "Optimizing the root position on the branch 1678 ... 0.0136606506"              
## [1681] "Optimizing the root position on the branch 1679 ... 0.0123588846"              
## [1682] "Optimizing the root position on the branch 1680 ... 0.0124366378"              
## [1683] "Optimizing the root position on the branch 1681 ... 0.0124558402"              
## [1684] "Optimizing the root position on the branch 1682 ... 0.0144799364"              
## [1685] "Optimizing the root position on the branch 1683 ... 0.0144918469"              
## [1686] "Optimizing the root position on the branch 1684 ... 0.0142062198"              
## [1687] "Optimizing the root position on the branch 1685 ... 0.0137060247"              
## [1688] "Optimizing the root position on the branch 1686 ... 0.0111814784"              
## [1689] "Optimizing the root position on the branch 1687 ... 0.0110384729"              
## [1690] "Optimizing the root position on the branch 1688 ... 0.0137109897"              
## [1691] "Optimizing the root position on the branch 1689 ... 0.0141339090"              
## [1692] "Optimizing the root position on the branch 1690 ... 0.0141270866"              
## [1693] "Optimizing the root position on the branch 1691 ... 0.0101593617"              
## [1694] "Optimizing the root position on the branch 1692 ... 0.0132208393"              
## [1695] "Optimizing the root position on the branch 1693 ... 0.0131311388"              
## [1696] "Optimizing the root position on the branch 1694 ... 0.0118907822"              
## [1697] "Optimizing the root position on the branch 1695 ... 0.0119702055"              
## [1698] "Optimizing the root position on the branch 1696 ... 0.0100203832"              
## [1699] "Optimizing the root position on the branch 1697 ... 0.0113456131"              
## [1700] "Optimizing the root position on the branch 1698 ... 0.0114626697"              
## [1701] "Optimizing the root position on the branch 1699 ... 0.0099162048"              
## [1702] "Optimizing the root position on the branch 1700 ... 0.0147626377"              
## [1703] "Optimizing the root position on the branch 1701 ... 0.0148195804"              
## [1704] "Optimizing the root position on the branch 1702 ... 0.0121069219"              
## [1705] "Optimizing the root position on the branch 1703 ... 0.0130842701"              
## [1706] "Optimizing the root position on the branch 1704 ... 0.0141631683"              
## [1707] "Optimizing the root position on the branch 1705 ... 0.0141505681"              
## [1708] "Optimizing the root position on the branch 1706 ... 0.0097278335"              
## [1709] "Optimizing the root position on the branch 1707 ... 0.0114796429"              
## [1710] "Optimizing the root position on the branch 1708 ... 0.0117611464"              
## [1711] "Optimizing the root position on the branch 1709 ... 0.0119257248"              
## [1712] "Optimizing the root position on the branch 1710 ... 0.0119266865"              
## [1713] "Optimizing the root position on the branch 1711 ... 0.0118025727"              
## [1714] "Optimizing the root position on the branch 1712 ... 0.0123385039"              
## [1715] "Optimizing the root position on the branch 1713 ... 0.0123548920"              
## [1716] "Optimizing the root position on the branch 1714 ... 0.0111779366"              
## [1717] "Optimizing the root position on the branch 1715 ... 0.0100515676"              
## [1718] "Optimizing the root position on the branch 1716 ... 0.0120806934"              
## [1719] "Optimizing the root position on the branch 1717 ... 0.0143651849"              
## [1720] "Optimizing the root position on the branch 1718 ... 0.0143662730"              
## [1721] "Optimizing the root position on the branch 1719 ... 0.0096960354"              
## [1722] "Optimizing the root position on the branch 1720 ... 0.0161591047"              
## [1723] "Optimizing the root position on the branch 1721 ... 0.0162493748"              
## [1724] "Optimizing the root position on the branch 1722 ... 0.0090052411"              
## [1725] "Optimizing the root position on the branch 1723 ... 0.0152500062"              
## [1726] "Optimizing the root position on the branch 1724 ... 0.0150470600"              
## [1727] "Optimizing the root position on the branch 1725 ... 0.0208071377"              
## [1728] "Optimizing the root position on the branch 1726 ... 0.0208553373"              
## [1729] "Optimizing the root position on the branch 1727 ... 0.0136783094"              
## [1730] "Optimizing the root position on the branch 1728 ... 0.0138173236"              
## [1731] "Optimizing the root position on the branch 1729 ... 0.0101009135"              
## [1732] "Optimizing the root position on the branch 1730 ... 0.0114249582"              
## [1733] "Optimizing the root position on the branch 1731 ... 0.0115595701"              
## [1734] "Optimizing the root position on the branch 1732 ... 0.0115823862"              
## [1735] "Optimizing the root position on the branch 1733 ... 0.0126614373"              
## [1736] "Optimizing the root position on the branch 1734 ... 0.0126877120"              
## [1737] "Optimizing the root position on the branch 1735 ... 0.0117071904"              
## [1738] "Optimizing the root position on the branch 1736 ... 0.0116929612"              
## [1739] "Optimizing the root position on the branch 1737 ... 0.0112777729"              
## [1740] "Optimizing the root position on the branch 1738 ... 0.0120428900"              
## [1741] "Optimizing the root position on the branch 1739 ... 0.0121998318"              
## [1742] "Optimizing the root position on the branch 1740 ... 0.0121996593"              
## [1743] "Optimizing the root position on the branch 1741 ... 0.0121892793"              
## [1744] "Optimizing the root position on the branch 1742 ... 0.0122374654"              
## [1745] "Optimizing the root position on the branch 1743 ... 0.0122385418"              
## [1746] "Optimizing the root position on the branch 1744 ... 0.0121120720"              
## [1747] "Optimizing the root position on the branch 1745 ... 0.0125129770"              
## [1748] "Optimizing the root position on the branch 1746 ... 0.0124902392"              
## [1749] "Optimizing the root position on the branch 1747 ... 0.0113776929"              
## [1750] "Optimizing the root position on the branch 1748 ... 0.0113951170"              
## [1751] "Optimizing the root position on the branch 1749 ... 0.0099269329"              
## [1752] "Optimizing the root position on the branch 1750 ... 0.0134674008"              
## [1753] "Optimizing the root position on the branch 1751 ... 0.0135494767"              
## [1754] "Optimizing the root position on the branch 1752 ... 0.0126410289"              
## [1755] "Optimizing the root position on the branch 1753 ... 0.0102225799"              
## [1756] "Optimizing the root position on the branch 1754 ... 0.0129799689"              
## [1757] "Optimizing the root position on the branch 1755 ... 0.0130039282"              
## [1758] "Optimizing the root position on the branch 1756 ... 0.0122333012"              
## [1759] "Optimizing the root position on the branch 1757 ... 0.0115332345"              
## [1760] "Optimizing the root position on the branch 1758 ... 0.0115474359"              
## [1761] "Optimizing the root position on the branch 1759 ... 0.0114760287"              
## [1762] "Optimizing the root position on the branch 1760 ... 0.0128469267"              
## [1763] "Optimizing the root position on the branch 1761 ... 0.0128779102"              
## [1764] "Optimizing the root position on the branch 1762 ... 0.0102780296"              
## [1765] "Optimizing the root position on the branch 1763 ... 0.0109165099"              
## [1766] "Optimizing the root position on the branch 1764 ... 0.0118680107"              
## [1767] "Optimizing the root position on the branch 1765 ... 0.0120648117"              
## [1768] "Optimizing the root position on the branch 1766 ... 0.0120565476"              
## [1769] "Optimizing the root position on the branch 1767 ... 0.0122095198"              
## [1770] "Optimizing the root position on the branch 1768 ... 0.0122014637"              
## [1771] "Optimizing the root position on the branch 1769 ... 0.0122429104"              
## [1772] "Optimizing the root position on the branch 1770 ... 0.0122380742"              
## [1773] "Optimizing the root position on the branch 1771 ... 0.0120216681"              
## [1774] "Optimizing the root position on the branch 1772 ... 0.0118509940"              
## [1775] "Optimizing the root position on the branch 1773 ... 0.0118536038"              
## [1776] "Optimizing the root position on the branch 1774 ... 0.0133535866"              
## [1777] "Optimizing the root position on the branch 1775 ... 0.0136782976"              
## [1778] "Optimizing the root position on the branch 1776 ... 0.0136493126"              
## [1779] "Optimizing the root position on the branch 1777 ... 0.0141825965"              
## [1780] "Optimizing the root position on the branch 1778 ... 0.0141646447"              
## [1781] "Optimizing the root position on the branch 1779 ... 0.0133239039"              
## [1782] "Optimizing the root position on the branch 1780 ... 0.0133425214"              
## [1783] "Optimizing the root position on the branch 1781 ... 0.0129446517"              
## [1784] "Optimizing the root position on the branch 1782 ... 0.0129423531"              
## [1785] "Optimizing the root position on the branch 1783 ... 0.0119912070"              
## [1786] "Optimizing the root position on the branch 1784 ... 0.0140920015"              
## [1787] "Optimizing the root position on the branch 1785 ... 0.0140362103"              
## [1788] "Optimizing the root position on the branch 1786 ... 0.0149412603"              
## [1789] "Optimizing the root position on the branch 1787 ... 0.0156529613"              
## [1790] "Optimizing the root position on the branch 1788 ... 0.0156551401"              
## [1791] "Optimizing the root position on the branch 1789 ... 0.0155213126"              
## [1792] "Optimizing the root position on the branch 1790 ... 0.0177123994"              
## [1793] "Optimizing the root position on the branch 1791 ... 0.0178367655"              
## [1794] "Optimizing the root position on the branch 1792 ... 0.0095636501"              
## [1795] "Optimizing the root position on the branch 1793 ... 0.0171160223"              
## [1796] "Optimizing the root position on the branch 1794 ... 0.0170258770"              
## [1797] "Optimizing the root position on the branch 1795 ... 0.0153806270"              
## [1798] "Optimizing the root position on the branch 1796 ... 0.0153628853"              
## [1799] "Optimizing the root position on the branch 1797 ... 0.0139911693"              
## [1800] "Optimizing the root position on the branch 1798 ... 0.0140581926"              
## [1801] "Optimizing the root position on the branch 1799 ... 0.0103061608"              
## [1802] "Optimizing the root position on the branch 1800 ... 0.0101060656"              
## [1803] "Optimizing the root position on the branch 1801 ... 0.0136804505"              
## [1804] "Optimizing the root position on the branch 1802 ... 0.0135903028"              
## [1805] "Optimizing the root position on the branch 1803 ... 0.0106093112"              
## [1806] "Optimizing the root position on the branch 1804 ... 0.0107837253"              
## [1807] "Optimizing the root position on the branch 1805 ... 0.0108191653"              
## [1808] "Optimizing the root position on the branch 1806 ... 0.0106901956"              
## [1809] "Optimizing the root position on the branch 1807 ... 0.0105309310"              
## [1810] "Optimizing the root position on the branch 1808 ... 0.0107612781"              
## [1811] "Optimizing the root position on the branch 1809 ... 0.0109865472"              
## [1812] "Optimizing the root position on the branch 1810 ... 0.0112932701"              
## [1813] "Optimizing the root position on the branch 1811 ... 0.0115353471"              
## [1814] "Optimizing the root position on the branch 1812 ... 0.0115866622"              
## [1815] "Optimizing the root position on the branch 1813 ... 0.0116053401"              
## [1816] "Optimizing the root position on the branch 1814 ... 0.0116085177"              
## [1817] "Optimizing the root position on the branch 1815 ... 0.0116941185"              
## [1818] "Optimizing the root position on the branch 1816 ... 0.0117032825"              
## [1819] "Optimizing the root position on the branch 1817 ... 0.0124528386"              
## [1820] "Optimizing the root position on the branch 1818 ... 0.0125471353"              
## [1821] "Optimizing the root position on the branch 1819 ... 0.0115452718"              
## [1822] "Optimizing the root position on the branch 1820 ... 0.0116740241"              
## [1823] "Optimizing the root position on the branch 1821 ... 0.0119563829"              
## [1824] "Optimizing the root position on the branch 1822 ... 0.0125201818"              
## [1825] "Optimizing the root position on the branch 1823 ... 0.0125214138"              
## [1826] "Optimizing the root position on the branch 1824 ... 0.0119684398"              
## [1827] "Optimizing the root position on the branch 1825 ... 0.0118908783"              
## [1828] "Optimizing the root position on the branch 1826 ... 0.0120151489"              
## [1829] "Optimizing the root position on the branch 1827 ... 0.0120285061"              
## [1830] "Optimizing the root position on the branch 1828 ... 0.0105278552"              
## [1831] "Optimizing the root position on the branch 1829 ... 0.0129587811"              
## [1832] "Optimizing the root position on the branch 1830 ... 0.0129452775"              
## [1833] "Optimizing the root position on the branch 1831 ... 0.0122746558"              
## [1834] "Optimizing the root position on the branch 1832 ... 0.0117985904"              
## [1835] "Optimizing the root position on the branch 1833 ... 0.0132077082"              
## [1836] "Optimizing the root position on the branch 1834 ... 0.0134719082"              
## [1837] "Optimizing the root position on the branch 1835 ... 0.0134950858"              
## [1838] "Optimizing the root position on the branch 1836 ... 0.0140573108"              
## [1839] "Optimizing the root position on the branch 1837 ... 0.0140589611"              
## [1840] "Optimizing the root position on the branch 1838 ... 0.0168456335"              
## [1841] "Optimizing the root position on the branch 1839 ... 0.0166872110"              
## [1842] "Optimizing the root position on the branch 1840 ... 0.0094634079"              
## [1843] "Optimizing the root position on the branch 1841 ... 0.0171638698"              
## [1844] "Optimizing the root position on the branch 1842 ... 0.0171262883"              
## [1845] "Optimizing the root position on the branch 1843 ... 0.0157822156"              
## [1846] "Optimizing the root position on the branch 1844 ... 0.0100677537"              
## [1847] "Optimizing the root position on the branch 1845 ... 0.0104267120"              
## [1848] "Optimizing the root position on the branch 1846 ... 0.0104208046"              
## [1849] "Optimizing the root position on the branch 1847 ... 0.0139324394"              
## [1850] "Optimizing the root position on the branch 1848 ... 0.0138876081"              
## [1851] "Optimizing the root position on the branch 1849 ... 0.0115421552"              
## [1852] "Optimizing the root position on the branch 1850 ... 0.0140690619"              
## [1853] "Optimizing the root position on the branch 1851 ... 0.0140728628"              
## [1854] "Optimizing the root position on the branch 1852 ... 0.0130864383"              
## [1855] "Optimizing the root position on the branch 1853 ... 0.0130721248"              
## [1856] "Optimizing the root position on the branch 1854 ... 0.0130076069"              
## [1857] "Optimizing the root position on the branch 1855 ... 0.0127542904"              
## [1858] "Optimizing the root position on the branch 1856 ... 0.0101853242"              
## [1859] "Optimizing the root position on the branch 1857 ... 0.0120288675"              
## [1860] "Optimizing the root position on the branch 1858 ... 0.0120631174"              
## [1861] "Optimizing the root position on the branch 1859 ... 0.0121089512"              
## [1862] "Optimizing the root position on the branch 1860 ... 0.0121632553"              
## [1863] "Optimizing the root position on the branch 1861 ... 0.0109302395"              
## [1864] "Optimizing the root position on the branch 1862 ... 0.0114745181"              
## [1865] "Optimizing the root position on the branch 1863 ... 0.0114768729"              
## [1866] "Optimizing the root position on the branch 1864 ... 0.0112787236"              
## [1867] "Optimizing the root position on the branch 1865 ... 0.0110838546"              
## [1868] "Optimizing the root position on the branch 1866 ... 0.0110884815"              
## [1869] "Optimizing the root position on the branch 1867 ... 0.0110962914"              
## [1870] "Optimizing the root position on the branch 1868 ... 0.0110966356"              
## [1871] "Optimizing the root position on the branch 1869 ... 0.0110958602"              
## [1872] "Optimizing the root position on the branch 1870 ... 0.0125313940"              
## [1873] "Optimizing the root position on the branch 1871 ... 0.0125004372"              
## [1874] "Optimizing the root position on the branch 1872 ... 0.0142229358"              
## [1875] "Optimizing the root position on the branch 1873 ... 0.0141641268"              
## [1876] "Optimizing the root position on the branch 1874 ... 0.0154987712"              
## [1877] "Optimizing the root position on the branch 1875 ... 0.0156269461"              
## [1878] "Optimizing the root position on the branch 1876 ... 0.0139328206"              
## [1879] "Optimizing the root position on the branch 1877 ... 0.0149954242"              
## [1880] "Optimizing the root position on the branch 1878 ... 0.0158361716"              
## [1881] "Optimizing the root position on the branch 1879 ... 0.0158183300"              
## [1882] "Optimizing the root position on the branch 1880 ... 0.0128794680"              
## [1883] "Optimizing the root position on the branch 1881 ... 0.0128864365"              
## [1884] "Optimizing the root position on the branch 1882 ... 0.0128397462"              
## [1885] "Optimizing the root position on the branch 1883 ... 0.0128565055"              
## [1886] "Optimizing the root position on the branch 1884 ... 0.0109358949"              
## [1887] "Optimizing the root position on the branch 1885 ... 0.0119389114"              
## [1888] "Optimizing the root position on the branch 1886 ... 0.0118939249"              
## [1889] "Optimizing the root position on the branch 1887 ... 0.0141459392"              
## [1890] "Optimizing the root position on the branch 1888 ... 0.0154828416"              
## [1891] "Optimizing the root position on the branch 1889 ... 0.0154744330"              
## [1892] "Optimizing the root position on the branch 1890 ... 0.0154476857"              
## [1893] "Optimizing the root position on the branch 1891 ... 0.0156066355"              
## [1894] "Optimizing the root position on the branch 1892 ... 0.0155944385"              
## [1895] "Optimizing the root position on the branch 1893 ... 0.0160251629"              
## [1896] "Optimizing the root position on the branch 1894 ... 0.0160260959"              
## [1897] "Optimizing the root position on the branch 1895 ... 0.0108353957"              
## [1898] "Optimizing the root position on the branch 1896 ... 0.0122936882"              
## [1899] "Optimizing the root position on the branch 1897 ... 0.0123116107"              
## [1900] "Optimizing the root position on the branch 1898 ... 0.0121303048"              
## [1901] "Optimizing the root position on the branch 1899 ... 0.0123251390"              
## [1902] "Optimizing the root position on the branch 1900 ... 0.0123227125"              
## [1903] "Optimizing the root position on the branch 1901 ... 0.0121471673"              
## [1904] "Optimizing the root position on the branch 1902 ... 0.0121130884"              
## [1905] "Optimizing the root position on the branch 1903 ... 0.0121207065"              
## [1906] "Optimizing the root position on the branch 1904 ... 0.0121659744"              
## [1907] "Optimizing the root position on the branch 1905 ... 0.0121620775"              
## [1908] "Optimizing the root position on the branch 1906 ... 0.0103861710"              
## [1909] "Optimizing the root position on the branch 1907 ... 0.0115292121"              
## [1910] "Optimizing the root position on the branch 1908 ... 0.0121001995"              
## [1911] "Optimizing the root position on the branch 1909 ... 0.0128221472"              
## [1912] "Optimizing the root position on the branch 1910 ... 0.0128190014"              
## [1913] "Optimizing the root position on the branch 1911 ... 0.0131968065"              
## [1914] "Optimizing the root position on the branch 1912 ... 0.0131364452"              
## [1915] "Optimizing the root position on the branch 1913 ... 0.0114774844"              
## [1916] "Optimizing the root position on the branch 1914 ... 0.0115053275"              
## [1917] "Optimizing the root position on the branch 1915 ... 0.0114251340"              
## [1918] "Optimizing the root position on the branch 1916 ... 0.0114159170"              
## [1919] "Optimizing the root position on the branch 1917 ... 0.0111086158"              
## [1920] "Optimizing the root position on the branch 1918 ... 0.0110862968"              
## [1921] "Optimizing the root position on the branch 1919 ... 0.0112686686"              
## [1922] "Optimizing the root position on the branch 1920 ... 0.0111990693"              
## [1923] "Optimizing the root position on the branch 1921 ... 0.0112926700"              
## [1924] "Optimizing the root position on the branch 1922 ... 0.0112559382"              
## [1925] "Optimizing the root position on the branch 1923 ... 0.0106364469"              
## [1926] "Optimizing the root position on the branch 1924 ... 0.0206062620"              
## [1927] "Optimizing the root position on the branch 1925 ... 0.0205367424"              
## [1928] "Optimizing the root position on the branch 1926 ... 0.0110600909"              
## [1929] "Optimizing the root position on the branch 1927 ... 0.0110539802"              
## [1930] "Optimizing the root position on the branch 1928 ... 0.0143289086"              
## [1931] "Optimizing the root position on the branch 1929 ... 0.0146846033"              
## [1932] "Optimizing the root position on the branch 1930 ... 0.0147782941"              
## [1933] "Optimizing the root position on the branch 1931 ... 0.0147905277"              
## [1934] "Optimizing the root position on the branch 1932 ... 0.0153595583"              
## [1935] "Optimizing the root position on the branch 1933 ... 0.0153552044"              
## [1936] "Optimizing the root position on the branch 1934 ... 0.0139998171"              
## [1937] "Optimizing the root position on the branch 1935 ... 0.0139771290"              
## [1938] "Optimizing the root position on the branch 1936 ... 0.0138723029"              
## [1939] "Optimizing the root position on the branch 1937 ... 0.0091662313"              
## [1940] "Optimizing the root position on the branch 1938 ... 0.0118463041"              
## [1941] "Optimizing the root position on the branch 1939 ... 0.0118505316"              
## [1942] "Optimizing the root position on the branch 1940 ... 0.0131747372"              
## [1943] "Optimizing the root position on the branch 1941 ... 0.0134460199"              
## [1944] "Optimizing the root position on the branch 1942 ... 0.0134567288"              
## [1945] "Optimizing the root position on the branch 1943 ... 0.0134603628"              
## [1946] "Optimizing the root position on the branch 1944 ... 0.0134597342"              
## [1947] "Optimizing the root position on the branch 1945 ... 0.0133406522"              
## [1948] "Optimizing the root position on the branch 1946 ... 0.0118196659"              
## [1949] "Optimizing the root position on the branch 1947 ... 0.0160964338"              
## [1950] "Optimizing the root position on the branch 1948 ... 0.0161824150"              
## [1951] "Optimizing the root position on the branch 1949 ... 0.0101233429"              
## [1952] "Optimizing the root position on the branch 1950 ... 0.0102799972"              
## [1953] "Optimizing the root position on the branch 1951 ... 0.0108135739"              
## [1954] "Optimizing the root position on the branch 1952 ... 0.0108848151"              
## [1955] "Optimizing the root position on the branch 1953 ... 0.0121296600"              
## [1956] "Optimizing the root position on the branch 1954 ... 0.0124074431"              
## [1957] "Optimizing the root position on the branch 1955 ... 0.0130925876"              
## [1958] "Optimizing the root position on the branch 1956 ... 0.0131060346"              
## [1959] "Optimizing the root position on the branch 1957 ... 0.0126627981"              
## [1960] "Optimizing the root position on the branch 1958 ... 0.0126577653"              
## [1961] "Optimizing the root position on the branch 1959 ... 0.0178627321"              
## [1962] "Optimizing the root position on the branch 1960 ... 0.0187977431"              
## [1963] "Optimizing the root position on the branch 1961 ... 0.0187986677"              
## [1964] "Optimizing the root position on the branch 1962 ... 0.0154289322"              
## [1965] "Optimizing the root position on the branch 1963 ... 0.0154413534"              
## [1966] "Optimizing the root position on the branch 1964 ... 0.0097891039"              
## [1967] "Optimizing the root position on the branch 1965 ... 0.0199917068"              
## [1968] "The pre-estimated root is on the branch 759, the branches are re-enumerated"   
## [1969] "Improving the root position with constraints around the pre-estimated root ..."
## [1970] "Optimizing the root position on the original branch 521 ... 0.0093212696"      
## [1971] "Optimizing the root position on the branch 1210 ... 0.0098405988"              
## [1972] "Optimizing the root position on the branch 761 ... 0.0098405988"               
## [1973] "Optimizing the root position on the branch 522 ... 0.0093032799"               
## [1974] "Optimizing the root position on the branch 719 ... 0.0093000863"               
## [1975] "Optimizing the root position on the branch 757 ... 0.0093060315"               
## [1976] "Optimizing the root position on the branch 720 ... 0.0093059911"               
## [1977] "Optimizing the root position on the branch 523 ... 0.0092910713"               
## [1978] "Optimizing the root position on the branch 569 ... 0.0093173899"               
## [1979] "Optimizing the root position on the branch 524 ... 0.0093173899"               
## [1980] "Optimizing the root position on the branch 520 ... 0.0093212000"               
## [1981] "Optimizing the root position on the branch 819 ... 0.0093212313"               
## [1982] "Optimizing the root position on the branch 31 ... 0.0093186229"                
## [1983] "Optimizing the root position on the branch 32 ... 0.0093302392"                
## [1984] "Optimizing the root position on the branch 30 ... 0.0093302392"                
## [1985] "The tree is rooted on the branch 523"                                          
## [1986] "Starting QPD on the new rooted tree ..."                                       
## [1987] "rate\t0.00249613\ttMRCA\t-0.379016\tobjective_function\t0.389403"              
## [1988] "Elapsed time: 6.85474 seconds"

The following R code extracts the evolutionary rate.

procresult <- function(fn){
  result <- readLines(fn)
  line <- result[grep("Tree 1 rate ",result)]
  line.split <- strsplit(line, " ")[[1]]
  list(rate=as.double(line.split[4]),tmrca=as.double(line.split[6]))
}
procresult(paste(stub,"_result.txt",sep=""))
## $rate
## [1] 0.002496
## 
## $tmrca
## [1] -0.379016

We can compare this with our result from rtt.

rate
## tr.rtt.tipdates 
##     0.002587992
tmrca
## [1] 0.732083

The time tree is saved with "_result_newick_date.txt" added at the end, as a Newick string.

lsd.tree <- read.tree(paste(stub,"_result_newick_date.txt",sep=""))

How close is the tree to the true tree? We can compare the patristic distances between tips using cophenetic.

true.tree <- read.tree("village_true.tre")
treedist(chronos.tree,true.tree)
##      symmetric.difference   branch.score.difference 
##                   96.0000                  211.7161 
##           path.difference quadratic.path.difference 
##                  877.5665                 5623.4920
treedist(lsd.tree,true.tree)
##      symmetric.difference   branch.score.difference 
##                  96.00000                  28.65891 
##           path.difference quadratic.path.difference 
##                 877.56652                1255.70485

Compare pairwise distances.

true.tree.d <- cophenetic.phylo(true.tree)
true.tree.ltri <- true.tree.d[lower.tri(true.tree.d)]
lsd.tree.d <- cophenetic.phylo(lsd.tree)
chronos.tree.d <- cophenetic.phylo(chronos.tree)
ids <- row.names(true.tree.d)
idx <- match(ids,row.names(lsd.tree.d))
lsd.tree.d2 <- lsd.tree.d[idx,idx]
lsd.tree.d2.ltri <- lsd.tree.d2[lower.tri(lsd.tree.d2)]
idx <- match(ids,row.names(chronos.tree.d))
chronos.tree.d2 <- chronos.tree.d[idx,idx]
chronos.tree.d2.ltri <- chronos.tree.d2[lower.tri(chronos.tree.d2)]

The LSD tree gives a very good correlation between the tip-to-tip distances.

summary(lm(lsd.tree.d2.ltri~true.tree.ltri))
## 
## Call:
## lm(formula = lsd.tree.d2.ltri ~ true.tree.ltri)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -6.3266 -0.6171 -0.1560  1.2715  5.3937 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    0.6949265  0.0148551   46.78   <2e-16 ***
## true.tree.ltri 1.0085462  0.0001893 5326.69   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.184 on 483634 degrees of freedom
## Multiple R-squared:  0.9832,    Adjusted R-squared:  0.9832 
## F-statistic: 2.837e+07 on 1 and 483634 DF,  p-value: < 2.2e-16
plot(lsd.tree.d2.ltri~true.tree.ltri,xlab="True distance (years)",ylab="Inferred distance (years)",main="LSD")

plot of chunk unnamed-chunk-22
Figure: plot of chunk unnamed-chunk-22

The chronos tree captures the overall rate, but the association is lower.

summary(lm(chronos.tree.d2.ltri~true.tree.ltri))
## 
## Call:
## lm(formula = chronos.tree.d2.ltri ~ true.tree.ltri)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -58.792  -1.510   0.886   3.125  35.868 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    4.849e+01  7.369e-02   658.0   <2e-16 ***
## true.tree.ltri 3.858e-01  9.393e-04   410.7   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.874 on 483634 degrees of freedom
## Multiple R-squared:  0.2586,    Adjusted R-squared:  0.2586 
## F-statistic: 1.687e+05 on 1 and 483634 DF,  p-value: < 2.2e-16
plot(chronos.tree.d2.ltri~true.tree.ltri,xlab="True distance (years)",ylab="Inferred distance (years)",main="chronos")

plot of chunk unnamed-chunk-23
Figure: plot of chunk unnamed-chunk-23