本節(jié)來介紹如何使用ggtree繪制動物之間的系統(tǒng)發(fā)育樹并添加動物輪廓。圖形的繪制過程比較簡單，主要是數(shù)據(jù)的整理。

加載R包

library(rotl)
library(ggtree)  
library(tidyverse)
library(ggimage) 
library(extrafont)

構(gòu)建分類信息

taxa <- tnrs_match_names(names = c("Ornithorhynchus anatinus", 
 "Tachyglossus aculeatus", "Phascolarctos cinereus",  
"Macropus giganteus", "Ailuropoda melanoleuca",  
"Carcharodon carcharias", "Megaptera novaeangliae", 
 "Eudyptula minor",  "Tiliqua scincoides", "Notechis scutatus",  
"Dromaius novaehollandiae", "Dacelo novaeguineae", 
"Myrmecia gulosa",  "Musca domestica"))

繪制分類系統(tǒng)發(fā)育樹

tree_data <- tol_induced_subtree(ott_ids = ott_id(taxa))

導(dǎo)入動物輪廓數(shù)據(jù)

phylopic_info <- read_tsv("phylopic.xls")

繪制系統(tǒng)發(fā)育樹

ggtree(tree_data) %<+% phylopic_info +
  geom_tiplab(aes(image=phylopic),geom="phylopic", alpha=.5, color='#3CB2EC', offset=.1) +
  geom_tiplab(aes(label=common_name),offset = .75, col="black") + 
  xlim(NA,10) +
  theme(plot.caption=element_text(size=10, face='italic'))

好了本節(jié)介紹到此結(jié)束，非常簡單的一個案例喜歡的觀眾老爺歡迎分享轉(zhuǎn)發(fā),更多精彩案例歡迎關(guān)注我的公粽號R語言數(shù)據(jù)分析指南