diff --git a/Disorder_prediction.ipynb b/Disorder_prediction.ipynb index 9823bd4..1b49085 100644 --- a/Disorder_prediction.ipynb +++ b/Disorder_prediction.ipynb @@ -423,59 +423,6 @@ " plt.savefig(\"fig/MobiDB_{}_disorder.svg\".format(prot))\n", " plt.show()" ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "## Tree preparation" - ] - }, - { - "cell_type": "code", - "execution_count": 13, - "metadata": {}, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "/home/malszyck/anaconda3/envs/py3/lib/python3.7/site-packages/rpy2/robjects/pandas2ri.py:17: FutureWarning: pandas.core.index is deprecated and will be removed in a future version. The public classes are available in the top-level namespace.\n", - " from pandas.core.index import Index as PandasIndex\n" - ] - } - ], - "source": [ - "%load_ext rpy2.ipython\n", - "\n", - "taxnames = list(dis_probs.index)" - ] - }, - { - "cell_type": "code", - "execution_count": 15, - "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAMAAABKCk6nAAADAFBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////isF19AAAgAElEQVR4nO2dCXwT1fb4J91Ll7RN972lG5Sl7EtZKluAylJqEYGCyCJIeVDZVzdQEJQCSsWKiPWxCYj+hfIaoc8F8Yc8K248DfAQRNRqBQLS0jb3P/feTJKmSTNJm8wknu/HT9rmzCShXzO9ufeecxgEODWM0C8AsC0g2MkBwU4OCHZyQLCTA4KdHBDs5IBgJwcEOzkg2MkBwU4OCHZyQLCTA4KdHBDs5IBgJwcEOzkg2MkBwU4OCHZyQLCTA4KdHBDsDFS6yGTSYb8ZC4FgZ6AyCqF7Y54yFgLBzgAW/NfQI6goLWVJw82R7dKPa0Mg2DE55Z+opRBVusbFeaZXnex0696wkpICdGG09kAQ7Jh8GndWy8/kHVy/ZtyS5xB6O/dM5PQ9N7UHgmCHoVShR1GUfqgyYnaYV0rIkvUIHchFN56Mjq3lQiDYUfjBtZsebT31Y//xyPnfX+O9lrCXaHnJqrU1V0O1b2EQ7Cj86Kv/0+kU/Z82M0GysJw3Ijf7hC9u6JAR4/O8dqAFgh2FxoKP+eVpKUIrpuP7alyuvDEGXYpuqJDrBlog2DF4ecuzHvo/31y2Q8tptHwGvq/G9YpKdmfjElQhPxPZfthMciAIdgzi+g93NR19J35WmFenxRFqNOFQ769YwWSgRUIg2DHofbTxJboxtb5Rp2+vc52K0PuDMhArmAy0SAgEOwbNCz4SNz3Eq0NhpFrZzStDvTymzZNZMT5zSAgEOwa9j170WGqCw7pBltLlFdX+gbWqpEsV8jJyIgh2CDq6rlaNMCV4v26QpWT/8M7cgdAfNSDYoegV/f+aC78T2cM35qn3ItTKdIRmlCB0QQWCHYr+qUebC6vcBv54KT2gBGHBe+X1VbFVINihMCP4QjA7yEpeViFXps8vvpHo4Yk/C4NgB6J/3PIDJvkCqadETC6+gj//ovnFsFzogPTv3HeISTawB1zfNyd0FxUMy4UOxSBi1ntNswd9+PL0YL8O7bDgWcXoxuE5sFzoMHRahq/CibubPegr90FXVOO8yjuiu0nFq9aiv2C50GHo8z6+7X+k2YPKZF2lIQ+9UpkfE+A74FonH1//gvnRfUgIBIsck4Jvfsht2fkVrcuj932Zq65L+bbM85ubvuvReHIXCBY5JgXP9eM23f0DvTAB39PQgM4WPeZ9umwwQnFT5sI72CEwKXjLZN335W3r2dvRu/d33fdN/9NlY1jB19FzJASCRQ4VnDJsliEZHXQHzQwvqK5d4nIwdg36xh8EOxRU8DOr1xsyuJfuoA/Dx0tlbSfUfJSanj9vCAh2JKhgI+hfotXxHyEU/1GFXL08NelFzYwWDYFgkWNScAHDkY/QigJ0Jk5dIefWgrHgcnIYCBY5JgWjqmoN9xA6H9GweAX7zuXWgrHg1eQoECxyTAvWQ8W4+bt0VW4I5daC8ZTlOhICwSKHp+AX+3ZdNGpDKF0LPkemLM+QEAgWOTwFX3N7sbz9hsBRvj6ppUVp/rEPjB2eSkIgWOQ0ERyQkGjA46zgOlQzM7/M85x66H6aRLp/4FlyOAgWNz9JVhrcIzlx1oBrrGCpNHDEz3iO8vFtNIl05o675HAQLG4uuzSapLz0tkLye9OjvmXfwQj5lI5BaOE2mkSKh1sYECxuLrs2ukQ/FNqNudr0qEaCK0gSKTvcIiEQLG4MBD85DzX3Dh7SfVLaxpFhnlGLT8oXJZMQCBY3RgSPyTNgK/qfDw76XMWJDXTPHf4oTADB4kZP8Pk1O3YMGYNWvrLDgFPoMhX8E05soHvuQLCDoCd4dXReXlxXo0fdlkwJ9uvjUoP3vdM9dxVyxWASAsHiRk9wyThyiTZGg1+S8obc8yoWTPfcnetY8xMJgWBxw09wWfjo4JAJr1QqZRkpU7OifTpEZk/sEk5CIFjcGAie081IduERw013Xt/W37df81cYBIubyxIpNyMZEIrQGw8ZEbzXyKa7x7eBYIfgsmQINyO5fKjJo95lPFk60k13W3v2HYznO0CwI3BZksd9iy/RJrjFPFZduz3x1kK86W5rz1FVINhR4CdYxYyXykZfQE97eCQ91rVnpvxmQnh8NxICweJGT/CWXibSR8+xgtkrdFv0RdZtdcEmRTYUQnMcLkviuW9fTzWRPvosWQ9meTamV69Oj2LBuhRSECxuLnss53GUipkW7Nfv1JpChP64hgWT6SwSAsHihp/gm0zeT3XvhO2Jq6ofuwcLJtNZJASCRcpuWhe61I2P4MPMbfZ2T+XaYA/pjLKs1IgxWYGeIFjM/C6hdaHTJXMumqNaO5NFR1l40RDPaZG7QLA4qWbU5Otl9zaB5pitncmioyy8aIjntMgDgGBxohXsx+doLn2UjrLwmhKe0yIhECxOLBP8BdOtunZ7ODOfjLKw4IVrUCUJgWBxYpngyiA/qWz02IjitfHJC0g9w69T0x8iIRAsTiwUnHpfJapNmFg8+lWE0j+jDbJySAgEi5NmBZ+SMPpMYQUXr0BHp00tPjwMfZeGaG6DihwLgsVJs4LPxlXrU8sK/jUNTTs2tbg2vPqJ9YjmNoBgMWNM8J4Nmm2USwMNjj7FeLp4SIZMLUYFuzpe3kBzG0CwmDEmuG8XzUbogV4GR59i6op7P4xYwZ+nDVfJaG7DHRICweLEmOCJ2zTfnE0wOJoV/KvrUdS/64Ao/7dy3HKLIjwCptNylSBYnDQrWOGl3+Wu2zK6Hux5oDTgxyttXlDJNFOW5FgQLE4MBa+fuXRp3FxNsGGvfp9KxSVuPbg052B0YR4rmE5ZkmNBsDgxFDxoyNKlsY+ZOvoXRiYLyFQujkborXxW8JpCpBi2j4RAsDgxFDxpq94lugm/MHeRetGoxd5/1I/cdScAKeKq/hW+goRAsDixVLCPj4+X92JpoFvc7dGxchQaEZUAdbLEjDHB2o9JBrzE1ehY7IprdNDZLEW2gpwIgsWJMcH/3GCYOEr5mKvRsTgSpzTQ2SwQLG44wf91ozsnpROaO/ocGUVLd5KUBjybdQ0EixxOcN1zdO9ztxXNHc0JJikNeDYLgWCRwwnmwJdo03CCM+WlKUlR05Ln5g0IjoHGWGKmieAt+U0qRnMcQ0qJlIVRVshLA36sHyrbxd6OKiIngmBx0kTwBmZGk6LvGnYjpQwfI8WCc9DBoL7sLdpLd1qCYHHSRPDzdDbSOOfITJbfd4/KSyfhySx6S0IgWJwYCu7Shrln+mj8N1i9yP2cv7w0Gk9m0VsSAsHixFDwTx80/w5mg+Uuw1zDSlO83dvf2eQb7Df3RRICweLEUDC6yew11X30a24mSyVDayPwQuEmySuqk53IeSBYnDQRXBNtsvvoWm4mixU83gcvFG7yRWRjFgLBYqWJ4GbRfA4+RxYKNbkNeGMWAsFixSrBX5GFQk1uQwVcosWMVYKVXdPSgqXJY0cOciuYH01LA4BgccJXcMNXFy/e0M5k0R3v0H3UAeAr+C23wMBHtTNZdMc76T56GVrbiZdFO3Zs5ilYQcrPagWTHe/Qu1Dk3GAeyBvDT/CpHsH4i2aqUkl3vINgkXOTHTLxvEQXh5JKHJqpSiVOK1zcAIJFjgWC93cll2iyL7q8PSrMSJmjLhtxf2R2DXcECBYflgjuIGGYqdxUpX41Ye4IECw++Ao+kZHXOeyX6uoabqpSv5owdxAIFh98Be9N3jGlPflOxRSE+/b7mFQTJh3eF25TDIY9WWKFr+B3+6H9Q8h3Kmbk5bry6Cm4mrBGcE0VCBYrlgs+wVSzt/teSE3vEdgmb7Sye7tuFXIQLFbMC679Ele429EdbfTv1q3bCrRlhCZA5ypxqTsQLF7MC17jgSvcebdBf/5ToVBcRFuyEcrx9Byoqc6BS92NjYV60WLFvODd9+Nb9hKtQcEEyQIyX85EdK6SrBbCO1i0WC5YxUypbpjkl4noXCUIFjfWCJ4ZLO0QP6goLUSW8sDw+9yyayrakllqECxC+Ap+LYSmjxZzM1ma5WAyk1WecZscCoLFB1/Bl+bQ1JUT3EyWbjn48W2rn6CHgmDx0VhwjQvTFIP6d5oaHQ8mKbJXh9KJjnWraQgEi4/Ggu8w16sN2T6y8RmaGh3l7d8dOXQzFXymA1S6EyuGgu82OYJeonVgwai8fUVajM+EYWUd01J6bFUP9iAhECw+eAqet06TPnqcG2SxH42yK+R0qLUxASrdiRWegts8qEkf3cUNsqhgOtR6pFsSORQEiw+egmU/aH/WNMbSCCazWTNKaAgEiw9OMB1QXTMi+AmfxMREyb/1ziCNsQ4QwXQ2a6+8noRAsPjQCL7I1XW/0+SIOxVnz571172DNY2xSuQLfSOia1d5+afVLw9yJyEQLD40gv/nT34ydokmSN/QFCP93khjrP0Da+FjkljhKbhthqac8BojjbFm7tAcBYLFB0/BehQzy9jbnkx27mBN8ig3xgLBIsRywV+4BVbXbvcNnx/7vSZ5lB1jVZEQCBaQ3LBEY8RbLLiybaiv7P7oicXTvZN7pMaFI2X3dn5pJASCBWRKwVljfGgg+Dsz3Udv6vomVcjx6CrpEh5nnSSng2ABmfWc0bsbX6JrPcx1H52l65tUIcejqz9q8DhrCTkdBAsIL8F8qExF933e9h4WjEdXF1SkSAcJgWABaTXB+5jlxb0ffpd5qEKOR1exVWRXFgmBYAFpPcGuqb+6Hs33ZgWrFyWnlSIsmAKCBaQZwb8zS2e5834gbffRCjlC84tvjmyXflwTAsEC0ozguhFLZzYvOFvb+32G/igaCy4pQBdGaw4DwQLSjGAWM5fox+ZzH5Nu6Y+iseAzkdP33NQcBoIFhJ/g94334ug9Q3fG/zGhHkHur1PBs7ZtuHGgdyzd0AGChYSf4KxORpvphA/VnfF/jKq49xjvCcXnOqK7SUXea1VBoZq3MAgWEH6Cpz9v9Kh5T+q+ZwWzo+ioUcWFQW3ajs1wC/GVZOA3MwLBgtIiwcMiue6jq9BdRoXQ6ZDbtEiHSobY/0Cw8JgXvG7F+vXtHjF6VOUhrvuokhUslfoyaxAt0qERHEsOA8ECYl5w50GzZsVMNPtAd5kpwX6d3e8sl+IiHe+5wztYHJgXPPyfJi/R+txhxv1U947Lp7mpuEhHtRTdCUDseAsDggVg/Uq6Yz19mtGw5YLfZX5nb8O3Fnv7tYkfVOHfbYTUM5u+8UGwAHQaRHesx4w3GjYQPLS3qYZJGl7Tbror88TtRyvkerXQQLAAyPfQr81fopV42yz7Hn9zicmeZ5Ry7aY7mjnKCtbVQgPBAsBPMPq1uvq+Ul4PWJ7waJhXp267aWIhK5h8JSEQLAA8BSN6ieZBnW/il7fmur4MgsVBqws+EpsnlY3eEXlsSPdJaR0Xysv8PlP3IoWGQbAQ8Bfcc4KZv7/0b/CK6fjgGpcreLMdbfG+fyCkjwoGf8GzHjQzgqaj6OVkZanG9QrebEdbvHOpDSBYAPgL5slexpMlIyibbNUhLd5nlECDaMFodcE3mAeu1pVHPk4Fkxbve+VHSQgEC0CrC1YxM0K8OpQqspWR7RJn1Ib3HhA1rTCKhECwlazobLJZpFl8CuljmBR8XGEJSlYwvkQfUGST/NE1oaTDO1yiW8SjY011ezVPiqbgvgnBd0O6WcQq2pQDIUU2zh+N85lPOryD4BaxZJn155q5RFuOiunhG/OUWpFNm4/SDu8guEWISvBvjKLhSuZrimzafJR2eAfBLUJUgr9ifHx8vNsrEjKC/ZIf7t43pE9QQNtZJASCraQ1BE9fbCYxlB8qpJ4SMbn4CqJbspQu06aQJFLyHCDYSlpD8BAfc5mhvMDTWNf3zQndRbdkKWNLp5MkUvIcINhKWkNwq/HhLvbmRA/aN4kr0XGBhECwlYhK8DfhY8O8QvosJH2TaImOfw2IISEQbCWiElyfFuAXnBX2dGp6/rwhWLB6UVIyrbMDgq1EVIKPJOBuo+9GqpXdJ02g3Ucr5DkkBIKtpCWCE3xbZXDFMavxejDdcVchh0p3LaIlgn9RtsrHI46bjdeDuZ13ILhFtERwq7OX9Gx4NUKN//5yG7NAcIsQlWC8Hnwvz6sE6QumRWpBsJWISjBZD46PQoXtPOaodwZMmjC+XYV8PgmBYCsRWHD7Trrlwicat3invUcRAsEtQmDBERt0y8vfNm7xjgdaILjFCCS47gStZxn+H/17NY2xdFOVpFYHCYFgKxFI8C4vWpHWZa/+vZrGWLqpSlKrg4RAsJUIJPh4D/o13uAdTBpjnfbxDBgzhGy8m9QlfCwJgWArEVhw0Djdxvdd3CCLFkCjjRs0Cf4g2FoEEvxMIv268Bld6spxbpBFC6DRxg0guIUIJDjGy9i9v5Bts9sUAw/Pia2lG+9AcAsRSHDfeGP3/kLL/i9fgf4KvUk33mlKdIBgaxGj4INBqbLQpBfXxocFscOs+0gIBFuJQIJ7BTVNHz1BMxvaIq5jAx1mHSQngGArEUhwVlzTqpXFXGYD17GBDrNofgwIthK7CL45yfDdGm/0abXdR2nHBjrMWk5CINhK7CL4S2/DZO+wBcaO03YfpR0b6DAL2uq0CLsI/j7Y8J6c14wdp+0+WhjqEfBm2cg014Au0X1ICARbiWCCG94yTB+9qC2ERtcLyzy/uem7HtEqayDYSgQTfMbDMH10hbYQGl0vxJuy4i4jmvgEgq1EMMHfhTU9LlHiwX5QSthN1wvxnp246yC4ZQgkuGv3wUaasah7TQqp3Z54i64XguDWQCDBb21a42vkwOsjJLLRF9AUT4+UeUMK/RJnxO4ZAd1HW4SNBW/bjDurPNmmaeS70GoDcFbDeXzhpmMsOo9V5vUtORwEW4mNBSdm4jmqoUYux1+7MAZMRRrBdIxF57HwSAsDgvnxi+HbpmCRTZ+v37v4tukl2hREMB1j0Xks/IcYA4J5USExfNsw7W36hBYK/jPXxa/fKTrGUoQNrR+7BwRbxNkEw3tsfIm2THBDrxEhde+EfUDSR9E07+QFahBsESIXXBaFB1p7KlFRWsqSBkU2bWFIQiCYF0IJduOTaPqodqryZKdb94aVgGDLEUqwxyYe6aM3tFOVS55D6O1cEGw5Agk+77qZ19HlbX/3WohG716yHqEDVLAiClJXLMCOgr//FKempB/C31/xeIzXOepeWSERLyfeqmAv0XL2En2uI7qRAIItwI6CO4fg1BS3peSH+3fzO+l6UBuXvhfUhQHeAVknu2ejtJSsUBBsAXYUPJ787aSXaJSezqM68Vr2Wu53Z/Yj6GD/2uvhG+kf4Aq5gjwCCOaFUIKP7uZRnfhr9tVMRf+W3i3ciNBMEGwNQgnmSX2Eu0TCrF3+PEKFrGBFdnYHEGwRIhd8zDv7ct0E71d719xI2Xiuo2JEJAi2DJELHtLmFkIX3N94KjnMP2ROfkZEZofNwSDYAgQTPHui+b5Ju9GWEZqTv8xduSHlW0U2nevAgGBeCCa402Tznc+OoS3ZCOV4eg5kX+jQvt6nQbDFWCJ4WB/r+7Gw+OfjB+EEZx3i8/I+ifjFayE6lKlZEQbBFmOJ4C6Lre/HwpLxNH6QCG7Z+R+8Xt+g3iGRn3bIRJoV4WwymUUAwbywRHD/Iy16KnqJvq3ZOJLJ6x2MboX6STp+0nnEQDe/ifNi2/qMyr4ZRfdzgWBe2F0wR5fZO8zyMTeTRevM0q13JQX0AUAwLwQT/OCIpumihrzEzWRpOruTrXdnIkmJYRDMD8EE8+MyI2EpoWVI6UAL3ThMQiCYFyIX/KZEjf7qIaOC6UBr1VoaAsG8MCe4zDdRi2vLNtRaI3iUG3vzX8nTuMP7tqUeHukFXQbL3EgIBPPCnOCTKWe1dDrcoqcaNd+i1qMsl9BVT3zm8J240Czdl4VTSMnDgWBemBNc0Un3fQsv0SNSLOs92q3bMo1g+eu40Czdl8UlNoBgfthRsDVc9Zwe7Jcp+2ynH/vCyL4sbls0COaHyAX/6DL+pxsjPa5iwXRfFgi2DJELLmWCZP7jtlfubDMqKntTolf7B3tnaEIgmBeWCO73lmGimo2p1W58L/M8px66Hw+1Dg9D35G7QDAvLBEc0yRPzcZM0W58p1NZeKhVG179BHk1IJgXlggWgPK29ezt6N10pgOXfEcFu+hyEgjmhcgFd3JxZTzd/A/30Qn+PG04CYFgXohcsKZGR5meYJSyj0RAsI67irOmKI00PFhcgmmGf9kAdhBdgwozUuao32qTTQIgWMfzbRJNEeFmeLAAgn/zNJ4+OhtxgskgmqwHFwW7f0NOAsE69g8xGRLFJfq6u/H00WrECdZbD4YiLE0Rs+Djuw4ceLXJZYRDxQS4BGQq9deDoYRDU8QsOKP9kCH9JKaiKuabMPWiUfrrwSC4KRYJfjhLf2vywyE2fWUIjSthL9EepqJYMCpvfyisXXrutq/DPX3TB20IpSEQrMMiwVtH6ycXDAuw6Sujgq+4GN9J/Rw0xuKHRYIbY/kl+ou55lMW9EhcxZ6z0fhO6i+hMRY/7Cp4bbz5pCM9QuY092Ca7qPoBtcY6/pKENwUuwp+fYz5Y/TAl2jTaLqPrlrLNca63xcyG5ri2ILZQdanfv5+W4vSgqXJSW65+ckkBIJ1iFlwn+FLTbKXG2Tp9typZOgkfUkgWIe9BH+Gpxy9Ai04A6H1s00LPsINsnR77ljBS6DiuyH2Enw25uLFi8/JLTjDHJpBFl5GonvusOD1JASCddhNMH4wCy/RzaMZZGHBdM/dnQDuJYFgHY4tmB1kKcPTUpZsTpP5dpkdLi+i3VlAsA4bC/6V2yR3Mo79aUuG+bxQnpziBll0gFWecftWh5Xs9+RZQbAO2wpeptsl58L+eKC7+bxQnmzlBll0gLX6CYSeWakZY4FgPWwruGQc3wezHM0giw6w1q1GaFSmZowFgvVwZMFkkEUHWGc6qP6Ky6yAS3QTHFswO8hCOR4ePevUA918ou5DWV4kBIL/XM4NVqammzzIrOB0s/kHW7TD5s+jTG7us4prjQdZL8TV/Ow16ctcNXkyEPyuPzdYyTDS9lGDOcHvNm270wTtkvFZf5Ob+6zi8caDrHmdktr3WYnOFpEnA8EfaxsgteASzQPdJbr1UTEF4b79PqaDrOXPK7KfWrm/60oSAsHOIXjk5bry6J1kkHWy51F555UL10AxUoozCD7B4K2z+94uSguTJc54Ktp/2ph2SfEkBILtJniUrdJH7+mqzdL9WIrs/V2qELyDKfYSvNpm6aP5umqzdD+WItJ9BwLBGuwl2JZ8EoFb2x3KRLTzqCL+vfBqEKzBGQSXuY2vVo/zy0R4P9bYPYpsNG0WCNbgDILRpKRgH4/3B+HNmskL1KzgXwNOgWCKUwj+M/JwfBlCO9slzqgtG3F/ZHYNKkojERDsKIKfCTQ1k4Vrhr/nUsCNoWlRYVgP5nAUwVuHmZqL/pWNHvIajzRjaJpECuvBHLwEnwlrUZn+1mDKINP/hjJp8KlO942mY2iaYwjrwRy8BJ8Lt7SAZKsTFmP636COTEc7XM7QMTQVDOvBHLwEi4DtD5qOvdAu4nB8fq+6HHf3DjXPhN0fmbgZ5dBkUxDsKIKXhptqxLMBwSCrGRxF8PebTf11/gLBIKsZHEVwc5wPDD3VaRedqLz8yBjFYBhk6XAGwd+5r0SVQRfIIOt1/zE1VTDI0uEMgpe61iL0VK/5Ee7SBTmukRXy3PRR/jS7DQTrBG+PN53BJ2YOoPMy/PppRyyVDFXIccsGX/KPAsE6wWceEFqVdRzSCKYV0KhgdqAVR+4DwTrBjst5/56+MU/tIxXQVP7ZrOAxIJjDGQRXSk42XMkcRiqg3fEFwY1wBsH/Ynx8fEJfS5UFBUSrk0M2hPTyLJjvSkJiFPyxZSWkWsj0YKH/vRZQtt5I+uhppJ4SMbn4imaUpcgmQywRfw5enGxRCakWMtxf6H+vBQzvYCR9FKcwXN83J3QXHWWxgvEQ6zI9Q4yCi/Lt+WwOdYme+7TRuz/cPjvMK7Et1+Cd/AW+TkMgWDSCO3Q0v2boN9DoqV95DLh4I89brmnwDoL1EY3gxCLzfSiHLTB66pGwrtKQh94MTY2LivfsNeJYcmpSIAjWIBrB7T4xf8yM2Ua27PyMVpBu3zUuV+hS4f6BtaqkSyfJGSDYoQTnhRrZdFeIls/AwRrXK3SpcOYOhP6oWULOAMH2EZzqZ7yhhh6S7VY/+l7GkyUjQk2368woQeiCCgRT7CO47yvGG2rokcjjHWyCG8wDV+vKI0sQFbxXXl8VW1VBQiDYPoIHvGP2ED6XaBOomBkhXh1Kb46M8TuOFuamBQamTe7VQEIgWDSC08qsTB+tYwXjS/QBWtKfVupQ93+JPCoIFo1gqbXpo5NJERYWWtKfVupASvGuB/9NBbcAjWBa0p9W6kCfSMldINgZBGsq3dGS/rQc2t2UMhICwbYUXKkdInc9aLtnwYJ9fHzcA3/OSk3dyv7+0lIWNyyaQEMg2IaCT7hpP+S6PG6zZ0GaS3R5e6x2SQNNIlUvTyUhEGxDwXo5ia11iW7wS2gyk7WgcaU7mt+wf2AtOQEEO5bge8xHTeairzeudEfzG/B0JQYEi0lw7S5zy0llzC0j5zUqJ0yTSGdo+vCAYDEJPuptbkG4Kyl5ZkijcsI0iXSvvJ6EQLCYBH+WbO6IeybewXrlhBFJIm3oJd70URBsmjomp8merM0wyGoO+wjuPIdPe41F5nd4Ptn0rM+EHWQ1nDO7UGbAyjybvyg97CP4gRF82mv0b2PdE51ilrK37zI5+oMse9XJ2upudqnbAM8Im78oPewjmB/mL9HGOcWksLf5sij9Qda/SMj2gt8ZYOkZzniJ5of1grMqUW3CxLy0wNCUOc+0CV2gTp8URUIguCzYZml/E/0sfC0ngqx5moOoMrV4BY7sPc0AABAdSURBVDo6bWoxzW04PAx9l2a3WpViF3xpdGt71TLR18LXcmWUNU/zNiv41zQ07djUYprbUBte/cR6EGx77FYf7xQjc/fzOj+1mOY2oIJdHa+BYNtjR8F1xb07jJpavJDkNqDP04bbryA4CLY9rOBfXde2j5n0darniHmJUV2kOaiIVsYDwTaEj+DdpS2uYvkFN5NFB1flGbe/91hmv0JoILhZorqYKmDHmw3cTBYdXK1+4mD0Q3bsPgqCm6XH8VZ4IhVDU80KSsKvIdx99IWJbe2WAA6CmyPe47lWeCJO8Ocp7gjh7qO9J/aGS7Tt4SHYL+ZIKzwRK1i9PDXpReTrmosKozxkQ6dnF/iQEAi2ITwEB7Tf1uLuo7+wGqtIyujpeBn6csizdcldns+ukJMnAME2hI/gxOgWdx+djwXj5cE3ov4pQ+jTDKlkmgIE24ET7mblSE63yjMxVTRlVCVDXKEOEGwHmu6ANMT3Yms8z1UmyNdt9DexVXcC0EJNoY4d9spN+jsLNk9AKwmuUi8IaFOKUFf516np+fOGKLJrxpMQCBaWVhLsyd7ciziNCjNS5qhpZyy4RIsCj1boPjqZCkbDd9LlYFqJBQSLgjut0j+YCpa/TpeDaSUWEOxEHGdkMunQsM/oCFrTOAkEOw/HJQj9mRCnGUGDYKcDv4P9Q3cp0/2l+Rl+UQnqsoQOshASEqPgx/xaPLnzt+JxVOkaF+eZXqV0eUWlN84iv0wxCq76oMXTs38rrqHKKITq14xTxiKkN84iv0wxCgYs5RQuo+ThrUxHeKayjw/+M6yIJSEQ7AycYupQw7IHsOCFa277kHFWzWMkBIKdgVPsICss5/pbvghNDvOTxPZNjp5QEU9CIFiMrM3kk6rGsYWrdIc/Gc0vVskQHm2J+GMSMHcon2RTjk+5Qmg6wbH0BwSCRchv1dWF8y07RcVUeS2kTmcVv+eO8B9jECxS/pQwDMkGtQAV81JYZAM61xHdTSr+PQgEi5lqRo1WLrTsHBXj4inx6lMaF+Adu+2Y+80sz/TjIFikYMHTe1qSXfgOQuf97sx+BJUG/Fg/qqhCTssKU0Cw2MCCd06yRPA+hJZMRf+W3i3NQWhvXoWclhWmgGCx8QPTYPE5/3XxkUpd2myahNBb+RUdp5OywjQEgsXGJSsEv+aiRmiZ+6boP+pH7qromEHKCtMQCBYb1gjOxtV5LjCr+0d4Bjx6suOQADfpZmV3EgLBYuMSs8/C9NFKpJRIWZhNOWStsFT6v7pRRUoX8nAgWDB+TzM68ZjNWJo+uh4pSYt36abpZK2wdCweaylpBRgQLBg/ehideNzMjqIt5WNGJvMb46+ke7JKwzPYsRae7EAgWECuGi9rV22NYAlCfy1wU9I9Wa9JfmHHWiBYaFpZMLoneft1n/T8zmldJdLkfqkJtO4lCBYMk4I/t3DLzm/aje/a5SRSd5Y8HAgWDBOCVd6Wbrqbp934rhVM6s6ShwPBgmFCsFUPRXOTPqPrhVgwLu5PQiBYMFpV8J9e/yjoX7ugY+26JCyYFPcnIRAsGK0pmPFxkeRcV8nyB3iMxYJJcywSAsGC0YqCEer5ZsgJlOOWy/5Hc0i5AAi2mEubWlybjrDdw/rXEKOfPprPrQezb132P5rbwB0Jgi1mbUiLa9MRMl2sfw39S/XSR2u59WAqmOY2cEeCYIvZfX/rPI5Vl+iFdD4z+pVG93bE+7iYEVSwpqSwBhBsMUIKvsOMIysSvs80uvtYwndhaFnOnQCkrcLChUCwxQgr+C75atCG6cHnzoehC15n20RM6BqX5hGUPy8e173DgGCLEYPg2CS9dnfL2DtYwQjp0kdJ3TtyJAi2GBsL/nduM/vrHtcI/rRMr2Hl/xAnWJs+OnPH3LnTyZG2EVz3m26MV9rXJk8hHDYW/HwyD8GG/Jnr4tfvlCZ9lB1izShBZUNIxDaCJ+t/TPO2yVMIh40Fvz6mmXPuGBfc0GtESN07YR+T9FE8xNorr99Pf++2Ebxutu57uESboBUFl0WdYy/Re97zexV97ZmYmRQY1SEpJo2EQLDFCCs4t+kurm1onabbI23aQBMcysgSIgi2nFYTLDFeuM6/uZMWG9nF9Ql6YQKONTTQpg00wQGX6UAg2ApaSzC6ZbRu3bZRlj9SecIjwX79+u2mHbFKSYJDGb0SgGCLaTXBxmn2Em2Cet+kb+886nqedsQqJQkOINhaRCi4LDxPKhu9uRIVekbMUS8OCnTtOC84iIRAsMXYWPDG9ust5Kh2kPVlduKdlG8Xu35z03c9smW9aBBsPf+Sz7KQN7WDrIOhjzzmfXpxJEJxlxHt1wOCLcbGgq1hD8lsGLZbV2027joI5s+LRfqfSiZ0Evr1NOGKS0H1jUFet3TVZkGwJcT1159X6Bgj9OtpwlWPB6Wy+0NPH/WSxs9LifFd0hAzuF0YCYFgHvQ6pv+T4JfoXzc33tz1lTazAa8XkqSGElymgxwMgnkgMsFbghpv7lqnzWzA64UkqSEXl+kgB4NgHohM8MH7DO85gwdZo8I+w+uFJKkhF904PIeEQDAP7CE4M5P3dsyUIMOTz7AaqzpLaQE0nNRQsmot+ouEQDAP7CE4+QneG6r/0cHwZPwO9h8bdvoD70kTUI6HR8+672UesJrEG3sI7nyS96Fv977YiD9hkNVCRCa4xCuwEbNhkNVC+AmOi21JZw3XF1rwArn0URhkWQc/wQlvtKSzRuIHLXiBx5nl6HZBe+ZpGGRZBT/ByZ+15DksuEQ35TjjKvMf90AEFkwyRxt+zkpNJSEQzAPRC65Mva8S1SZMLKY1HG6ObJd+HNEMUhDMA6OC95QafHoJ/ndLniPQ+u6jU1nBxSvQ0WlTNYJpOeH95HFBMA+MCo7tYjD/4P7PljzHTeu7j9awgn9NQ9OOcYJpOeGZ5HFBMA86vK//k0Zwn/cNjmrZJbpFVKb+6TK57T0qeFYxujG1W2ztDBICwTyQzNT/SZSCt/t7TEVTi2nXhlVrNy4NvbmXhEAwD7w26/8kSsE9X2KeRVPz2/lG5I4dEukvlW3d2Y6EQDAPGgve3I5sdQvZYXBUQLtuwrCcq9HxRdZtdcGmg5k1fySsZL8nrwoE86Cx4PJhZKubtMjgqOPlCmG4xNXoeDamV69OjxZuRGjxSvZ78qpAMA8aC9bQ5BItIJoaHWsKETrcb1RPhJ5aib/HgGAeiF6wpkaHIq6qfvSO5SF3VZ1Xst+TEAjmgegFczU6onz8HpCXpnq7Jy9D2e4kBIJ5YFRwfD9LOoTakpeRXo2OCjlt2XCy0y3yMkEwD4wK3lxkJJNTED5EejU6WMGkZcMSui0aBPPBqGBxcUayEGl6UpZOxNmjeFkYA4J54ACC17hENnCCI6tweyzBL9GrsoT+28Ubl5U2+S21Jp2kISfQD8GpSY/Ju3bsnRx0ShlDKwU0K7jOWMEAPshH6B7ElOCHRrXqHyJb4roOv+BnVlua12knyriZLFL+7BIt1MEOuMwL/t7dyv/lEzPMCy54qrmnFhX0Ep0yzNK8TjvxBjeTNXMHQn/U0EIdpCgaMiP4sp+VvxE+l2iHE9z/iNCvwyQ3GIlEwrw0owShCypcqGP1yt0+NASCeeAajPc9ui0S+nWY5BBzD6l7he+V11fFVuFCHRtXfkvLZIFgPsSS/ZKdDgr9OkySy9Qh9IbkD7mHe1at+nG3lBkr3/LlsScLBFOS/g/fdl3c8jL+NuFrpGIF18zM1ywX9ki4nrCyQs5jTxYIplDBWX1aXsbfJjzDCpZKA0f8TJcL75ceQotZwTz2ZIFgChUsYlTMtGC/fqfwEuEf15Y/j5cLK+Q89mSBYIroBd9k8n6qeydsT1xV/dg9J3vi5cIKOY89WbYS3OCK16e7Wvng9kf0gg8zeH/OnspCT7/U0a+uTfG6f0xWYBQJCSK4nrlcXT3jSSsf3P4kHL4oZm5oC6Hh9UJdxVlyl0CCbzjUJTpUGihmZmkLoeHpK13FWQIIdgbK2+K36+jdpOu3tuIsAQQ7A18wXkHSdnG3iGBtxVkSAsHOQGX4g1JZ+HyyIIxQyj5WcP8+ySQEgp2BSnbE/NfQI8rukyZUyK8lFhSjnDZdqAMQbCPcrU8HtZSHUaVrXJxnehXZdJcRfWh+cXnc1Fu0Fg8Itg1q5gu7fUxSkXdw/ZpxdNMdTiBd/QRCtL2hcIInOc6ODit4ham28ndnDZUR04P9ukTQPVk4gXTdaoRoVRfbC95u7KryG0JLBrTqvimxwfxs5e/OGv7jMf6n2nFeuDEWTSA900H1V28Ssr1gdKNJSnoVfgc7N2q7voNfZoJkYTnb30svzEhJuS+vS3hi7+TEYBKyg+Cm1P8NBEsy7JY+ukrXs4G0Hv0i6/aXfTa9RxuMg2Absd9+6aNK7VQlbT36bEyvHlJpbAJ5Hc0LtqYJOQYE25fytr97LUSjac+GrclkWZhP36SPrN0WD4Lti7pXVkjEy4m0Z8PWXmRZmI9guEQ7CteD2rj0vVA6zNu97WNds9NkyQvSi7JJxCkFlw4Vei+6XSnhMhto21FFNl0SLiG/C6cUvKar0NkkduUDLrOBth1VZNMl4TPkdyGQ4MvW13XjwbOTrXzdjkp9uJtEwqyjbUcV2XRJmCKI4AYXG0+/h1j5uh2Vo57jf6ob57mFtB1lBZMl4VUkJIhgW1OUL+jT258B/jUIXXBf2t3XL7ku1z9lTsquke14lBN2VMF/s7/B5dxMVmmeyysqPJuVGF3Mq2eDowreM0zoga1d2cnNZL2Zh3NGzxbJJc/x69ngqIL/duxkcNk9z+fxeiHtQMqvZwMIdhC+c23z39rtrl9iwbQDKb+eDSDYQTgf3dFdNjpYqYxsl5ibEheV4DEt2s2bhJxE8Gfv2G/xRnT8gM7H1bb/AIUpaf5omSffFu+OIzgzxm6rr+JjNSsYfZRWG6ak+aNlg/m2eHccwdM22fkJxcV/GJnMIy5YuUa7UBi3f/B0ErKN4OnRds6BDswx/6KcmP8wdeh3f1elQrtQGPdjlS0Fn11r502MXR+38pU6B1gwWsooi0LcA+c3POib9GLcIbktBdudKcvttgtZfNzianTQxqO0HFqFXEF+M04iOMtH6BxOAZnJ1eigjUdpOTQnE/w3B7+DEdd4lJZDA8HOhIqp8lrINR6l5dBAsDOhYl4Ki2zQNB5VL0pOK0Ug2Lno+WbICVQ2clR4wfzoPlWleQOiptWSQLOCL7gIPX4AeDCT23THzVEW0d135gWjH4T+BADw4Ba36Y6bo3yC7r7jIRhwDOojfaRSl4NkCus6Fkx235EQCHYGjiWoEVqWoxNMdt+REAh2Bh7EK0cXvPZoBffvkzyXfDQGwU5J6XTttyDYGQHBfx9AsJMDgp0cEOzkgGAnBwQ7OSDYyQHBTg4IdnJAsJMDgp0cEOzkgGAnBwQ7OSDYyQHBTg4IdnJAsJMDgp0cEOzkgGAnBwQ7OSDYyfn/ITnDH/rW4VAAAAAASUVORK5CYII=\n" - }, - "metadata": {}, - "output_type": "display_data" - } - ], - "source": [ - "%%R -i taxnames\n", - "library(ggtree)\n", - "tree <- read.tree(\"data/TimeTree.nwk\")\n", - "tree_matching_labels <- tree$tip.label[match(taxnames, tree$tip.label)]\n", - "\n", - "library(ape)\n", - "pruned_tree <- drop.tip(tree, tree$tip.label[-match(taxnames, tree$tip.label)])\n", - "\n", - "write.tree(pruned_tree, file = \"data/TimeTree_pruned.nwk\")\n", - "ggtree(pruned_tree) + geom_tiplab()" - ] } ], "metadata": {