Module 7 Exercises
Overview
In todays workshop we will be looking to:
- Explore various tools for predicting protein secondary structure, function, and localization from primary sequence.
We will also see how to:
- Use R to perform repetitive structure prediction of proteins
Distance-based trees
There are many ways to build a phylogenetic tree, and each one has its pros and cons. Here, we will cover two of these. UPGMA and NJ Trees can be built from any sequences, whether DNA, RNA, or protein. What you are tracking is the history of mutations. Today we will build trees using an available data set for mammalian DNA.
Download the data file GGH.dna (https://bit.ly/3n0ql1n). Upload this data to your RStudio server account or store this data in your RStudio/R working directory on your personal computer. Go ahead and click on this file. You should see this open a new text tab. You will notice that this file contains aligned gamma-glutamyl hydrolase genes from 40 mammalian species. (If you are unfamiliar with the Latin names of mammals, they are listed here: http://evolution.gs.washington.edu/gs570/2016/data/species)
The file starts with two numbers, which are the dimensions of the data matrix. The first number is the number of sequences (lines in this case) and the second number is the number of aligned base pairs.
To perform our exercise today, we will need to load a number of libraries using the following commands:
library(ape)
library(phangorn)
library(seqinr)##
## Attaching package: 'seqinr'## The following objects are masked from 'package:ape':
##
## as.alignment, consensusThis loads the programs that we will be using and sets your working
directory to the relevant folder. You may get an error message regarding
seqinr and ape but you can ignore it. However,
always make note of such conflicts in your own work. Which function,
from which library, you use can have large (often unintentional) effects
on your analysis.
Now execute:
mammals<-read.dna("GGH.dna",format="sequential")
mammals_phyDat<-phyDat(mammals,type="DNA",levels=NULL)The first line of code reads the data file into the program and gives
it a handy name (mammals). The second line of code converts
the data file into a format that the program can manipulate. In other
words, you have renamed your data file twice, first calling it
mammals and then calling it
mammals_phyDat.
Next type:
dna_dist<-dist.ml(mammals_phyDat)This converts the original data matrix to a matrix of distances. To see what this looks like, type:
dna_dist## Dasypus Echinops Pongo Ictidomys Loxodonta
## Echinops 0.292734977
## Pongo 0.202036851 0.242208060
## Ictidomys 0.202149558 0.263548411 0.138162071
## Loxodonta 0.215647557 0.218454552 0.197325452 0.211277630
## Canis 0.211008250 0.296077362 0.181132821 0.194594452 0.223778310
## Macaca 0.207406803 0.244066797 0.032021939 0.143236086 0.204196778
## Myotis 0.215647557 0.301572723 0.206213919 0.213304262 0.220346567
## Felis 0.212008981 0.281664453 0.191186094 0.187130807 0.228885559
## Macropus 0.419932722 0.453185956 0.368837072 0.415875515 0.420353141
## Callithrix 0.226710654 0.262651822 0.073439912 0.162381089 0.221060294
## Erinaceus 0.251809411 0.359679810 0.221378793 0.260640871 0.280574587
## Sorex 0.333407422 0.368640004 0.249168106 0.283237814 0.299767731
## Choloepus 0.149477029 0.287969177 0.200656474 0.200761054 0.219559031
## Nomascus 0.241055104 0.272937079 0.038786551 0.166639617 0.213362474
## Ornithorhy 0.421638950 0.435200874 0.405268478 0.421662818 0.398871008
## Vicugna 0.197134441 0.282224679 0.168984723 0.170339775 0.199277374
## Monodelphi 0.434221089 0.455741150 0.409138573 0.431316060 0.435111135
## Sarcophilu 0.439186648 0.465738550 0.403795432 0.427582077 0.437028046
## Cavia 0.271115497 0.337071172 0.244971614 0.242430634 0.290423016
## Equus 0.211008250 0.292098572 0.164034316 0.191873361 0.213923640
## Rattus 0.310797600 0.339180878 0.281639745 0.264048248 0.326918846
## Ochotona 0.199741145 0.281553781 0.179568939 0.165696750 0.233043121
## Gorilla 0.197361664 0.238885006 0.015666366 0.147350977 0.197860943
## Procavia 0.220349941 0.254397775 0.198684723 0.219273726 0.146071694
## Dipodomys 0.268532897 0.339180878 0.225196750 0.193262907 0.271586000
## Otolemur 0.275037075 0.337071172 0.178476776 0.220949452 0.251198391
## Microcebus 0.150503022 0.230991353 0.078638471 0.106637770 0.168725710
## Mus 0.286355518 0.347679552 0.272539532 0.271586000 0.320442821
## Pteropus 0.225903066 0.300077371 0.183768825 0.202909342 0.210505751
## Tursiops 0.200255378 0.290117065 0.182464379 0.191873361 0.209739529
## Tarsius 0.219865756 0.271620242 0.113339584 0.171915372 0.242406714
## Mustela 0.207564744 0.291002457 0.183768825 0.178350532 0.219740343
## Tupaia 0.243406132 0.370154570 0.238412338 0.250975270 0.298006951
## Ailuropoda 0.236189759 0.334619474 0.191249952 0.203411687 0.232224449
## Pan 0.202036851 0.245930152 0.016908078 0.139427360 0.195958709
## Bos 0.209205365 0.320403819 0.199993267 0.191873361 0.238084895
## Homo 0.202036851 0.245930152 0.016908078 0.141964361 0.197325452
## Oryctolagu 0.211008250 0.266738549 0.157560388 0.157369267 0.230897490
## Sus 0.209205365 0.298074700 0.195911624 0.204196778 0.215323547
## Canis Macaca Myotis Felis Macropus
## Echinops
## Pongo
## Ictidomys
## Loxodonta
## Canis
## Macaca 0.186473298
## Myotis 0.165917747 0.206213919
## Felis 0.120256988 0.188485821 0.169134214
## Macropus 0.418281907 0.368837072 0.438553411 0.436256333
## Callithrix 0.178847830 0.087180754 0.223690833 0.218165946 0.401529692
## Erinaceus 0.224666470 0.216474143 0.257765811 0.230613525 0.449219974
## Sorex 0.270843305 0.267186991 0.277647726 0.266591015 0.414170522
## Choloepus 0.226876372 0.209915851 0.243495807 0.220372588 0.393233479
## Nomascus 0.198120341 0.053403335 0.220180678 0.212511900 0.376971642
## Ornithorhy 0.411424464 0.397138318 0.442985923 0.434415977 0.389557025
## Vicugna 0.147096010 0.174931235 0.166029054 0.147674567 0.420507428
## Monodelphi 0.435147372 0.410966728 0.452393500 0.442753265 0.169836211
## Sarcophilu 0.438943664 0.396601132 0.460133667 0.435147372 0.147117966
## Cavia 0.265530904 0.247874375 0.298615778 0.290851733 0.443133578
## Equus 0.154986383 0.162735052 0.181777369 0.167226648 0.413850580
## Rattus 0.304885407 0.286231596 0.315622237 0.310832438 0.483404611
## Ochotona 0.185501732 0.189483162 0.222119174 0.225428994 0.422875068
## Gorilla 0.175080008 0.043308650 0.204367069 0.191348420 0.381424780
## Procavia 0.220320839 0.210627651 0.234113452 0.216828736 0.410958273
## Dipodomys 0.238084895 0.222362547 0.262549746 0.250668310 0.459460039
## Otolemur 0.223551786 0.183798306 0.253193221 0.244429733 0.409445281
## Microcebus 0.132919878 0.076459820 0.183869131 0.190181695 0.363238536
## Mus 0.289308523 0.284697853 0.315622237 0.301343646 0.473734856
## Pteropus 0.154779709 0.185113513 0.147578585 0.168413258 0.402743684
## Tursiops 0.151141863 0.189157864 0.173805637 0.163306778 0.398543926
## Tarsius 0.171123928 0.121326880 0.217128914 0.200882791 0.385848298
## Mustela 0.092786568 0.181086658 0.167817847 0.102714714 0.431735183
## Tupaia 0.218152243 0.252229655 0.266306314 0.249858335 0.486623975
## Ailuropoda 0.102553787 0.196808965 0.177483179 0.104677792 0.401007630
## Pan 0.175830104 0.041901336 0.207595135 0.189834742 0.377205360
## Bos 0.170564615 0.205470243 0.204835242 0.185795235 0.443133578
## Homo 0.178476776 0.041901336 0.210365221 0.192539885 0.379312103
## Oryctolagu 0.174510262 0.167945668 0.211754111 0.188485821 0.431735183
## Sus 0.161438035 0.197269705 0.163306778 0.156818848 0.438553411
## Callithrix Erinaceus Sorex Choloepus Nomascus
## Echinops
## Pongo
## Ictidomys
## Loxodonta
## Canis
## Macaca
## Myotis
## Felis
## Macropus
## Callithrix
## Erinaceus 0.253200679
## Sorex 0.250015485 0.309926366
## Choloepus 0.195993787 0.304098831 0.279588900
## Nomascus 0.106431711 0.243831353 0.245615701 0.234360890
## Ornithorhy 0.428594827 0.410685176 0.454476738 0.378008466 0.410509271
## Vicugna 0.209047292 0.212473106 0.201999700 0.215761554 0.190175009
## Monodelphi 0.414159738 0.463483228 0.489049362 0.394184204 0.442933510
## Sarcophilu 0.418940630 0.449781982 0.486381143 0.399220436 0.430495007
## Cavia 0.251942761 0.310074958 0.318117930 0.294651540 0.265177965
## Equus 0.184089552 0.236288174 0.269012920 0.223074083 0.184812944
## Rattus 0.285303396 0.349567297 0.391715300 0.323985342 0.298088496
## Ochotona 0.223375833 0.242361867 0.218932806 0.213218563 0.210580442
## Gorilla 0.076471750 0.221996548 0.231121234 0.209251611 0.049197962
## Procavia 0.191729771 0.315824865 0.322038969 0.203855211 0.196318791
## Dipodomys 0.232240625 0.304098831 0.313325234 0.263325118 0.243069327
## Otolemur 0.196464612 0.281103342 0.276361417 0.221180142 0.208636287
## Microcebus 0.085046755 0.206608293 0.192976853 0.165449611 0.119643424
## Mus 0.273358883 0.328384888 0.355843484 0.319640652 0.284528804
## Pteropus 0.201831361 0.219975915 0.262330140 0.237647824 0.202707804
## Tursiops 0.194684196 0.234616864 0.274517532 0.226876372 0.198120341
## Tarsius 0.133854160 0.223942941 0.235549229 0.201623100 0.140454046
## Mustela 0.192907996 0.224940081 0.278946204 0.216587091 0.205736582
## Tupaia 0.287876137 0.295527309 0.251574552 0.298613232 0.241624180
## Ailuropoda 0.212055228 0.239692648 0.301787581 0.228404391 0.180254031
## Pan 0.071928576 0.218105465 0.250950638 0.206198371 0.047278352
## Bos 0.210862175 0.254922673 0.289397127 0.248140683 0.219301779
## Homo 0.071928576 0.219740343 0.256323795 0.208054807 0.047278352
## Oryctolagu 0.187604497 0.241324567 0.259927406 0.230698037 0.173179167
## Sus 0.198249265 0.234616864 0.261735725 0.228784770 0.214712270
## Ornithorhy Vicugna Monodelphi Sarcophilu Cavia
## Echinops
## Pongo
## Ictidomys
## Loxodonta
## Canis
## Macaca
## Myotis
## Felis
## Macropus
## Callithrix
## Erinaceus
## Sorex
## Choloepus
## Nomascus
## Ornithorhy
## Vicugna 0.417631396
## Monodelphi 0.384913477 0.430239115
## Sarcophilu 0.367411837 0.419859610 0.167688424
## Cavia 0.458181812 0.260196454 0.467992285 0.438943664
## Equus 0.436566483 0.164555576 0.427623797 0.425801901 0.261077650
## Rattus 0.478320244 0.297455741 0.543038658 0.498218001 0.340337764
## Ochotona 0.432436549 0.193485840 0.435253381 0.452086695 0.296578521
## Gorilla 0.401192382 0.167505432 0.413077411 0.411230890 0.240815488
## Procavia 0.439455596 0.179268368 0.477404741 0.439731761 0.294613610
## Dipodomys 0.458244696 0.234884261 0.498599118 0.470064475 0.285790989
## Otolemur 0.464789784 0.243657747 0.458181812 0.460071822 0.283613146
## Microcebus 0.380143023 0.164049645 0.370706042 0.374330055 0.235837966
## Mus 0.453855510 0.297455741 0.527864730 0.508570132 0.312796873
## Pteropus 0.415467600 0.152020306 0.435111135 0.435129311 0.276611884
## Tursiops 0.415556706 0.117486901 0.429497627 0.422089022 0.277537164
## Tarsius 0.409535292 0.199099199 0.374966961 0.384030516 0.239352763
## Mustela 0.421798073 0.157231299 0.438925510 0.431346201 0.264545259
## Tupaia 0.505260628 0.271653930 0.480890581 0.504758170 0.308313641
## Ailuropoda 0.398744783 0.152205633 0.443464673 0.459326841 0.277225566
## Pan 0.399162610 0.166029054 0.418324236 0.414717955 0.244971614
## Bos 0.403227661 0.151422945 0.464052683 0.465938866 0.292820423
## Homo 0.397138318 0.164555576 0.423890125 0.416553941 0.240638429
## Oryctolagu 0.442985923 0.167505432 0.454293274 0.454250318 0.259599089
## Sus 0.415556706 0.118870666 0.446576164 0.440839999 0.279051551
## Equus Rattus Ochotona Gorilla Procavia
## Echinops
## Pongo
## Ictidomys
## Loxodonta
## Canis
## Macaca
## Myotis
## Felis
## Macropus
## Callithrix
## Erinaceus
## Sorex
## Choloepus
## Nomascus
## Ornithorhy
## Vicugna
## Monodelphi
## Sarcophilu
## Cavia
## Equus
## Rattus 0.303313079
## Ochotona 0.230505006 0.305267966
## Gorilla 0.163973836 0.273082442 0.185501732
## Procavia 0.213039228 0.350505384 0.209220315 0.186759698
## Dipodomys 0.255621009 0.300940932 0.236022390 0.225078444 0.259518082
## Otolemur 0.229204191 0.324015031 0.230505006 0.179287593 0.191609199
## Microcebus 0.142361387 0.262301273 0.160227631 0.085212655 0.156288302
## Mus 0.304885407 0.147829015 0.324226322 0.265129027 0.324837480
## Pteropus 0.183768825 0.310832438 0.198566618 0.185634107 0.231619684
## Tursiops 0.177152272 0.297056526 0.193485840 0.176479915 0.220320839
## Tarsius 0.195544690 0.273974295 0.194563650 0.116524303 0.212980496
## Mustela 0.163882551 0.300940932 0.205623664 0.180643544 0.228117682
## Tupaia 0.255724128 0.362569987 0.273582010 0.240125646 0.244066796
## Ailuropoda 0.162278336 0.346207378 0.217270379 0.189483162 0.228739381
## Pan 0.162735052 0.278594055 0.185501732 0.004442806 0.196318791
## Bos 0.213761437 0.300178279 0.193485840 0.203587003 0.247641740
## Homo 0.164034316 0.280115354 0.185501732 0.005557632 0.198684723
## Oryctolagu 0.199993267 0.283167240 0.110013808 0.151672877 0.220320839
## Sus 0.161438035 0.309622265 0.197509994 0.193485840 0.230140946
## Dipodomys Otolemur Microcebus Mus Pteropus
## Echinops
## Pongo
## Ictidomys
## Loxodonta
## Canis
## Macaca
## Myotis
## Felis
## Macropus
## Callithrix
## Erinaceus
## Sorex
## Choloepus
## Nomascus
## Ornithorhy
## Vicugna
## Monodelphi
## Sarcophilu
## Cavia
## Equus
## Rattus
## Ochotona
## Gorilla
## Procavia
## Dipodomys
## Otolemur 0.251198391
## Microcebus 0.213459426 0.147127090
## Mus 0.283806230 0.320792684 0.279315235
## Pteropus 0.250668310 0.226617878 0.150503021 0.304493257
## Tursiops 0.242430634 0.243524436 0.154332910 0.298615778 0.175750931
## Tarsius 0.216333324 0.179754633 0.102956675 0.283881969 0.192006948
## Mustela 0.250135241 0.239530679 0.154332910 0.300940932 0.159790150
## Tupaia 0.291087352 0.284282010 0.255795537 0.342568802 0.256867732
## Ailuropoda 0.282850951 0.275158025 0.179153122 0.334894931 0.169253998
## Pan 0.225196750 0.179803623 0.080823469 0.268030500 0.185113513
## Bos 0.255621009 0.270010758 0.183869131 0.303313079 0.182426542
## Homo 0.225196750 0.179803623 0.083014851 0.266533494 0.186460617
## Oryctolagu 0.200066433 0.212373191 0.125915936 0.277075836 0.190516508
## Sus 0.239530679 0.242080044 0.178864942 0.309622265 0.161270446
## Tursiops Tarsius Mustela Tupaia Ailuropoda
## Echinops
## Pongo
## Ictidomys
## Loxodonta
## Canis
## Macaca
## Myotis
## Felis
## Macropus
## Callithrix
## Erinaceus
## Sorex
## Choloepus
## Nomascus
## Ornithorhy
## Vicugna
## Monodelphi
## Sarcophilu
## Cavia
## Equus
## Rattus
## Ochotona
## Gorilla
## Procavia
## Dipodomys
## Otolemur
## Microcebus
## Mus
## Pteropus
## Tursiops
## Tarsius 0.172842141
## Mustela 0.159967797 0.192006948
## Tupaia 0.257477489 0.269530501 0.245289177
## Ailuropoda 0.169416997 0.224603911 0.071085914 0.187024545
## Pan 0.175830104 0.114930253 0.182426542 0.238412338 0.189406067
## Bos 0.128479627 0.206259158 0.170452896 0.252229655 0.176624260
## Homo 0.181132821 0.114930253 0.185113513 0.240125646 0.193098382
## Oryctolagu 0.173192738 0.181492756 0.182426542 0.240125646 0.206166541
## Sus 0.116179228 0.179754633 0.165192023 0.257477489 0.176624260
## Pan Bos Homo Oryctolagu
## Echinops
## Pongo
## Ictidomys
## Loxodonta
## Canis
## Macaca
## Myotis
## Felis
## Macropus
## Callithrix
## Erinaceus
## Sorex
## Choloepus
## Nomascus
## Ornithorhy
## Vicugna
## Monodelphi
## Sarcophilu
## Cavia
## Equus
## Rattus
## Ochotona
## Gorilla
## Procavia
## Dipodomys
## Otolemur
## Microcebus
## Mus
## Pteropus
## Tursiops
## Tarsius
## Mustela
## Tupaia
## Ailuropoda
## Pan
## Bos 0.197269705
## Homo 0.004191419 0.201358766
## Oryctolagu 0.152421181 0.190503760 0.154986383
## Sus 0.194555998 0.153702685 0.194555998 0.190503760Which writes the matrix to the screen. The matrix is large so it wraps around and is a little hard to read. Scroll to the top of the output so you can see the beginning of the matrix. The first column of numbers is labeled Dasypus, the second is Echinops, the third is Pongo, etc. The species names are in the same order down the left side of the matrix. This shows you that the distance between Dasypus (armadillo) and Echinops (tenrec) is 0.292605364, that between Dasypus (armadillo) and Pongo (orangutan) is 0.202628039, etc.
Question 1
What sequences are most distant from Dasypus? What do those animals have in common?
Now use the distance matrix to produce a tree. We will compare two trees, one produced by the Unweighted Pair Group Method with Averaging (UPGMA), and one produced by the Neighbor Joining (NJ) method.
mammals_UPGMA<-upgma(dna_dist)To see this tree:
plot(mammals_UPGMA)
To save the tree, do it two ways.
First, just save as a pdf, using the Export menu in your Plots menu.
Second, write it to a file:
write.tree(mammals_UPGMA,file="UPGMA_tree")You can download this file to your desktop via the Files panel (same
window that you can visualize your plots).
Files -> More -> Export. You can use this tree file
to view your tree in tree viewing software such as FigTree (http://tree.bio.ed.ac.uk/software/figtree/).
Now produce a neighbor-joining tree:
mammals_NJ<-NJ(dna_dist)
plot(mammals_NJ)
Save the plot as a pdf and
write.tree(mammals_NJ,file="NJ_tree")Now, let’s compare the two trees.
Question 2: What is different? What is the same?
Sometimes, it is useful to redefine the root of a tree (or reroot a tree). Let’s take a look at how we can do this in R.
First, plot your tree as we’ve done before. Except, now we will use
an additional command from the ape package called
nodelabels(). This will label all tree nodes with their
numerical node ID in the tree. We can then use this node ID to define
the root of the tree at this node.
plot(mammals_NJ)
nodelabels()
Now, let’s reroot the tree around node 77 (the clade with marsupials). It is often best to export these labeled trees to a PDF to see them more clearly.
marsupial_rooted_NJ <- root(mammals_NJ, node = 77)
plot(marsupial_rooted_NJ)
See if you can do the same for your UPGMA tree!
Maximum likelihood trees, maximum parsimony trees, and bootstrapping
We just built phylogenetic trees using the clustering algorithms UPGMA and neighbor joining. Now, we will use the same data set, GGH.dna, to build a tree using maximum likelihood and maximum parsimony. We will also do a bootstrap analysis.
Maximum likelihood
To get started, let’s do a preliminary neighbor-joining tree. This is exactly as we did above. We need to make this tree as a starting point for a maximum likelihood optimization of the tree.
dna_dist<-dist.ml(mammals_phyDat)
mammals_NJ<-NJ(dna_dist)Now we can calculate the likelihood of the NJ tree.
fit<-pml(mammals_NJ, mammals_phyDat)
print(fit)##
## loglikelihood: -15999.05
##
## unconstrained loglikelihood: -6087.627
##
## Rate matrix:
## a c g t
## a 0 1 1 1
## c 1 0 1 1
## g 1 1 0 1
## t 1 1 1 0
##
## Base frequencies:
## 0.25 0.25 0.25 0.25The output will give you a loglikelihood. Write down this number somewhere.
Now see if you can improve the likelihood of the NJ tree by adjusting the tree topology and branch lengths. This, in effect, is the process of maximum likelihood optimization. Let’s do this first using a Jukes-Cantor model of evolution.
fitJC<-optim.pml(fit,model="JC",rearrangement="stochastic")## optimize edge weights: -15999.05 --> -15909.49
## optimize edge weights: -15909.49 --> -15909.49
## optimize topology: -15909.49 --> -15881.9
## optimize topology: -15881.9 --> -15876.25
## optimize topology: -15876.25 --> -15866.6
## optimize topology: -15866.6 --> -15854.15
## optimize topology: -15854.15 --> -15848.86
## NNI moves: 12
## optimize edge weights: -15848.86 --> -15848.86
## optimize topology: -15848.86 --> -15841.21
## optimize topology: -15841.21 --> -15837.22
## optimize topology: -15837.22 --> -15833.67
## optimize topology: -15833.67 --> -15830.68
## optimize topology: -15830.68 --> -15821.98
## NNI moves: 7
## optimize edge weights: -15821.98 --> -15821.98
## optimize topology: -15821.98 --> -15818.52
## optimize topology: -15818.52 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 2
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 0
## -16590.2 -> -16577.1
## -16588.41 -> -16575.01
## -16588.94 -> -16575.01
## -16533.94 -> -16531.02
## -16572.74 -> -16557.79
## -16578.21 -> -16563.22
## -16578.49 -> -16563.64
## -16562.43 -> -16558.54
## -16578.45 -> -16562.31
## -16578.61 -> -16573.85
## -16581.57 -> -16575.03
## -16566.58 -> -16564.67
## -16570.2 -> -16565.8
## -16580.69 -> -16572.1
## -16606.65 -> -16540.69
## -16610.04 -> -16540.69
## -16507.55 -> -16414.98
## -16565.45 -> -16482.9
## -16592.66 -> -16506.71
## -16566.7 -> -16500.78
## -16579.03 -> -16575.13
## -16596.05 -> -16574.06
## -16532.41 -> -16525.79
## -16579.53 -> -16563.8
## -16580.77 -> -16575.57
## -16583.63 -> -16576.98
## -16555.55 -> -16555.54
## -16575.9 -> -16571.67
## -16560.72 -> -16552.08
## -16577.96 -> -16572.53
## -16585.14 -> -16573.91
## -16588.14 -> -16575.75
## -16601.13 -> -16539.03
## -16597.74 -> -16537.31
## -16405.57 -> -16405.46
## -16563.86 -> -16486.67
## -16578.77 -> -16565.12
## -16573.24 -> -16559.46
## optimize edge weights: -16097.9 --> -16066.45
## optimize topology: -16578.44 --> -16066.45
## -16064.71 -> -16059.95
## -16065.17 -> -16061
## -16066.45 -> -15938.49
## -16066.45 -> -16065.84
## -16066.45 -> -16063.58
## -16066.45 -> -16033
## -16064.05 -> -16063.94
## -16057.5 -> -16052.25
## -16063.99 -> -16062.14
## optimize edge weights: -15884.25 --> -15880.19
## optimize topology: -16066.45 --> -15880.19
## -15880.19 -> -15880.19
## -15880.19 -> -15880.19
## -15880.19 -> -15877.56
## -15880.19 -> -15879.51
## -15878.09 -> -15877.81
## -15878.18 -> -15877.7
## -15870.81 -> -15867.49
## -15877.62 -> -15875.68
## -15877.28 -> -15877.13
## optimize edge weights: -15864.87 --> -15864.69
## optimize topology: -15880.19 --> -15864.69
## -15864.69 -> -15864.69
## -15864.69 -> -15864.69
## -15864.69 -> -15839.99
## -15864.69 -> -15863.93
## -15861.77 -> -15861.63
## -15864.68 -> -15864.6
## -15862.04 -> -15859.95
## optimize edge weights: -15836.82 --> -15836.49
## optimize topology: -15864.69 --> -15836.49
## -15827.94 -> -15815.62
## optimize edge weights: -15815.62 --> -15815.55
## optimize topology: -15836.49 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 16
## [1] "Ratchet iteration 1 , best pscore so far: -15815.5475996227"
## -16093.06 -> -16090.17
## -16084.37 -> -16084.32
## -16092.38 -> -16089.83
## -16094.35 -> -16070.53
## -16088.97 -> -16066.89
## -16084.14 -> -16059.97
## -16019.2 -> -16014.44
## -16091.4 -> -16090.24
## -16099.04 -> -16080.58
## -16099.69 -> -16080.58
## -16091.62 -> -16091.5
## -16091.61 -> -16091.5
## -16058.45 -> -16057.75
## -16080.87 -> -16071.61
## -16074.78 -> -16072.81
## -16088.63 -> -16082.79
## -16090.44 -> -16087.42
## -16094.6 -> -16090.81
## -16098.02 -> -16090.7
## -16094.11 -> -16089.38
## -16082.53 -> -16075.28
## -16091.64 -> -16083.2
## -16088.64 -> -16085.81
## -16088.72 -> -16085.48
## -16088.31 -> -16083.37
## -16095.19 -> -16087.94
## -16096.26 -> -16082.04
## -16087.79 -> -16075.5
## -16086.64 -> -16074.84
## -16095.3 -> -16083.95
## -16066.03 -> -16065.28
## -16097.15 -> -16084.15
## -16081.05 -> -16077.21
## -16093.71 -> -16085.59
## -16073.35 -> -16069.69
## -16092.4 -> -16083.4
## optimize edge weights: -15940.11 --> -15912.33
## optimize topology: -16091.5 --> -15912.33
## -15903.73 -> -15897.2
## -15910.08 -> -15909.39
## -15904.49 -> -15900.72
## -15910.71 -> -15909.47
## -15912.33 -> -15907.1
## -15912.33 -> -15912.33
## -15900.77 -> -15894.8
## -15909.93 -> -15907.03
## -15912.05 -> -15911.4
## -15905.6 -> -15902.65
## -15905.84 -> -15894.73
## -15912.15 -> -15912.13
## optimize edge weights: -15850.11 --> -15849.63
## optimize topology: -15912.33 --> -15849.63
## -15850.14 -> -15848.75
## -15849.63 -> -15843.21
## -15847.78 -> -15844.59
## -15849.69 -> -15849
## -15850.06 -> -15848.96
## -15840.55 -> -15834.89
## optimize edge weights: -15828.6 --> -15828.38
## optimize topology: -15849.63 --> -15828.38
## -15828.23 -> -15827.81
## -15827.95 -> -15827
## -15827.08 -> -15826.23
## -15827.68 -> -15826.65
## optimize edge weights: -15826.23 --> -15826.18
## optimize topology: -15828.38 --> -15826.18
## optimize edge weights: -15826.18 --> -15826.18
## optimize topology: -15826.18 --> -15826.18
## NNI moves: 13
## [1] "Ratchet iteration 2 , best pscore so far: -15815.5475996227"
## -16685.55 -> -16683.44
## -16690.07 -> -16687.22
## -16696.67 -> -16644.03
## -16663.64 -> -16631.77
## -16543.23 -> -16460.61
## -16808.53 -> -16522.22
## -16563.68 -> -16473.34
## -16749.27 -> -16538.95
## -16510.28 -> -16299.09
## -16699.27 -> -16549.91
## -16705.79 -> -16653.88
## -16691.71 -> -16642.81
## -16686.39 -> -16686.05
## -16688.49 -> -16686.8
## -16693.17 -> -16693.08
## -16693.35 -> -16693.23
## -16678.71 -> -16678.71
## -16693.74 -> -16688.26
## -16694.29 -> -16694.09
## -16694.24 -> -16694.09
## -16670.34 -> -16670.26
## -16688.64 -> -16684.45
## -16692.61 -> -16678.19
## -16700 -> -16682.76
## -16653.23 -> -16653.13
## -16683.4 -> -16673.49
## -16694.03 -> -16691.26
## -16676.37 -> -16676.29
## -16699.36 -> -16693.17
## -16692.03 -> -16691.5
## -16692.28 -> -16692
## optimize edge weights: -16203.56 --> -16154.7
## optimize topology: -16694.14 --> -16154.7
## -16154.29 -> -16154.2
## -16151.28 -> -16142.14
## -16154.7 -> -15968.56
## -16151.86 -> -16151.31
## -16153.94 -> -16153.92
## -16143.88 -> -16139.25
## -16152.83 -> -16151.89
## optimize edge weights: -15946.23 --> -15920.58
## optimize topology: -16154.7 --> -15920.58
## -15913.07 -> -15898.49
## -15913.07 -> -15844.31
## -15919.02 -> -15918.59
## optimize edge weights: -15842.29 --> -15841.19
## optimize topology: -15920.58 --> -15841.19
## -15831.71 -> -15816.5
## -15839.44 -> -15837.22
## -15840.95 -> -15840.31
## optimize edge weights: -15815.81 --> -15815.55
## optimize topology: -15841.19 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 14
## [1] "Ratchet iteration 3 , best pscore so far: -15815.547426795"
## -16276.06 -> -16268.49
## -16276.15 -> -16267.11
## -16272.73 -> -16267.02
## -16273.66 -> -16268.68
## -16273.93 -> -16271.09
## -16259.1 -> -16217.82
## -16270.64 -> -16226.8
## -16352.01 -> -16259.5
## -16354.55 -> -16259.5
## -16021.27 -> -16020.85
## -16250.8 -> -16151.81
## -16322.19 -> -16259.83
## -16327.36 -> -16259.04
## -16256.82 -> -16250.82
## -16271.21 -> -16266.06
## -16271.57 -> -16266.82
## -16276.54 -> -16269.28
## -16268.42 -> -16258.19
## -16271.18 -> -16261.86
## -16272.43 -> -16263.98
## -16280.65 -> -16264.36
## -16256.24 -> -16256.21
## -16271.55 -> -16268.13
## -16269.39 -> -16266.23
## -16273.41 -> -16269.54
## -16272.64 -> -16270.1
## -16263.67 -> -16256.15
## -16264.61 -> -16257.54
## -16253.1 -> -16253.04
## -16275.94 -> -16269.03
## -16271.12 -> -16257.33
## -16277.25 -> -16263.95
## -16273.07 -> -16269.52
## -16271.02 -> -16266.98
## -16249.46 -> -16249.44
## -16271.72 -> -16266.27
## -16268.49 -> -16268.48
## -16270.14 -> -16269.5
## optimize edge weights: -15945.96 --> -15919.98
## optimize topology: -16271.18 --> -15919.98
## -15917.33 -> -15915.39
## -15919.24 -> -15919.17
## -15918.45 -> -15918.11
## -15916.39 -> -15915.7
## -15919.79 -> -15919.73
## -15901.77 -> -15885.11
## -15911.87 -> -15907.8
## -15918.35 -> -15917.25
## -15910.91 -> -15906.19
## -15917.3 -> -15916.95
## optimize edge weights: -15854.55 --> -15854.11
## optimize topology: -15919.98 --> -15854.11
## -15851.33 -> -15849.68
## -15850.67 -> -15832
## -15853.3 -> -15852.55
## -15854.11 -> -15854.11
## -15854.11 -> -15854.11
## -15850.79 -> -15850.05
## optimize edge weights: -15823.56 --> -15823.45
## optimize topology: -15854.11 --> -15823.45
## -15821.96 -> -15815.57
## optimize edge weights: -15815.57 --> -15815.55
## optimize topology: -15823.45 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 15
## [1] "Ratchet iteration 4 , best pscore so far: -15815.547426795"
## -16422.12 -> -16352.45
## -16405.83 -> -16352.26
## -16392.54 -> -16373.53
## -16394.93 -> -16376.54
## -16238.43 -> -16238.21
## -16406.73 -> -16316.99
## -16395.54 -> -16381.66
## -16387.69 -> -16376.21
## -16425.12 -> -16386.98
## -16427.52 -> -16387.53
## -16388.46 -> -16387.76
## -16391.83 -> -16390.01
## -16392.18 -> -16372.16
## -16403.76 -> -16377.8
## -16367.85 -> -16354.33
## -16377.23 -> -16362.47
## -16396.54 -> -16363.84
## -16391.76 -> -16358.94
## -16370.13 -> -16346.01
## -16394.91 -> -16370.03
## -16423.33 -> -16391.45
## -16421.41 -> -16391.45
## -16370.4 -> -16364.59
## -16396.7 -> -16383.54
## -16316.27 -> -16311.51
## -16379.65 -> -16354.1
## -16344.68 -> -16324.78
## -16390.46 -> -16355.1
## -16397.56 -> -16364.33
## -16408.47 -> -16365.6
## -16391.42 -> -16390.38
## -16393.87 -> -16392.92
## -16396.13 -> -16390.05
## -16393.96 -> -16389.21
## -16378.29 -> -16363.12
## -16402.61 -> -16370.35
## -16399.66 -> -16394.57
## -16394.77 -> -16394.54
## -16394.76 -> -16394.54
## -16331.51 -> -16325.18
## -16393.18 -> -16368
## optimize edge weights: -15975.15 --> -15954.28
## optimize topology: -16394.6 --> -15954.28
## -15951.15 -> -15940.18
## -15951.07 -> -15949.89
## -15954.28 -> -15952.4
## -15954.28 -> -15953.77
## -15945.65 -> -15938.05
## -15954.41 -> -15953.81
## -15952.29 -> -15951.12
## -15953.45 -> -15951.77
## -15953.22 -> -15945.11
## -15953.8 -> -15952.84
## -15950.81 -> -15947.73
## -15948.77 -> -15935.82
## -15945.68 -> -15940.98
## optimize edge weights: -15887.23 --> -15884.69
## optimize topology: -15954.28 --> -15884.69
## -15877.97 -> -15863.46
## -15884.69 -> -15883.7
## -15884.69 -> -15884.14
## -15884.32 -> -15883.76
## -15882.45 -> -15880.89
## -15884.43 -> -15884.06
## -15868.39 -> -15851.9
## optimize edge weights: -15830.74 --> -15830.3
## optimize topology: -15884.69 --> -15830.3
## -15829.71 -> -15829.34
## -15829.85 -> -15829.7
## -15830.11 -> -15829.43
## -15827.82 -> -15826.38
## -15828.39 -> -15827.78
## optimize edge weights: -15823.58 --> -15823.45
## optimize topology: -15830.3 --> -15823.45
## -15821.96 -> -15815.57
## optimize edge weights: -15815.57 --> -15815.55
## optimize topology: -15823.45 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 16
## [1] "Ratchet iteration 5 , best pscore so far: -15815.5472187106"
## -16954.85 -> -16683.25
## -16802.96 -> -16668.82
## -16824.54 -> -16602.23
## -16345.68 -> -16267.03
## -16933.79 -> -16778.61
## -16932.78 -> -16778.61
## -16867.68 -> -16628.13
## -16870.09 -> -16628.14
## -16873.14 -> -16802.35
## -16865 -> -16802.35
## -16807.12 -> -16785.07
## -16836.4 -> -16801.31
## -16808.24 -> -16796.13
## -16836.46 -> -16818.43
## -16825.45 -> -16825.14
## -16811.47 -> -16806.99
## -16838.42 -> -16821.26
## -16831.85 -> -16809.6
## -16839.11 -> -16809.6
## -16834.37 -> -16823.58
## -16834.4 -> -16823.58
## -16789.66 -> -16775.48
## -16825.71 -> -16796.33
## -16825.84 -> -16804.88
## -16786.46 -> -16781.18
## -16814.99 -> -16814.98
## -16816.59 -> -16815.36
## -16825.96 -> -16823.28
## -16821.4 -> -16819.21
## -16823.92 -> -16821.51
## -16764.03 -> -16757.74
## -16825.22 -> -16809.52
## -16840.61 -> -16821.23
## -16847.64 -> -16821.27
## -16801.87 -> -16800.96
## -16825.68 -> -16822.45
## -16820.51 -> -16817.39
## -16832.59 -> -16822.35
## -16825.08 -> -16820.61
## -16821.65 -> -16816.21
## optimize edge weights: -16137.42 --> -16068.3
## optimize topology: -16826 --> -16068.3
## -16065.81 -> -15994.65
## -16060.29 -> -16041.08
## -16065.71 -> -16048.09
## -16067.89 -> -16067.6
## -16060.67 -> -16054.24
## -16067.34 -> -16067.16
## -16068.3 -> -16031.04
## -16068.3 -> -16068.29
## -16065.39 -> -16065.19
## -16068.04 -> -16067.99
## -16058.5 -> -16053.31
## -16068.09 -> -16068.02
## -16066.36 -> -16066.1
## -16065.5 -> -16064.8
## -16068.97 -> -16066.9
## optimize edge weights: -15918.74 --> -15917.91
## optimize topology: -16068.3 --> -15917.91
## -15912.38 -> -15884.47
## -15909.07 -> -15900.44
## -15918.23 -> -15917.58
## -15915.44 -> -15915.37
## -15916.7 -> -15916.57
## -15912.17 -> -15911.45
## -15915.4 -> -15913.43
## optimize edge weights: -15860.44 --> -15859.67
## optimize topology: -15917.91 --> -15859.67
## -15857.33 -> -15857.27
## -15859.24 -> -15859.17
## -15859.69 -> -15859.65
## -15859.29 -> -15857.15
## optimize edge weights: -15854.73 --> -15854.43
## optimize topology: -15859.67 --> -15854.43
## -15845.93 -> -15833.15
## optimize edge weights: -15833.15 --> -15833.04
## optimize topology: -15854.43 --> -15833.04
## -15832.84 -> -15832.82
## optimize edge weights: -15832.82 --> -15832.79
## optimize topology: -15833.04 --> -15832.79
## -15821.25 -> -15815.69
## optimize edge weights: -15815.69 --> -15815.55
## optimize topology: -15832.79 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 17
## [1] "Ratchet iteration 6 , best pscore so far: -15815.5454303595"
## -16108.28 -> -16089.39
## -16099.2 -> -16087.63
## -16033.24 -> -16032.56
## -16098.73 -> -16068.39
## -16110.72 -> -16087.07
## -16119.05 -> -16087.02
## -16097.61 -> -16090.07
## -16082.86 -> -16082.67
## -16086.44 -> -16084.62
## -16106.08 -> -16091.84
## -16099.34 -> -16089.12
## -16083.17 -> -16079.88
## -16087.92 -> -16083.87
## -16107.42 -> -16062.44
## -16095.35 -> -16060.52
## -16077.97 -> -16043.83
## -16092.54 -> -16055.82
## -16093.96 -> -16092.54
## -16093.96 -> -16092.54
## -16081.21 -> -16060.99
## -16094.99 -> -16063.75
## -16086.34 -> -16063.81
## -16066.87 -> -16051.99
## -16047.81 -> -16031.5
## -16093.12 -> -16057.61
## -16098.76 -> -16085.98
## -16098.3 -> -16085.98
## -16097.79 -> -16092.04
## -16086.05 -> -16086.03
## -16091.59 -> -16091.57
## optimize edge weights: -15963.05 --> -15942.86
## optimize topology: -16092.59 --> -15942.86
## -15938.45 -> -15933.78
## -15938.38 -> -15936.75
## -15940.88 -> -15938.35
## -15941.37 -> -15938.44
## -15938.89 -> -15936.09
## -15938.66 -> -15936.76
## -15936.86 -> -15935.01
## -15933.54 -> -15929.09
## -15941.45 -> -15941.07
## -15942.86 -> -15937.67
## -15942.86 -> -15942.86
## -15942.86 -> -15942.6
## optimize edge weights: -15920.21 --> -15919.89
## optimize topology: -15942.86 --> -15919.89
## -15915.55 -> -15910.62
## -15916.95 -> -15915.88
## -15917.72 -> -15914.4
## -15918.08 -> -15914
## -15919.89 -> -15919.1
## -15919.89 -> -15896.81
## optimize edge weights: -15881.85 --> -15881.24
## optimize topology: -15919.89 --> -15881.24
## -15869.95 -> -15848.99
## -15882.04 -> -15880.8
## -15881.18 -> -15879.66
## -15881.25 -> -15881.23
## -15873.44 -> -15853.28
## -15879.82 -> -15878.65
## optimize edge weights: -15818.43 --> -15817.51
## optimize topology: -15881.24 --> -15817.51
## -15816.45 -> -15815.59
## optimize edge weights: -15815.59 --> -15815.55
## optimize topology: -15817.51 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 17
## [1] "Ratchet iteration 7 , best pscore so far: -15815.5452452212"
## -16284.2 -> -16283.24
## -16260.6 -> -16260.28
## -16264.61 -> -16263.48
## -16317.81 -> -16275.57
## -16318.88 -> -16275.57
## -16229.33 -> -16228.46
## -16277.06 -> -16253.65
## -16299.74 -> -16281.99
## -16297.22 -> -16279.31
## -16279.24 -> -16273.6
## -16272.25 -> -16257.02
## -16287.92 -> -16278.65
## -16280.94 -> -16280.51
## -16281.02 -> -16280.32
## -16273.62 -> -16267.21
## -16289.72 -> -16279.86
## -16261.82 -> -16257.81
## -16285.73 -> -16278.01
## -16265.81 -> -16261.68
## -16287.26 -> -16276.76
## -16278.53 -> -16269.84
## -16283.36 -> -16272.04
## -16270.7 -> -16261.7
## -16287.78 -> -16264.2
## -16203.16 -> -16172.47
## -16279.66 -> -16239.27
## -16363.62 -> -16275.55
## -16364.16 -> -16275.55
## -16219.54 -> -16212.84
## -16277.65 -> -16263.23
## -16283.55 -> -16268.21
## -16279.79 -> -16266.22
## -16284.01 -> -16268.11
## -16283.62 -> -16268.11
## -16223.61 -> -16222.03
## -16282.56 -> -16261.05
## -16308.07 -> -16281.28
## -16308.92 -> -16281.28
## -16298.63 -> -16262.16
## -16302.73 -> -16262.58
## -16279.83 -> -16276.42
## -16258.18 -> -16239.38
## -16272.59 -> -16245.27
## -16286.18 -> -16280.92
## -16275.52 -> -16274.03
## optimize edge weights: -16008.49 --> -15965.65
## optimize topology: -16283.27 --> -15965.65
## -15965.65 -> -15956.03
## -15965.65 -> -15957.86
## -15964.01 -> -15963.73
## -15964.78 -> -15964.7
## -15963.61 -> -15963.15
## -15959.03 -> -15950.14
## -15964.53 -> -15964.1
## -15965.65 -> -15961.25
## -15965.65 -> -15952.64
## -15963.44 -> -15962.58
## -15960.44 -> -15954.68
## -15964.1 -> -15962.03
## -15964.12 -> -15962.7
## -15965.65 -> -15941.05
## -15965.65 -> -15965.19
## -15965.65 -> -15965.65
## -15964.52 -> -15964.25
## optimize edge weights: -15903.44 --> -15902.73
## optimize topology: -15965.65 --> -15902.73
## -15899.56 -> -15898.5
## -15902.03 -> -15901.86
## -15902.73 -> -15898.15
## -15902.64 -> -15902.25
## -15902.11 -> -15891.43
## -15895.33 -> -15887.91
## -15901.34 -> -15900.62
## -15901.95 -> -15900.95
## -15894.08 -> -15874.6
## -15901.4 -> -15901.1
## optimize edge weights: -15854.6 --> -15854.16
## optimize topology: -15902.73 --> -15854.16
## -15848.53 -> -15847.99
## -15853.23 -> -15853.13
## -15854.16 -> -15848.42
## -15854.16 -> -15853.27
## -15854.16 -> -15843.6
## optimize edge weights: -15837.42 --> -15836.53
## optimize topology: -15854.16 --> -15836.53
## -15835.01 -> -15834.67
## -15836.53 -> -15818.71
## optimize edge weights: -15816.67 --> -15816.31
## optimize topology: -15836.53 --> -15816.31
## -15815.83 -> -15815.56
## optimize edge weights: -15815.56 --> -15815.55
## optimize topology: -15816.31 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 20
## [1] "Ratchet iteration 8 , best pscore so far: -15815.5452452212"
## -16153.09 -> -16141.6
## -16171.06 -> -16152.89
## -16140.37 -> -16126.19
## -16162.11 -> -16133.96
## -16126.07 -> -16114.81
## -16143.9 -> -16129.23
## -16158.8 -> -16130.94
## -16152.55 -> -16131.06
## -16159.89 -> -16143.38
## -16165.75 -> -16144.42
## -16089.28 -> -16089.21
## -16147.35 -> -16129.3
## -16138.55 -> -16120.9
## -16161.67 -> -16140.72
## -16149.24 -> -16144.39
## -16154.49 -> -16149.69
## -16155.56 -> -16136.13
## -16169.81 -> -16138.54
## -16139.43 -> -16116.82
## -16162.75 -> -16123.48
## -16130.17 -> -16102.45
## -16142.48 -> -16123.43
## -16155.19 -> -16135.02
## -16151.43 -> -16134.2
## -16155.16 -> -16147.3
## -16151.82 -> -16144
## -16147.75 -> -16145.55
## -16146.48 -> -16145.06
## -16150.46 -> -16150.31
## -16156.57 -> -16152.39
## optimize edge weights: -15987.6 --> -15957.53
## optimize topology: -16153.54 --> -15957.53
## -15957 -> -15956.47
## -15954.58 -> -15953.67
## -15957.68 -> -15957.38
## -15957.02 -> -15956.62
## -15956.48 -> -15955.61
## -15953.67 -> -15948.85
## -15954.42 -> -15952.52
## -15949.7 -> -15919.08
## -15956.04 -> -15955.12
## -15955.95 -> -15950.67
## -15956.88 -> -15955.33
## -15957.53 -> -15945.51
## -15957.53 -> -15951.81
## optimize edge weights: -15898.05 --> -15895.74
## optimize topology: -15957.53 --> -15895.74
## -15878.49 -> -15870.43
## -15894.21 -> -15893.76
## -15894.21 -> -15888.96
## -15895.11 -> -15893.54
## -15895.18 -> -15894.73
## -15895.62 -> -15895.54
## -15888.51 -> -15863.26
## -15895.21 -> -15894.7
## optimize edge weights: -15855.58 --> -15854.84
## optimize topology: -15895.74 --> -15854.84
## -15835.43 -> -15825.9
## -15853.07 -> -15852.55
## -15851.75 -> -15848
## -15854.19 -> -15853.49
## -15852.67 -> -15851.58
## -15854.86 -> -15854.73
## optimize edge weights: -15825.9 --> -15825.55
## optimize topology: -15854.84 --> -15825.55
## -15825.5 -> -15825.39
## -15823.3 -> -15821.26
## -15825.02 -> -15823.87
## -15822.97 -> -15821.27
## optimize edge weights: -15821.26 --> -15821.08
## optimize topology: -15825.55 --> -15821.08
## -15817.71 -> -15815.65
## -15819.24 -> -15818.46
## optimize edge weights: -15815.65 --> -15815.55
## optimize topology: -15821.08 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 14
## [1] "Ratchet iteration 9 , best pscore so far: -15815.5450292405"
## -16299.11 -> -16292.42
## -16302.24 -> -16297.64
## -16309.35 -> -16301.63
## -16302.9 -> -16298.68
## -16304.72 -> -16300.84
## -16214.54 -> -16192.36
## -16367.46 -> -16243.18
## -16247.43 -> -16166.07
## -16316.23 -> -16236.29
## -16304.99 -> -16295.94
## -16310.49 -> -16299.07
## -16276.56 -> -16268.11
## -16313.93 -> -16292.15
## -16307.25 -> -16301.7
## -16305.31 -> -16300.64
## -16307.47 -> -16305.61
## -16322.64 -> -16304.66
## -16306.44 -> -16296.9
## -16302.65 -> -16294.87
## -16302.48 -> -16290.78
## -16316.25 -> -16294.67
## -16315.22 -> -16294.73
## -16324.59 -> -16298.87
## -16276.34 -> -16272.37
## -16307.97 -> -16297.4
## -16232.61 -> -16232.08
## -16297.96 -> -16276.09
## -16290.01 -> -16288
## -16307.19 -> -16299.11
## optimize edge weights: -16045.38 --> -16020.5
## optimize topology: -16305.62 --> -16020.5
## -16010.65 -> -16006.48
## -16015.46 -> -16013.29
## -16007.69 -> -15991.31
## -16018.63 -> -16017.11
## -16011.17 -> -16005.18
## -16015.65 -> -16011.7
## -16018.74 -> -16015.07
## -16019.14 -> -16014.8
## -15999.26 -> -15893.67
## optimize edge weights: -15836.94 --> -15834.98
## optimize topology: -16020.5 --> -15834.98
## -15825.9 -> -15822.91
## -15830.32 -> -15828.4
## -15832.22 -> -15829.68
## -15834.8 -> -15833.97
## optimize edge weights: -15822.91 --> -15822.8
## optimize topology: -15834.98 --> -15822.8
## -15821.99 -> -15821.06
## -15823.8 -> -15822.57
## -15820.45 -> -15818.27
## optimize edge weights: -15815.98 --> -15815.54
## optimize topology: -15822.8 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 12
## [1] "Ratchet iteration 10 , best pscore so far: -15815.544866991"
## -16595.6 -> -16538.85
## -16607.3 -> -16539.46
## -16570.45 -> -16567.05
## -16507.74 -> -16505.33
## -16573.1 -> -16554.35
## -16518.43 -> -16513.92
## -16569.67 -> -16555.61
## -16572.22 -> -16557.02
## -16572.09 -> -16557.02
## -16577.83 -> -16564.25
## -16573.55 -> -16562.66
## -16525.28 -> -16523.37
## -16574.8 -> -16554.55
## -16578.33 -> -16568.34
## -16580.05 -> -16568.6
## -16558.69 -> -16557.25
## -16560.41 -> -16557.72
## -16581.44 -> -16547.94
## -16591.02 -> -16553.21
## -16625.12 -> -16557.55
## -16635.83 -> -16557.55
## -16349.81 -> -16349.77
## -16586.58 -> -16481.83
## -16572.46 -> -16565.05
## -16570.33 -> -16565.02
## -16536.97 -> -16529.54
## -16577.88 -> -16559.29
## -16560.55 -> -16544.3
## -16573.46 -> -16560.31
## -16574.52 -> -16552.69
## -16571.25 -> -16550.07
## -16518.53 -> -16511.93
## -16592.34 -> -16551.92
## -16563.83 -> -16485.24
## -16411.75 -> -16410.7
## -16570.8 -> -16555.43
## -16564.37 -> -16548.96
## optimize edge weights: -15962.48 --> -15930.35
## optimize topology: -16569.67 --> -15930.35
## -15921.87 -> -15917.49
## -15930.35 -> -15930.35
## -15930.35 -> -15921.28
## -15925.22 -> -15910.6
## -15930.35 -> -15908.77
## -15930.35 -> -15930.35
## -15927.63 -> -15921.78
## -15930.35 -> -15930.34
## -15929.56 -> -15928.74
## -15930.7 -> -15929.67
## optimize edge weights: -15877.79 --> -15874.41
## optimize topology: -15930.35 --> -15874.41
## -15866.46 -> -15862.34
## -15874.03 -> -15873.9
## -15874.5 -> -15874.32
## -15866.22 -> -15862.06
## -15862.66 -> -15842.28
## optimize edge weights: -15816.59 --> -15815.54
## optimize topology: -15874.41 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 14
## [1] "Ratchet iteration 11 , best pscore so far: -15815.5447450837"
## -16124.27 -> -16124.22
## -16138.46 -> -16128.48
## -16134.22 -> -16125.66
## -16128.25 -> -16127.89
## -16130.6 -> -16129.8
## -16125.09 -> -16099.71
## -16125.97 -> -16099.66
## -16092.48 -> -16075.96
## -16130.54 -> -16089.74
## -16117.26 -> -16080.66
## -16101.68 -> -16076.35
## -16140.92 -> -16111.47
## -16139.72 -> -16111.17
## -16112 -> -16107.18
## -16130.35 -> -16120.64
## -16113.28 -> -16106.54
## -16137.1 -> -16124.18
## -16143.49 -> -16111.7
## -16139.69 -> -16109.56
## -16130.01 -> -16125.61
## -16133.04 -> -16127.26
## -16087.71 -> -16068.72
## -16130.02 -> -16109.66
## -16132.98 -> -16117.32
## -16131.09 -> -16117.32
## -16070.22 -> -16067.72
## -16125.91 -> -16108.05
## -16139.58 -> -16123.5
## -16139.52 -> -16123.5
## -16140.57 -> -16115.01
## -16128.74 -> -16112.43
## -16121.96 -> -16121.45
## -16121.78 -> -16121.43
## -16127.15 -> -16127.15
## -16129.5 -> -16128.86
## optimize edge weights: -15926.13 --> -15903.58
## optimize topology: -16130.17 --> -15903.58
## -15901.81 -> -15901.03
## -15903.01 -> -15902.07
## -15902.26 -> -15902.03
## -15903.58 -> -15899.76
## -15903.58 -> -15900.72
## -15902.95 -> -15901.34
## -15902.82 -> -15895.84
## -15899.61 -> -15886.87
## -15902.71 -> -15901.73
## -15901.3 -> -15900.11
## -15900.75 -> -15898.96
## -15901.81 -> -15896.79
## -15902.6 -> -15900.53
## -15903.58 -> -15903.58
## optimize edge weights: -15869.8 --> -15869.51
## optimize topology: -15903.58 --> -15869.51
## -15868.26 -> -15868.08
## -15869.51 -> -15866.39
## -15869.51 -> -15860.64
## -15867.28 -> -15865.59
## -15867.33 -> -15865.29
## optimize edge weights: -15857.61 --> -15857.14
## optimize topology: -15869.51 --> -15857.14
## -15845.05 -> -15840.89
## -15853.26 -> -15848
## -15855.75 -> -15855.55
## optimize edge weights: -15840.89 --> -15839.92
## optimize topology: -15857.14 --> -15839.92
## -15834.03 -> -15826.65
## -15838.24 -> -15837.71
## -15837.26 -> -15831.09
## -15838.61 -> -15838.44
## optimize edge weights: -15825.09 --> -15824.81
## optimize topology: -15839.92 --> -15824.81
## -15824.22 -> -15816.04
## -15826.05 -> -15824.76
## optimize edge weights: -15816.04 --> -15815.84
## optimize topology: -15824.81 --> -15815.84
## -15815.71 -> -15815.61
## optimize edge weights: -15815.61 --> -15815.55
## optimize topology: -15815.84 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 17
## [1] "Ratchet iteration 12 , best pscore so far: -15815.5447450837"
## -16444.6 -> -16355.05
## -16287.76 -> -16283.62
## -16431.62 -> -16409.29
## -16428.53 -> -16405.78
## -16485.02 -> -16393.1
## -16482.56 -> -16393.1
## -16413.76 -> -16400.69
## -16436.74 -> -16411.69
## -16436.16 -> -16427.11
## -16441.76 -> -16426.7
## -16422.19 -> -16421.76
## -16446.54 -> -16431.31
## -16436.72 -> -16427.89
## -16442.12 -> -16429.03
## -16441.33 -> -16426.67
## -16437.53 -> -16424.36
## -16398.85 -> -16398.79
## -16427.01 -> -16416.45
## -16423.09 -> -16422.21
## -16435 -> -16429.56
## -16431.96 -> -16431.85
## -16427.12 -> -16426.46
## -16430.03 -> -16428.82
## -16395.61 -> -16386.19
## -16432.36 -> -16398.89
## -16442.91 -> -16424.33
## -16444.85 -> -16424.24
## -16426.62 -> -16421.23
## -16431.63 -> -16425.16
## -16419.19 -> -16416.78
## -16433.47 -> -16423.35
## -16417.64 -> -16417.63
## -16432.71 -> -16427.43
## -16430.8 -> -16430.75
## -16433.4 -> -16430.75
## -16440.4 -> -16432.2
## -16437.15 -> -16430.86
## -16449.6 -> -16405.2
## -16453.86 -> -16406
## -16508.56 -> -16416.56
## -16514.18 -> -16416.61
## -16178.97 -> -16178.88
## -16416.96 -> -16317.48
## -16414.59 -> -16380.78
## -16440.84 -> -16393.8
## optimize edge weights: -15954.39 --> -15912.04
## optimize topology: -16433.1 --> -15912.04
## -15908.9 -> -15901.41
## -15909.93 -> -15904.85
## -15912.04 -> -15903.05
## -15912.04 -> -15911.94
## -15901.95 -> -15879.95
## -15912.29 -> -15911.87
## -15903.5 -> -15893.37
## -15909.05 -> -15908.39
## -15912.04 -> -15912.04
## optimize edge weights: -15867.18 --> -15865.27
## optimize topology: -15912.04 --> -15865.27
## -15861.77 -> -15859.56
## -15858.98 -> -15848.62
## -15857.16 -> -15843.58
## -15864.86 -> -15864.85
## -15863.48 -> -15857.63
## optimize edge weights: -15843.58 --> -15842.96
## optimize topology: -15865.27 --> -15842.96
## -15841.15 -> -15834.66
## -15842.96 -> -15824.54
## optimize edge weights: -15824.54 --> -15824.21
## optimize topology: -15842.96 --> -15824.21
## -15822.71 -> -15816.33
## -15823.73 -> -15823.46
## optimize edge weights: -15815.58 --> -15815.55
## optimize topology: -15824.21 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 15
## [1] "Ratchet iteration 13 , best pscore so far: -15815.5447450837"
## -16404.01 -> -16365.75
## -16403.54 -> -16365.88
## -16370.55 -> -16369.93
## -16371.09 -> -16370.37
## -16389.83 -> -16368.27
## -16383.86 -> -16367.21
## -16337.81 -> -16329.68
## -16364.09 -> -16356.42
## -16374.18 -> -16354.28
## -16375.6 -> -16354.85
## -16371.95 -> -16363.75
## -16358.84 -> -16350.01
## -16379.82 -> -16361.8
## -16377.53 -> -16361.01
## -16409.6 -> -16367.47
## -16411.41 -> -16367.47
## -16348.9 -> -16340.86
## -16310.84 -> -16301.27
## -16404.54 -> -16368.64
## -16397.77 -> -16367.92
## -16369.98 -> -16363.97
## -16370.32 -> -16363.81
## -16377.55 -> -16359.91
## -16377.94 -> -16359.91
## -16313.37 -> -16313.25
## -16374.12 -> -16353.58
## -16328.43 -> -16325.51
## -16376.66 -> -16354.45
## -16373.5 -> -16346.47
## -16366.01 -> -16339.35
## -16362.06 -> -16358
## -16360.61 -> -16356.64
## -16381.89 -> -16366.16
## -16382.98 -> -16366.28
## -16394.99 -> -16374.32
## -16387.73 -> -16373.52
## -16375.33 -> -16374.05
## -16376.94 -> -16311.31
## -16372.05 -> -16306.3
## -16196.67 -> -16196.16
## -16381.37 -> -16280.68
## optimize edge weights: -15964.04 --> -15932.72
## optimize topology: -16374.34 --> -15932.72
## -15932.72 -> -15931.5
## -15932.72 -> -15932.33
## -15924.66 -> -15919.62
## -15926.19 -> -15920.48
## -15929.33 -> -15920.54
## -15932.21 -> -15931.92
## -15932 -> -15931.36
## -15932.29 -> -15930.41
## -15932.72 -> -15932.72
## -15932.72 -> -15932.33
## -15931.26 -> -15925.01
## optimize edge weights: -15899.33 --> -15898.29
## optimize topology: -15932.72 --> -15898.29
## -15898.29 -> -15896.74
## -15898.29 -> -15897.32
## -15891.99 -> -15885.64
## -15895.81 -> -15893.51
## -15897.92 -> -15897.7
## -15897.91 -> -15897
## -15896.38 -> -15895.49
## -15898.29 -> -15897.86
## -15898.29 -> -15894.9
## optimize edge weights: -15885.64 --> -15885.5
## optimize topology: -15898.29 --> -15885.5
## -15885.5 -> -15882.96
## -15885.5 -> -15884.64
## -15883.33 -> -15881.3
## -15885.03 -> -15884.47
## -15884.98 -> -15884.05
## -15883.34 -> -15882.83
## -15885.5 -> -15885.05
## -15885.5 -> -15881.88
## optimize edge weights: -15875.16 --> -15874.83
## optimize topology: -15885.5 --> -15874.83
## -15869.4 -> -15841.13
## -15874.22 -> -15873.57
## -15864.94 -> -15852.03
## optimize edge weights: -15818.62 --> -15818.3
## optimize topology: -15874.83 --> -15818.3
## -15816.84 -> -15816.54
## -15819.03 -> -15817.83
## optimize edge weights: -15815.7 --> -15815.54
## optimize topology: -15818.3 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 19
## [1] "Ratchet iteration 14 , best pscore so far: -15815.5445479794"
## -16359.11 -> -16358.68
## -16358.03 -> -16357.9
## -16360.16 -> -16358.85
## -16356.53 -> -16352.26
## -16355.02 -> -16350.69
## -16298.16 -> -16298.09
## -16359.12 -> -16338.8
## -16363.2 -> -16343.96
## -16357.28 -> -16340.75
## -16360.89 -> -16350.2
## -16359 -> -16349.52
## -16333.23 -> -16331.81
## -16360.71 -> -16343.26
## -16353.22 -> -16326.32
## -16332.08 -> -16313.08
## -16379.22 -> -16351.06
## -16381.1 -> -16351.89
## -16351.07 -> -16322.42
## -16333.25 -> -16304.34
## -16335.3 -> -16313.1
## -16353.02 -> -16323.02
## -16409.22 -> -16333.78
## -16420.79 -> -16335.3
## -16087.45 -> -16086.75
## -16368.66 -> -16247.54
## -16367.69 -> -16344.86
## -16357.18 -> -16341.23
## -16350.56 -> -16345.55
## -16360.27 -> -16352.48
## -16355.43 -> -16354.87
## -16360.5 -> -16356.81
## -16335.9 -> -16335.86
## -16355.12 -> -16350.66
## -16340.18 -> -16338.41
## -16336.92 -> -16336.74
## -16359.83 -> -16358.84
## optimize edge weights: -15923.96 --> -15892.39
## optimize topology: -16358.95 --> -15892.39
## -15892.41 -> -15892.25
## -15891.54 -> -15890.7
## -15890.98 -> -15881.33
## -15884.44 -> -15878.25
## -15889.33 -> -15886.81
## -15891.46 -> -15890.82
## -15892.39 -> -15872.17
## -15892.39 -> -15892.39
## -15889.05 -> -15886.7
## -15892.39 -> -15890.01
## optimize edge weights: -15856.59 --> -15856.4
## optimize topology: -15892.39 --> -15856.4
## -15847.21 -> -15838.85
## -15853.23 -> -15851.08
## -15855.47 -> -15854.92
## -15856.98 -> -15856.17
## optimize edge weights: -15838.85 --> -15838.04
## optimize topology: -15856.4 --> -15838.04
## -15831.76 -> -15830.58
## -15837.76 -> -15837.67
## -15836.76 -> -15836.65
## -15834.86 -> -15833.46
## -15837.87 -> -15837.71
## -15838.7 -> -15837.91
## optimize edge weights: -15830.58 --> -15830.36
## optimize topology: -15838.04 --> -15830.36
## -15829.99 -> -15829.88
## -15827.32 -> -15825.98
## -15830.27 -> -15830.12
## -15830.11 -> -15829.29
## optimize edge weights: -15825.49 --> -15825.41
## optimize topology: -15830.36 --> -15825.41
## -15821.81 -> -15820.31
## -15825.26 -> -15824.26
## optimize edge weights: -15818.87 --> -15818.3
## optimize topology: -15825.41 --> -15818.3
## -15816.84 -> -15816.54
## -15819.03 -> -15817.83
## optimize edge weights: -15815.7 --> -15815.54
## optimize topology: -15818.3 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 18
## [1] "Ratchet iteration 15 , best pscore so far: -15815.5443150919"
## -16283.14 -> -16283.12
## -16297.69 -> -16292.26
## -16301.98 -> -16298.25
## -16301.65 -> -16298.25
## -16305.01 -> -16296.9
## -16290.72 -> -16289.03
## -16296.25 -> -16293.76
## -16294.18 -> -16291.12
## -16285.3 -> -16283.46
## -16299.71 -> -16295.48
## -16300.9 -> -16293.25
## -16302.68 -> -16293.76
## -16283.99 -> -16271.73
## -16299.61 -> -16281.82
## -16296.89 -> -16279.79
## -16295.38 -> -16277.69
## -16318.13 -> -16292.99
## -16308.37 -> -16289.08
## -16298.71 -> -16296.55
## -16290.49 -> -16275.79
## -16301.33 -> -16287.13
## -16288.96 -> -16277.01
## -16286.45 -> -16276.32
## -16301.84 -> -16285.6
## -16305.26 -> -16284.93
## -16305.62 -> -16293.73
## -16307.74 -> -16294.37
## -16314.39 -> -16277.27
## -16326.33 -> -16278.6
## -16230.78 -> -16210.19
## -16292.45 -> -16265.37
## -16293.55 -> -16266.49
## -16303.5 -> -16271
## -16296.33 -> -16281.19
## -16296.72 -> -16282.89
## -16297.18 -> -16296.95
## -16304.06 -> -16293.8
## -16291.55 -> -16288.29
## -16326.59 -> -16255.56
## -16329.52 -> -16255.92
## -16296.17 -> -16273.84
## -16294.13 -> -16271.38
## -16142.11 -> -16140.97
## -16302.88 -> -16217.11
## optimize edge weights: -16001.95 --> -15959.15
## optimize topology: -16298.29 --> -15959.15
## -15959.15 -> -15958.42
## -15959.15 -> -15950.66
## -15952.39 -> -15951.74
## -15958.45 -> -15957.82
## -15957.97 -> -15957.75
## -15957.28 -> -15955.29
## -15957.26 -> -15956.39
## -15959.15 -> -15954.03
## -15959.15 -> -15952.25
## -15953.47 -> -15937.67
## -15958.34 -> -15957.91
## -15959.15 -> -15943.94
## -15959.15 -> -15949.8
## -15957.96 -> -15956.89
## -15959.02 -> -15958.9
## -15951.16 -> -15946.86
## optimize edge weights: -15913.06 --> -15911.66
## optimize topology: -15959.15 --> -15911.66
## -15911.22 -> -15910.99
## -15902.97 -> -15902.14
## -15909.13 -> -15907.95
## -15911.56 -> -15911.35
## -15898.75 -> -15876.46
## -15911.09 -> -15906.75
## -15910.72 -> -15905.67
## -15911.66 -> -15895.69
## -15911.66 -> -15902.11
## -15910.46 -> -15909.64
## optimize edge weights: -15851.21 --> -15850.01
## optimize topology: -15911.66 --> -15850.01
## -15848.79 -> -15845.77
## -15847 -> -15843.07
## -15849.18 -> -15848.15
## -15850.07 -> -15849.98
## -15849.08 -> -15848.88
## -15838.82 -> -15835.78
## -15851.13 -> -15849.94
## -15849.57 -> -15849.32
## optimize edge weights: -15828.5 --> -15828.22
## optimize topology: -15850.01 --> -15828.22
## -15827.52 -> -15826.81
## -15827.97 -> -15827.85
## optimize edge weights: -15826.81 --> -15826.74
## optimize topology: -15828.22 --> -15826.74
## -15826.4 -> -15825.81
## -15824.58 -> -15823.51
## optimize edge weights: -15823.51 --> -15823.45
## optimize topology: -15826.74 --> -15823.45
## -15821.95 -> -15815.57
## optimize edge weights: -15815.57 --> -15815.54
## optimize topology: -15823.45 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 19
## [1] "Ratchet iteration 16 , best pscore so far: -15815.5440825489"
## -16502.1 -> -16483.25
## -16522.08 -> -16498.04
## -16469.9 -> -16459.43
## -16530.46 -> -16492.15
## -16587.31 -> -16504.7
## -16564.11 -> -16500.97
## -16527.82 -> -16460.4
## -16523.03 -> -16455.83
## -16353.82 -> -16353.62
## -16529.61 -> -16428.76
## -16526.24 -> -16520.96
## -16529.01 -> -16522.63
## -16501.29 -> -16501.29
## -16523.7 -> -16518.47
## -16500.41 -> -16500.38
## -16508.06 -> -16507.1
## -16526.48 -> -16514.7
## -16528.29 -> -16513.14
## -16561.41 -> -16503.21
## -16561.24 -> -16503.12
## -16599.07 -> -16519.64
## -16604.07 -> -16519.64
## -16289.92 -> -16289.62
## -16541.32 -> -16425.02
## -16505.84 -> -16475.02
## -16533.3 -> -16490.58
## -16522.48 -> -16515.52
## -16521.09 -> -16515.1
## -16509.62 -> -16507.56
## -16519.54 -> -16515.49
## -16483.88 -> -16483.72
## -16521.73 -> -16511.6
## -16526.33 -> -16519.8
## -16527.55 -> -16520.12
## -16522.6 -> -16507.98
## -16522.69 -> -16507.98
## -16518.47 -> -16503.15
## -16518.36 -> -16500.38
## -16525.25 -> -16523.58
## -16527.39 -> -16524.35
## -16518.76 -> -16518.03
## -16519.64 -> -16519.08
## optimize edge weights: -15961.37 --> -15929.93
## optimize topology: -16524.6 --> -15929.93
## -15929.93 -> -15928.12
## -15929.93 -> -15929.54
## -15926.54 -> -15925.59
## -15929.47 -> -15929.43
## -15927.51 -> -15926.71
## -15929.93 -> -15929.01
## -15929.93 -> -15928.85
## -15929.93 -> -15922.28
## -15929.93 -> -15927.47
## -15929.45 -> -15929.42
## -15929.93 -> -15927.87
## -15929.93 -> -15927.25
## -15929.93 -> -15902.88
## optimize edge weights: -15893.19 --> -15892.15
## optimize topology: -15929.93 --> -15892.15
## -15890.19 -> -15889.46
## -15884.75 -> -15877.15
## -15892.15 -> -15889.96
## -15892.15 -> -15891.46
## -15888.88 -> -15879.07
## -15891.04 -> -15890.11
## -15889.02 -> -15888.14
## -15892.15 -> -15889.68
## -15892.15 -> -15891.74
## -15891.78 -> -15891.74
## -15873.79 -> -15857.06
## optimize edge weights: -15837.93 --> -15837.09
## optimize topology: -15892.15 --> -15837.09
## -15831.68 -> -15824.49
## -15835.47 -> -15834.98
## -15837.78 -> -15836.54
## -15837.09 -> -15834.45
## -15837.09 -> -15837.09
## -15837.09 -> -15837.09
## -15837.09 -> -15837.09
## -15837.09 -> -15834.47
## -15837.09 -> -15836.52
## optimize edge weights: -15823.72 --> -15823.51
## optimize topology: -15837.09 --> -15823.51
## -15823.51 -> -15821.09
## -15823.51 -> -15819.96
## -15823.51 -> -15822.75
## optimize edge weights: -15817.56 --> -15817.46
## optimize topology: -15823.51 --> -15817.46
## -15817.66 -> -15816.63
## optimize edge weights: -15816.63 --> -15816.48
## optimize topology: -15817.46 --> -15816.48
## -15816.9 -> -15815.64
## optimize edge weights: -15815.64 --> -15815.55
## optimize topology: -15816.48 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 20
## [1] "Ratchet iteration 17 , best pscore so far: -15815.5440825489"
## -16575.35 -> -16532.77
## -16573.23 -> -16532.66
## -16518.03 -> -16517.58
## -16519.8 -> -16518.71
## -16541.23 -> -16540.23
## -16538.54 -> -16511.14
## -16461.08 -> -16456.07
## -16542.37 -> -16518.42
## -16545.11 -> -16519.43
## -16543.04 -> -16532.11
## -16551.24 -> -16533.38
## -16535.59 -> -16535.25
## -16518.52 -> -16492.81
## -16543.61 -> -16518.48
## -16550.51 -> -16534.57
## -16545.17 -> -16531.8
## -16530.15 -> -16526.68
## -16530.35 -> -16521.83
## -16539.93 -> -16530.59
## -16535.85 -> -16528.19
## -16534.87 -> -16526.94
## -16478.81 -> -16477.92
## -16551.42 -> -16518.48
## -16553.01 -> -16512.08
## -16557.63 -> -16513.93
## -16565.47 -> -16451.3
## -16571 -> -16451.39
## -16361.91 -> -16350.19
## -16510.96 -> -16421.34
## -16548.76 -> -16445.36
## -16536.23 -> -16440.66
## -16540.68 -> -16530.01
## -16530.58 -> -16516.43
## -16540.2 -> -16535.58
## -16543.16 -> -16537.11
## -16545.78 -> -16526.97
## -16546.27 -> -16526.97
## -16479.59 -> -16479.59
## -16540.35 -> -16520.61
## optimize edge weights: -16135.99 --> -16081.3
## optimize topology: -16540.43 --> -16081.3
## -16081.3 -> -16080.6
## -16081.3 -> -16079.05
## -16081.3 -> -16073.25
## -16081.3 -> -16073.87
## -16076.51 -> -16076
## -16064.91 -> -16048.71
## -16070.67 -> -16063.88
## -16078.72 -> -16078.14
## -16081.3 -> -16075.79
## -16081.3 -> -16071.91
## -16081.3 -> -16079.54
## -16066.33 -> -16045.81
## -16081.3 -> -15953.97
## -16081.3 -> -16080.7
## optimize edge weights: -15871.93 --> -15869.31
## optimize topology: -16081.3 --> -15869.31
## -15869.31 -> -15868.45
## -15869.31 -> -15866.79
## -15868.82 -> -15868.12
## -15869.68 -> -15868.73
## -15863.93 -> -15863.42
## -15868.58 -> -15868.54
## -15858.07 -> -15852.11
## -15868.68 -> -15868.56
## -15869.31 -> -15846.85
## -15869.31 -> -15868.56
## optimize edge weights: -15818.87 --> -15817.72
## optimize topology: -15869.31 --> -15817.72
## -15816.05 -> -15815.59
## optimize edge weights: -15815.59 --> -15815.54
## optimize topology: -15817.72 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 17
## [1] "Ratchet iteration 18 , best pscore so far: -15815.5439318904"
## -16185.72 -> -16182.55
## -16192.24 -> -16182.79
## -16182.5 -> -16181.44
## -16172.94 -> -16171.66
## -16175.5 -> -16169.3
## -16182.44 -> -16174.73
## -16191.85 -> -16177.63
## -16202.6 -> -16182.18
## -16176.25 -> -16168.1
## -16188.55 -> -16180.45
## -16188.23 -> -16182.5
## -16188.43 -> -16182.5
## -16182.62 -> -16182.42
## -16189.58 -> -16156.94
## -16179.46 -> -16156.72
## -16188.79 -> -16141.1
## -16102.75 -> -16094.64
## -16197.79 -> -16142.38
## -16190.02 -> -16142.47
## -16143.65 -> -16139.09
## -16176.73 -> -16164.75
## -16167.99 -> -16156.03
## -16175.81 -> -16165.25
## -16189.59 -> -16177.67
## -16177.19 -> -16172.22
## -16179.01 -> -16168.92
## -16177.33 -> -16166.61
## -16170.25 -> -16163.53
## -16178.27 -> -16171.02
## -16182.88 -> -16180.89
## -16197.3 -> -16180.1
## -16201.99 -> -16180.76
## optimize edge weights: -16034.82 --> -15995.88
## optimize topology: -16183.11 --> -15995.88
## -15995.88 -> -15994.91
## -15995.88 -> -15995.77
## -15992.5 -> -15991.59
## -15995.88 -> -15991.83
## -15990.55 -> -15978.19
## -15989.71 -> -15953.26
## -15989.43 -> -15978.98
## -15989.84 -> -15979.09
## -15995.88 -> -15994.52
## -15995.88 -> -15993.61
## -15995.88 -> -15982.91
## -15995.88 -> -15994.52
## -15994.83 -> -15992.65
## optimize edge weights: -15911 --> -15910.52
## optimize topology: -15995.88 --> -15910.52
## -15910.52 -> -15909.04
## -15910.52 -> -15910.24
## -15910.52 -> -15907.98
## -15910.52 -> -15909.83
## -15909.29 -> -15908.9
## -15905.34 -> -15898.02
## -15910.68 -> -15910.06
## -15910.52 -> -15910.52
## -15910.52 -> -15909.14
## -15910.52 -> -15903.68
## -15910.52 -> -15908.89
## optimize edge weights: -15888.98 --> -15888.65
## optimize topology: -15910.52 --> -15888.65
## -15888.65 -> -15886.59
## -15888.65 -> -15888.32
## -15886.88 -> -15886.16
## -15888.65 -> -15888.12
## -15888.65 -> -15870.78
## -15888.65 -> -15843.47
## -15888.65 -> -15888.65
## optimize edge weights: -15825.43 --> -15824.53
## optimize topology: -15888.65 --> -15824.53
## -15824.53 -> -15823.09
## -15824.53 -> -15824.07
## -15822.73 -> -15822.27
## -15823.89 -> -15823.82
## -15822.91 -> -15822.36
## -15823.2 -> -15820.97
## optimize edge weights: -15817.74 --> -15817.51
## optimize topology: -15824.53 --> -15817.51
## -15815.93 -> -15815.6
## optimize edge weights: -15815.6 --> -15815.54
## optimize topology: -15817.51 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 18
## [1] "Ratchet iteration 19 , best pscore so far: -15815.543758178"
## -16516.54 -> -16510.76
## -16550.3 -> -16532.83
## -16408.97 -> -16187.19
## -16539.56 -> -16450.93
## -16640.32 -> -16493.92
## -16670.57 -> -16495.39
## -16529.59 -> -16526.11
## -16538.68 -> -16530.68
## -16521.6 -> -16516.44
## -16544.43 -> -16529.63
## -16544.93 -> -16534.99
## -16513.65 -> -16504.56
## -16539.5 -> -16526.02
## -16534.69 -> -16524.55
## -16531.8 -> -16524.19
## -16540.87 -> -16518
## -16543.96 -> -16519.35
## -16495.16 -> -16488.82
## -16551.74 -> -16531.25
## -16533.38 -> -16529.17
## -16533.25 -> -16528.67
## -16475.73 -> -16475.64
## -16536.5 -> -16516.15
## -16546.99 -> -16521.53
## -16540.23 -> -16518.31
## -16536.89 -> -16535.58
## -16534.36 -> -16533.31
## -16556.64 -> -16496.67
## -16530.64 -> -16483.2
## -16522.77 -> -16496.22
## -16531.15 -> -16501.95
## optimize edge weights: -16005.19 --> -15988.66
## optimize topology: -16536.5 --> -15988.66
## -15986.48 -> -15983.71
## -15979.2 -> -15970.58
## -15983.08 -> -15980.76
## -15988.66 -> -15973.86
## -15988.66 -> -15979.33
## -15988.66 -> -15987.8
## -15988.66 -> -15984.26
## -15987.63 -> -15986.89
## optimize edge weights: -15954.01 --> -15953.47
## optimize topology: -15988.66 --> -15953.47
## -15953.23 -> -15951.74
## -15951.09 -> -15948.31
## -15947.21 -> -15944.43
## -15953.47 -> -15952.79
## -15953.47 -> -15947.83
## optimize edge weights: -15939.06 --> -15938.73
## optimize topology: -15953.47 --> -15938.73
## -15938.73 -> -15938
## -15938.73 -> -15933.11
## optimize edge weights: -15933.11 --> -15933
## optimize topology: -15938.73 --> -15933
## -15932.7 -> -15932.51
## optimize edge weights: -15932.51 --> -15932.38
## optimize topology: -15933 --> -15932.38
## optimize edge weights: -15932.38 --> -15932.38
## optimize topology: -15932.38 --> -15932.38
## NNI moves: 13
## [1] "Ratchet iteration 20 , best pscore so far: -15815.543758178"
## -16367.24 -> -16332.74
## -16314.8 -> -16286.05
## -16360.81 -> -16360.69
## -16365.74 -> -16363.08
## -16383.13 -> -16335.38
## -16345.07 -> -16313.44
## -16385.33 -> -16356.09
## -16372.64 -> -16350.04
## -16374 -> -16363.99
## -16373.16 -> -16364.69
## -16368.22 -> -16362.53
## -16358.67 -> -16351.09
## -16355.54 -> -16348.66
## -16379.18 -> -16274
## -16385.25 -> -16274
## -16022.75 -> -16016.69
## -16352.93 -> -16218.74
## -16365.35 -> -16364.38
## -16361.17 -> -16361.14
## -16366.6 -> -16360.58
## -16372.37 -> -16363
## -16340.99 -> -16340.98
## -16361.06 -> -16357.23
## optimize edge weights: -15902 --> -15892.84
## optimize topology: -16365.27 --> -15892.84
## -15881.85 -> -15861.93
## -15892.84 -> -15892.1
## -15892.84 -> -15884.43
## -15891.64 -> -15891.37
## optimize edge weights: -15852.12 --> -15851.05
## optimize topology: -15892.84 --> -15851.05
## -15846.86 -> -15846.4
## -15849.93 -> -15849.11
## -15846.71 -> -15843.27
## -15846.68 -> -15844.48
## -15844.56 -> -15831.36
## optimize edge weights: -15823.59 --> -15823.37
## optimize topology: -15851.05 --> -15823.37
## -15822.44 -> -15822.28
## -15818.22 -> -15817.74
## -15822.17 -> -15821.5
## optimize edge weights: -15817.74 --> -15817.51
## optimize topology: -15823.37 --> -15817.51
## -15815.93 -> -15815.6
## optimize edge weights: -15815.6 --> -15815.54
## optimize topology: -15817.51 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 12
## [1] "Ratchet iteration 21 , best pscore so far: -15815.543758178"
## -16192.61 -> -16185.01
## -16205.77 -> -16195.7
## -16156.19 -> -16154.21
## -16213.33 -> -16189.65
## -16221.67 -> -16201.26
## -16219.23 -> -16198.62
## -16236.32 -> -16210.3
## -16220.47 -> -16202.19
## -16211.5 -> -16201.76
## -16223.27 -> -16209.55
## -16183.84 -> -16173.34
## -16216.87 -> -16200.09
## -16212.75 -> -16202.53
## -16193.33 -> -16192.04
## -16204.09 -> -16173.3
## -16214.86 -> -16181.03
## -16215.01 -> -16208.79
## -16208.09 -> -16206.46
## -16222.27 -> -16212.55
## -16216.03 -> -16176.99
## -16180.64 -> -16163.8
## -16212.55 -> -16181.3
## -16188.69 -> -16170.92
## -16211.07 -> -16202.46
## -16214.01 -> -16203.77
## -16213.9 -> -16211.44
## -16211.09 -> -16210.21
## -16213.75 -> -16210.86
## -16218.83 -> -16211.94
## -16201.98 -> -16199.62
## -16213.59 -> -16208.54
## optimize edge weights: -16053.79 --> -16010.92
## optimize topology: -16213.02 --> -16010.92
## -16009.47 -> -16009.24
## -16009.49 -> -16009.41
## -16009.86 -> -16009.55
## -16002.94 -> -15986.41
## -15999.37 -> -15993.43
## -16009.46 -> -16008.84
## -16006.33 -> -15997.53
## -16010.92 -> -16004.49
## -16010.92 -> -16010.92
## -16010.92 -> -16010.67
## -16010.92 -> -16010.92
## -16010.92 -> -15983.85
## -16010.92 -> -16009.52
## -16000.85 -> -15997.74
## optimize edge weights: -15915.42 --> -15912.35
## optimize topology: -16010.92 --> -15912.35
## -15908.89 -> -15907.07
## -15899.68 -> -15872.06
## -15912.29 -> -15912.21
## -15907.34 -> -15882.52
## -15912.35 -> -15901.46
## -15912.35 -> -15910.39
## -15912.35 -> -15912.33
## -15911.47 -> -15910.94
## optimize edge weights: -15842.02 --> -15841.26
## optimize topology: -15912.35 --> -15841.26
## -15838.59 -> -15836.03
## -15841.26 -> -15830.71
## -15841.26 -> -15839.15
## -15841.26 -> -15840.23
## -15840.94 -> -15840.84
## optimize edge weights: -15825.55 --> -15825.09
## optimize topology: -15841.26 --> -15825.09
## -15825.09 -> -15815.55
## -15825.09 -> -15825.09
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15825.09 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 16
## [1] "Ratchet iteration 22 , best pscore so far: -15815.543758178"
## -16445.87 -> -16367.28
## -16435.6 -> -16367.12
## -16236 -> -16226.96
## -16372.86 -> -16294.08
## -16383.65 -> -16371.2
## -16377.61 -> -16364.53
## -16394.17 -> -16367.11
## -16392.34 -> -16367.11
## -16393.27 -> -16373.91
## -16394.45 -> -16373.91
## -16333.26 -> -16328.72
## -16385.08 -> -16370.81
## -16340.55 -> -16340.09
## -16385.35 -> -16371.88
## -16389.7 -> -16374.73
## -16372.93 -> -16368.81
## -16386.36 -> -16382.57
## -16385.28 -> -16381.75
## -16388.91 -> -16349.97
## -16394.57 -> -16353.78
## -16383.49 -> -16342.35
## -16388.9 -> -16343.88
## -16350.09 -> -16294.52
## -16391.77 -> -16340.61
## -16318.81 -> -16311.6
## -16376.96 -> -16350.84
## -16399.16 -> -16384.08
## -16385.39 -> -16384.21
## -16377.61 -> -16377.5
## -16380.93 -> -16347.29
## -16397.02 -> -16352.61
## -16381 -> -16380.79
## -16383.17 -> -16381.92
## -16383.75 -> -16383.58
## -16385.09 -> -16384.54
## optimize edge weights: -16019.45 --> -15961.64
## optimize topology: -16385.08 --> -15961.64
## -15961.64 -> -15961.64
## -15961.64 -> -15952.17
## -15955.28 -> -15950.31
## -15956.53 -> -15951.8
## -15957.91 -> -15954.87
## -15949.97 -> -15923.4
## -15961.71 -> -15961.63
## -15961.64 -> -15958.5
## optimize edge weights: -15903.07 --> -15902.45
## optimize topology: -15961.64 --> -15902.45
## -15901.47 -> -15901.16
## -15900.02 -> -15898.89
## -15897.96 -> -15892.35
## -15897.43 -> -15891.97
## -15902.09 -> -15898.25
## optimize edge weights: -15891.97 --> -15891.06
## optimize topology: -15902.45 --> -15891.06
## -15888.7 -> -15886.92
## -15891.32 -> -15891.02
## -15881.33 -> -15876.51
## -15891.06 -> -15871.74
## -15891.06 -> -15891.02
## optimize edge weights: -15871.74 --> -15871.54
## optimize topology: -15891.06 --> -15871.54
## -15869.45 -> -15867.84
## -15863.22 -> -15859.9
## optimize edge weights: -15859.9 --> -15859.54
## optimize topology: -15871.54 --> -15859.54
## -15852.01 -> -15833.2
## -15859.56 -> -15859.43
## optimize edge weights: -15833.2 --> -15832.12
## optimize topology: -15859.54 --> -15832.12
## -15830.55 -> -15829.42
## -15826.02 -> -15816.37
## optimize edge weights: -15816.37 --> -15815.84
## optimize topology: -15832.12 --> -15815.84
## -15815.71 -> -15815.61
## optimize edge weights: -15815.61 --> -15815.54
## optimize topology: -15815.84 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 15
## [1] "Ratchet iteration 23 , best pscore so far: -15815.543758178"
## -16505.63 -> -16501.17
## -16506.62 -> -16500.33
## -16507.72 -> -16504.49
## -16516.64 -> -16506.7
## -16503.13 -> -16489.44
## -16504.81 -> -16493.74
## -16513.44 -> -16496.39
## -16499.56 -> -16491.38
## -16490.64 -> -16490.58
## -16513.4 -> -16506.47
## -16508.02 -> -16495
## -16518.82 -> -16502.43
## -16509.36 -> -16506.57
## -16505.04 -> -16504.96
## -16508.34 -> -16507.64
## -16476.15 -> -16471.3
## -16511.84 -> -16492.75
## -16525.09 -> -16506.14
## -16548.3 -> -16507.7
## -16504.65 -> -16504.01
## -16510.85 -> -16507.33
## -16493.14 -> -16484.52
## -16510.91 -> -16493.5
## -16505.12 -> -16492.91
## -16515.12 -> -16496.85
## -16503.25 -> -16500.58
## -16507.58 -> -16505.45
## -16529.23 -> -16494.6
## -16531.42 -> -16495.64
## -16441.04 -> -16439.67
## -16497.44 -> -16475.43
## -16492.85 -> -16434.59
## -16512.33 -> -16440.66
## -16428.31 -> -16318.27
## -16505.04 -> -16414.95
## -16534.61 -> -16437.57
## -16479.95 -> -16425.5
## -16524.96 -> -16504.14
## -16520.85 -> -16503.24
## -16475.01 -> -16461.14
## -16508.45 -> -16484.97
## -16508.86 -> -16486.77
## -16481.76 -> -16476.21
## optimize edge weights: -16181.71 --> -16122.65
## optimize topology: -16508.51 --> -16122.65
## -16122.13 -> -16121.48
## -16121.19 -> -16120.08
## -16119.33 -> -16101.52
## -16122 -> -16121.4
## -16122.56 -> -16121.1
## -16121.09 -> -16119.94
## -16120.63 -> -16119.7
## -16122.21 -> -16121.15
## -16122.55 -> -16121.63
## -16121.2 -> -16120.5
## -16110.64 -> -16083.58
## -16122.22 -> -16121.96
## -16122.65 -> -16120.02
## -16122.65 -> -16109.89
## -16122.65 -> -15988.3
## -16122.65 -> -16122.19
## -16122.17 -> -16122.15
## -16122.64 -> -16082.38
## optimize edge weights: -15879.19 --> -15878.55
## optimize topology: -16122.65 --> -15878.55
## -15877.95 -> -15877.24
## -15877.29 -> -15876.42
## -15878.4 -> -15876.86
## -15879.01 -> -15878.27
## -15876.71 -> -15876.06
## -15877.91 -> -15876.24
## -15877.14 -> -15871.82
## optimize edge weights: -15869.8 --> -15869.51
## optimize topology: -15878.55 --> -15869.51
## -15869.66 -> -15868.11
## -15870.32 -> -15869.36
## -15866.11 -> -15865.38
## optimize edge weights: -15865.38 --> -15865.1
## optimize topology: -15869.51 --> -15865.1
## -15865.11 -> -15863.63
## optimize edge weights: -15863.63 --> -15863.56
## optimize topology: -15865.1 --> -15863.56
## -15862.09 -> -15855.46
## optimize edge weights: -15855.46 --> -15855.44
## optimize topology: -15863.56 --> -15855.44
## optimize edge weights: -15855.44 --> -15855.44
## optimize topology: -15855.44 --> -15855.44
## NNI moves: 17
## [1] "Ratchet iteration 24 , best pscore so far: -15815.543758178"
## -16012.54 -> -16001.58
## -15985.73 -> -15984.85
## -16017.32 -> -16004.67
## -16015.16 -> -16001.43
## -15968.31 -> -15968.29
## -16007.85 -> -15997.65
## -16009.28 -> -16007.13
## -16012.58 -> -16009.91
## -15999.52 -> -15999.5
## -16003.67 -> -16003.05
## -16001.15 -> -16000.7
## -16007.92 -> -16004.86
## -16012.14 -> -16008.9
## -16009.73 -> -16005.47
## -16013.34 -> -15992.13
## -16016.3 -> -15993.82
## -15943.45 -> -15933.52
## -16025.38 -> -15989.62
## optimize edge weights: -15878.78 --> -15864.25
## optimize topology: -16011.46 --> -15864.25
## -15864.25 -> -15847.17
## -15862.46 -> -15861.95
## -15861.16 -> -15857.53
## -15862.28 -> -15860.49
## -15864.25 -> -15847.6
## -15864.25 -> -15855.28
## optimize edge weights: -15830.49 --> -15830.25
## optimize topology: -15864.25 --> -15830.25
## -15823.18 -> -15822.51
## -15828.97 -> -15828.35
## -15827.14 -> -15823.36
## -15828.31 -> -15826.57
## -15830.11 -> -15830.01
## -15829.71 -> -15829.25
## optimize edge weights: -15822.51 --> -15822.19
## optimize topology: -15830.25 --> -15822.19
## -15818.85 -> -15816.79
## -15820.35 -> -15819.57
## -15821.99 -> -15821.86
## -15821.6 -> -15821.08
## optimize edge weights: -15815.65 --> -15815.54
## optimize topology: -15822.19 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 8
## [1] "Ratchet iteration 25 , best pscore so far: -15815.543758178"
## -16036.79 -> -16036.09
## -16047.72 -> -16037.92
## -16026.45 -> -16025.69
## -16042.54 -> -16033.14
## -16044.66 -> -16033.7
## -16047.33 -> -16036.29
## -16048.32 -> -16036.29
## -16030.66 -> -16022.77
## -16026.25 -> -16019.34
## -16040.22 -> -16019.73
## -16028.72 -> -16014.67
## -16044.04 -> -16038.28
## -16042.05 -> -16038.23
## -15996.69 -> -15995.65
## -16036.86 -> -16021.62
## -16044.96 -> -16030.34
## -16064.17 -> -16033.21
## -16029.14 -> -16025.21
## -16044.28 -> -16034.85
## -16040.95 -> -16031.38
## -16047.98 -> -16035.42
## -16035.15 -> -16035.04
## -16044.06 -> -16031.13
## -16033.75 -> -16025.79
## -16019.28 -> -15997.1
## -16038.98 -> -16025.53
## -16034.27 -> -16029.41
## -16034.72 -> -16030.47
## -16041.52 -> -16036.65
## -16045.12 -> -16037.78
## -16015.28 -> -16015.21
## -16038.59 -> -16033.77
## -16034.1 -> -16032.31
## -16031.34 -> -16031.17
## optimize edge weights: -15916.45 --> -15886.2
## optimize topology: -16038.52 --> -15886.2
## -15888.29 -> -15885.4
## -15885.22 -> -15882.9
## -15884.01 -> -15883.3
## -15886.2 -> -15884.33
## -15886.2 -> -15885.11
## -15886.31 -> -15885.47
## -15884.2 -> -15878.81
## -15884.25 -> -15884.02
## -15884.15 -> -15883.2
## -15883.56 -> -15880.98
## -15886.2 -> -15877.71
## -15886.2 -> -15885.86
## optimize edge weights: -15873.65 --> -15873.21
## optimize topology: -15886.2 --> -15873.21
## -15869.65 -> -15864.62
## -15870.72 -> -15865.1
## -15872.83 -> -15872.7
## -15873.14 -> -15868.88
## -15873.17 -> -15871.75
## -15872.26 -> -15872.16
## optimize edge weights: -15862.04 --> -15860.64
## optimize topology: -15873.21 --> -15860.64
## -15857.41 -> -15856.68
## -15851.92 -> -15849.37
## -15851.65 -> -15843.34
## -15858.08 -> -15852.18
## -15860.16 -> -15860.02
## -15855.66 -> -15855.18
## optimize edge weights: -15837.96 --> -15837.29
## optimize topology: -15860.64 --> -15837.29
## -15835.15 -> -15834.77
## -15828.26 -> -15826.03
## -15835 -> -15828.98
## -15835.7 -> -15835.31
## optimize edge weights: -15823.9 --> -15823.74
## optimize topology: -15837.29 --> -15823.74
## -15823.62 -> -15823.5
## -15822.22 -> -15815.87
## optimize edge weights: -15815.87 --> -15815.84
## optimize topology: -15823.74 --> -15815.84
## -15815.71 -> -15815.61
## optimize edge weights: -15815.61 --> -15815.54
## optimize topology: -15815.84 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 15
## [1] "Ratchet iteration 26 , best pscore so far: -15815.543758178"
## -16365.95 -> -16365.87
## -16383.72 -> -16379.94
## -16385.46 -> -16381.53
## -16381.4 -> -16378.18
## -16383.51 -> -16382.62
## -16371.59 -> -16371.57
## -16376.42 -> -16375.54
## -16383.05 -> -16363.69
## -16385.62 -> -16365.27
## -16359.33 -> -16353.97
## -16382.54 -> -16367.09
## -16407.29 -> -16349.75
## -16425 -> -16351.07
## -16222.95 -> -16221.46
## -16368.22 -> -16295.77
## -16383.87 -> -16370.2
## -16378.07 -> -16364.11
## -16392.58 -> -16382.14
## -16370.7 -> -16351.42
## -16387.58 -> -16357.31
## -16397.8 -> -16360.73
## -16392.39 -> -16356.54
## -16265.3 -> -16246.12
## -16405.09 -> -16340.76
## -16400.66 -> -16374.84
## -16405.52 -> -16375.58
## -16382.39 -> -16368.52
## -16381.79 -> -16370.25
## -16313.26 -> -16311.19
## -16382.31 -> -16355.58
## -16413.95 -> -16382.02
## -16409.97 -> -16381.75
## -16354.65 -> -16342.22
## -16372.98 -> -16356.82
## -16397.65 -> -16368.63
## -16399.97 -> -16370.01
## -16329.78 -> -16316.3
## -16383.35 -> -16368.48
## optimize edge weights: -15992.3 --> -15950.65
## optimize topology: -16383.35 --> -15950.65
## -15950.64 -> -15950.6
## -15938.16 -> -15918.38
## -15950.66 -> -15950.57
## -15949.27 -> -15948.44
## -15947.81 -> -15943.13
## -15939.5 -> -15914.69
## -15950.01 -> -15949.97
## -15940.74 -> -15933.63
## -15947.05 -> -15945.95
## -15945.08 -> -15944.36
## -15950.65 -> -15948.55
## -15950.65 -> -15950.37
## optimize edge weights: -15878.87 --> -15876.49
## optimize topology: -15950.65 --> -15876.49
## -15873.8 -> -15871.82
## -15876.08 -> -15875.78
## -15870.53 -> -15869.9
## -15865.17 -> -15859.24
## optimize edge weights: -15847.63 --> -15847.43
## optimize topology: -15876.49 --> -15847.43
## -15847.79 -> -15847.39
## optimize edge weights: -15847.39 --> -15847.25
## optimize topology: -15847.43 --> -15847.25
## -15841.03 -> -15835.88
## optimize edge weights: -15835.88 --> -15835.72
## optimize topology: -15847.25 --> -15835.72
## -15828.83 -> -15819
## optimize edge weights: -15819 --> -15818.51
## optimize topology: -15835.72 --> -15818.51
## -15817.93 -> -15815.74
## optimize edge weights: -15815.74 --> -15815.54
## optimize topology: -15818.51 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 16
## [1] "Ratchet iteration 27 , best pscore so far: -15815.543758178"
## -16577.96 -> -16516.64
## -16579.57 -> -16514.97
## -16555.97 -> -16536.62
## -16553.78 -> -16533.84
## -16399.52 -> -16399.35
## -16571.87 -> -16479.53
## -16563.84 -> -16557.29
## -16551.54 -> -16549.69
## -16559.9 -> -16550.4
## -16563.58 -> -16551.09
## -16555.94 -> -16543.12
## -16502.14 -> -16498.04
## -16548.3 -> -16545.67
## -16558.92 -> -16553.04
## -16547.27 -> -16542.56
## -16523.68 -> -16522.93
## -16558.86 -> -16557.72
## -16567.33 -> -16550.69
## -16556.05 -> -16543.92
## -16587.16 -> -16556.97
## -16591.24 -> -16555.65
## -16560.46 -> -16524.24
## -16521.99 -> -16493.53
## -16598.52 -> -16534.25
## -16605.84 -> -16534.25
## -16490.32 -> -16388.67
## -16544.84 -> -16449.67
## -16574.06 -> -16484.88
## -16529.51 -> -16474.85
## -16551.22 -> -16544.23
## -16552 -> -16540.27
## -16561.85 -> -16551.49
## -16557.69 -> -16549.77
## -16542.06 -> -16541.67
## -16568.34 -> -16553.18
## -16562.36 -> -16552.18
## -16565.49 -> -16553.15
## -16566.36 -> -16550.96
## -16564.45 -> -16550.94
## -16522.09 -> -16522.03
## -16552.33 -> -16541.43
## -16550.92 -> -16550.06
## -16554.44 -> -16552.17
## optimize edge weights: -16116.91 --> -16057.74
## optimize topology: -16558.45 --> -16057.74
## -16054.53 -> -16052.25
## -16057.32 -> -16057.23
## -16057.74 -> -16056.15
## -16057.74 -> -16056.04
## -16057.74 -> -16054.76
## -16056 -> -16054.2
## -16049.11 -> -16031.14
## -16056.74 -> -16056.52
## -16057.8 -> -16057.14
## -16056.81 -> -16056.35
## -16052 -> -16027.27
## -16057.74 -> -15938.94
## -16057.74 -> -16057.06
## -16055.79 -> -16049.37
## -16057.74 -> -16053.68
## -16057.74 -> -16055.8
## optimize edge weights: -15904.74 --> -15902.59
## optimize topology: -16057.74 --> -15902.59
## -15901.92 -> -15901.77
## -15901.95 -> -15901.84
## -15902.59 -> -15900.73
## -15902.59 -> -15901.09
## -15896.19 -> -15895.56
## -15902.21 -> -15900.77
## -15901.2 -> -15899.42
## -15898.13 -> -15867.04
## -15903.29 -> -15902.49
## -15899.97 -> -15893.27
## optimize edge weights: -15855.16 --> -15854.72
## optimize topology: -15902.59 --> -15854.72
## -15847.44 -> -15841.98
## -15854.62 -> -15854.45
## -15847.16 -> -15846.5
## -15853.03 -> -15851.19
## -15853.63 -> -15852.69
## -15853.96 -> -15853.22
## -15854.19 -> -15853.22
## optimize edge weights: -15837.82 --> -15837.31
## optimize topology: -15854.72 --> -15837.31
## -15830.24 -> -15829.69
## -15835.38 -> -15834.9
## -15827.21 -> -15822.4
## optimize edge weights: -15822.4 --> -15822.27
## optimize topology: -15837.31 --> -15822.27
## -15816.2 -> -15815.67
## -15821.24 -> -15820.59
## optimize edge weights: -15815.67 --> -15815.54
## optimize topology: -15822.27 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 18
## [1] "Ratchet iteration 28 , best pscore so far: -15815.543758178"
## -16531.54 -> -16522.48
## -16542.89 -> -16533.42
## -16552.35 -> -16545.04
## -16548.1 -> -16545.03
## -16545.93 -> -16528.18
## -16547.88 -> -16528.81
## -16483.06 -> -16473.95
## -16544.31 -> -16513.31
## -16564.18 -> -16544.5
## -16567.58 -> -16545.37
## -16561.59 -> -16516.59
## -16554.77 -> -16515.41
## -16550.68 -> -16538.58
## -16544.55 -> -16536.88
## -16565.34 -> -16529.98
## -16539.45 -> -16516.39
## -16548.27 -> -16481.15
## -16543.78 -> -16476.52
## -16377.77 -> -16377.73
## -16557.43 -> -16452.36
## -16543.16 -> -16538.81
## -16543.24 -> -16538.67
## -16555.34 -> -16531.78
## -16556.13 -> -16531.78
## -16484.42 -> -16484.34
## -16545.98 -> -16525.41
## -16515.58 -> -16515.03
## -16547.57 -> -16534.66
## -16553.13 -> -16541.29
## -16554.33 -> -16541.29
## -16552.82 -> -16542.31
## -16550.19 -> -16541.81
## -16559.69 -> -16543.39
## -16576.68 -> -16511.09
## -16581.26 -> -16511.1
## -16474.71 -> -16381.21
## -16534.04 -> -16449.23
## -16562.84 -> -16474.94
## -16536.32 -> -16469.41
## optimize edge weights: -16046.11 --> -16013.2
## optimize topology: -16545.98 --> -16013.2
## -16007.6 -> -16001.72
## -16012.15 -> -16011.74
## -16011.04 -> -16006.93
## -16012.41 -> -16010.42
## -16012.49 -> -16012.28
## -16011.14 -> -16008.94
## -16013.2 -> -15885.65
## -16013.2 -> -16012.61
## -16012.16 -> -16011.34
## optimize edge weights: -15871.04 --> -15869.14
## optimize topology: -16013.2 --> -15869.14
## -15867.88 -> -15867.48
## -15858.57 -> -15845.46
## -15867.17 -> -15864.99
## -15869.14 -> -15866.81
## -15869.14 -> -15845.3
## optimize edge weights: -15821.32 --> -15821
## optimize topology: -15869.14 --> -15821
## -15819.41 -> -15819.08
## -15819.95 -> -15819.69
## -15819.02 -> -15817.57
## optimize edge weights: -15815.64 --> -15815.54
## optimize topology: -15821 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 14
## [1] "Ratchet iteration 29 , best pscore so far: -15815.543758178"
## -16556.35 -> -16536.95
## -16568.84 -> -16541.82
## -16550.79 -> -16507.69
## -16566.27 -> -16505.88
## -16554.62 -> -16522.03
## -16554.06 -> -16520.65
## -16527.51 -> -16459.14
## -16551.43 -> -16487.31
## -16589.41 -> -16509.03
## -16553.58 -> -16495.43
## -16557.58 -> -16546.18
## -16565.95 -> -16548.21
## -16558.85 -> -16546.6
## -16569.21 -> -16553.15
## -16491.45 -> -16490.9
## -16568.32 -> -16536.36
## -16530.35 -> -16527.17
## -16553.78 -> -16545.57
## -16555.01 -> -16555.01
## -16553.65 -> -16550.29
## -16560.96 -> -16552.5
## -16572.41 -> -16536.36
## -16582.78 -> -16540.86
## -16612.88 -> -16543.17
## -16623.44 -> -16543.17
## -16336.02 -> -16336
## -16572.91 -> -16468.17
## -16550.17 -> -16545.02
## -16554.34 -> -16549.94
## -16514.93 -> -16514.66
## -16549.55 -> -16541.15
## -16558.48 -> -16549.81
## -16559.73 -> -16550.11
## -16494.67 -> -16485.42
## -16557.72 -> -16527.36
## -16551.81 -> -16542.07
## -16557.16 -> -16545.26
## -16557.59 -> -16555.96
## -16559.39 -> -16556.51
## -16566.17 -> -16524.79
## -16519.14 -> -16492.85
## -16572.6 -> -16541.56
## -16563.24 -> -16539.99
## optimize edge weights: -16039.46 --> -16004.44
## optimize topology: -16556.83 --> -16004.44
## -15996.08 -> -15921.29
## -16003.07 -> -16001.45
## -16000.61 -> -15988.94
## -16003.55 -> -16002.49
## -16003.08 -> -16002.92
## -16004.44 -> -15992.89
## -16004.44 -> -16004.44
## -16004.44 -> -16004.44
## -16000.89 -> -15992.9
## -16003.45 -> -16003.17
## -16000.96 -> -15995.95
## optimize edge weights: -15881.77 --> -15880.46
## optimize topology: -16004.44 --> -15880.46
## -15878.19 -> -15875.23
## -15870.43 -> -15861.01
## -15874.11 -> -15867.83
## -15880.46 -> -15878.82
## -15880.46 -> -15879.63
## optimize edge weights: -15859.91 --> -15859.5
## optimize topology: -15880.46 --> -15859.5
## -15853.25 -> -15838.54
## -15856.89 -> -15850.27
## -15854.03 -> -15850.16
## optimize edge weights: -15838.54 --> -15838.4
## optimize topology: -15859.5 --> -15838.4
## -15837.28 -> -15837.09
## -15836.55 -> -15828.85
## -15838.38 -> -15838.35
## -15833.2 -> -15829.53
## optimize edge weights: -15828.85 --> -15828.78
## optimize topology: -15838.4 --> -15828.78
## -15827.71 -> -15827.54
## -15822.65 -> -15818.26
## optimize edge weights: -15818.26 --> -15818.1
## optimize topology: -15828.78 --> -15818.1
## -15816.32 -> -15815.93
## -15817.72 -> -15817.57
## optimize edge weights: -15815.64 --> -15815.55
## optimize topology: -15818.1 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 20
## [1] "Ratchet iteration 30 , best pscore so far: -15815.543758178"
## -16660.27 -> -16609.06
## -16637 -> -16603.21
## -16651.91 -> -16599.07
## -16636.41 -> -16591.7
## -16451.9 -> -16451.55
## -16622.98 -> -16527.11
## -16618.49 -> -16558.13
## -16618.69 -> -16557.26
## -16609.59 -> -16598.83
## -16602.51 -> -16591.16
## -16615.39 -> -16599.19
## -16609.35 -> -16596.89
## -16607.94 -> -16590.76
## -16612.56 -> -16593.83
## -16624.35 -> -16602.08
## -16625.13 -> -16602.08
## -16562.54 -> -16562.48
## -16615.16 -> -16595.73
## -16616.13 -> -16598.58
## -16606.43 -> -16597.45
## -16598.71 -> -16595.91
## -16614.7 -> -16605.65
## -16617.57 -> -16614.01
## -16613.38 -> -16610.4
## -16607.53 -> -16575.09
## -16559.64 -> -16535.52
## -16622.7 -> -16582.9
## -16611.87 -> -16576.88
## -16625.55 -> -16613.19
## -16619.54 -> -16611.04
## -16628.95 -> -16611.34
## -16639.18 -> -16593.05
## -16644.36 -> -16593.99
## -16677.22 -> -16592.92
## -16689.89 -> -16592.92
## -16378.45 -> -16378.44
## -16599.05 -> -16508.17
## -16635.44 -> -16613.87
## -16643.15 -> -16614.93
## -16574.01 -> -16573.6
## -16617.06 -> -16593.44
## -16626.55 -> -16614.4
## -16603.38 -> -16590.05
## -16605.09 -> -16591.19
## -16585.35 -> -16582.83
## -16618.45 -> -16596.25
## -16615.86 -> -16600.44
## -16613.96 -> -16600.15
## -16617.17 -> -16611.5
## -16622.84 -> -16613.77
## -16590.35 -> -16590.33
## -16612.05 -> -16607.36
## optimize edge weights: -15979.51 --> -15932.68
## optimize topology: -16615.16 --> -15932.68
## -15929.15 -> -15927.18
## -15932.38 -> -15931.59
## -15933.09 -> -15932.07
## -15932.68 -> -15932.34
## -15932.68 -> -15926.91
## -15932.58 -> -15932.28
## -15932.68 -> -15928.3
## -15932.68 -> -15931.08
## -15924.9 -> -15921.22
## -15930.97 -> -15929.67
## -15929.92 -> -15928.71
## -15933.58 -> -15932.54
## -15929.38 -> -15925.36
## -15930.66 -> -15928.63
## optimize edge weights: -15907.77 --> -15907.34
## optimize topology: -15932.68 --> -15907.34
## -15901.22 -> -15887.66
## -15907.96 -> -15907.26
## -15907.34 -> -15894.77
## -15907.34 -> -15905.8
## -15907.34 -> -15903.03
## -15907.34 -> -15907.29
## -15907.34 -> -15905.87
## -15905.03 -> -15903.76
## -15908.23 -> -15907.16
## optimize edge weights: -15882.84 --> -15882.55
## optimize topology: -15907.34 --> -15882.55
## -15882.55 -> -15864.31
## -15882.55 -> -15880.92
## -15878.45 -> -15877.77
## -15881.48 -> -15880.86
## -15882.55 -> -15872.15
## -15882.55 -> -15881.05
## optimize edge weights: -15852.28 --> -15852.02
## optimize topology: -15882.55 --> -15852.02
## -15843.04 -> -15821.36
## -15852.56 -> -15851.89
## optimize edge weights: -15821.04 --> -15820.53
## optimize topology: -15852.02 --> -15820.53
## optimize edge weights: -15820.53 --> -15820.53
## optimize topology: -15820.53 --> -15820.53
## NNI moves: 20
## [1] "Ratchet iteration 31 , best pscore so far: -15815.543758178"
## -16497.39 -> -16477.3
## -16487.8 -> -16475.32
## -16476.93 -> -16476.64
## -16505.3 -> -16486.78
## -16489.77 -> -16481.19
## -16496.8 -> -16485.1
## -16486.95 -> -16484.68
## -16492.69 -> -16489.3
## -16489.51 -> -16474.27
## -16497.15 -> -16479.59
## -16423.31 -> -16421.12
## -16488.41 -> -16465.68
## -16481.36 -> -16479.8
## -16493.41 -> -16479.3
## -16470.99 -> -16467.86
## -16531.71 -> -16467.05
## -16528.99 -> -16466.96
## -16315.51 -> -16315.17
## -16484.62 -> -16385.33
## -16494.44 -> -16425.12
## -16497.03 -> -16426.11
## -16496.04 -> -16470.99
## -16518.13 -> -16476.86
## -16451.13 -> -16451.02
## -16489.66 -> -16480.16
## -16494.93 -> -16486.99
## -16496.15 -> -16487.29
## -16483.58 -> -16469.23
## -16493.81 -> -16469.73
## -16478.64 -> -16464.76
## -16493.9 -> -16468.72
## -16493.57 -> -16492.2
## -16495.79 -> -16492.93
## -16493.19 -> -16492.75
## -16495.16 -> -16490.47
## -16494.43 -> -16491.82
## -16480.94 -> -16443.59
## -16500.19 -> -16455.78
## -16263.27 -> -16262.79
## -16513.95 -> -16394.27
## -16546.8 -> -16480.5
## -16537.54 -> -16481
## optimize edge weights: -15928.33 --> -15915.04
## optimize topology: -16493.02 --> -15915.04
## -15913.16 -> -15911.8
## -15913.51 -> -15912.07
## -15914.68 -> -15914.63
## -15915.04 -> -15912.16
## -15915.04 -> -15876.91
## -15906.61 -> -15901.78
## -15911.85 -> -15911.24
## -15915.47 -> -15914.59
## -15915.04 -> -15911.44
## -15911.89 -> -15911.63
## optimize edge weights: -15857.2 --> -15856.82
## optimize topology: -15915.04 --> -15856.82
## -15855.84 -> -15850.16
## -15856.62 -> -15856.4
## -15855.31 -> -15855.18
## -15848.64 -> -15840.97
## -15856.82 -> -15856.82
## -15856.82 -> -15856.82
## -15853.91 -> -15853.34
## optimize edge weights: -15835.28 --> -15834.37
## optimize topology: -15856.82 --> -15834.37
## -15833.06 -> -15830.88
## -15833.09 -> -15827.12
## optimize edge weights: -15823.96 --> -15823.77
## optimize topology: -15834.37 --> -15823.77
## -15823.32 -> -15823.2
## optimize edge weights: -15823.2 --> -15823.14
## optimize topology: -15823.77 --> -15823.14
## -15819.6 -> -15818.01
## -15822.92 -> -15822.83
## optimize edge weights: -15818.01 --> -15817.72
## optimize topology: -15823.14 --> -15817.72
## -15816.06 -> -15815.59
## optimize edge weights: -15815.59 --> -15815.55
## optimize topology: -15817.72 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 19
## [1] "Ratchet iteration 32 , best pscore so far: -15815.543758178"
## -16311.67 -> -16278.99
## -16246.91 -> -16239.72
## -16356.11 -> -16290.07
## -16369.24 -> -16290.31
## -16309.68 -> -16242.46
## -16305.68 -> -16238.69
## -16138.95 -> -16138.71
## -16310.99 -> -16211.26
## -16323.24 -> -16301.23
## -16311.27 -> -16298.13
## -16312.27 -> -16297.75
## -16293.48 -> -16289.24
## -16288.98 -> -16288.91
## -16312.22 -> -16304.77
## -16309.72 -> -16288.46
## -16300.95 -> -16283.09
## -16289.97 -> -16279.28
## -16320.98 -> -16299.08
## -16312.53 -> -16284.6
## -16298 -> -16279.68
## -16310.47 -> -16305.96
## -16309.28 -> -16306.38
## -16312.26 -> -16280.11
## -16325.11 -> -16280.68
## -16229.79 -> -16228.9
## -16294.44 -> -16269.86
## -16313.72 -> -16286.83
## -16320.83 -> -16286.66
## -16285.49 -> -16284.5
## -16303.4 -> -16295.96
## -16298.71 -> -16290.82
## -16312.93 -> -16301.45
## -16313.46 -> -16305.91
## -16301.69 -> -16292.63
## -16305.93 -> -16304.28
## -16289.84 -> -16289.77
## -16302.38 -> -16299.33
## -16297.72 -> -16297.72
## -16307.93 -> -16297.72
## -16310.93 -> -16304.5
## -16313.2 -> -16304.84
## optimize edge weights: -15989.47 --> -15944.26
## optimize topology: -16306.98 --> -15944.26
## -15936.18 -> -15907.35
## -15942.04 -> -15941.77
## -15941.92 -> -15938.97
## -15938.3 -> -15937.14
## -15940.88 -> -15933.31
## -15938.75 -> -15934.67
## -15935.14 -> -15933.18
## -15943.38 -> -15943.21
## optimize edge weights: -15889.9 --> -15888.17
## optimize topology: -15944.26 --> -15888.17
## -15886.43 -> -15885.68
## -15886.18 -> -15886.05
## -15885.24 -> -15884.73
## -15881.01 -> -15879.11
## -15876.24 -> -15854.65
## -15886.9 -> -15886.78
## -15886.81 -> -15884.54
## optimize edge weights: -15841.4 --> -15840.69
## optimize topology: -15888.17 --> -15840.69
## -15840.69 -> -15840.31
## -15840.69 -> -15815.77
## -15838.4 -> -15836.28
## optimize edge weights: -15815.77 --> -15815.54
## optimize topology: -15840.69 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 14
## [1] "Ratchet iteration 33 , best pscore so far: -15815.543758178"
## -16187.88 -> -16175.4
## -16160.17 -> -16155.01
## -16180.08 -> -16168.53
## -16187.97 -> -16172.35
## -16110.63 -> -16109.97
## -16188.49 -> -16156.38
## -16149.31 -> -16146.08
## -16173.02 -> -16164.78
## -16178.09 -> -16174.8
## -16176.49 -> -16173
## -16154.08 -> -16154.07
## -16178.16 -> -16172.26
## -16172.98 -> -16168.79
## -16181.61 -> -16171.7
## -16171.96 -> -16143.23
## -16175.9 -> -16148.22
## -16183.35 -> -16155.49
## -16195.95 -> -16159.21
## -16188.84 -> -16173.51
## -16184.39 -> -16170.69
## -16174.99 -> -16174.76
## -16156.9 -> -16135.55
## -16170.83 -> -16151.34
## -16155.99 -> -16146.68
## -16180.74 -> -16164.22
## -16188.56 -> -16173.73
## -16170.97 -> -16164.89
## -16177.51 -> -16175.29
## -16176.51 -> -16174.99
## -16134.62 -> -16125.25
## -16154.28 -> -16136.89
## -16173.09 -> -16173.07
## -16175.65 -> -16175.02
## -16175.33 -> -16158.61
## -16166.12 -> -16150.7
## optimize edge weights: -15986.82 --> -15955.86
## optimize topology: -16175.92 --> -15955.86
## -15945.5 -> -15941.62
## -15952 -> -15939.92
## -15954.98 -> -15953.89
## -15953.71 -> -15949.72
## -15955 -> -15954.51
## -15956.29 -> -15954.7
## -15955.91 -> -15955.81
## -15952.48 -> -15942.22
## -15949.62 -> -15941.41
## -15954.82 -> -15954.28
## -15953.33 -> -15952.37
## -15954.9 -> -15953.88
## optimize edge weights: -15907.37 --> -15906.56
## optimize topology: -15955.86 --> -15906.56
## -15895.73 -> -15891.99
## -15903.15 -> -15897.99
## -15906.12 -> -15906.03
## -15898.7 -> -15895.19
## optimize edge weights: -15880.53 --> -15879.12
## optimize topology: -15906.56 --> -15879.12
## -15877.33 -> -15876.65
## -15873.54 -> -15865.29
## -15879.12 -> -15870.01
## -15879.12 -> -15878.68
## optimize edge weights: -15865.29 --> -15865.13
## optimize topology: -15879.12 --> -15865.13
## -15865.13 -> -15856.46
## -15865.13 -> -15864.7
## optimize edge weights: -15856.46 --> -15856.26
## optimize topology: -15865.13 --> -15856.26
## -15854.7 -> -15854.09
## -15841.57 -> -15815.94
## optimize edge weights: -15815.94 --> -15815.84
## optimize topology: -15856.26 --> -15815.84
## -15815.71 -> -15815.61
## optimize edge weights: -15815.61 --> -15815.54
## optimize topology: -15815.84 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 16
## [1] "Ratchet iteration 34 , best pscore so far: -15815.543758178"
## -16303.28 -> -16268.49
## -16300.67 -> -16265.25
## -16303.19 -> -16292.32
## -16307.77 -> -16297.25
## -16302.68 -> -16294.42
## -16304.1 -> -16294.71
## -16305.22 -> -16300.78
## -16304.28 -> -16299.79
## -16280.29 -> -16280.28
## -16303.25 -> -16297.35
## -16291.41 -> -16286.13
## -16311.02 -> -16301.11
## -16291.6 -> -16271.84
## -16311.62 -> -16278.9
## -16304.19 -> -16276.6
## -16243.7 -> -16227.53
## -16288 -> -16247.8
## -16304.83 -> -16261.38
## -16286.4 -> -16272.68
## -16294.97 -> -16274.84
## -16312.04 -> -16300.48
## -16287.97 -> -16277.21
## -16300.64 -> -16271.53
## -16286.01 -> -16269.88
## -16302.91 -> -16287.41
## -16297.81 -> -16278.7
## -16325.13 -> -16299.6
## -16316.76 -> -16296.09
## -16305.49 -> -16300.42
## -16313.93 -> -16297.09
## -16317.32 -> -16297.09
## -16239.2 -> -16238.35
## -16307.98 -> -16277.99
## -16320.13 -> -16295.88
## -16324.31 -> -16298.55
## optimize edge weights: -16102.04 --> -16055.95
## optimize topology: -16303.71 --> -16055.95
## -16055.95 -> -16055.95
## -16055.95 -> -16055.86
## -16054.19 -> -16052.4
## -16052.56 -> -16050.55
## -16052.53 -> -16050.8
## -16054.13 -> -16053.92
## -16055.95 -> -16046.82
## -16055.95 -> -16055.56
## -16055.95 -> -16010.42
## -16055.95 -> -16048.22
## -16055.95 -> -16053.44
## -16048.21 -> -16043.16
## -16054 -> -16053.8
## optimize edge weights: -16002.04 --> -16000.52
## optimize topology: -16055.95 --> -16000.52
## -15996.71 -> -15990.23
## -16000.25 -> -15999.74
## -16001.22 -> -16000.11
## -15998.21 -> -15996.76
## -16000.52 -> -15985.62
## -16000.52 -> -15995.4
## -16000.52 -> -15991.94
## -16000.52 -> -15999.18
## -15995.04 -> -15992.03
## -15999.77 -> -15999.73
## -15992.9 -> -15980.61
## optimize edge weights: -15955.55 --> -15955.24
## optimize topology: -16000.52 --> -15955.24
## -15954.03 -> -15953.92
## -15955.24 -> -15955.06
## -15955.49 -> -15954.39
## -15952.09 -> -15949.86
## -15955.24 -> -15946.13
## -15955.24 -> -15955.24
## -15949.05 -> -15945.47
## optimize edge weights: -15934.98 --> -15934.74
## optimize topology: -15955.24 --> -15934.74
## -15934.74 -> -15934.74
## -15933.63 -> -15932.25
## -15929.88 -> -15927.28
## -15933.6 -> -15931.08
## -15932.57 -> -15928.57
## -15934.4 -> -15934.19
## optimize edge weights: -15927.28 --> -15926.87
## optimize topology: -15934.74 --> -15926.87
## -15926.17 -> -15923.44
## -15926.5 -> -15922.71
## -15926.87 -> -15925.46
## -15926.87 -> -15926.66
## -15925.52 -> -15923.39
## -15924.25 -> -15921.88
## -15926.53 -> -15926.35
## optimize edge weights: -15921.88 --> -15921.8
## optimize topology: -15926.87 --> -15921.8
## -15921.07 -> -15918.29
## -15921.43 -> -15917.7
## -15921.8 -> -15920.94
## -15921.8 -> -15921.6
## -15921.8 -> -15913.59
## optimize edge weights: -15909.5 --> -15909.44
## optimize topology: -15921.8 --> -15909.44
## -15909.44 -> -15908.59
## -15909.44 -> -15909.29
## -15908.93 -> -15908.67
## optimize edge weights: -15908.59 --> -15908.57
## optimize topology: -15909.44 --> -15908.57
## -15908.06 -> -15907.79
## -15896.03 -> -15892.38
## optimize edge weights: -15892.38 --> -15892.05
## optimize topology: -15908.57 --> -15892.05
## -15892.32 -> -15890.3
## -15886.7 -> -15878.84
## -15892.16 -> -15892.03
## -15891.52 -> -15891.24
## optimize edge weights: -15878.04 --> -15876.78
## optimize topology: -15892.05 --> -15876.78
## -15865.29 -> -15859.27
## -15868.09 -> -15863.85
## optimize edge weights: -15859.27 --> -15859.05
## optimize topology: -15876.78 --> -15859.05
## -15847.71 -> -15828.23
## -15857.96 -> -15856.82
## optimize edge weights: -15828.23 --> -15827.69
## optimize topology: -15859.05 --> -15827.69
## -15825.15 -> -15823.79
## optimize edge weights: -15823.79 --> -15823.74
## optimize topology: -15827.69 --> -15823.74
## -15823.62 -> -15823.5
## -15822.22 -> -15815.87
## optimize edge weights: -15815.87 --> -15815.84
## optimize topology: -15823.74 --> -15815.84
## -15815.71 -> -15815.61
## optimize edge weights: -15815.61 --> -15815.54
## optimize topology: -15815.84 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 26
## [1] "Ratchet iteration 35 , best pscore so far: -15815.543758178"
## -16364.68 -> -16364.52
## -16372.5 -> -16361.64
## -16389.12 -> -16364.06
## -16332.19 -> -16317.56
## -16366.39 -> -16337.55
## -16354.01 -> -16331.92
## -16362.38 -> -16335
## -16365.19 -> -16345.6
## -16316.05 -> -16314.56
## -16364.83 -> -16364.44
## -16364.17 -> -16363.84
## -16356.91 -> -16353.08
## -16368.64 -> -16361.49
## -16370.11 -> -16363.25
## -16368.48 -> -16361.89
## -16374.77 -> -16362.43
## -16379.66 -> -16362.45
## -16337.62 -> -16321.23
## -16379.24 -> -16347.41
## -16377.98 -> -16359.95
## -16386.46 -> -16364.74
## -16366.54 -> -16335.98
## -16343.08 -> -16324.55
## -16304.03 -> -16279.99
## -16364.19 -> -16324.37
## -16361.5 -> -16361.37
## -16340.8 -> -16340.49
## -16366.39 -> -16357.38
## -16389.36 -> -16350.32
## -16399.1 -> -16354.68
## -16368.25 -> -16301.58
## -16362.65 -> -16296.13
## -16188.72 -> -16188.43
## -16374.75 -> -16271.58
## optimize edge weights: -16032.29 --> -15998.06
## optimize topology: -16365.25 --> -15998.06
## -15991.18 -> -15988.33
## -15998.06 -> -15998
## -15998.06 -> -15994.8
## -15994.16 -> -15991.29
## -15993.67 -> -15991.51
## -15997.19 -> -15997.11
## -15995.14 -> -15962.42
## -15989.62 -> -15985.1
## optimize edge weights: -15932.79 --> -15928.08
## optimize topology: -15998.06 --> -15928.08
## -15927.83 -> -15927.58
## -15928.08 -> -15928.02
## -15928.08 -> -15924.57
## -15912.41 -> -15883.53
## -15917.12 -> -15908.19
## -15926.98 -> -15926.89
## optimize edge weights: -15879.47 --> -15878.45
## optimize topology: -15928.08 --> -15878.45
## -15878.45 -> -15878.41
## -15878.45 -> -15841.15
## -15874.8 -> -15862.75
## -15877.57 -> -15876.4
## optimize edge weights: -15825.45 --> -15825.09
## optimize topology: -15878.45 --> -15825.09
## -15825.09 -> -15815.55
## -15825.09 -> -15825.09
## -15825.09 -> -15825.09
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15825.09 --> -15815.55
## -15815.55 -> -15815.55
## -15815.55 -> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## optimize edge weights: -15815.55 --> -15815.55
## optimize topology: -15815.55 --> -15815.55
## NNI moves: 17
## [1] "Ratchet iteration 36 , best pscore so far: -15815.543758178"
## -16107.48 -> -16105.95
## -16113.01 -> -16103.42
## -16097.64 -> -16095.38
## -16114.97 -> -16100.07
## -16114.23 -> -16100.07
## -16076.12 -> -16076.09
## -16096.84 -> -16090.64
## -16084.92 -> -16084.17
## -16117.01 -> -16103.12
## -16122.9 -> -16094.6
## -16115.99 -> -16093.34
## -16120.13 -> -16105.65
## -16120.19 -> -16105.95
## -16045.08 -> -16044.74
## -16107.41 -> -16088.59
## -16114.69 -> -16100.94
## -16107.12 -> -16094.86
## -16114.69 -> -16107.61
## -16112.92 -> -16101.41
## -16118.79 -> -16103.09
## -16108.46 -> -16106.17
## -16111.29 -> -16086.79
## -16106.59 -> -16084.21
## -16101.17 -> -16077.6
## -16040.17 -> -16035.52
## -16133.32 -> -16104.91
## -16127.71 -> -16105.49
## -16050.78 -> -16050.61
## -16114.5 -> -16086.62
## -16120.95 -> -16103.48
## -16114.93 -> -16102.12
## -16106.69 -> -16106.69
## -16107.44 -> -16106.73
## optimize edge weights: -15880.81 --> -15866.94
## optimize topology: -16108.59 --> -15866.94
## -15863.11 -> -15859.6
## -15858.27 -> -15849.46
## -15865.59 -> -15865.38
## -15861.06 -> -15851.42
## -15866.31 -> -15866.08
## optimize edge weights: -15840.3 --> -15839.71
## optimize topology: -15866.94 --> -15839.71
## -15839.24 -> -15838.99
## -15833.24 -> -15832.79
## -15839.56 -> -15838.77
## -15840.54 -> -15839.34
## -15835.12 -> -15823.15
## -15838.95 -> -15838.19
## optimize edge weights: -15822.42 --> -15822.27
## optimize topology: -15839.71 --> -15822.27
## -15816.2 -> -15815.67
## -15821.24 -> -15820.59
## optimize edge weights: -15815.67 --> -15815.54
## optimize topology: -15822.27 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 11
## [1] "Ratchet iteration 37 , best pscore so far: -15815.543758178"
## -16605.62 -> -16518.4
## -16598.82 -> -16518.16
## -16540.23 -> -16515.43
## -16541.76 -> -16517.64
## -16400.23 -> -16387.06
## -16543.74 -> -16453.55
## -16529.17 -> -16502.73
## -16542.69 -> -16510.52
## -16541.05 -> -16532.49
## -16545.58 -> -16534.77
## -16502 -> -16501.32
## -16532.2 -> -16520.22
## -16539.06 -> -16504.31
## -16491.19 -> -16475.71
## -16533.05 -> -16519.85
## -16513.34 -> -16497.96
## -16537.27 -> -16533.79
## -16539.15 -> -16535.24
## -16493.7 -> -16493.55
## -16531.05 -> -16521.94
## -16553.77 -> -16529.75
## -16550.23 -> -16529.57
## -16500.69 -> -16486.31
## -16536.44 -> -16511.13
## -16537.02 -> -16514.31
## -16510.62 -> -16505.08
## -16537.36 -> -16536.12
## -16528.27 -> -16527.77
## -16514.3 -> -16511.54
## -16539.57 -> -16533.82
## -16537.92 -> -16533.82
## -16537.55 -> -16536.5
## -16553.4 -> -16514.63
## -16562.84 -> -16517
## -16594.79 -> -16510.53
## -16607.77 -> -16510.53
## -16300.31 -> -16300.31
## -16518.86 -> -16429.75
## optimize edge weights: -16040.66 --> -15996.15
## optimize topology: -16536.55 --> -15996.15
## -15995.22 -> -15994.46
## -15993.7 -> -15991.46
## -15996.15 -> -15995.33
## -15996.15 -> -15947.76
## -15989.56 -> -15981.43
## -15996.15 -> -15968.19
## -15996.15 -> -15995.85
## -15996.09 -> -15996.04
## -15996.15 -> -15955.35
## -15994.67 -> -15993.96
## optimize edge weights: -15878.22 --> -15876
## optimize topology: -15996.15 --> -15876
## -15874.9 -> -15873.96
## -15870.99 -> -15861.87
## -15876 -> -15869.88
## -15866.89 -> -15852.77
## -15874.15 -> -15873.53
## optimize edge weights: -15852.77 --> -15851.01
## optimize topology: -15876 --> -15851.01
## -15849.45 -> -15848.43
## -15844.85 -> -15835.78
## -15850.59 -> -15837.11
## -15850.83 -> -15848.79
## -15850.38 -> -15850.22
## -15848.61 -> -15846.12
## optimize edge weights: -15835.36 --> -15834.59
## optimize topology: -15851.01 --> -15834.59
## -15830.11 -> -15820.6
## -15833.11 -> -15832.25
## -15832.01 -> -15829.59
## optimize edge weights: -15820.6 --> -15820.48
## optimize topology: -15834.59 --> -15820.48
## -15818.09 -> -15816.08
## optimize edge weights: -15816.08 --> -15815.84
## optimize topology: -15820.48 --> -15815.84
## -15815.71 -> -15815.61
## optimize edge weights: -15815.61 --> -15815.54
## optimize topology: -15815.84 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 15
## [1] "Ratchet iteration 38 , best pscore so far: -15815.543758178"
## -16584.35 -> -16575.07
## -16588.15 -> -16577.93
## -16590.39 -> -16587.03
## -16586.66 -> -16583.47
## -16563.2 -> -16542.37
## -16597.94 -> -16556.46
## -16587.76 -> -16573.12
## -16595.53 -> -16576.86
## -16597.1 -> -16575.24
## -16597.89 -> -16575.24
## -16529.82 -> -16529.82
## -16587.93 -> -16568.89
## -16496.06 -> -16479.09
## -16585.08 -> -16546.77
## -16601.64 -> -16579.41
## -16600.59 -> -16579.41
## -16605.42 -> -16578.33
## -16593.88 -> -16572.66
## -16609.19 -> -16570.89
## -16617.7 -> -16573.1
## -16646.08 -> -16579.68
## -16654.88 -> -16579.68
## -16371.71 -> -16371.71
## -16604.43 -> -16499.82
## -16575.03 -> -16570.06
## -16584.22 -> -16578.43
## -16586.05 -> -16585.76
## -16588.34 -> -16587.58
## -16597.46 -> -16577.14
## -16584.87 -> -16569.61
## -16588.85 -> -16565.37
## -16596.31 -> -16569.31
## -16561.66 -> -16556.71
## -16590.32 -> -16578.46
## -16582.68 -> -16554.91
## -16556.02 -> -16543.99
## -16550.36 -> -16531.62
## -16593.96 -> -16561.58
## -16608.04 -> -16575.74
## -16594.02 -> -16569.11
## -16591.16 -> -16522.2
## -16585.79 -> -16517.08
## -16415.13 -> -16415.11
## -16598.68 -> -16492.79
## optimize edge weights: -15953.6 --> -15936.25
## optimize topology: -16587.93 --> -15936.25
## -15934.15 -> -15933.23
## -15924.66 -> -15903.46
## -15934.97 -> -15934.8
## -15936.22 -> -15935.63
## -15924.63 -> -15917.79
## -15936.25 -> -15890.68
## -15936.25 -> -15934.83
## optimize edge weights: -15857.26 --> -15855.58
## optimize topology: -15936.25 --> -15855.58
## -15854.51 -> -15854.22
## -15847.3 -> -15842.71
## -15855.57 -> -15855.36
## -15854.01 -> -15852.9
## -15855.63 -> -15855.33
## -15854.73 -> -15853.43
## -15855.06 -> -15854.02
## -15854.62 -> -15853.9
## -15846.75 -> -15841.82
## optimize edge weights: -15826.98 --> -15826.74
## optimize topology: -15855.58 --> -15826.74
## -15826.4 -> -15825.81
## -15824.58 -> -15823.51
## optimize edge weights: -15823.51 --> -15823.45
## optimize topology: -15826.74 --> -15823.45
## -15821.95 -> -15815.57
## optimize edge weights: -15815.57 --> -15815.54
## optimize topology: -15823.45 --> -15815.54
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 14
## [1] "Ratchet iteration 39 , best pscore so far: -15815.543758178"
## optimize edge weights: -15815.54 --> -15815.54
## optimize topology: -15815.54 --> -15815.54
## NNI moves: 0Is this analysis faster or slower than the NJ analysis you did earlier?
Get the log likelihood by typing:
logLik(fitJC)## 'log Lik.' -15815.54 (df=77)Question 3: What is this log likelihood and is it higher or lower than the log likelihood of the NJ analysis?
Look at the tree topology by typing:
plot(fitJC)
How does this tree compare to the neighbor joining tree you got earlier?
Save the tree as a pdf if you want. Write the tree to a file by typing:
write.tree(fitJC$tree,file="fitJC")Now change the model of evolution to the General Time Reversible model (GTR). Type:
fitGTR<-optim.pml(fit,model="GTR",rearrangement="stochastic")## optimize edge weights: -15999.05 --> -15909.49
## optimize base frequencies: -15909.49 --> -15899.83
## optimize rate matrix: -15899.83 --> -15473.94
## optimize edge weights: -15473.94 --> -15471
## optimize topology: -15471 --> -15439.66
## optimize topology: -15439.66 --> -15432.77
## optimize topology: -15432.77 --> -15419.57
## optimize topology: -15419.57 --> -15405.04
## optimize topology: -15405.04 --> -15396.54
## NNI moves: 13
## optimize base frequencies: -15396.54 --> -15386.1
## optimize rate matrix: -15386.1 --> -15381.44
## optimize edge weights: -15381.44 --> -15381.41
## optimize topology: -15381.41 --> -15379.91
## optimize topology: -15379.91 --> -15379.58
## optimize topology: -15379.58 --> -15378.68
## optimize topology: -15378.68 --> -15367.25
## optimize topology: -15367.25 --> -15361.79
## NNI moves: 5
## optimize base frequencies: -15361.79 --> -15360.06
## optimize rate matrix: -15360.06 --> -15359.47
## optimize edge weights: -15359.47 --> -15359.43
## optimize topology: -15359.43 --> -15355.59
## optimize topology: -15355.59 --> -15354.16
## optimize topology: -15354.16 --> -15354.16
## NNI moves: 2
## optimize base frequencies: -15354.16 --> -15353.97
## optimize rate matrix: -15353.97 --> -15353.87
## optimize edge weights: -15353.87 --> -15353.86
## optimize topology: -15353.86 --> -15353.86
## NNI moves: 0
## -15895 -> -15885.96
## -15889.58 -> -15877.71
## -15886.52 -> -15876.63
## -15888.25 -> -15885.57
## -15891.01 -> -15885.57
## -15855.09 -> -15854.83
## -15886.71 -> -15873.31
## -15911.29 -> -15842.68
## -15910.67 -> -15841.73
## -15898.53 -> -15873.43
## -15898.69 -> -15874.61
## -15712.78 -> -15712.65
## -15902.99 -> -15801.22
## -15898.22 -> -15867.21
## -15900.23 -> -15867.12
## -15865.09 -> -15832.38
## -15887.01 -> -15845.94
## -15882.05 -> -15856.21
## -15888.76 -> -15857.46
## -15897.9 -> -15880.72
## -15898.85 -> -15881.68
## -15886.37 -> -15874.84
## -15883.36 -> -15873.73
## -15882.76 -> -15846.68
## -15886.37 -> -15848.21
## -15877.94 -> -15872.11
## -15884.5 -> -15873.82
## -15810.43 -> -15805.46
## -15887.51 -> -15852.14
## -15893.48 -> -15880.21
## -15893.24 -> -15880.21
## -15886.11 -> -15886.06
## -15886.18 -> -15886.04
## -15860.14 -> -15857.08
## -15883.97 -> -15868.76
## -15861.19 -> -15838.86
## -15883.9 -> -15869.06
## -15901.06 -> -15886.34
## -15861.1 -> -15856.68
## -15888.09 -> -15872.79
## -15884.9 -> -15881.31
## -15882.94 -> -15881.26
## optimize edge weights: -15571.68 --> -15524.09
## optimize topology: -15886.42 --> -15524.09
## -15521.39 -> -15520.84
## -15523.34 -> -15523.05
## -15515.99 -> -15489.47
## -15517.07 -> -15508.03
## -15522.68 -> -15522.18
## -15524.18 -> -15523.79
## -15519.58 -> -15517.14
## -15518.83 -> -15510.28
## -15516.21 -> -15508.19
## -15524.09 -> -15524.09
## -15524.09 -> -15522.25
## -15524.09 -> -15524.09
## -15524.09 -> -15524.09
## -15524.09 -> -15524.09
## optimize edge weights: -15471.64 --> -15470.82
## optimize topology: -15524.09 --> -15470.82
## -15468.79 -> -15468.21
## -15464.7 -> -15453.88
## -15469.43 -> -15468.86
## -15470.79 -> -15470.22
## -15466.13 -> -15465.2
## -15465.71 -> -15457.85
## -15469.83 -> -15469.74
## -15470.89 -> -15470.69
## -15470.82 -> -15455.57
## -15470.82 -> -15470.73
## -15470.82 -> -15470.82
## optimize edge weights: -15425.01 --> -15424.4
## optimize topology: -15470.82 --> -15424.4
## -15422.33 -> -15420.87
## -15421.95 -> -15421.11
## -15423.77 -> -15423.37
## -15422.04 -> -15420.55
## -15426.57 -> -15423.83
## -15419.17 -> -15382.55
## -15422.71 -> -15421.1
## optimize edge weights: -15378.9 --> -15378.59
## optimize topology: -15424.4 --> -15378.59
## -15370.53 -> -15356.72
## -15378.58 -> -15378.54
## -15380.82 -> -15378.16
## -15376.59 -> -15376.02
## -15378.33 -> -15378.02
## optimize edge weights: -15356.72 --> -15356.54
## optimize topology: -15378.59 --> -15356.54
## -15354.55 -> -15353.99
## -15356.79 -> -15356.38
## optimize edge weights: -15353.99 --> -15353.87
## optimize topology: -15356.54 --> -15353.87
## optimize edge weights: -15353.87 --> -15353.87
## optimize topology: -15353.87 --> -15353.87
## NNI moves: 16
## [1] "Ratchet iteration 1 , best pscore so far: -15353.864132096"
## -15859.31 -> -15844.92
## -15856.67 -> -15844.89
## -15851.85 -> -15850.59
## -15853.17 -> -15851.27
## -15820.58 -> -15771.9
## -15853.27 -> -15784.8
## -15816.57 -> -15770.44
## -15853.61 -> -15787.24
## -15897.37 -> -15840.51
## -15895.91 -> -15839.84
## -15898.34 -> -15798.87
## -15869.2 -> -15789.19
## -15872.84 -> -15844.87
## -15873.88 -> -15846.35
## -15859.3 -> -15847.34
## -15866.77 -> -15851.05
## -15869.52 -> -15849.37
## -15868.83 -> -15849.52
## -15787.14 -> -15786.92
## -15868.66 -> -15830.72
## -15867.22 -> -15840.5
## -15874.8 -> -15843.18
## -15844.52 -> -15844.29
## -15848.08 -> -15846.27
## -15831.47 -> -15831.46
## -15860.55 -> -15847.76
## -15855.14 -> -15849.5
## -15856.08 -> -15850.09
## -15854.78 -> -15854.78
## -15855.12 -> -15854.78
## -15864.93 -> -15833.9
## -15865.84 -> -15833.9
## -15830.5 -> -15819.1
## -15847.82 -> -15832.6
## -15785.94 -> -15785.61
## -15856.62 -> -15827.58
## -15859.51 -> -15846.58
## -15863.33 -> -15848.39
## -15829.01 -> -15828.92
## -15841.34 -> -15839.77
## optimize edge weights: -15658.8 --> -15589.08
## optimize topology: -15855.09 --> -15589.08
## -15589.08 -> -15588.28
## -15589.08 -> -15589.08
## -15581.02 -> -15568.99
## -15584.51 -> -15583.39
## -15579.95 -> -15577.35
## -15585.86 -> -15583.45
## -15586.02 -> -15584.51
## -15578.58 -> -15565.26
## -15581.66 -> -15577.36
## -15581.7 -> -15576.96
## -15589.16 -> -15588.94
## -15589.08 -> -15572.92
## -15589.08 -> -15582.49
## -15578.98 -> -15572.78
## optimize edge weights: -15520 --> -15518.2
## optimize topology: -15589.08 --> -15518.2
## -15518.2 -> -15496.77
## -15504.71 -> -15478.99
## -15514.8 -> -15511.95
## -15514.96 -> -15513.26
## -15518.2 -> -15425.29
## -15518.2 -> -15518.2
## -15511.6 -> -15508.08
## -15518.38 -> -15518.12
## optimize edge weights: -15385.84 --> -15384.94
## optimize topology: -15518.2 --> -15384.94
## -15378.09 -> -15363.21
## -15381.67 -> -15379.28
## -15381.48 -> -15379.8
## -15386.61 -> -15384.23
## -15382.47 -> -15381.63
## -15385.15 -> -15384.88
## optimize edge weights: -15354.09 --> -15353.86
## optimize topology: -15384.94 --> -15353.86
## optimize edge weights: -15353.86 --> -15353.86
## optimize topology: -15353.86 --> -15353.86
## NNI moves: 16
## [1] "Ratchet iteration 2 , best pscore so far: -15353.863179606"
## -15714.58 -> -15709.15
## -15713.33 -> -15709.13
## -15683.68 -> -15682.1
## -15720.82 -> -15702.59
## -15682.45 -> -15679.78
## -15687.42 -> -15683.6
## -15707.18 -> -15707.16
## -15729.09 -> -15699.15
## -15728.45 -> -15699.04
## -15698.83 -> -15665.95
## -15708.88 -> -15672.28
## -15736.85 -> -15706.58
## -15733.9 -> -15706.19
## -15706.02 -> -15705.86
## -15719.41 -> -15688.08
## -15716.55 -> -15686.75
## -15725.39 -> -15707.65
## -15735.57 -> -15707.73
## -15709.3 -> -15708.76
## -15712.9 -> -15706.8
## -15707.69 -> -15703.83
## -15711.32 -> -15700.53
## -15711.07 -> -15700.53
## -15669.69 -> -15664.45
## -15710.52 -> -15697.57
## -15689.45 -> -15686.09
## -15698.46 -> -15692.63
## -15715.82 -> -15708.55
## -15710.62 -> -15705.45
## -15732.79 -> -15705.1
## -15737.67 -> -15705.13
## -15686.06 -> -15685.29
## -15712.17 -> -15705.55
## -15665.97 -> -15657.57
## -15714.5 -> -15697.15
## -15702.29 -> -15702.13
## -15708.25 -> -15706.43
## optimize edge weights: -15556.68 --> -15503.01
## optimize topology: -15710.49 --> -15503.01
## -15491.47 -> -15482.07
## -15503.01 -> -15492.95
## -15503.01 -> -15492.79
## -15501.12 -> -15497.3
## -15499.72 -> -15499.17
## -15503.01 -> -15502.58
## -15503.01 -> -15494.65
## -15503.01 -> -15491.36
## -15503.01 -> -15501.55
## -15500.13 -> -15499.27
## -15496.52 -> -15484.64
## -15494.52 -> -15488.32
## -15503.01 -> -15498.86
## -15503.01 -> -15501.73
## optimize edge weights: -15447.87 --> -15447.43
## optimize topology: -15503.01 --> -15447.43
## -15444.87 -> -15439.72
## -15441.98 -> -15440.14
## -15447.43 -> -15447.43
## -15447.43 -> -15438.71
## -15447.43 -> -15436.04
## -15447.43 -> -15440.82
## -15439.01 -> -15432.49
## -15447.43 -> -15444.97
## -15447.43 -> -15445.51
## optimize edge weights: -15403.28 --> -15402.53
## optimize topology: -15447.43 --> -15402.53
## -15389.21 -> -15364.18
## -15397.51 -> -15392.43
## -15402.53 -> -15402.47
## optimize edge weights: -15364.18 --> -15364.04
## optimize topology: -15402.53 --> -15364.04
## -15359.05 -> -15354
## -15363.86 -> -15363.8
## optimize edge weights: -15354 --> -15353.86
## optimize topology: -15364.04 --> -15353.86
## optimize edge weights: -15353.86 --> -15353.86
## optimize topology: -15353.86 --> -15353.86
## NNI moves: 16
## [1] "Ratchet iteration 3 , best pscore so far: -15353.8621617233"
## -15797.21 -> -15792.7
## -15808.36 -> -15797.04
## -15812.46 -> -15786.97
## -15812.01 -> -15787.1
## -15827.85 -> -15799.19
## -15816.52 -> -15797.48
## -15815.67 -> -15795.11
## -15813.27 -> -15795.11
## -15762.91 -> -15749.95
## -15799.5 -> -15776.28
## -15799.95 -> -15775.59
## -15770.02 -> -15764.28
## -15805.24 -> -15791.91
## -15809.99 -> -15795.23
## -15740.79 -> -15734.12
## -15796.36 -> -15777.89
## -15811.26 -> -15795.33
## -15806.33 -> -15793.54
## -15798.02 -> -15797.01
## -15806.87 -> -15766.29
## -15810.61 -> -15765.7
## -15738.22 -> -15724.14
## -15796.52 -> -15768.84
## -15804.56 -> -15783.63
## -15806.16 -> -15782.2
## -15784.72 -> -15780.22
## -15802.94 -> -15782.06
## -15812.26 -> -15786.12
## -15818.86 -> -15786.18
## -15774.37 -> -15772.68
## -15785.32 -> -15780.44
## -15812.98 -> -15795.55
## -15820.74 -> -15796.68
## -15735.53 -> -15733.5
## -15813.55 -> -15774.86
## -15821.52 -> -15793.85
## -15809.97 -> -15787.67
## -15796.88 -> -15796.35
## -15799.53 -> -15798.38
## -15800.54 -> -15799.18
## -15773.56 -> -15773.4
## -15783.15 -> -15782.69
## optimize edge weights: -15616.96 --> -15540.24
## optimize topology: -15799.75 --> -15540.24
## -15539.28 -> -15539.2
## -15540.49 -> -15539.36
## -15539.82 -> -15539.81
## -15540.24 -> -15506.78
## -15540.24 -> -15540.24
## -15545.03 -> -15537.03
## -15540.24 -> -15503.39
## -15540.24 -> -15540.18
## -15538.92 -> -15538.5
## -15527.96 -> -15499.68
## -15535.75 -> -15524.07
## -15540.24 -> -15540.24
## -15540.24 -> -15540.1
## optimize edge weights: -15463.16 --> -15462.54
## optimize topology: -15540.24 --> -15462.54
## -15463.26 -> -15462.26
## -15462.54 -> -15429.76
## -15462.54 -> -15462.45
## -15454.6 -> -15434.2
## optimize edge weights: -15429.63 --> -15428.93
## optimize topology: -15462.54 --> -15428.93
## -15426.24 -> -15424.51
## -15428.77 -> -15428.46
## -15429.19 -> -15428.91
## -15428.96 -> -15428.91
## -15423.23 -> -15420.64
## -15424.39 -> -15421.43
## -15421.13 -> -15400.33
## -15427.3 -> -15427.12
## optimize edge weights: -15393.84 --> -15393.21
## optimize topology: -15428.93 --> -15393.21
## -15392.49 -> -15387.6
## -15393.12 -> -15392.69
## -15390.04 -> -15389.26
## -15390.52 -> -15389.28
## -15387.46 -> -15371.31
## -15387.57 -> -15386.67
## optimize edge weights: -15359.2 --> -15358.91
## optimize topology: -15393.21 --> -15358.91
## -15359.12 -> -15358.86
## -15354.61 -> -15353.75
## -15355.97 -> -15354.24
## -15358.54 -> -15358.47
## -15358.34 -> -15357.59
## optimize edge weights: -15353.75 --> -15353.22
## optimize topology: -15358.91 --> -15353.22
## -15353.43 -> -15353.21
## -15350.88 -> -15348.21
## -15353.22 -> -15353.04
## -15352.65 -> -15351.83
## -15352.57 -> -15352.34
## optimize edge weights: -15347.95 --> -15347.76
## optimize topology: -15353.22 --> -15347.76
## -15347.76 -> -15347.1
## -15347.76 -> -15347.68
## -15347.2 -> -15346.4
## -15347.31 -> -15346.85
## optimize edge weights: -15345.49 --> -15345.48
## optimize topology: -15347.76 --> -15345.48
## -15345.48 -> -15344.81
## -15345.48 -> -15345.31
## optimize edge weights: -15344.81 --> -15344.79
## optimize topology: -15345.48 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 24
## [1] "Ratchet iteration 4 , best pscore so far: -15344.7896495107"
## -15747.38 -> -15733.92
## -15752.4 -> -15735.6
## -15736.47 -> -15735.92
## -15734.17 -> -15734.02
## -15735.42 -> -15735.35
## -15738.15 -> -15738.01
## -15737.13 -> -15733.73
## -15742.15 -> -15737.34
## -15736.86 -> -15726.04
## -15728.37 -> -15722.64
## -15751.1 -> -15738.1
## -15692.6 -> -15692.32
## -15738.42 -> -15718.82
## -15741.3 -> -15720.07
## -15748.19 -> -15723.41
## -15736.27 -> -15723.43
## -15742.53 -> -15724.28
## -15695.71 -> -15692.45
## -15734.67 -> -15716.25
## -15744.11 -> -15724.04
## -15744.82 -> -15724.74
## -15749.22 -> -15728.32
## -15745.57 -> -15725.99
## -15748.8 -> -15734.08
## -15742.11 -> -15731.91
## -15721.27 -> -15720.57
## -15734.65 -> -15728.89
## -15728.48 -> -15721.5
## -15744.4 -> -15734.35
## -15749.9 -> -15725.57
## -15772.19 -> -15733.53
## -15691.62 -> -15691.42
## -15740.74 -> -15714.6
## -15760.34 -> -15687.71
## -15764.67 -> -15687.63
## -15751.11 -> -15723.97
## -15751.36 -> -15725.19
## -15558.15 -> -15558.14
## -15746.71 -> -15647.74
## optimize edge weights: -15432.53 --> -15394.38
## optimize topology: -15738.41 --> -15394.38
## -15393.32 -> -15393.12
## -15393.6 -> -15391.72
## -15394.38 -> -15392.96
## -15394.38 -> -15394.04
## -15392.08 -> -15371.46
## -15394.06 -> -15393.98
## -15391.67 -> -15390.59
## -15390.55 -> -15388.83
## optimize edge weights: -15367.33 --> -15366.92
## optimize topology: -15394.38 --> -15366.92
## -15365.55 -> -15364.68
## -15366.71 -> -15366.65
## -15363.99 -> -15362.86
## -15362.57 -> -15361.24
## optimize edge weights: -15359.65 --> -15359.41
## optimize topology: -15366.92 --> -15359.41
## -15358.93 -> -15358.73
## -15356.84 -> -15355.41
## optimize edge weights: -15355.41 --> -15355.4
## optimize topology: -15359.41 --> -15355.4
## -15356.25 -> -15355.11
## optimize edge weights: -15355.11 --> -15355.06
## optimize topology: -15355.4 --> -15355.06
## optimize edge weights: -15355.06 --> -15355.06
## optimize topology: -15355.06 --> -15355.06
## NNI moves: 13
## [1] "Ratchet iteration 5 , best pscore so far: -15344.7896495107"
## -16123.65 -> -16103.94
## -16147.48 -> -16130.29
## -16114.1 -> -16112.72
## -16143.67 -> -16124.37
## -16167.94 -> -16126.42
## -16174.91 -> -16126.04
## -16076.53 -> -16059.94
## -16141.7 -> -16112.81
## -16137.41 -> -16112.43
## -16145.28 -> -16115.63
## -16148.09 -> -16142.43
## -16142.14 -> -16116.27
## -16144.41 -> -16116.74
## -16134.19 -> -16090.12
## -16220.2 -> -16094.75
## -16082.13 -> -15982.8
## -16153.32 -> -16072.65
## -16137.11 -> -16136.28
## -16147.67 -> -16141.82
## -16130.26 -> -16130.17
## -16143.16 -> -16139.34
## -16152.39 -> -16130.04
## -16151.65 -> -16130.04
## -16094.05 -> -16093.98
## -16143.08 -> -16123.71
## -16142.44 -> -16138.84
## -16146.82 -> -16140.3
## -16146.53 -> -16142.55
## -16147.87 -> -16132.12
## -16111.97 -> -16073.01
## -16168.62 -> -16137.33
## -16154.91 -> -16132.3
## -16003.46 -> -16003.42
## -16159.58 -> -16070.68
## -16152.12 -> -16096.21
## -16151.77 -> -16095.57
## -16157.26 -> -16136.92
## -16155.45 -> -16136.68
## optimize edge weights: -15587.27 --> -15567.4
## optimize topology: -16143.08 --> -15567.4
## -15554.13 -> -15534.96
## -15563.06 -> -15551.35
## -15566.96 -> -15566.74
## -15567.4 -> -15552.43
## -15567.4 -> -15561.15
## -15567.4 -> -15564.34
## -15567.4 -> -15555.54
## -15554.87 -> -15424.72
## -15565.3 -> -15564.3
## -15567.4 -> -15565.75
## -15567.4 -> -15567.26
## -15567.4 -> -15566.27
## -15567.4 -> -15567.4
## optimize edge weights: -15370.05 --> -15369.22
## optimize topology: -15567.4 --> -15369.22
## -15369.22 -> -15350.46
## -15369.22 -> -15367.87
## -15369.22 -> -15367.42
## -15369.22 -> -15363.64
## optimize edge weights: -15344.86 --> -15344.79
## optimize topology: -15369.22 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 15
## [1] "Ratchet iteration 6 , best pscore so far: -15344.7896495107"
## -16046.69 -> -16041.72
## -16052.72 -> -16017.66
## -16062.91 -> -16017.66
## -15977.25 -> -15976.57
## -16037.42 -> -16012.41
## -16054.58 -> -16002.9
## -16052.07 -> -16003.1
## -16070.33 -> -15988.8
## -16075.77 -> -15989.51
## -15733.36 -> -15730.6
## -16061.24 -> -15909.08
## -16043.38 -> -15954.2
## -16051.87 -> -15958.16
## -16046.24 -> -16042.02
## -16044.46 -> -16030.17
## -16027.28 -> -16019.87
## -16042.31 -> -16038.63
## -16047.29 -> -16039.92
## -16083.03 -> -16040.9
## -16087.27 -> -16040.9
## -16027.07 -> -16022.6
## -16033.96 -> -16024.93
## -16047.56 -> -16023.69
## -16039.98 -> -16020.16
## -16057.11 -> -16040.66
## -16067.15 -> -16040.66
## -16034.62 -> -16034.6
## -16044.06 -> -16034.62
## -16043.25 -> -16029.74
## -16043.8 -> -16029.74
## -16034.94 -> -16027.81
## -16055.1 -> -16036.5
## -16041.28 -> -16038.34
## -16040.49 -> -16025.63
## -16040.37 -> -16025.15
## -16021.24 -> -15990.26
## -16041.36 -> -16025.13
## -16017.38 -> -16009.74
## -16054.25 -> -16022.85
## -16040.91 -> -16029.63
## -16036.91 -> -16026.99
## -16031.71 -> -16031.09
## -16047.91 -> -16037.15
## optimize edge weights: -15609.23 --> -15571.28
## optimize topology: -16042.52 --> -15571.28
## -15567.52 -> -15566.34
## -15565.39 -> -15547.95
## -15570.57 -> -15570.37
## -15568.13 -> -15565.88
## -15568.22 -> -15566.84
## -15564.89 -> -15561.98
## -15571.28 -> -15567.56
## -15571.28 -> -15555.15
## -15571.28 -> -15570.86
## -15571.28 -> -15570.82
## -15571.28 -> -15571.28
## -15570.15 -> -15564.12
## -15564.08 -> -15550.57
## -15569.12 -> -15568.38
## -15570.76 -> -15570.54
## optimize edge weights: -15511.19 --> -15507.16
## optimize topology: -15571.28 --> -15507.16
## -15504.08 -> -15501.96
## -15503.68 -> -15502.01
## -15507.16 -> -15487.78
## -15507.16 -> -15507.16
## -15506.57 -> -15506.08
## -15502.14 -> -15497.26
## -15507.16 -> -15482.08
## -15507.16 -> -15507.16
## -15507.16 -> -15506.75
## -15507.16 -> -15506.91
## -15507.16 -> -15507.16
## -15499.44 -> -15474.09
## -15506.66 -> -15506.51
## -15504.64 -> -15503.68
## optimize edge weights: -15423.95 --> -15423.51
## optimize topology: -15507.16 --> -15423.51
## -15423.3 -> -15422.23
## -15423.41 -> -15423.23
## -15417.42 -> -15414.32
## -15423.51 -> -15396.93
## -15423.51 -> -15423.43
## -15423.51 -> -15421.56
## -15423.51 -> -15423.51
## -15423.51 -> -15423.51
## -15420.57 -> -15419.25
## optimize edge weights: -15390.99 --> -15390.57
## optimize topology: -15423.51 --> -15390.57
## -15390.77 -> -15390.53
## -15383.46 -> -15378.51
## -15390.57 -> -15388.91
## -15390.57 -> -15390.57
## optimize edge weights: -15376.71 --> -15376.37
## optimize topology: -15390.57 --> -15376.37
## -15368.23 -> -15365.93
## -15373.12 -> -15371.82
## -15376.24 -> -15376.19
## optimize edge weights: -15365.93 --> -15365.69
## optimize topology: -15376.37 --> -15365.69
## -15365.69 -> -15348.84
## -15365.69 -> -15364.82
## optimize edge weights: -15348.84 --> -15348.45
## optimize topology: -15365.69 --> -15348.45
## -15347.56 -> -15344.92
## optimize edge weights: -15344.92 --> -15344.8
## optimize topology: -15348.45 --> -15344.8
## optimize edge weights: -15344.8 --> -15344.8
## optimize topology: -15344.8 --> -15344.8
## NNI moves: 23
## [1] "Ratchet iteration 7 , best pscore so far: -15344.7896495107"
## -15737.78 -> -15737.69
## -15745.21 -> -15743.08
## -15745.52 -> -15745.01
## -15743.25 -> -15742.86
## -15745.66 -> -15745
## -15751.98 -> -15735.05
## -15761.04 -> -15735.05
## -15688.14 -> -15688.04
## -15758.92 -> -15727.24
## -15756.92 -> -15735.74
## -15771.71 -> -15740.19
## -15735.53 -> -15730.12
## -15741.83 -> -15732.69
## -15732.09 -> -15721.82
## -15742.08 -> -15730.61
## -15749.48 -> -15732.33
## -15750.16 -> -15732.33
## -15735.6 -> -15722.77
## -15746.25 -> -15725.76
## -15742.89 -> -15727.62
## -15752.7 -> -15727.92
## -15719.57 -> -15719.4
## -15739.35 -> -15732.84
## -15726.98 -> -15701.28
## -15743.44 -> -15720.31
## -15733.06 -> -15726.93
## -15746.76 -> -15726.9
## -15762.88 -> -15739.91
## -15755.53 -> -15739.59
## -15769.58 -> -15740.78
## -15757.37 -> -15735.83
## -15609.05 -> -15609.02
## -15763.82 -> -15675.13
## -15754.86 -> -15699.17
## -15754.43 -> -15698.52
## optimize edge weights: -15459.43 --> -15435.26
## optimize topology: -15745.81 --> -15435.26
## -15432.09 -> -15431.05
## -15434.69 -> -15430.31
## -15434.92 -> -15434.09
## -15435.26 -> -15423.89
## -15435.26 -> -15433.4
## -15435.2 -> -15435.17
## -15435.26 -> -15401.55
## -15450.12 -> -15432.06
## -15425.78 -> -15401.83
## optimize edge weights: -15356.75 --> -15356.11
## optimize topology: -15435.26 --> -15356.11
## -15352.59 -> -15351.02
## -15356.11 -> -15353.98
## -15356.11 -> -15350.31
## optimize edge weights: -15345.32 --> -15345.16
## optimize topology: -15356.11 --> -15345.16
## -15345.18 -> -15344.86
## optimize edge weights: -15344.86 --> -15344.79
## optimize topology: -15345.16 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 12
## [1] "Ratchet iteration 8 , best pscore so far: -15344.7896495107"
## -15614.57 -> -15614.53
## -15617.5 -> -15615.79
## -15618.91 -> -15617.03
## -15607.15 -> -15606.8
## -15615.21 -> -15612.76
## -15622.09 -> -15613.13
## -15627.32 -> -15613.73
## -15631.28 -> -15615.18
## -15611.62 -> -15611.45
## -15620.89 -> -15601.89
## -15622.06 -> -15603.31
## -15584.5 -> -15582.24
## -15612.82 -> -15596.27
## -15574.44 -> -15574.3
## -15615.91 -> -15604.06
## -15618.56 -> -15616.87
## -15660.72 -> -15604.15
## -15640.63 -> -15599.34
## -15630.45 -> -15565.31
## -15626.68 -> -15565.45
## -15462.33 -> -15462.23
## -15598.06 -> -15525.64
## -15588.86 -> -15588.61
## -15615.67 -> -15607.52
## optimize edge weights: -15404.6 --> -15379.19
## optimize topology: -15617.06 --> -15379.19
## -15377.15 -> -15377.03
## -15379.19 -> -15364.04
## -15379.19 -> -15378.51
## -15372.41 -> -15363.48
## -15378.06 -> -15377.28
## optimize edge weights: -15345.99 --> -15345.31
## optimize topology: -15379.19 --> -15345.31
## -15344.89 -> -15344.82
## optimize edge weights: -15344.82 --> -15344.79
## optimize topology: -15345.31 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 8
## [1] "Ratchet iteration 9 , best pscore so far: -15344.7896495107"
## -16055.33 -> -16055.18
## -16058.27 -> -16056.69
## -16058.47 -> -16053.84
## -16061.77 -> -16054.97
## -16058.59 -> -16048.44
## -16058.33 -> -16048.44
## -16004.18 -> -15995.43
## -16057.34 -> -16029.83
## -16052.5 -> -16050.61
## -16058.86 -> -16057.9
## -16051.72 -> -16051.65
## -16065.59 -> -16056.81
## -16060.22 -> -16054.35
## -16064.68 -> -16053.96
## -16062.89 -> -16053.87
## -16043.95 -> -16029.05
## -16061.52 -> -16049.98
## -16055.53 -> -16055.4
## -16057.71 -> -16057.61
## -16060.86 -> -16051.05
## -16057.58 -> -16048.11
## -16044.65 -> -16034.52
## -16065.92 -> -16048.45
## -16062.77 -> -15995.05
## -16063.71 -> -15996.82
## -16140.43 -> -15958.37
## -16157.31 -> -15958.7
## -15695.9 -> -15637.67
## -16039.56 -> -15852.55
## -16177.63 -> -15960.65
## -16081.22 -> -15956.54
## optimize edge weights: -15539.8 --> -15494.16
## optimize topology: -16058.27 --> -15494.16
## -15493.09 -> -15492.68
## -15494.54 -> -15494.14
## -15494.16 -> -15494.16
## -15492.02 -> -15490.05
## -15494.16 -> -15488.48
## -15494.16 -> -15493.71
## -15490.43 -> -15489.36
## -15486.46 -> -15484.21
## -15481.41 -> -15426.27
## -15492.07 -> -15491.92
## -15490.56 -> -15488.99
## optimize edge weights: -15412.1 --> -15411.44
## optimize topology: -15494.16 --> -15411.44
## -15410.43 -> -15410.01
## -15408.18 -> -15399.94
## -15404.78 -> -15389.71
## -15412.37 -> -15411.23
## optimize edge weights: -15388.07 --> -15387.7
## optimize topology: -15411.44 --> -15387.7
## -15385.14 -> -15381.48
## -15379.04 -> -15372.08
## optimize edge weights: -15372.08 --> -15371.49
## optimize topology: -15387.7 --> -15371.49
## -15370.66 -> -15370.18
## -15367.87 -> -15346.21
## optimize edge weights: -15346.21 --> -15346.1
## optimize topology: -15371.49 --> -15346.1
## -15345.6 -> -15345.18
## -15346.13 -> -15345.81
## optimize edge weights: -15344.88 --> -15344.8
## optimize topology: -15346.1 --> -15344.8
## optimize edge weights: -15344.8 --> -15344.8
## optimize topology: -15344.8 --> -15344.8
## NNI moves: 17
## [1] "Ratchet iteration 10 , best pscore so far: -15344.7896495107"
## -16015.96 -> -15991.43
## -16018.88 -> -15994.2
## -16014.5 -> -16003.58
## -16017.45 -> -16004.51
## -15965.21 -> -15964.63
## -16013.74 -> -15997.62
## -16014.2 -> -15981.74
## -15949.97 -> -15942.48
## -16023.35 -> -15984.97
## -15997.84 -> -15977.28
## -16007.32 -> -15987.07
## -16008.6 -> -15984.39
## -15963.56 -> -15958.46
## -16002.6 -> -15985.65
## -16013.64 -> -15999.83
## -16024.22 -> -15995.75
## -16036.8 -> -15996.1
## -16011.57 -> -15995.07
## -16013.64 -> -15995.54
## -16014.17 -> -15992.38
## -16013.45 -> -15992.38
## -15955.18 -> -15955.17
## -16004.88 -> -15986.04
## -16004.44 -> -15994.97
## -16001.83 -> -15994.36
## -15977.92 -> -15966.22
## -16003.07 -> -15979.89
## -15995.21 -> -15985.17
## -16003.75 -> -15991.97
## -16008.59 -> -16003.4
## -16008.17 -> -16003.37
## -16013.7 -> -15989.41
## -16026.6 -> -15990.46
## -15980.92 -> -15974.85
## -16016.61 -> -15988.06
## -16010.62 -> -15987.61
## -15995.2 -> -15968.36
## -16008.6 -> -15993.43
## -16008.37 -> -15993.94
## -16023.16 -> -15994.47
## -16022.05 -> -15995.8
## -15943.8 -> -15942.89
## -16013.57 -> -15979.5
## -16020.97 -> -16001.22
## -16016.37 -> -16001.22
## -16026.7 -> -15989.45
## -16016 -> -15985.11
## -15991.77 -> -15921.35
## -16001.08 -> -15939.59
## -15989.53 -> -15930.28
## -16009.5 -> -15941.12
## optimize edge weights: -15724.71 --> -15627.97
## optimize topology: -16004.88 --> -15627.97
## -15627.97 -> -15624.82
## -15605.39 -> -15590.04
## -15616.81 -> -15590.31
## -15627.97 -> -15561.69
## -15627.97 -> -15625.68
## -15627.97 -> -15627.9
## -15627.97 -> -15627.41
## -15627.48 -> -15618.43
## -15627.97 -> -15620.67
## -15627.97 -> -15627.35
## -15627.75 -> -15623.14
## -15627.51 -> -15589.58
## -15627.97 -> -15622.43
## -15627.97 -> -15626.18
## -15623.14 -> -15619.5
## -15613.36 -> -15591.97
## optimize edge weights: -15413.95 --> -15412.09
## optimize topology: -15627.97 --> -15412.09
## -15412.09 -> -15411.56
## -15412.09 -> -15407.7
## -15411.44 -> -15402.06
## -15412.09 -> -15405.31
## -15412.09 -> -15412.02
## -15412.06 -> -15406.68
## -15412.08 -> -15410.5
## -15408.58 -> -15403.34
## -15409.69 -> -15408.67
## optimize edge weights: -15392.96 --> -15392.34
## optimize topology: -15412.09 --> -15392.34
## -15392.34 -> -15391.54
## -15392.34 -> -15388.22
## -15392.34 -> -15370.11
## optimize edge weights: -15370.11 --> -15369.86
## optimize topology: -15392.34 --> -15369.86
## -15369.86 -> -15368.98
## -15369.86 -> -15367
## -15369.52 -> -15368.97
## optimize edge weights: -15366.12 --> -15366.03
## optimize topology: -15369.86 --> -15366.03
## optimize edge weights: -15366.03 --> -15366.03
## optimize topology: -15366.03 --> -15366.03
## NNI moves: 17
## [1] "Ratchet iteration 11 , best pscore so far: -15344.7896495107"
## -15723.75 -> -15722.89
## -15726.78 -> -15725.28
## -15733.6 -> -15728.98
## -15733.09 -> -15727.89
## -15732.72 -> -15726.59
## -15729.98 -> -15725.05
## -15745.83 -> -15720.95
## -15730.73 -> -15714.3
## -15705.45 -> -15699.17
## -15749.87 -> -15716.7
## -15702.37 -> -15688.36
## -15733.61 -> -15716.79
## -15722.79 -> -15721.28
## -15720.66 -> -15718.02
## -15736.99 -> -15729.98
## -15731.21 -> -15727.81
## -15735.09 -> -15729.96
## -15734.03 -> -15729.86
## -15723.44 -> -15715.4
## -15740.38 -> -15723.94
## -15731.78 -> -15729.57
## -15712.44 -> -15707.86
## -15726.34 -> -15718.74
## -15729.95 -> -15729.58
## -15731.94 -> -15730.85
## -15716.46 -> -15716.46
## -15731.24 -> -15725.46
## -15775.96 -> -15711.69
## -15756.44 -> -15707.38
## -15721.95 -> -15649.47
## -15734.9 -> -15668.97
## -15709.58 -> -15661.1
## -15734.02 -> -15669.38
## -15701.77 -> -15701.34
## -15725.15 -> -15718.63
## optimize edge weights: -15557.91 --> -15511.38
## optimize topology: -15731.63 --> -15511.38
## -15510.5 -> -15510.43
## -15498.12 -> -15478.36
## -15511.38 -> -15509.8
## -15511.38 -> -15506.56
## -15511.38 -> -15510.81
## -15511.38 -> -15511.3
## -15501.92 -> -15497.97
## -15510.18 -> -15509.77
## -15511.28 -> -15511.11
## -15511.38 -> -15444.34
## -15511.38 -> -15508.97
## -15504.38 -> -15492.01
## -15510.84 -> -15510.44
## optimize edge weights: -15392.54 --> -15390.28
## optimize topology: -15511.38 --> -15390.28
## -15380.5 -> -15375.06
## -15389.61 -> -15389.13
## -15390.27 -> -15389.97
## -15390.35 -> -15389.44
## -15390.38 -> -15389.64
## -15386.05 -> -15373.92
## -15389.18 -> -15387.09
## -15389.16 -> -15388.75
## -15390.34 -> -15390.1
## optimize edge weights: -15357.23 --> -15356.7
## optimize topology: -15390.28 --> -15356.7
## -15354.66 -> -15354.51
## -15356.05 -> -15355.53
## -15356.32 -> -15355.55
## -15356.74 -> -15355.77
## -15357.07 -> -15356.33
## optimize edge weights: -15354.51 --> -15354.51
## optimize topology: -15356.7 --> -15354.51
## -15353.79 -> -15353.22
## -15354.09 -> -15353.25
## -15354.48 -> -15353.54
## -15354.81 -> -15354.08
## optimize edge weights: -15353.22 --> -15353.18
## optimize topology: -15354.51 --> -15353.18
## optimize edge weights: -15353.18 --> -15353.18
## optimize topology: -15353.18 --> -15353.18
## NNI moves: 15
## [1] "Ratchet iteration 12 , best pscore so far: -15344.7896495107"
## -16735.35 -> -16730.84
## -16733.31 -> -16728.02
## -16711.54 -> -16711.46
## -16730.4 -> -16725.74
## -16730.04 -> -16728.52
## -16732 -> -16730.13
## -16711.31 -> -16702.08
## -16731.4 -> -16720.19
## -16729.6 -> -16721.09
## -16728.42 -> -16720.49
## -16729.14 -> -16722.84
## -16714.17 -> -16712.83
## -16734.76 -> -16731.32
## -16734.02 -> -16731.14
## -16735.09 -> -16733.02
## -16751.47 -> -16706.12
## -16740.66 -> -16702.59
## -16767.41 -> -16657.21
## -16773.16 -> -16657.21
## -16395.83 -> -16392.66
## -16717.95 -> -16576.43
## -16740.83 -> -16634.9
## -16745.7 -> -16634.98
## -16675.78 -> -16672.08
## -16733.33 -> -16714.36
## -16727.05 -> -16720.52
## -16729.65 -> -16721.83
## -16742.54 -> -16719.16
## -16741.79 -> -16719.16
## -16679.97 -> -16679.79
## -16733.15 -> -16712.82
## -16734.38 -> -16718.27
## -16738.02 -> -16719.31
## -16771.57 -> -16722.58
## -16767.08 -> -16719.9
## -16621.83 -> -16548.24
## -16818.68 -> -16592.51
## -16603.63 -> -16544.43
## -16766.78 -> -16595.56
## -16480.71 -> -16251.02
## -16733.43 -> -16530.58
## -16750.75 -> -16562.24
## -16808.82 -> -16562.69
## -16736.95 -> -16727.39
## -16734.6 -> -16727.06
## optimize edge weights: -15744.74 --> -15691.72
## optimize topology: -16733.15 --> -15691.72
## -15691.72 -> -15676.73
## -15691.72 -> -15684.75
## -15691.25 -> -15690.78
## -15686.33 -> -15679.65
## -15691.66 -> -15691.49
## -15688.3 -> -15685.47
## -15688.28 -> -15686.5
## -15686.53 -> -15682.1
## -15691.41 -> -15690.95
## -15691.72 -> -15689.65
## -15692.03 -> -15691.22
## -15693.43 -> -15691.29
## -15691.72 -> -15513.52
## -15691.72 -> -15691.72
## optimize edge weights: -15477.23 --> -15449.48
## optimize topology: -15691.72 --> -15449.48
## -15445.46 -> -15443.79
## -15449.48 -> -15369.25
## -15449.44 -> -15449.25
## -15449.39 -> -15449.16
## -15447.94 -> -15447.75
## optimize edge weights: -15367.52 --> -15366.51
## optimize topology: -15449.48 --> -15366.51
## -15357.92 -> -15353.8
## -15366.58 -> -15366.47
## -15361.78 -> -15361.03
## optimize edge weights: -15348.26 --> -15347.32
## optimize topology: -15366.51 --> -15347.32
## -15347.32 -> -15346.31
## -15346.78 -> -15345.95
## optimize edge weights: -15344.96 --> -15344.79
## optimize topology: -15347.32 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 18
## [1] "Ratchet iteration 13 , best pscore so far: -15344.7896495107"
## -15835.19 -> -15832.98
## -15837.29 -> -15825.38
## -15843.13 -> -15828.97
## -15764.78 -> -15759.31
## -15837.52 -> -15801.41
## -15865.87 -> -15829.54
## -15863.38 -> -15829.43
## -15813.71 -> -15802.48
## -15827.27 -> -15813.41
## -15860.41 -> -15822.37
## -15860.12 -> -15822.37
## -15753.75 -> -15750.45
## -15835.3 -> -15805
## -15805.72 -> -15802.45
## -15836.24 -> -15826.82
## -15864.28 -> -15834.39
## -15868.56 -> -15834.76
## -15834.66 -> -15815.49
## -15798.45 -> -15793.22
## -15839 -> -15834.68
## -15836.29 -> -15834.68
## -15835.87 -> -15822.23
## -15838.19 -> -15822.23
## -15833.08 -> -15827.4
## -15845.96 -> -15834.84
## -15833.42 -> -15825.44
## -15842.53 -> -15826.4
## -15837.56 -> -15831.1
## -15839.38 -> -15831.59
## -15833.54 -> -15833.23
## -15814.29 -> -15784.56
## -15836.28 -> -15796.77
## -15847.56 -> -15790.22
## -15749.95 -> -15739.08
## -15851.6 -> -15834.43
## -15848.72 -> -15834.4
## -15806.19 -> -15793.17
## -15835.31 -> -15811.68
## -15835.51 -> -15813.69
## -15811.41 -> -15805.57
## optimize edge weights: -15572.05 --> -15524.27
## optimize topology: -15835.65 --> -15524.27
## -15515.07 -> -15488.71
## -15524.11 -> -15524.11
## -15516.15 -> -15510.68
## -15521.57 -> -15520.88
## -15517.22 -> -15509.47
## -15524.2 -> -15524.13
## -15519.13 -> -15510.41
## -15524.19 -> -15524.15
## -15524.5 -> -15523.7
## -15521.97 -> -15521.57
## -15524.27 -> -15523.94
## -15524.27 -> -15509.98
## -15524.27 -> -15489.73
## -15524.27 -> -15524.27
## optimize edge weights: -15425.34 --> -15424.03
## optimize topology: -15524.27 --> -15424.03
## -15415.37 -> -15414.56
## -15413.73 -> -15408
## -15417.12 -> -15409.86
## -15424.22 -> -15423.94
## -15424.02 -> -15423.55
## -15423.56 -> -15421.77
## -15424.03 -> -15414.88
## -15424.03 -> -15424.03
## -15420.92 -> -15420.2
## optimize edge weights: -15370.62 --> -15370.18
## optimize topology: -15424.03 --> -15370.18
## -15370.24 -> -15369.92
## -15370.17 -> -15369.8
## -15370.18 -> -15370.18
## -15370.18 -> -15368.52
## -15370.2 -> -15370.01
## optimize edge weights: -15368.45 --> -15368.29
## optimize topology: -15370.18 --> -15368.29
## -15368.29 -> -15345.6
## -15368.29 -> -15368.24
## optimize edge weights: -15345.6 --> -15345.16
## optimize topology: -15368.29 --> -15345.16
## -15345.18 -> -15344.86
## optimize edge weights: -15344.86 --> -15344.79
## optimize topology: -15345.16 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 17
## [1] "Ratchet iteration 14 , best pscore so far: -15344.7896495107"
## -15719.91 -> -15715.1
## -15719.95 -> -15714.37
## -15714.3 -> -15712.63
## -15713.96 -> -15712.27
## -15704.42 -> -15699.67
## -15700.7 -> -15696.12
## -15714.7 -> -15714.68
## -15717.67 -> -15716.26
## -15716.88 -> -15712.79
## -15719.92 -> -15713.9
## -15726.02 -> -15703.96
## -15725.28 -> -15703.96
## -15663.36 -> -15663.36
## -15716.63 -> -15697.61
## -15699.7 -> -15698.74
## -15712.26 -> -15704.84
## -15719.17 -> -15712.4
## -15714.76 -> -15710.3
## -15717.77 -> -15715.79
## -15692.25 -> -15682.93
## -15713.69 -> -15701.8
## -15711.56 -> -15702.41
## -15711.44 -> -15702.56
## -15716.37 -> -15714.37
## -15735.49 -> -15705.13
## -15731.36 -> -15705.81
## -15648.26 -> -15647.86
## -15722.87 -> -15688.26
## -15730.35 -> -15710.07
## -15723.05 -> -15709.79
## -15740.45 -> -15709.91
## -15728.3 -> -15706
## -15581.1 -> -15581.07
## -15734.92 -> -15646.21
## -15725.53 -> -15670.11
## -15725.22 -> -15669.46
## optimize edge weights: -15401.53 --> -15396.26
## optimize topology: -15716.63 --> -15396.26
## -15394.15 -> -15393.13
## -15390.96 -> -15389.1
## -15387.09 -> -15384.47
## -15395.09 -> -15394.58
## -15396.25 -> -15396.02
## -15396.26 -> -15379.57
## -15396.26 -> -15389.69
## -15395.49 -> -15394.6
## optimize edge weights: -15367.77 --> -15367.39
## optimize topology: -15396.26 --> -15367.39
## -15365.33 -> -15364.45
## -15366.92 -> -15366.16
## -15366.32 -> -15365.86
## -15367.61 -> -15367.29
## optimize edge weights: -15361.92 --> -15361.69
## optimize topology: -15367.39 --> -15361.69
## -15359.08 -> -15356.17
## optimize edge weights: -15356.17 --> -15356.03
## optimize topology: -15361.69 --> -15356.03
## -15352.95 -> -15348.8
## -15354.69 -> -15352.66
## -15353.83 -> -15353.45
## optimize edge weights: -15348.8 --> -15348.57
## optimize topology: -15356.03 --> -15348.57
## -15345.47 -> -15344.83
## optimize edge weights: -15344.83 --> -15344.79
## optimize topology: -15348.57 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 14
## [1] "Ratchet iteration 15 , best pscore so far: -15344.7896495107"
## -15839.83 -> -15830.69
## -15841.55 -> -15830.35
## -15828.33 -> -15806.82
## -15831.61 -> -15809.29
## -15788.3 -> -15775.14
## -15825.37 -> -15804.46
## -15819.42 -> -15798.72
## -15847.95 -> -15805.76
## -15810.83 -> -15805.48
## -15835.54 -> -15820.5
## -15813.55 -> -15813.54
## -15830.94 -> -15826.56
## -15835.57 -> -15834.16
## -15835.71 -> -15834.27
## -15831.06 -> -15830.89
## -15836.2 -> -15831.52
## -15834.47 -> -15831.34
## -15851.87 -> -15828.94
## -15859.77 -> -15828.94
## -15838.63 -> -15830.83
## -15838.65 -> -15830.83
## -15824.74 -> -15824.74
## -15833.5 -> -15824.75
## -15809.39 -> -15798.92
## -15836.61 -> -15823.91
## -15828.29 -> -15798.92
## -15848.55 -> -15808.58
## -15861.37 -> -15781.98
## -15928.7 -> -15793.3
## -15801.87 -> -15724.63
## -15836.31 -> -15755.59
## -15819.38 -> -15782.17
## -15833.96 -> -15794.79
## -15845.07 -> -15796.38
## -15871.81 -> -15797.57
## -15832.32 -> -15831.39
## -15835.23 -> -15832.88
## optimize edge weights: -15593.35 --> -15556.49
## optimize topology: -15834.51 --> -15556.49
## -15552.1 -> -15533.84
## -15554.44 -> -15547.9
## -15550.13 -> -15539.97
## -15556.47 -> -15556.42
## -15554.53 -> -15527.39
## -15551.84 -> -15549.87
## -15556.28 -> -15555.58
## -15549.91 -> -15546.56
## -15556.49 -> -15468.92
## optimize edge weights: -15417.17 --> -15414.38
## optimize topology: -15556.49 --> -15414.38
## -15401.6 -> -15381.64
## -15406.89 -> -15397.99
## -15405.3 -> -15400.69
## -15414.07 -> -15413.51
## -15411.43 -> -15410.51
## optimize edge weights: -15363.95 --> -15363.08
## optimize topology: -15414.38 --> -15363.08
## -15355.63 -> -15345.81
## -15362.59 -> -15362.13
## optimize edge weights: -15344.88 --> -15344.79
## optimize topology: -15363.08 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 14
## [1] "Ratchet iteration 16 , best pscore so far: -15344.7896495107"
## -15720.91 -> -15698.59
## -15696.76 -> -15687.12
## -15700.29 -> -15692.83
## -15713.06 -> -15693.79
## -15733.33 -> -15712.34
## -15676.87 -> -15663.08
## -15713.04 -> -15692.37
## -15724.82 -> -15704.39
## -15725.07 -> -15704.39
## -15720.56 -> -15713.81
## -15715.64 -> -15714.62
## -15697.91 -> -15689.67
## -15717.33 -> -15709.68
## -15705.2 -> -15682.95
## -15721.44 -> -15701.41
## -15738.51 -> -15703.89
## -15754.19 -> -15705.64
## -15711.57 -> -15698.76
## -15729.01 -> -15702.59
## -15710.54 -> -15710.29
## -15717.43 -> -15703.5
## -15725.07 -> -15703.83
## -15705.27 -> -15695.51
## -15706.01 -> -15689.9
## -15723.33 -> -15712.53
## -15739.14 -> -15713.27
## -15717.52 -> -15709.84
## -15723.17 -> -15710.76
## -15722.69 -> -15708.82
## -15718.48 -> -15707.08
## -15684.54 -> -15684.27
## -15708.78 -> -15699.81
## -15707.99 -> -15707.82
## -15708.85 -> -15707.27
## -15713.11 -> -15710.09
## -15717.55 -> -15713.63
## -15694.57 -> -15694.45
## -15718.57 -> -15711.86
## optimize edge weights: -15554.51 --> -15499.81
## optimize topology: -15714.77 --> -15499.81
## -15493.35 -> -15477.12
## -15499.81 -> -15468.95
## -15514.6 -> -15496.59
## -15497.39 -> -15495.61
## -15499.84 -> -15497.1
## -15499.4 -> -15499.06
## -15499.66 -> -15499.55
## -15499.81 -> -15498.98
## -15499.81 -> -15491.86
## -15492.06 -> -15491.56
## -15499.81 -> -15495.94
## -15499.81 -> -15498.64
## -15499.12 -> -15499.02
## -15497.59 -> -15497.42
## optimize edge weights: -15431.95 --> -15429.8
## optimize topology: -15499.81 --> -15429.8
## -15429.8 -> -15428.06
## -15429.8 -> -15429.65
## -15429.8 -> -15429.13
## -15429.8 -> -15429.78
## -15429.8 -> -15428.8
## -15429.8 -> -15422.24
## -15430.07 -> -15429.73
## -15429.8 -> -15425.38
## -15429.8 -> -15428.45
## -15426.79 -> -15425.08
## optimize edge weights: -15416.57 --> -15416.4
## optimize topology: -15429.8 --> -15416.4
## -15416.4 -> -15390.52
## -15416.4 -> -15416.4
## -15417.27 -> -15412.23
## -15416.4 -> -15411.91
## -15416.4 -> -15415.73
## optimize edge weights: -15387.02 --> -15386.46
## optimize topology: -15416.4 --> -15386.46
## -15386.7 -> -15381.28
## -15386.44 -> -15386.26
## -15386.46 -> -15350.51
## -15386.46 -> -15385.09
## optimize edge weights: -15350.51 --> -15349.95
## optimize topology: -15386.46 --> -15349.95
## -15350 -> -15344.81
## optimize edge weights: -15344.81 --> -15344.79
## optimize topology: -15349.95 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 17
## [1] "Ratchet iteration 17 , best pscore so far: -15344.7896495107"
## -16559.69 -> -16555.63
## -16562.47 -> -16557.98
## -16563.02 -> -16559.23
## -16561.69 -> -16557.81
## -16538.15 -> -16538.15
## -16559.65 -> -16554.28
## -16563.4 -> -16553.2
## -16570.43 -> -16556.26
## -16501.09 -> -16500.4
## -16571.57 -> -16542.06
## -16531.65 -> -16529
## -16556.84 -> -16548.69
## -16553.81 -> -16549.84
## -16554.89 -> -16547.02
## -16564.05 -> -16558.99
## -16570.38 -> -16559.18
## -16554.71 -> -16551.34
## -16555.31 -> -16553.23
## -16571.89 -> -16557.14
## -16573.88 -> -16557.62
## -16579.61 -> -16544.42
## -16577.47 -> -16544.49
## -16599.49 -> -16508.53
## -16607.57 -> -16508.53
## -16308.74 -> -16308.47
## -16566.17 -> -16450.66
## -16561.81 -> -16558.34
## -16560.12 -> -16557.87
## -16510.93 -> -16502.08
## -16560.92 -> -16532.05
## -16554.4 -> -16544.79
## -16559.96 -> -16548.11
## -16559.02 -> -16552.25
## -16563.6 -> -16553.49
## -16521.01 -> -16520.93
## -16563.25 -> -16549.55
## -16595.07 -> -16549.99
## -16595.95 -> -16549.99
## -16609.64 -> -16478.31
## -16608.34 -> -16478.96
## -16380.03 -> -16364.47
## -16568.7 -> -16434.83
## -16640.86 -> -16476.6
## -16666.39 -> -16476.64
## -16060.24 -> -16015.93
## -16556.27 -> -16317.9
## -16616.07 -> -16463.58
## -16626.19 -> -16463.58
## -16556.42 -> -16556.06
## -16558.86 -> -16557.72
## optimize edge weights: -15607.49 --> -15528.72
## optimize topology: -16559.63 --> -15528.72
## -15527.15 -> -15522.17
## -15528.72 -> -15513.71
## -15528.72 -> -15527.92
## -15528.72 -> -15528.28
## -15528.72 -> -15510.33
## -15528.72 -> -15527.57
## -15528.72 -> -15528.72
## -15528.72 -> -15523.52
## -15528.72 -> -15527.11
## -15529.38 -> -15528.63
## -15527.16 -> -15526.65
## -15528.72 -> -15528.72
## -15528.72 -> -15527.12
## -15528.72 -> -15528.72
## -15518.39 -> -15500.48
## optimize edge weights: -15460.45 --> -15459.86
## optimize topology: -15528.72 --> -15459.86
## -15453.72 -> -15372.72
## -15459.86 -> -15458.61
## -15459.86 -> -15459.86
## -15459.86 -> -15454.99
## -15459.86 -> -15458.26
## -15459.72 -> -15459.55
## -15458.62 -> -15458.15
## -15459.86 -> -15458.27
## -15459.86 -> -15459.86
## optimize edge weights: -15365.02 --> -15363.52
## optimize topology: -15459.86 --> -15363.52
## -15363.64 -> -15363.36
## -15362.39 -> -15361.91
## -15354.98 -> -15351.13
## optimize edge weights: -15351.13 --> -15350.87
## optimize topology: -15363.52 --> -15350.87
## -15350.92 -> -15350.68
## -15349.75 -> -15349.29
## -15351.69 -> -15350.56
## -15348.43 -> -15347.11
## optimize edge weights: -15347.11 --> -15346.96
## optimize topology: -15350.87 --> -15346.96
## -15347 -> -15346.77
## -15345.85 -> -15345.4
## -15347.61 -> -15346.42
## optimize edge weights: -15344.85 --> -15344.8
## optimize topology: -15346.96 --> -15344.8
## optimize edge weights: -15344.8 --> -15344.8
## optimize topology: -15344.8 --> -15344.8
## NNI moves: 21
## [1] "Ratchet iteration 18 , best pscore so far: -15344.7896495107"
## -15839.96 -> -15826.74
## -15840.41 -> -15826.74
## -15798.3 -> -15797.06
## -15817.45 -> -15810.59
## -15841.96 -> -15812.88
## -15844.03 -> -15812.88
## -15760.14 -> -15759.38
## -15820.05 -> -15802.36
## -15807.61 -> -15803.33
## -15832.49 -> -15820.64
## -15836.7 -> -15832.4
## -15833.35 -> -15828.97
## -15814.87 -> -15814.83
## -15836.85 -> -15830.7
## -15831.12 -> -15824.63
## -15832.73 -> -15825.03
## -15833.29 -> -15832.82
## -15824.58 -> -15824.58
## -15826.09 -> -15824.98
## -15831.79 -> -15831.45
## -15833.4 -> -15832.35
## -15802.02 -> -15800.76
## -15832.66 -> -15818.99
## -15832.35 -> -15820.5
## -15834.35 -> -15821.27
## -15862.95 -> -15823.35
## -15858.43 -> -15823.35
## -15827.13 -> -15759.09
## -15836.38 -> -15768.45
## -15679.5 -> -15678.76
## -15838.21 -> -15740.8
## -15806.74 -> -15786.77
## -15845.03 -> -15798.19
## -15837.12 -> -15805.45
## -15847.94 -> -15805.63
## -15773.89 -> -15762.73
## -15823.87 -> -15798.03
## -15843.1 -> -15814.13
## -15844.54 -> -15814.11
## optimize edge weights: -15485.4 --> -15443.48
## optimize topology: -15833.31 --> -15443.48
## -15442.51 -> -15442.1
## -15443.48 -> -15424.31
## -15443.48 -> -15441.71
## -15433.32 -> -15426.9
## -15443.48 -> -15428.74
## -15443.48 -> -15442.75
## -15442.43 -> -15442.02
## -15443.74 -> -15443.4
## -15443.48 -> -15441.84
## -15440.72 -> -15439.79
## -15439.36 -> -15429.77
## -15437.83 -> -15421.59
## optimize edge weights: -15364.51 --> -15363.83
## optimize topology: -15443.48 --> -15363.83
## -15362.45 -> -15362.02
## -15362.95 -> -15362.94
## -15362.98 -> -15362.59
## -15362.39 -> -15360.5
## -15361.59 -> -15359.17
## optimize edge weights: -15356.96 --> -15356.55
## optimize topology: -15363.83 --> -15356.55
## -15356.37 -> -15351.09
## -15354.77 -> -15354.46
## -15355.16 -> -15353.95
## -15355.57 -> -15352.89
## optimize edge weights: -15347.63 --> -15347.5
## optimize topology: -15356.55 --> -15347.5
## -15345.94 -> -15345.21
## -15346.5 -> -15344.86
## optimize edge weights: -15344.86 --> -15344.79
## optimize topology: -15347.5 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 18
## [1] "Ratchet iteration 19 , best pscore so far: -15344.7896495107"
## -15838.21 -> -15836.57
## -15838.78 -> -15837.03
## -15836.31 -> -15826.34
## -15853.35 -> -15829.39
## -15862.53 -> -15809.25
## -15869.39 -> -15812.85
## -15888.44 -> -15809.89
## -15899.85 -> -15809.89
## -15600.83 -> -15600.82
## -15860.48 -> -15741.58
## -15816.11 -> -15810.45
## -15825.46 -> -15816.94
## -15843.68 -> -15834.72
## -15844.03 -> -15834.86
## -15815.67 -> -15807.41
## -15840.97 -> -15821.82
## -15847.45 -> -15825.04
## -15846.72 -> -15825.04
## -15788.92 -> -15788.85
## -15838.15 -> -15818.7
## -15837.54 -> -15834.03
## -15841.35 -> -15835.29
## -15812.28 -> -15799.53
## -15842.94 -> -15821.49
## -15828.33 -> -15820.7
## -15829.75 -> -15824.35
## -15839.15 -> -15837.92
## -15838.19 -> -15837.19
## -15815.45 -> -15815.44
## -15834.34 -> -15829.86
## -15837.72 -> -15828.18
## -15842.28 -> -15830.81
## -15808.13 -> -15806.78
## -15833.5 -> -15826.49
## -15818.81 -> -15818.59
## -15846.05 -> -15834.8
## -15832.96 -> -15832.94
## optimize edge weights: -15436.5 --> -15413.79
## optimize topology: -15838.15 --> -15413.79
## -15413.7 -> -15413.62
## -15414.85 -> -15413.73
## -15413.79 -> -15399.77
## -15413.79 -> -15411.04
## -15413.77 -> -15413.76
## -15407.89 -> -15388.51
## -15413.79 -> -15405.16
## -15413.79 -> -15410.14
## -15412.63 -> -15411.96
## -15412.31 -> -15410.31
## -15413.79 -> -15412.38
## -15413.79 -> -15413.45
## optimize edge weights: -15368.83 --> -15368.09
## optimize topology: -15413.79 --> -15368.09
## -15367.64 -> -15367.38
## -15359.47 -> -15354.86
## -15364.28 -> -15360.95
## -15368.09 -> -15359.24
## -15368.09 -> -15363.92
## -15368.15 -> -15367.03
## optimize edge weights: -15354.06 --> -15353.35
## optimize topology: -15368.09 --> -15353.35
## -15351.28 -> -15344.81
## -15350.8 -> -15349.98
## optimize edge weights: -15344.81 --> -15344.79
## optimize topology: -15353.35 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 14
## [1] "Ratchet iteration 20 , best pscore so far: -15344.7896495107"
## -15636.97 -> -15634.23
## -15641.55 -> -15635.78
## -15641.34 -> -15632.48
## -15642.52 -> -15632.55
## -15618.74 -> -15618.67
## -15635.17 -> -15630.03
## -15638.4 -> -15637.84
## -15639.69 -> -15639.06
## -15642.05 -> -15638.43
## -15650.7 -> -15622.79
## -15649.39 -> -15622.79
## -15599.96 -> -15591.42
## -15639.68 -> -15622.68
## -15617.86 -> -15606.14
## -15628 -> -15617.19
## -15633.57 -> -15625.19
## -15642.15 -> -15634.98
## -15638.1 -> -15623.08
## -15632.1 -> -15621.52
## -15638.83 -> -15631.93
## -15630.48 -> -15630.46
## -15634.18 -> -15633.14
## -15641.04 -> -15626.61
## -15629.95 -> -15619.06
## -15646.94 -> -15628.91
## -15647.33 -> -15627.2
## -15615.25 -> -15611.91
## -15634.91 -> -15627.5
## -15601.7 -> -15595.55
## -15652.28 -> -15631.5
## optimize edge weights: -15503.49 --> -15469.37
## optimize topology: -15639.68 --> -15469.37
## -15469.37 -> -15468.13
## -15469.27 -> -15468.86
## -15469.19 -> -15468.79
## -15468.9 -> -15468.24
## -15465.93 -> -15460.99
## -15469.37 -> -15454.11
## -15469.37 -> -15468.33
## -15469.09 -> -15468.94
## -15458.09 -> -15426.48
## -15469.37 -> -15463.96
## -15469.37 -> -15466.68
## -15469.37 -> -15464.73
## -15459.17 -> -15453.26
## optimize edge weights: -15395.2 --> -15394.45
## optimize topology: -15469.37 --> -15394.45
## -15394.45 -> -15392.88
## -15394.36 -> -15394.32
## -15394.01 -> -15393.7
## -15393.89 -> -15393.09
## -15394.45 -> -15385.45
## -15390.55 -> -15370.76
## -15393.96 -> -15393.02
## -15394.45 -> -15384.23
## -15394.45 -> -15393.18
## -15394.45 -> -15391.66
## -15394.45 -> -15389.75
## optimize edge weights: -15360.41 --> -15360.28
## optimize topology: -15394.45 --> -15360.28
## -15360.28 -> -15359.01
## -15360.28 -> -15350.16
## -15360.28 -> -15358.94
## -15360.28 -> -15357.16
## -15360.28 -> -15355.75
## -15360.11 -> -15360.01
## -15359.87 -> -15359.74
## optimize edge weights: -15349.24 --> -15348.93
## optimize topology: -15360.28 --> -15348.93
## -15348.93 -> -15346.65
## -15348.93 -> -15345.57
## -15348.93 -> -15348.66
## -15348.93 -> -15348.3
## optimize edge weights: -15344.91 --> -15344.79
## optimize topology: -15348.93 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 15
## [1] "Ratchet iteration 21 , best pscore so far: -15344.7896495107"
## -15731.71 -> -15726.85
## -15740.65 -> -15729.81
## -15716.68 -> -15716.67
## -15739.71 -> -15725.76
## -15736.36 -> -15728.34
## -15734.75 -> -15728.86
## -15734.04 -> -15731.7
## -15736.44 -> -15732.87
## -15734.81 -> -15731.06
## -15711.22 -> -15711.22
## -15733.64 -> -15727.93
## -15738.32 -> -15728.28
## -15745.15 -> -15730.6
## -15678.87 -> -15678.86
## -15743.9 -> -15718.06
## -15730.48 -> -15728.06
## -15735.02 -> -15732.97
## -15719.72 -> -15715.5
## -15730.74 -> -15728.93
## -15732.23 -> -15731.11
## -15723.2 -> -15721.41
## -15736.28 -> -15728.94
## -15731.19 -> -15724.77
## -15736.63 -> -15728
## -15720.12 -> -15715.21
## -15739.86 -> -15723.68
## -15739.72 -> -15722.74
## -15739.83 -> -15722.74
## -15691.68 -> -15691.27
## -15733.12 -> -15716.39
## -15737.24 -> -15731.17
## -15734.47 -> -15730.25
## -15713.01 -> -15696.71
## -15718.49 -> -15695.75
## -15724.56 -> -15662.66
## -15721.47 -> -15663.63
## -15730.5 -> -15668.24
## -15731.2 -> -15676.63
## -15729.53 -> -15729.31
## -15729.91 -> -15729.71
## optimize edge weights: -15504.89 --> -15462.45
## optimize topology: -15733.14 --> -15462.45
## -15461.39 -> -15459.41
## -15461.39 -> -15458.79
## -15461.62 -> -15460.6
## -15460.89 -> -15460.46
## -15460.87 -> -15460.25
## -15462.45 -> -15458.87
## -15458.34 -> -15450.56
## -15462.45 -> -15455.67
## -15462.45 -> -15457.13
## -15460.35 -> -15398.62
## -15462.92 -> -15460.98
## optimize edge weights: -15382.66 --> -15382.11
## optimize topology: -15462.45 --> -15382.11
## -15382.99 -> -15382.1
## -15379.55 -> -15373.42
## -15380.86 -> -15379.35
## -15379.48 -> -15377.7
## -15380.58 -> -15379.99
## -15382.11 -> -15380.11
## -15382.11 -> -15382.11
## -15382.11 -> -15371.25
## -15382.11 -> -15375.66
## -15379.81 -> -15378.97
## optimize edge weights: -15360.92 --> -15360.45
## optimize topology: -15382.11 --> -15360.45
## -15355.07 -> -15350.64
## -15360.46 -> -15360.39
## -15360.45 -> -15358.38
## -15360.45 -> -15360.45
## -15357.85 -> -15356.99
## optimize edge weights: -15347.17 --> -15346.95
## optimize topology: -15360.45 --> -15346.95
## -15346.56 -> -15346.5
## -15346.95 -> -15345.32
## -15346.95 -> -15346.81
## optimize edge weights: -15344.83 --> -15344.79
## optimize topology: -15346.95 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 18
## [1] "Ratchet iteration 22 , best pscore so far: -15344.7896495107"
## -16218.65 -> -16212.64
## -16218.67 -> -16213.57
## -16207.68 -> -16198.97
## -16183.12 -> -16178.9
## -16242.09 -> -16201.28
## -16242.74 -> -16201.28
## -16232.23 -> -16218.72
## -16222.02 -> -16193.57
## -16201.12 -> -16189.25
## -16212.95 -> -16199.63
## -16215.66 -> -16204.09
## -16172.98 -> -16164.07
## -16216.97 -> -16191.06
## -16233.66 -> -16211.72
## -16238.38 -> -16211.66
## -16239.86 -> -16211.55
## -16235.17 -> -16210.24
## -16196.07 -> -16180.54
## -16216.68 -> -16204.07
## -16225.59 -> -16203.81
## -16225.53 -> -16203.81
## -16187.08 -> -16180.79
## -16219.11 -> -16197.35
## -16186.05 -> -16170.99
## -16218.94 -> -16197.3
## -16225.41 -> -16208.62
## -16224.53 -> -16208.26
## -16203.82 -> -16202.93
## -16208.25 -> -16204.74
## -16235.46 -> -16192.17
## -16241.71 -> -16194.65
## -16272.24 -> -16198.31
## -16282.63 -> -16198.31
## -15981.47 -> -15981.46
## -16239.6 -> -16120.57
## -16248.34 -> -16187.19
## -16240.25 -> -16185.62
## -16217.35 -> -16115.17
## -16213.53 -> -16128.25
## -16246.68 -> -16150.53
## -16231.68 -> -16142.99
## -16227.12 -> -16167.28
## -16228.23 -> -16167.28
## -16191.82 -> -16114.34
## -16236.13 -> -16148.76
## -16231.68 -> -16189.47
## -16248.43 -> -16191.87
## -16239.48 -> -16204.48
## -16236.34 -> -16203.72
## optimize edge weights: -15758.02 --> -15645.19
## optimize topology: -16218.85 --> -15645.19
## -15644.21 -> -15643.98
## -15644.83 -> -15644.76
## -15633.1 -> -15616.82
## -15635.74 -> -15613.87
## -15643.55 -> -15639.25
## -15644.16 -> -15640.32
## -15644.23 -> -15643.59
## -15645.19 -> -15624.75
## -15645.19 -> -15644.85
## -15645.19 -> -15613.03
## -15639.63 -> -15636.9
## -15642.19 -> -15630.24
## -15644.21 -> -15643.29
## optimize edge weights: -15560.22 --> -15559.08
## optimize topology: -15645.19 --> -15559.08
## -15557.36 -> -15556.91
## -15542.89 -> -15527.09
## -15558.29 -> -15557.84
## -15558.08 -> -15557.38
## -15559.08 -> -15536.31
## -15559.08 -> -15558.89
## -15555.76 -> -15551.19
## -15558.11 -> -15556.41
## -15559.46 -> -15557.42
## optimize edge weights: -15502.95 --> -15500.24
## optimize topology: -15559.08 --> -15500.24
## -15495.87 -> -15494.01
## -15500.24 -> -15486.5
## -15500.24 -> -15497.02
## -15500.82 -> -15499.51
## -15499.25 -> -15498.07
## -15500.31 -> -15500.13
## -15500.22 -> -15500.18
## -15498.21 -> -15469.79
## optimize edge weights: -15456.29 --> -15455.91
## optimize topology: -15500.24 --> -15455.91
## -15451.46 -> -15449.54
## -15455.91 -> -15357.39
## -15455.37 -> -15454.54
## optimize edge weights: -15349.93 --> -15348.44
## optimize topology: -15455.91 --> -15348.44
## -15347.55 -> -15344.91
## optimize edge weights: -15344.91 --> -15344.79
## optimize topology: -15348.44 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 18
## [1] "Ratchet iteration 23 , best pscore so far: -15344.7896495107"
## -15762.13 -> -15752.85
## -15749.67 -> -15744.45
## -15779.41 -> -15763.9
## -15779.36 -> -15763.9
## -15769.46 -> -15755.31
## -15767.03 -> -15755.05
## -15769.4 -> -15751.46
## -15759.75 -> -15740.14
## -15765.91 -> -15741
## -15759.07 -> -15735.87
## -15744.93 -> -15728.25
## -15765.91 -> -15740.47
## -15766.9 -> -15749.96
## -15757.51 -> -15739.35
## -15765.76 -> -15746.25
## -15749.5 -> -15729.76
## -15720.45 -> -15715.65
## -15767.88 -> -15744.08
## -15742.29 -> -15733.34
## -15767.3 -> -15747
## -15780.62 -> -15762.54
## -15782.54 -> -15762.58
## -15776.36 -> -15757.54
## -15780.08 -> -15757.64
## -15740.55 -> -15716.65
## -15767.82 -> -15744.27
## -15779.81 -> -15765.14
## -15796.7 -> -15767.41
## -15781.04 -> -15765.9
## -15782.87 -> -15764.58
## -15763.47 -> -15758.97
## -15764.72 -> -15756.84
## -15761.83 -> -15757.22
## -15766.91 -> -15762.24
## -15755.82 -> -15754.58
## -15769.69 -> -15767.21
## -15769.1 -> -15766.97
## -15778.55 -> -15764.91
## -15774.91 -> -15760.04
## -15708.43 -> -15707.99
## -15771.09 -> -15752.9
## -15763.25 -> -15763.25
## optimize edge weights: -15598.78 --> -15550.45
## optimize topology: -15768.03 --> -15550.45
## -15550.45 -> -15545.38
## -15550.45 -> -15549.37
## -15550.32 -> -15549.87
## -15547.94 -> -15544.78
## -15543.18 -> -15538.98
## -15540.71 -> -15536.09
## -15550.45 -> -15545.03
## -15550.45 -> -15550.45
## -15550.45 -> -15534.34
## -15550.45 -> -15550.45
## -15549.69 -> -15548.02
## -15549.06 -> -15546.62
## -15550.02 -> -15549.11
## -15549.78 -> -15549.54
## -15549.24 -> -15548.96
## -15548.3 -> -15548.24
## -15544.28 -> -15543.81
## optimize edge weights: -15499.43 --> -15493.08
## optimize topology: -15550.45 --> -15493.08
## -15491.69 -> -15484.85
## -15480.85 -> -15457.67
## -15492.22 -> -15490.19
## -15491.45 -> -15489.1
## -15494.39 -> -15492.87
## -15483.34 -> -15477.08
## -15493.06 -> -15492.95
## -15492.33 -> -15492.31
## -15493.08 -> -15456.79
## optimize edge weights: -15405.36 --> -15404.67
## optimize topology: -15493.08 --> -15404.67
## -15404.57 -> -15397.14
## -15395.91 -> -15393.73
## -15403.91 -> -15401.64
## -15403.24 -> -15400.39
## -15405.85 -> -15403.5
## -15402.57 -> -15402.42
## -15404.41 -> -15404.29
## -15404.67 -> -15395.6
## optimize edge weights: -15380.96 --> -15380.73
## optimize topology: -15404.67 --> -15380.73
## -15378.72 -> -15372.51
## -15379.53 -> -15379.22
## -15378.58 -> -15371.96
## -15378.3 -> -15377.7
## -15378.71 -> -15378.58
## optimize edge weights: -15368.38 --> -15368.24
## optimize topology: -15380.73 --> -15368.24
## -15360.99 -> -15345.25
## optimize edge weights: -15345.25 --> -15344.79
## optimize topology: -15368.24 --> -15344.79
## optimize edge weights: -15344.79 --> -15344.79
## optimize topology: -15344.79 --> -15344.79
## NNI moves: 21
## [1] "Ratchet iteration 24 , best pscore so far: -15344.7896495107"
## optimize base frequencies: -15344.79 --> -15344.73
## optimize rate matrix: -15344.73 --> -15344.7
## optimize edge weights: -15344.7 --> -15344.7
## optimize topology: -15344.7 --> -15344.7
## NNI moves: 0
## optimize base frequencies: -15344.7 --> -15344.69
## optimize rate matrix: -15344.69 --> -15344.69
## optimize edge weights: -15344.69 --> -15344.69
## optimize base frequencies: -15344.69 --> -15344.69
## optimize rate matrix: -15344.69 --> -15344.69
## optimize edge weights: -15344.69 --> -15344.69
## optimize base frequencies: -15344.69 --> -15344.69
## optimize rate matrix: -15344.69 --> -15344.69
## optimize edge weights: -15344.69 --> -15344.69
## optimize base frequencies: -15344.69 --> -15344.69
## optimize rate matrix: -15344.69 --> -15344.69
## optimize edge weights: -15344.69 --> -15344.69
## optimize base frequencies: -15344.69 --> -15344.69
## optimize rate matrix: -15344.69 --> -15344.69
## optimize edge weights: -15344.69 --> -15344.69When the analysis finishes, get the loglikelihood of the tree by typing:
logLik(fitGTR)## 'log Lik.' -15344.69 (df=85)Is this log likelihood higher or lower than the log likelihood of the NJ analysis? How does it compare to the log likelihood of the JC analysis?
Look at the tree topology by typing:
plot(fitGTR)
Is this tree different from the one using the JC model?
Save the tree as a pdf if you want. Write the tree to a file by typing:
write.tree(fitGTR$tree,file="fitGTR")Maximum parsimony
Another method for constructing trees is maximum parsimony. In this
case, the MP approach looks to minimize the number of changes in
character states between branches. This contrasts to ML that looks to
find the tree with the highest probability of each character state in
their respective positions in the tree. Let’s see how we can implement
MP in R using phangorn and ape.
Like before, it is best to start with an existing model to optimize.
We can use the NJ tree we made above (mammals_NJ). Now,
let’s determine the parsimony score for this tree.
parsimony(mammals_NJ, mammals_phyDat)## [1] 3165Similar to ML, we can perform an optimization of the tree except optimizing for parsimony.
fitNJ_MP <- optim.parsimony(mammals_NJ, mammals_phyDat, perturbation="stochastic")## Final p-score 3119 after 5 nni operationsLet’s see what we’ve got!
plot(fitNJ_MP)
Boostrapping
Finally, let’s assess the support for a tree using the bootstrap.
Here, we’ll examine the ML tree we produced, above. Run the bootstrap
analysis using the command below. We are setting the number of
iterations to 10. Typically, you would set this value much higher to
produce more accurate branch support estimates. If you are using the
RStudio server, be aware this process will take some time to execute
when you choose a much larger bootstrap value (bs).
bs<-bootstrap.pml(fitGTR, bs=10, optNni=TRUE, multicore=TRUE, control=pml.control(trace=0))Then summarize all the bootstrap trees with:
plotBS((fitGTR$tree),bs,p=50,type="p")
Save this tree as a pdf.
Question 4: Which branches are poorly supported? Which ones are strongly supported? How do the poorly supported branches relate to branches that vary among the NJ and ML trees?
Write the set of bootstrap trees to a file using write.tree:
write.tree(bs,file="bootstrap.GTR.tree")