专栏首页生物信息学、python、R、linux不同的peak calling软件比较

不同的peak calling软件比较

我一般用MACS2做peak calling,但是不知道效果是不是最好的,去搜了一下,发现14年有一篇文章用DNase-seq的数据比较了主流的几个peak caller的效果。虽然发表在Plos one上,不过还是很有参考意义(A Comparison of Peak Callers Used for DNase-Seq Data)。这篇文章比较了如下四个软件:

image

用 ENCODE中的K562, GM12878 和 HelaS3的DNase-seq的数据从以下几方面去比较这几种软件的效果:

  1. sensitivity 和 specificity

作者先从ENCODE中下载了K562, GM12878 和 HelaS3的几十套转录因子结合(TFBS)的narrow peaks的数据,用BEDOPS取了这些的并集作为比较这几种软件的"reference set"。以下是这四种数据的灵敏度(TPR)和特异度(1-FDR)的情况(其中ZINBA分为两种,一种是ZINBA_N,即输出为narrow peak模式,另一种是ZINBA_B,输出为broad peak):

image

ZINBA_N具有最小的TPR,即可能会错失更多正确的DHS,不过它同时也有最小的FDR。但是ZINBA_B的TPR和FDR都比较低,所以ZINBA_B与reference set和其他的方法差距较大。

2. 找到的peaks数量和peaks的长度以及coverage

image

image

image

3. 调整默认的信号阈值

image

在调整了阈值之后。F-seq的效果得到了提升,但是Hotspot与MACS表现则与之前差不多。

总体来看,对于DNase-seq的数据而言。F-seq的效果来说最好,MACS和Hotsplot也相对不错。但是ZINBA对于DNase-seq的数据无论是在特异度、灵敏度还是在计算时间上都不是很适合。

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • celltalker配体受体互作

    对单细胞RNAseq数据常规的分析经常会有对细胞间的相互作用 的分析(cell-cell communication)。这里介绍一个简单的分析包celltalk...

    生信编程日常
  • SCmut||分析单细胞数据突变

    单细胞RNA测序(scRNA-seq)和DNA测序(scDNA-seq)都可以应用于细胞水平基因组分析。对于突变分析,scDNA-seq似乎更常见。然而,这项任...

    生信编程日常
  • 推荐一些配色的网站

    可以选择不同色系的网站: https://www.toodoo.com/db/color.html

    生信编程日常
  • 【Rust日报】2020-06-08 - Rust/WinRT快速入门

    Apache Benchmark style http bench tool written in async rust.

    MikeLoveRust
  • ChIP‐Atlas(逆向收费读文献2019-21)

    本次分享的文献发表了一个网页数据库,把其它数据库(GEO, ArrayExpress, DDBJ, ENCODE等等)的表观数据(主要是ChIP-seq and...

    生信技能树
  • MySQL之父解读数据库最新亮点

    本文是MySQL创始人Monty在5月30日"腾讯云CDB/CynosDB技术揭秘"系列直播中的分享实录。 ---- 大家好,我是MariaDB的 Micha...

    腾讯云数据库 TencentDB
  • 系统编程语言Rust特点介绍(1)

    最近在学习Rust语言,感觉这门语言有点意思,因此写一篇文章分享。我不会去介绍Rust的基本语法,什么变量声明,if..else..,循环等等。这些东西,文档介...

    glinuxer
  • php自带的缓存扩展-APCu

    在很多业务中,咱们都会碰到需要缓存数据的时候,而缓存方案有很多,也可以借助redis,memcache等第三方来进行缓存

    仙士可
  • Win10下的 Anaconda的安装

    首先我们先从Anaconda官网(https://anaconda.org/)上下载对应自己系统版本的Anaconda。因为我的电脑是win10,64位的, ...

    学到老
  • haskell 求助

    findBonding :: Eq a => (a -> a -> Bool) -> [a] -> Maybe [(a,a)]

    用户6797589

扫码关注云+社区

领取腾讯云代金券