能看大秀直播平台app

淌若生物信息只学习一个R包，那一定是ggplot2。

淌若ggplot2里只可选拔一种，那一定是ggtree。

在生物信息学范围，进化树是照顾物种之间关系的蹙迫器具，是咱们分析责任中一定不可清寒的一环，而ggtree便是一款绘画tree的大大大神器。

这款由可视化大神余晖创陶冶研发的R神器，可谓是生物信息研发东说念主员东说念主手一份。

然则可惜的是，ggtree现在是全代码操作，关于号令行操作不熟习的陶冶使用起来较为苦处。为了浮浅大家的使用，特将笔者我方梳理的ggtree的使用纪录共享下~

一、ggtree的装置（仅需一次）

装置 ggtree前先确保您依然装置了 R 言语和 RStudio。接下来在 R 终端台中输入以下代码以装置 `ggtree` 偏激联系依赖包：

install.packages("BiocManager")BiocManager::install("ggtree")# install.packages("treeio")

二、读入文献

 library("ggtree") library("treeio") library("dplyr")# 读入tree文献 tree<-read.newick(nwkfile,node.label="label") # 读入列表文献（可用于后续分组、标记神色等） meta<-read.delim(info,header=TRUE)

这里再啰嗦两句，诠释一下这里读入的nwk文献是啥？nwk也便是Newick方法，是一种最常用于暗示进化树（系统发育树）的文献方法。同期亦然mega、fasttree等多种构树软件的输出文献方法。

nwk文献示例：

((A:0.1, B:0.2):0.3, C:0.4);

这便是一个粗心的Newick方法的树，暗示有一个四个节点的树，其中A和B是子节点，C是另一个子节点，且各个分支具有特定的长度。

三、绘画tree（最要道的来了）

最粗心，最径直，但是不好意思瞻念的tree

ggtree(tree)

图片

为了好意思化，咱们不错把柄需要将其他联所有据（如物种分组信息）澌灭到树数据框中。

# 读入列表文献（可用于后续分组、标记神色等）meta<-read.delim(info,header=TRUE)tree_df <- fortify(tree)  # 澌灭tree和分组信息  tree_df <- tree_df %>%    left_join(meta, by = c("label" = "ID"))  # 绘画tree，并添加分组神色p1 <- ggtree(tree_df,layout="roundrect",branch.length="none") +      geom_tippoint(aes(label = label, color = Group), size = 3) +      theme(legend.position = "right") p1 

图片

这么看，是不是自得了好多？分组更为了了明了一些啦。

具体讲解下ggtree示例中几个小参数的具体含义：

geom_tippoint:用于在叶子节点（tip）上添加点，频繁用于标示具体的样本或物种。不错通过size，color的具体诞生来作念节点的调遣。

scale_color_manual: 用于自界说神色的映射。当你思要为不同的组或分类指定特定的神色时很有用。

theme：用于自界说图形的举座主题和方法，比如字体、布景、网格线等。不错调遣图形的好意思不雅性和可读性。

annotate：用于在图形上添加防范，如文本、箭头或标记。这不错匡助讲解图上的某些特征或数据。

layout 指定树体式的布局计谋，比如树的分支地方、叶节点的位置等。

除示例中的roundrect外，泰国按摩群常用的另外两种布局类型：

"rectangular": 直角布局，默许的基本布局，树呈现为高下或独揽结构。

"circular": 圆形布局，与径向布局雷同，但叶节点沿圆的周围均匀散播。

四、效果保存

这里同ggplot的其他包，齐是用ggsave号令完告成果的保存。一般提议保存为pdf方法，浮浅后续修改。

ggsave("tree_plot.pdf", plot = p1, width = 10, height = 6)

五、批量伊始

说真话淌若是单个tree的绘画和好意思化，我并不提议用ggtree，网页版的iTol是更好的选拔（步步图解iTol-给进化树作念个好意思颜）

ggtree的优点是更适于批量化土产货伊始。比如咱们有一百个文献齐需要作念tree，咱们不错如下操作

for (file in file_list) {  # 轮回读取树文献    tree <- read.tree(file)  # 从文献名索要物种名  # 创成立数据框    tree_df <- fortify(tree)  # 澌灭分组信息   tree_df <- tree_df %>%      left_join(meta, by = c("label" = "ID"))# 绘画树    p1 <- ggtree(tree_df, layout = "roundrect") +       geom_tree() +      geom_tippoint(aes(label = label, color = Group), size = 1) +      theme(legend.position = "none") +       annotate("text", x = max(tree$edge.length) * 0.5, y = -3,                label = file, size = 2, vjust = -1)  # 保存每个图    ggsave(paste0(species_name, ".pdf"), plot = p1, width =11, height = 5)  # 将图添加到列表中    plots[[length(plots) + 1]] <- p1  }  # 使用cowplot将图澌灭，每行知道5个图  combined_plot <- plot_grid(plotlist = plots, ncol = 5)  # 保存澌灭后的图形  ggsave("combined_tree_plot.pdf", plot = combined_plot)

此时，咱们不仅得到了100个pdf方法的tree图，还得到了一张整图（涵盖了所灵验果）。这是iTol饱胀作念不了的事情。

诚然，ggtree还有更多其他的优点，比如，有更大的变化，支撑多样变形。就像乐高通常，ggtree其实莫得固定的套路，支撑每个玩家我方的创造。

其官方教程如下：

https://yulab-smu.top/treedata-book/

临了再次欷歔一下，学习真不是一件容易的事情。

长按怜惜

风作伴，梦作念马，追啊迎啊最强烈的年华

-少小的你啊

公众堪称号：微微悦明

科学的乐趣是得到新常识的旺盛~

高通量测序、大数据病原微生物检测和监测健康大数据行业资讯纪录与共享能看大秀直播平台app

本站仅提供存储处事，通盘本色均由用户发布，如发现存害或侵权本色，请点击举报。

• 能看大秀直播平台app