生物系统学中的进化树构建和分析R工具包V.PhyloMaker2的介绍和详细使用

V.PhyloMaker2是一个R语言的工具包,专门用于构建和分析生物系统学中的进化树(也称为系统发育树或phylogenetic tree)。以下是对V.PhyloMaker2的一些基本介绍和使用说明:

论文介绍:V.PhyloMaker2: An updated and enlarged R package that can generate very large phylogenies for vascular plants - ScienceDirect

 github仓库代码:jinyizju/V.PhyloMaker2: This package (an updated version of 'V.PhyloMaker') can generate a phylogenetic tree for vascular plants based on three different botanical nomenclature systems. (github.com)

介绍:

V.PhyloMaker2提供了一系列的函数和方法,帮助用户处理和分析分子序列数据,包括但不限于:

  1. 数据预处理:对分子序列数据进行质量控制、格式转换和多重比对。
  2. 进化树构建:支持多种流行的进化树构建方法,如最大似然法(Maximum Likelihood)、贝叶斯推断法(Bayesian Inference)等。
  3. 进化树优化:通过搜索最优的树形结构和参数组合来提高进化树的准确性。
  4. 进化树可视化:提供丰富的图形选项来定制和美化进化树的显示。
  5. 树形数据分析:包括节点支持度评估、分支长度分析、祖先状态重建等。

详细使用:

由于V.PhyloMaker2的具体使用会涉及到具体的代码操作和数据分析过程,以下是一些基本的使用步骤:

  1. 安装V.PhyloMaker2: 在R环境中,使用install.packages("V.PhyloMaker2")命令来安装这个包。

    #BioManager安装
    if (!requireNamespace("BiocManager", quietly = TRUE))install.packages("BiocManager")
    BiocManager::install("V.PhyloMaker2")#github 安装
    install.packages("devtools")library(devtools)
    install_github("JinYongJiang/V.PhyloMaker")
  2. 加载V.PhyloMaker2: 安装后,使用library(V.PhyloMaker2)命令来加载这个包。

  3. 数据预处理: 根据你的数据类型和格式,使用相应的函数进行数据导入和预处理。例如,如果你的数据是fasta格式的序列文件,可以使用read.FASTA()函数将其读入R。

    # 导入数据:首先,你需要将你的序列数据导入到R中。这通常是以fasta或 nexus格式存储的。
    library(ape)
    sequences <- read.fasta("your_file.fasta")#数据清理:检查并处理缺失数据、异质性(例如,核苷酸替换)、和错误。
    # 查看是否存在任何缺失数据
    sum(is.na(sequences))# 如果存在缺失数据,可以考虑删除含有缺失数据的行
    sequences <- sequences[!apply(sequences, 1, function(x) any(is.na(x))), ]# 或者用某种方法填补缺失数据(例如,通过平均或中位数)
    sequences[is.na(sequences)] <- median(sequences, na.rm = TRUE)
  4. 多重比对: 使用muscle()或其他比对函数对序列进行比对。

    #序列对齐:对于DNA或蛋白质序列,你需要进行序列对齐。
    aligned_sequences <- muscle(sequences)#转换为距离矩阵:将对齐后的序列转换为距离矩阵,这通常是后续构建系统发育树的步骤。
    dist_matrix <- dist.dna(aligned_sequences)
  5. 进化树构建: 使用build.tree()或其他相关函数,根据你的数据和研究目标选择合适的树构建方法。

    # 假设您已经有了一个包含序列数据的数据框df,并且列名是物种名称
    # df <- data.frame(sequence1, sequence2, ..., sequenceN)
    # 或前面的 data_matrix# 使用build.tree()函数构建进化树
    # 这里的参数是假设的,实际参数需要参考V.PhyloMaker包的文档
    tree <- build.tree(data = df(或data_matrix), seq_type = "dna",   # 数据类型,可以是"dna"、"rna"或"protein"method = "neighbor_joining",  # 构建树的方法,例如"neighbor_joining"(邻接法)或"maximum_likelihood"(最大似然法)distance_method = "kimura")  # 距离计算方法,例如"kimura"(金氏距离)
  6. 进化树优化: 对构建的初步树进行优化,例如使用optimize.tree()函数。

    # 假设你已经使用 build.tree() 建立了一个决策树模型
    # 假设 tree_model 是你建立的模型# 查看建立的树的概况
    summary(tree_model)# 根据交叉验证选择最佳的剪枝参数
    prune_model <- prune.tree(tree_model)# 查看剪枝后的树的概况
    summary(prune_model)# 如果需要,你可以根据需要进一步调整剪枝参数
    
  7. 进化树可视化: 使用plot.tree()函数将进化树可视化,并通过调整各种参数来定制图形。

    # 可视化决策树并调整参数
    plot(tree_model, type = "uniform", fsize = 0.8, cex = 0.8, label = "all")# 添加各种参数以定制图形
    plot(my_tree,type = "fan",       # 树的类型,可以是"phylogram"(分支长度代表进化时间)、"cladogram"(所有分支长度相等)或"fan"(扇形树)show.tip.label = TRUE,  # 是否显示叶节点的标签edge.width = 2,      # 分支线的宽度edge.color = "black",   # 分支线的颜色tip.color = "blue",    # 叶节点的颜色no.margin = TRUE,    # 是否移除图形边框cex = 0.8,           # 标签的字体大小font = 2,            # 标签的字体类型main = "My Evolutionary Tree",  # 图形的标题sub = "Customized with plot() function")  # 图形的副标题
  8. 树形数据分析: 根据你的研究问题,选择相应的函数进行树形数据分析,如节点支持度评估、分支长度分析等。

    # 安装并加载相关包
    install.packages("ape")
    install.packages("phytools")
    library(ape)
    library(phytools)# 假设 tree 是你的树形数据# 计算节点支持度
    bootstrap_tree <- bootstrap.phylo(tree, FUN = your_function_for_tree, B = 100)  # your_function_for_tree 是用于估计树的函数# 生成共识树
    consensus_tree <- consensus(bootstrap_tree)# 计算树的相似性矩阵
    coph_matrix <- cophenetic(tree)# 绘制共演化历史图
    cophyloplot(tree1, tree2)
    

补充分析示例:

树形数据分析可以使用R中的多个包来实现,例如apephangornggtree等。下面是一个简单的示例代码,使用了ape包来进行树形数据分析。

首先,我们需要安装并加载ape包:

install.packages("ape")
library(ape)

接下来,我们可以根据需求读取树形数据。假设我们有一棵简单的进化树,包含5个物种,并且我们想要计算节点的支持度值:

# 创建一个简单的进化树
tree <- rtree(5)# 计算节点的支持度值
supports <- node.depths(tree)

接下来,我们可以绘制树形图,并标记节点的支持度值:

# 绘制树形图
plot(tree, show.node.label = TRUE)# 标记节点支持度值
nodelabels(round(supports, 2), bg = "white")

要分析分支长度,我们可以使用cophenetic.phylo()函数计算树的协同形态矩阵,然后使用plot()函数绘制分支长度图:

# 计算协同形态矩阵
cophenetic_matrix <- cophenetic(tree)# 绘制分支长度图
plot(cophenetic_matrix, main = "Branch Lengths", xlab = "Pairwise Distances")

相似工具包S.PhyloMaker

S.PhyloMaker的介绍和使用看这里:种系进化树分析和构建工具R工具包S.phyloMaker的介绍和详细使用方法-CSDN博客

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/297088.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

互联网加大赛历届作品_匠心筑梦 ——家具漆服务系统2020年度涂装大赛成功举办...

时进九月&#xff0c;秋分将至&#xff0c;天气逐渐微凉&#xff0c;展辰新材家具漆2020年度涂装大赛进入最终评审阶段。大赛历时3个月&#xff0c;国内5大区域一百多位技术能手提交了参赛作品&#xff0c;经过层层筛选&#xff0c;最终入围作品175件&#xff0c;是历届参赛作品…

面对女人的喋喋不休时,男人内心的真实想法

1 笑而不语2 击了掌你就是我的狗了。。。。3 每次看猫咪背对着自己舔毛&#xff0c;总有一种饿的感觉。。4 你看对面的那个小朋友可不可爱&#xff0c;快去抱抱他~5 猫咪三明治&#xff0c;感觉很好吃的样子&#xff01;6 看剧的时候真的不能开弹幕&#xff01;你点的每个赞&am…

Dapr牵手.NET学习笔记:绑定

绑定有点像订阅发布&#xff0c;但又不一样&#xff0c;绑定更简单&#xff0c;绑定输出&#xff08;调用方&#xff09;-绑定输入&#xff08;被调用方&#xff09;。本例是用docker compose编排&#xff0c;并且用rabbitMQ来支持&#xff0c;因为rabbitMQ支持输入和输出绑定。…

pytorch 加载模型_福利,PyTorch中文版官方教程来啦(附下载)

PyTorch 中文版官方教程来了。PyTorch 是近期最为火爆的深度学习框架之一&#xff0c;然而其中文版官方教程久久不来。近日&#xff0c;一款完整的 PyTorch 中文版官方教程出炉&#xff0c;读者朋友从中可以更好的学习了解 PyTorch 的相关细节了。教程作者来自 pytorchchina.co…

为什么睡觉时身体突然抖一下?答案吓到我了!

全世界只有3.14 % 的人关注了爆炸吧知识不知你有没有过这样的经历&#xff1a;即将进入甜美梦乡&#xff0c;突然就像触电般地抖了一下&#xff0c;或者不由自主猛地踢一下&#xff0c;瞬间被惊醒&#xff01;这到底是啥情况&#xff1f;网上的答案五花八门&#xff1a;有人说缺…

Dapper防sql注入,同一条SQL支持多种数据库

前言防SQL注入&#xff0c;常用的方案是使用Dapper执行SQL的参数化查询。例如&#xff1a;using (IDbConnection conn CreateConnection()) {string sqlCommandText "SELECT * FROM USERS WHERE IDID";Users user conn.Query<Users>(sqlCommandText, new { …

基于css3的鼠标滑动按钮动画之CSS--续

2019独角兽企业重金招聘Python工程师标准>>> btn2 /*--按钮1--*/ .container1 {width:200px;display:block;margin:20px auto;position:relative;font-family:droid arabic kufi; } .con_down1 {display:block;cursor:pointer;background-color:#F6EB96;width:190px…

rabbitmq入门_Rabbit MQ 入门

Rabbit MQ是一个通用的消息中间件&#xff0c;支持AMQP&#xff0c;STOMP&#xff0c;MQTT等多种协议安装#在OSX下可以使用如下命令来安装 rabbitmqbrew install rabbitmq基本命令#ls -al ~/rabbitmq/3.7.14/sbin/total 1104drwxr-xr-x 10 jet admin 320 May 19 14:35 .d…

94年出生,她们如今都是985高校博士生导师!

全世界只有3.14 % 的人关注了爆炸吧知识鱼羊 萧萧 发自 凹非寺量子位 报道 | 公众号 QbitAI94年出生、博士研究方向与材料相关、目前都成了985重点高校的博士生导师。拥有相同经历的两个女生&#xff0c;概率有多大&#xff1f;就在今年&#xff0c;26岁的夏娟和李晟曼&#xf…

IBM沃森为存储系统开发人员带来的启发

前一段时间&#xff0c;IBM 沃森参加了CBS的益智节目《危险边缘》(Jeopardy)&#xff0c;这是他在全国观众面前首次亮相。确切地说&#xff0c;站在中央舞台选手答题台后面的IBM沃森实际上 是不断闪烁的虚拟头像。尽管如此&#xff0c;摆在沃森面前的答题按钮却是如假包换的&am…

01Prism WPF 入门实战 - 项目准备

1.概要这一系列将进行PrismWPF技术的实战讲解。实战项目内容选型为Email邮件收发的客户端&#xff08;WeMail&#xff09;&#xff0c;项目结构简单方便大家理解。相关技术&#xff1a;C#、WPF、Prism软件开发环境&#xff1a;VS2019 、 .NET5 、 windows11需掌握技能&#xf…

redis查看key的过期时间_面试官:你在Redis中设置过带过期时间的Key吗?

点击上方小伟后端笔记关注公众号每天阅读Java干货文章熟悉Redis的同学应该知道&#xff0c;Redis的每个Key都可以设置一个过期时间&#xff0c;当达到过期时间的时候&#xff0c;这个key就会被自动删除。在为key设置过期时间需要注意的事项1、 DEL/SET/GETSET等命令会清除过期时…

Hadoop学习系列之PageRank

昨晚上不想做其他的事&#xff0c;突然想起来好久都没更新博客了&#xff0c;shell也差不多学完了&#xff0c;只不过学习的时候都是只带着书出去了&#xff0c;改天总结总结。Hadoop么&#xff0c;黄宜华老师讲完了&#xff0c;自己也马马虎虎快学完了&#xff0c;也是没总结&…

您的屁股发热严重,请降温后使用。

▲ 点击查看不是坐在办公椅上的屁股都渴望自由&#xff0c;而是——最近天越来越热&#xff0c;屁股捂在椅子上&#xff0c;既不散热也不排汗&#xff0c;比戴口罩闷出痱子还难受&#xff0c;实在是坐不住。。最尴尬的就是站起来裤子时常黏在屁股缝里&#xff0c;难不成每次还要…

C# 正则表达式编写及验证方法

01—前言正则表达式应用很广泛&#xff0c;应该大多人都接触过了&#xff0c;这个语法规则既多又凌乱&#xff0c;每次用的时候都得重新看一遍语法&#xff0c;真的是让人头疼啊&#xff01;但是实际上我们并不要掌握很多的符号用法规则&#xff0c;牢记最常用的几个就能应付很…

Domino Web开发规则之二:DOMINO与开发相关的管理规范

1.服务器HTTP优化设置 调整活动线程数&#xff0c;HTTP服务器可以同时处理的请求数&#xff0c;而非连接数、会话数 单CPU服务器<64 多CPU服务器<80 并发运行Web代理 确保Web代理是线程安全的情况下&#xff0c;可以启用来提高性能。 服务器文档-> Internet协议 ->…

老师一定没有教!9个数学速算技巧!看到第1个就跪了!

全世界只有3.14 % 的人关注了寄语 与世界上的其他学科一样&#xff0c;数学也有自己的奇妙之处。由于各种各样的原因&#xff0c;老师没有将这些秘密告诉我们。下面就是一些有趣的数学技巧&#xff0c;看完后你也许会爱上这门学科。▼不管是几个1的平方&#xff0c;都是有规律…

愚蠢的领导才会用程序员祭天!!

“人非圣贤&#xff0c;孰能无过网络上频频流出程序员误删线上数据库&#xff0c;造成XX公司XX万损失的新闻&#xff0c;进而牵扯出“杀一个程序员祭天”的妙语&#xff0c;虽然杀一个程序员比杀一个进程要难很多&#xff0c;但也充满着讽刺的味道。程序员每天都往返于代码和服…

一般将来时语法课教案_速看,如何在考场写出一篇脱颖而出的教案

面试最重要的考试环节就是试讲&#xff0c;而想要进行一次优秀的表演&#xff08;试讲&#xff09;&#xff0c;则需要一个过硬的剧本&#xff08;教案&#xff09;&#xff0c;那我们如何写出一篇完美教案呢&#xff1f;今天上元教师就教给大家10分钟速写一份教案的方法&#…

你的感情路不顺么?

1 你的感情路不顺么&#xff1f;2 扎心了&#xff01;3 哈哈哈哈没毛病&#xff0c;差一点就信了&#xff01;4 这个反杀&#xff0c;满分&#xff01;5 一只认认真真学揉面的猫咪6 当你妈说你该出门运动下的时候。。 你点的每个赞&#xff0c;我都认真当成了喜欢