OmicCircos的一些简单的使用心得

对于circos图想必大家并不陌生,但是一般自己安装circos软件着实有点儿麻烦。而在R中的OmicCircos包却能让你很好的解决问题,可以去画一些简单的而又不落俗套的circo

1,首先是下载

## try http:// if https:// URLs are not supported

source("https://bioconductor.org/biocLite.R")

biocLite("OmicCircos")2,下载后之后在R中载入包

library(OmicCircos)

3,对OmicCircos的一些函数的了解

3.1函数的认识

OmicCircos主要有三个函数,分别是sim.circos,segAnglePo和circos。sim.circos是用来模拟一些数据的,如果我们用自己的数据,这个函数一般我们用不到。segAnglePo是将线性的数据结构转化为环形,circos则主要是画图。大多数情况下,我们只用得着后两种函数。

3.2数据的结构

一般画图我们只需要给予一个table就可以,但是这个table前三列必须是位置。如果你用染色体做外圈,那么应该为

Chr Start End 。。。。

chr1 500 1000

chr2 500 1000

chr3 500 1000

如果你用基因作为外圈,那么结构应该为

Genename Start End。。。

Tp53 5000 10000

EGFR 5000 10000

在OmicCircos的给出的例子,seg.f

chrom chromStart chromEnd name gieStain
chr1 2300000 5300000 p36.32 gpos25
chr1 5300000 7100000 p36.31 gneg
chr1 7100000 9200000 p36.23 gpos25
chr1 9200000 12600000 p36.22 gneg
chr1 12600000 16100000 p36.21 gpos50

给出的seg.v,也就是mapping数据表

chr start end value 。。。

chr1 10000 10090 2.5 。。。

3.3数据的转换

seg.name ← paste ( " chr" , 1:seg.num , sep="")###seg.num就是你染色体的数目
db← segAnglePo ( seg.f , seg=seg.name )###seg.f为线性的数据,如果是染色体画圈,前三列必须是‘chr,start,end‘;如果是基因,则为‘genename,start,end‘。而且seg.f必须是5列以上,除了前三列,后面可以随便加至少2列数据。

3.4数据画图

使circos这个函数,如

circos (R=400 , type=" chr " , cir=db , col=colors, print.chr.lab=TRUE, W=4, scale=TRUE)####画出最外圈的结构

circos(R=360,cir=db ,W=40,mapping=seg.v,col.v=8, type="s" ,B=FALSE, col=colors[9] , lwd =2, scale=TRUE) ###cir表示用说y图的线宽;scale就是你mapping数值的标尺,显示或者不显示。
circos(R=320,cir=db ,W=40,mapping=seg.v,col.v=3, type="sv",B=TRUE, col=colors[7] ,scale=TRUE)####

4.关于OmicCircos的一些信息

自带人hg19,hg18的染色体数据,还有mm10等。OmicCircos不仅可以画一些简单的图,也能画SV的,还需要自己去paper上自己看更多的细节。

时间: 2024-10-11 16:20:47

OmicCircos的一些简单的使用心得的相关文章

简单的学习心得:网易云课堂Android开发第三章自定义控件

这一章分三部分: (1)自定义控件:老师先简单讲解了一些细节,如为什么不用px,而要用dp,只因机型的屏幕分辨率不同,用px会导致差异太大.然后演示了制作自定义控件的步骤,先在xml文件中添加对应的自定义View类,并给出画布大小,然后在.java中onDraw,值得注意的是onDraw中的坐标都是相对于画布(即xml中定义的)而言的,所以用具体数值去绘制view,很容易出现问题. (2)Fragment:第一点是两种设置fragment的方法.第一种是在.xml文件中,这是静态设置,不能交互:

简单的学习心得:网易云课堂Android开发第四章服务、广播与酷特性

1.Service (1)含义:Service是在后台长时间运行的应用组件,所在线程是main线程. (2)使用方式:startService能够执行onStartCommand中操作,bindService能够将Service的实例对象返回给Activity,从而能够在Activity中使用Service的方法.而用new方法创建Service对象,没有报错,但无法成功调用Service的方法. (3)注意:在onStartCommand中,老师建议我们将返回值改为START_NOT_STIC

简单的学习心得:网易云课堂Android开发第五章SharedPreferences与文件管理

一.SharedPreferences (1)SharedPreferences能够用来保存一些属于基本数据类型的数据. (2)保存数据,删除数据都是由SharedPreferences的内部接口Editor负责. (3)通过这种方式保存的数据在root->data->包名->shared-prefs文件夹下的xml文件中. 二.文件管理 (1)创建文件的方式:new File(...)和openFileOutput,通过new File(...)的方式创建文件,还需要调用createN

简单的学习心得:网易云课堂Android开发第六章SQLite与ContentProvider

一.SQLite 1.基本操作: (1)创建数据库:在SQLiteOpenHelper的子类构造器中创建. (2)创建表:在SQLiteOpenHelper的子类onCreate方法中,调用execSQL方法实现. (3)增.删.改.查:分别调用SQLiteDatabase的insert.delete.update.query方法,又或者在execSQL方法中使用SQL语句来实现相应功能. 2.利用事务对数据库进行操作 (1)先开启事务db.beginTransaction(),接着执行操作 (

枚举进程再来两弹

看了刚出几个博友的博客,感觉人家的量大,详细,干货量实足啊, 所以我就把另外两种常见的枚举进程的方法简单说下心得, 一个是EnumProcesses和CreateToolhelp32Snapshot系列的Tool help API的 Process32First和Process32Next函数完成列举进程. 这两种都是比较简单实用的 https://github.com/Arsense/WindowsCode 蛮简单的  需要编译好的源码的 支持vs2015,低版本的VS想编译 简单右键工程 属

【转】七年阿里老人谈新人成长

先做事再做人 新同学常有这些困惑,到底是做技术还是做业务?是每天加班加点写代码还是空点时间来看书学习?是先追求广度还是先钻研一门技术?是两耳不闻窗外事一心埋头做事还是和同事打成一片参加各种活动?这些问题也曾困扰着我. 我觉得新人到一个团队最关键的是在前3年建立个人品牌,为什么?当别人对你完全不了解,你做的第一件事非常重要.你给别人留下的印象会蔓延到整个公司,所以第一件事一定要做好.怎么才叫做好?在职场中,并不需要时刻比拼智商,大部分情下,做事认真一点.考虑更深入一点.积极主动一点,已经好于70%

游戏网络编程(三)——WebSocket入门及实现自己的WebSocket协议

(一)WebSocket简介 短连接:在传统的Http协议中,客户端和服务器端的通信方式是短连接的方式,也就是服务器端并不会保持一个和客户端的连接,在消息发送后,会断开这个连接,客户端下次通信时,必须再建立和服务器的新连接,这就是短连接.在短链接的情况下,客户端必须不停的主动发起请求,而服务器始终被动的响应请求,来推送回数据.这种方式用到游戏开发中,显然是不适合的. 长连接:那么与之相对的就是长连接了.在长连接的情况下,客户端和服务器端始终保持一条有效的连接,那么客户端并不需要不停的主动发送消息

Alpha第五天

Alpha第五天 听说 031502543 周龙荣(队长) 031502615 李家鹏 031502632 伍晨薇 031502637 张柽 031502639 郑秦 1.前言 任务分配是VV.ZQ.ZC负责前端开发,由JP和LL负责建库和服务器. 界面开发的教辅材料是<第一行代码>,利用Android Studio进行开发. 2.燃尽图 3.项目进展及安排 日期 已完成 未完成 10月31日 数据库完成+数据库测试:登录界面+注册界面+密码找回 11月 1日 服务器初步搭建:写日记界面+我的

七年阿里老人谈新程序员成长:先做事,在做人!

新程序员常有这些困惑,到底是做技术还是做业务?是每天加班加点写代码还是空点时间来看书学习?是先追求广度还是先钻研一门技术?是两耳不闻窗外事一心埋头做事还是和同事打成一片参加各种活动?这些问题也曾困扰着我. 我觉得新人到一个团队最关键的是在前3年建立个人品牌,为什么?当别人对你完全不了解,你做的第一件事非常重要.你给别人留下的印象会蔓延到整个公 司,所以第一件事一定要做好.怎么才叫做好?在职场中,并不需要时刻比拼智商,大部分情下,做事认真一点.考虑更深入一点.积极主动一点,已经好于70% 的同学了