使用Y叔神包ggtree进行基因家族基因进化树构建 |
您所在的位置:网站首页 › 进化树分组 › 使用Y叔神包ggtree进行基因家族基因进化树构建 |
![]() 系列目录 基于全基因组的基因家族分析(1):数据准备 基于全基因组的基因家族分析(2):SlNRAMP家族基因成员鉴定 基于全基因组的基因家族分析(3):SlNRAMP家族基因CDS和Genomic DNA序列获取 ![]() 继续之前还没有完成的内容。因为最近在学习Y叔的R包--ggtree,所以就顺便拿这个内容来进行展示,作为一个例子来记录。nwk树文件和R代码文件我已经放在github:https://github.com/Lxmic/ggtree_note 1. ggtree安装关于这个包,y叔已经写了相当详细的说明,有一个完整的电子书《Data Integration, Manipulation and Visualization of Phylogenetic Trees》。这个包是在Bioconductor上面,有非常详细的介绍,可以去查找相关的内容。 #安装相关的包,包括ggtree以及ggplot2#对于R版本在3.6及以上的,需要使用BiocManager包来安装bioconductor上的包if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager")BiocManager::install("ggtree")#ggplot2安装install.packages("ggplot2")#载入两个包library(ggplot2)library(ggtree) 2. 读取及可视化树结构关于用什么算法以及什么软件来构建你自己需要的树,完全看个人的需要,y叔在电子文档《Data Integration, Manipulation and Visualization of Phylogenetic Trees》中有很详细的介绍各种树以及各种算法,有兴趣了解一下(反正我是看不太懂)。我就用最简单,最常用的方法来获得进化树——MEGA软件,可以输出newick格式的树,非常常用的进化树文件(我们需要保存其bootstrap值以及branch.length值)。 # 读取newick树,在当前工作目录中的nramp.nwk文件,并赋值给treetree |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |