r语言remarkdown展示图_使用R语言包circlize可视化展示blast双序列比对结果
circlize这个包还挺强大的，R语言里用来画圈图还挺方便的。 今天这篇文章记录用circlize这个包画圈图展示blast双序列比对结果的代码
植物线粒体基因组类的文章通常会分析细胞器基因组间基因转移情况，基本的分析方法就是blast比对。可视化展示可以选择用这个圈图来做
首先是使用blast建库比对
makeblastdb -in mt.fasta -dbtype nucl -out mt
blastn -query cp.fasta -db mt -outfmt 6 > output.txt
首先是准备数据用来画最外圈用来展示两条序列的部分
df<-data.frame(chr=c(rep("chloroplast",2),rep("mitochondrial",2)),x=c(1,131478,1,444567),y=c(0,1,0,1))> dfchr x y
1 chloroplast 1 0
2 chloroplast 131478 1
3 mitochondrial 1 0
4 mitochondrial 444567 1
然后是读入blast的输出结果
df1<-read.csv("output6.txt",stringsAsFactors = F,header=F,sep="t")
作图用到的代码
library(circlize)
library(RColorBrewer)
library(ComplexHeatmap)
col<-RColorBrewer::brewer.pal(6,"Paired")
circos.par("start.degree" = 130)
circos.initialize(factors = df$chr,x=df$x)
circos.trackPlotRegion(factors = df$chr,y=df$y,panel.fun = function(x,y){circos.axis()},track.height = 0.1)
highlight.sector(sector.index = "chloroplast",col=col[1])
highlight.sector(sector.index = "mitochondrial",col=col[2])
circos.text(x=70000,y=0.5,labels = "chloroplast",sector.index = "chloroplast")
circos.text(x=220000,y=0.5,labels = "mitochondrial",sector.index = "mitochondrial",facing = "outside")
col_fun = colorRamp2(c(70,90,100),c("green", "yellow", "red"))
for (i in 1:13){x<-sort(c(df1[i,8],df1[i,7]))y<-sort(c(df1[i,10],df1[i,9]))z<-df1[i,3]circos.link("chloroplast",x,"mitochondrial",y,col=add_transparency(col_fun(z)))
}
circos.clear()
lgd_links = Legend(at = c(70, 80, 90, 100), col_fun = col_fun, title_position = "topleft",title = "identity(%)")
lgd_list_vertical = packLegend(lgd_links)draw(lgd_list_vertical, x = unit(10, "mm"), y = unit(10, "mm"), just = c("left", "bottom"))
新学到的两个知识点
调整整体的角度
circos.par("start.degree" = 130)
调整用来表示染色体的外圈粗细
circos.trackPlotRegion(factors = df$chr,y=df$y,panel.fun = function(x,y){circos.axis()},track.height = 0.1)
画图的时候可以加一个track.height参数
遇到的问题是
调整外圈的刻度，现在展示的有点多，我想增大间隔，减少展示的数字，暂时不知道如何实现。
添加图例参考了文章
https://jokergoo.github.io/circlize_book/book/legends.html
欢迎大家关注我的公众号 小明的数据分析笔记本
相关文章：
吴恩达说“将引领下一波机器学习技术”的迁移学习到底好在哪？
AI技术年度盛会即将开启！11月8-9日，来自Google、Amazon、微软、Facebook、LinkedIn、阿里巴巴、百度、腾讯、美团、京东、小米、字节跳动、滴滴、商汤、旷视、思必驰、第四范式、云知声等企业的技术大咖将带来工业界AI应用的最新思维。 如果你是某个AI技…
获取系统相关属性
uname -a uname -r cat /etc/issue cat /etc/redhat-release #进入centos后忘记自己安装的是6还是7，可以从这里查看 转载于:https://www.cnblogs.com/biaopei/p/7730684.html
双水泵轮换工作原理图_「物业管理工作」水泵维护保养规程
1 目的为规范水泵维护保养程序，特制定本规程。2 适用范围适用于公司管辖范围内水泵(包括生活水泵、消防水泵及集水坑潜水泵)及其相关设施的维护保养。3 职责3.1 给排水维修工负责水泵的维护保养。3.2 值班锅炉工或电梯工负责晚间协助对水泵房设备运转情况的巡查。3.…
Dropbox如何使用机器学习从数十亿图片中自动提取文字
参加 2018 AI开发者大会，请点击 ↑↑↑作者 | Leonard Fink译者 | 刘旭坤整理 | Jane出品 | AI科技大本营【导读】提到 Dropbox，大家可能都知道这是一个文件同步、备份、共享的云存储软件。其实 Dropbox 可以实现的功能远不止这些。今天就为大家介绍 Dr…
001_性能压测工具
一、 Jmeter
SharePoint Server 2007 页面模型
虽然SharePoint Server 2007使用了ASP.NET 2.0的基础页面模型，SharePoint页面基本上也是基于标准的aspx技术来构建，但SharePoint Server 2007的页面模型仍然要比普通的ASP.NET应用复杂很多。对于一个SharePoint开发人员（和设计人员࿰…
oracle存储过程的常用语法
记录一下，主要方便自己以后查找： createorreplaceproceduresp_test(aa invarchar,bb innumberdefault0) isv_pos varchar(50);v_num number;begin---字符串查找和替换selectSUBSTR(PARAMETER, 1, INSTR(PARAMETER,branchId)-2) ||SUBSTR(PARAMETER, INST…
机器学习大神迈克尔 · 乔丹：我讨厌将机器学习称为AI
AI技术年度盛会即将开启！11月8-9日，来自Google、Amazon、微软、Facebook、LinkedIn、阿里巴巴、百度、腾讯、美团、京东、小米、字节跳动、滴滴、商汤、旷视、思必驰、第四范式、云知声等企业的技术大咖将带来工业界AI应用的最新思维。 如果你是某个AI技…
es6箭头函数(=)与展开特性运算符(...)的使用
2019独角兽企业重金招聘Python工程师标准>>> 箭头函数的几个写法 //当含有固定个参数的时候，es5写法 function test(a,b){ ..}//es6写法(es6参数可以设置默认值) let test (a,b) > {..}当函数体中仅有一行的时候，可以省略掉return与外部的…
华为atn980传输设备_在头发丝中实现每秒1000张高清DVD传输
这是近日央视《对话》栏目“超高清视频带来的产业诱惑”节目中，提到的随着4K、8K超高清视频发展，顺应时代发生的变化之一。超高清视频将为消费者带来全新的、超乎想象的体验。在这其中，黑科技的涌现更是对产业快速发展起着不可磨灭的推动作用…
李飞飞、吴恩达、Bengio等人的15大顶级深度学习课程，你收集全了吗？
AI技术年度盛会即将开启！11月8-9日，来自Google、Amazon、微软、Facebook、LinkedIn、阿里巴巴、百度、腾讯、美团、京东、小米、字节跳动、滴滴、商汤、旷视、思必驰、第四范式、云知声等企业的技术大咖将带来工业界AI应用的最新思维。 如果你是某个AI技…
2 拖拽模型_3Dmax插件 | 螺丝钉、挖洞、网格模型等可以分分钟搞定
今天给大家带来3Dmax插件Mesh Insert v1.14Mesh Insert是一款3dsmax建模插件，灵感来自Polystein for Modo。它通过替换库中模型资源快速插入到模型表面。非常适合硬表面和有机建模!1大纲1、安装说明2、插件功能3、视频讲解2过程步骤1安装说明1、关闭max，…
谈谈职业规划——CSDN对我的采访
职业规划就像软件工程电信、银行等行业一直是许多人非常向往的工作单位，清差厚禄，旱涝保收，陈皓却不以此为然。所以当记者采访他的时候，他连用了两个“最”字来形容他离开银行的成就感。陈皓毕业后的前两年就职于云南省工商银行&a…
最实用的机器学习算法优缺点分析，没有比这篇说得更好了
AI技术年度盛会即将开启！11月8-9日，来自Google、Amazon、微软、Facebook、LinkedIn、阿里巴巴、百度、腾讯、美团、京东、小米、字节跳动、滴滴、商汤、旷视、思必驰、第四范式、云知声等企业的技术大咖将带来工业界AI应用的最新思维。 如果你是某个AI技…
南大电子机器人入驻云南_云南大学外语协会 英文歌曲演唱比赛
云南大学外语协会英文歌曲演唱比赛顺利落幕2020年11月8日14：30，云南大学外语协会英文歌曲演唱比赛在文汇楼2116顺利开展。比赛流程简洁明了：参赛同学完成抽签之后，在主持人的引导下，按照流程有序进入指定位置等待表演&…
Spring BeanFactory与FactoryBean的区别及其各自的详细介绍于用法
1. BeanFactory BeanFactory，以Factory结尾，表示它是一个工厂类(接口)，用于管理Bean的一个工厂。在Spring中，BeanFactory是IOC容器的核心接口，它的职责包括：实例化、定位、配置应用程序中的对象及建立这些对…
智“逗”冰雪寒
和小K帮客户做完直播准备，他说一会儿要接冰雪寒冰－雪－寒多伟大的名字——我来51看的第一篇博客就是他的《网管日志》那曾经是51最火的文章之一<?xml:namespace prefix o ns "urn:schemas-microsoft-com:office:office" />…
装饰模式（Decorator Pattern）
装饰模式（Decorator Pattern）一句话 继承一个抽象类，加上自己的特点，然后再使自己也变成抽象类，然后让子类继承的模式，叫装饰模式。意图…
python 区域和检索_304. 二维区域和检索（Python）
题目难度：★★★☆☆类型：二维数组方法：动态规划给定一个二维矩阵，计算其子矩形范围内元素的总和，该子矩阵的左上角为 (row1, col1) ，右下角为 (row2, col2)。示例:给定 matrix [[3, 0, 1, 4, 2],[5, 6, 3…
观察者模式(Observer Pattern)（二）：HeadFirst中的气象站的实现
1 观察者模式的原理，首先由一个主题，当主题发送变化的时候，通知该主题的订阅者 按照上面的分析我们来进行设计 1.抽象主题Subject public interface Subject {public void registerObserver(Observer o);public void removeObserver(Observer…
专访图灵奖得主John Hopcroft：中国必须提升本科教育水平，才能在AI领域赶上美国
AI技术年度盛会即将开启！11月8-9日，来自Google、Amazon、微软、Facebook、LinkedIn、阿里巴巴、百度、腾讯、美团、京东、小米、字节跳动、滴滴、商汤、旷视、思必驰、第四范式、云知声等企业的技术大咖将带来工业界AI应用的最新思维。 如果你是某个AI技…
完整复现何恺明ICCV获奖论文结果并开源 !（附论文开源代码）
ICCV 作为计算机视觉的顶级会议，2017年共收到2143篇论文投稿，比上一届ICCV2015的1698篇增加了26.2%。共621篇被选为大会论文，录用比例28.9%；poster、spotlight、oral 比例分别为24.61%、2.61%和2.09%。组委会根据作者署名统计了不…
电子合同的履行_什么是电子合同履行？怎么履行电子合同？
随着互联网产业的发展，许多传统产业都与互联网接轨，人们的传统生活方式也在转变。现在许多人的各种消费都喜欢在方便、快捷的互联网上进行，就连严肃的合同签订也是如此。不过还有很多人不知道什么是电子合同履行，下面律图小编就为…
10月机器学习开源项目Top10
作者 | Mybridge 译者 | 林春眄 整理 | Jane 出品 | AI科技大本营 【导读】过去一个月里，我们对近 250 个机器学习开源项目进行了排名，并挑选出热度前 10 的项目。这份清单的平均 github star 数量高达 1345，涵盖了包括深度学习, Tensorfl…
用SDM架构Cisco IOS ***图文详解全攻略（一）——easy ***
在测试***的过程中发现网上资料少少，于是自己花点力气写几篇普及教程，希望转贴的朋友给与支持，不要忘记署上原创是水煮豆豆的大名，嘿嘿～Cisco给我们提供了管理路由器的很好的安全工具SDM，同时也给我们提供了…
HP 服务器使用 SmartStart CD 引导安装 windows 2008 操作系统
1. 使用SmartStart CD引导安装之前的准备工作：首先，先根据自己的需要配置好RAID，可以使用开机自检检测到阵列卡时按F8进入ORCA的方法配置，也可以使用SmartStart CD中附带的ACU工具配置。其次，需要确认BIOS中启动顺…
parameter缩略语_缩略语
缩略语AccAccelerate(加速)BFBreakerfailure(断路器失灵)BRCBBuferReportControlBloc(有缓存报告控制块)CIDConfiguredIEDDescription(IED实例配置文件)CTCurrenttransformer(电流互感器)DevDevice(设备)ErrError(错误)FstFirst(第一个)GOOSE,GoGenericobjectorientedsubstatio…
[译]怎样用VisualStudio查看非托管代码
(译者:这篇文章作者是一位美国的MVP,这是他的系列文章"Under the cover"的第一篇,文章的本意从最底层的角度来优化代码的性能,并作为阅读作者其他文章的技术基础,这种通过这样的做法虽然初看起来有些过分,但是对读者了解.Net许多底层运作是十分有益的) 我们从使用vi…
盘点互联网大厂AI战略变迁，开发者将怎样pick前进路线？
随着各大企业相继试水“全面 AI”，人工智能在技术落地层面也开始持续深入，泛人工智能时代正在逼近。越来越多的发展趋势表明，未来的人工智能将逐步迈入广泛普及阶段，继而深入影响人类日常的生产生活方式，重塑传统生产结…