CHIP-SEQ 芯片分析时,对于来自重复实验的数据,怎样进行MACS peaks calling 分析?
解决方法:
将bowtie得到的bam文件,转换为bed文件,并使用cat指令,将其所有重复实验文件合并,将合并文件作为输入文件,进行peaks calling。
参考链接:MACS官网
http://liulab.dfci.harvard.edu/MACS/00README.html
Notes:
3) For the experiment with several replicates, it is recommended to concatenate several ChIP-seq treatment files into a single file. To do this, under Unix/Mac or Cygwin (for windows OS), type:
$ cat replicate1.bed replicate2.bed replicate3.bed > all_replicates.bed
附录:
我的数据
样本 | 数据 |
重复样本1 | SRR6168965 |
重复样本2 | SRR6168971 |
重复样本3 | SRR6168972 |
过程代码
(base) s45@HP45:~/chip-seq_analysis/1-bowtie/2-lung$ bamToBed -i SRR6168965.bam >SRR6168965.bed
(base) s45@HP45:~/chip-seq_analysis/1-bowtie/2-lung$ bamToBed -i SRR6168971.bam >SRR6168971.bed
(base) s45@HP45:~/chip-seq_analysis/1-bowtie/2-lung$ bamToBed -i SRR6168972.bam >SRR6168972.bed
(base) s45@HP45:~/chip-seq_analysis/1-bowtie/2-lung$ cat SRR6168965.bed SRR6168971.bed SRR6168972.bed >lung_all.bed
lung_all.bed是最终三个样本数据合并的文件。
MACS分析指令:
(base) s45@HP45:~/chip-seq_analysis/2-macs2_change$ macs2 callpeak -t /home/s45/chip-seq_analysis/1-bowtie/2-lung/lung_all.bed -f BED -g mm -n lung_all -B -p 0.0001 --outdir /home/s45/chip-seq_analysis/2-macs2_change/
但是,后来实践发现,以上方法分析得到的peaks序列与原文献数量相比,大大增加,原因不明。
觉得不对。
后来发现一篇博文,解决了我的疑惑,请有相同疑惑的同学可以参考:
参考链接(https://www.jianshu.com/p/d8a7056b4294)使用bedtools取共有的overlap区,猜得到与文献相似的结果。猜想这是合适的对于重复样本取overlap的方式。
CHIP-SEQ 芯片分析时,对于来自重复实验的数据,怎样进行MACS peaks calling 分析?相关推荐
- 生物信息学入门 GEO芯片数据差异表达分析时需要log2处理的原因
首先借用一张图,通常使用limma处理时,需要经过log2后的矩阵作为表达矩阵输入.根据log2FC的定义,这个数字表示变化倍数经过log2后的一个值,比如log2FC=1,则变化为2倍:log2FC ...
- ICC 图文学习——LAB6:Chip Finishing 芯片完成
布线工作完成后,主体工作已经完成,chip finishing 阶段主要是做检查设计.提高良率和物理验证工作,主要内容如下: ·DRC 和 LVS检查 ·进行连线扩展和加宽连线来减少关键区域 ·修 ...
- 新的微芯片MCU增加了来自外部闪存的安全引导保护
新的微芯片MCU增加了来自外部闪存的安全引导保护 New Microchip MCU Adds Secure Boot Protection from External Flash 对于从外部SPI闪 ...
- Openstack数据库初始化时为何要重复授权
问题; 为什么Openstack在创建数据库时要同时赋予用户在%和localhost登陆的权限而在MySQL中的%已经包含了localhost? 要回答这个问题我们可以先看看不这样做会怎样; open ...
- c语言程序执行时无法输入字符串,C语言程序设计中键盘输入数据的方法分析
1. 引言 C语言程序设计是一门实践性很强的语言课程,在程序设计过程中经常需要从键盘输入所需要的数据.C语言没有专门的输入语句,可以调用C语言编译系统提供的函数库中的库函数来实现.有些程序虽然调试成功 ...
- 白话Elasticsearch44-深入聚合数据分析之案例实战_颜色+品牌下钻分析时按最深层metric进行排序
文章目录 概述 案例 概述 继续跟中华石杉老师学习ES,第44篇 课程地址: https://www.roncoo.com/view/55 案例 需求: 颜色+品牌下钻分析时按最深层metric进行排 ...
- 全球与中国光电通信芯片市场数据专项调研及竞争格局分析报告2022-2028年版
全球与中国光电通信芯片市场数据专项调研及竞争格局分析报告2022-2028年版 [报告目录]: 2021年全球光电通信芯片市场销售额达到了 亿美元,预计2028年将达到 亿美元,年复合增长率( ...
- r roc函数_一棵树专栏 | ROC分析时一定要告诉R分析谁
两个知名R包产生的矛盾嗨!大家好,我是一棵树,这是我第一次在解螺旋发文,还是蛮激动的.下面就开始吧!今天的主题是:ROC分析时一定要告诉R分析谁 用到的软件是:R语言 用到的R包是:pROC和ROCi ...
- 在MAC下安装一些软件时提示来自身份不明开发者
在MAC下安装一些软件时提示"来自身份不明开发者",其实这是MAC新系统启用了新的安全机制. 默认只信任 Mac App Store 下载的软件和拥有开发者 ID 签名的应用程序. ...
- mysql vba日期空值_使用SQL(VBA/ADO)命令填充来自重复记录的相应非空值的NULL字段...
我有一个数据库,其中有数百个记录已被复制.但是,重复的信息在所有字段中都不相同.对于任何两行,第一行将包含一些字段中的信息,而重复行的字段为空;但对于其他字段,重复(第二)行将包含信息,而第一行的字段 ...
最新文章
- pat1094. The Largest Generation (25)
- 用python写脚本看什么书-终于知晓python编写脚本入门教程
- maven安装_如何从官网下载Maven与安装Maven
- 67、Fragment实现Tab
- vue配置二级目录vue-axios跨域办法谷歌浏览器设置跨域
- Scrum 项目7.0
- 蓝桥杯 BEGIN-3 入门训练 圆的面积
- ansys16.0安装教程
- 【axios】get和post请求用法
- 字母组合---排列组合问题
- windows pe安装系统
- socket多线程图形化界面聊天室实例
- 不动产登记工作统计难 要靠三维不动产管理系统
- 腾讯云TVP大数据之约活动归来
- android 模拟器声音设置,逍遥安卓模拟器运行游戏过程中没有声音怎么解决
- 生活随笔与读书笔记20140302
- RHEL8.0快速入门系列笔记--理论知识储备(一)
- 如何开发一个企业级的 Go 应用?
- 华为USG6000系统进入诊断视图调节风扇转速配置
- 百度、高德离线地图SDK开发工具,局域网内离线地图开发环境