methylKit 进行差异甲基化分析

本文介绍了R包methylKit在差异甲基化分析中的应用,包括读取原始数据、合并样本数据和执行差异分析的步骤。通过methylKit可以处理多种格式的甲基化测序数据,并提供了方便的接口进行差异甲基化位点和区域的注释。
摘要由CSDN通过智能技术生成

methylKit 是一个用于分析甲基化测序数据的R包,不仅支持WGBSRRBS和目的区域甲基化测序,还支持oxBS-sq, TAB-seq等分析5hmc的数据。 其核心功能是差异甲基化分析和差异甲基化位点和区域的注释。

安装过程如下:

source(“http://bioconductor.org/biocLite.R“)
biocLite(“methylKit”)

推荐使用最新版本的R进行安装,这样可以使用最新版本的methylKist
利用methylKit 做差异分析包括3步

1. 读取原始数据

每个样本一个原始数据,methylKit支持两种格式的methylation calling文件

  • 纯文本格式

    内容如下

每一行是一个甲基化位点,coverage 代表覆盖这个位点的reads数,freqC 代表甲基化C的比例,freqT 代表非甲基化C的比例。这种纯文本格式内容非常直观,文件大小相比bam 文件小很多,读取的速度更快。
纯文本格式的读取过程如下

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值