看到朋友圈有人转发的推文 Nature|重大发现！人类的近亲肺鱼基因组被解析，点进去看到里面有一张进化树的图
正好自己最近在学习R语言的ggtree，之前也在ggtree的帮助文档看到过类似的图片，所以决定重复一下这张图
论文是
Giant lungfish genome elucidates the conquest of land by vertebrates
首先是将右侧的动物图片截图保存好
用拉丁名来命名，属和种之间用下划线分隔
接下来是模拟一个进化树文件
(((((((Anolis_carolinensis:0.4,Gallus_gallus:0.32)94:0.4,Homo_sapiens:0.3)95:0.4,(Ambystoma_mexicanum:0.4,Xenopus_laevis:0.41)93:0.6)90:0.41,Neoceratodus_forsteri:0.3)80:0.3,Latimeria_chalumnae:0.6)99:0.3,(Danio_rerio:0.3,Lepisosteus_oculatus:0.4)95:0.5)100:0.4,Callorhinchus_milii:0.3);
最基本的进化树展示
加载用到的包
library(stringr)
library(ggtree)
library(treeio)
读入进化树
tree1<-read.tree("Nature/Nature_tree_1.nwk")
展示
ggtree(tree1)+
geom_tiplab()+
xlim(NA,4.5)
接下来简单美化
-
去掉拉丁名中的下划线 -
拉丁名改为斜体 -
加粗线
ggtree(tree1,size=2)+
geom_tiplab(aes(label=str_replace(label,"_"," ")),
offset = 0.05,
font="italic")+
xlim(NA,4.5)
最后就是添加图片了
ggtree(tree1,size=2)+
geom_tiplab(aes(label=str_replace(label,"_"," ")),
offset = 0.05,
font="italic")+
xlim(NA,4.5)+
geom_tiplab(aes(image=paste0("Nature/",label,".png")),
geom = "image",size=0.2,offset = 1.25)
出图以后再手动编辑一下图片的位置就可以了 最终的结果如下
哈哈哈，距离发Nature又近了一步，还差9999999999999999步。
如果需要示例文件可以给本篇推文点赞和点击在看，然后在文末留言就好了。
最后是福利时间
前段时间有机械工业出版社的编辑找到我，借助公众号的平台开展一个送书活动，书是 R语言数据可视化实战，总共送出4本，我自己也收到了一本。我自己的这本也送出，所以总共是5本。送书规则如下
-
首先是给本篇推文点赞和点击在看 -
接下来是留言对自己的寄语集赞
留言点赞数第一名直接获得我的这一本，由我寄出
留言点赞数第2-21名获得抽奖资格，瓜分剩下的4本，到时候由出版社的编辑寄出。本周日在B站直播抽奖
书的封面如下
欢迎大家关注我的公众号
小明的数据分析笔记本