How to get the pooled-replicate narrowPeak file for IDR peak-list · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

How to get the pooled-replicate narrowPeak file for IDR peak-list

Hi,

I am newbie for peak calling. When I am trying to follow ENCODE's pipeline to merge narrowPeak from replicates, I stuck with this pooled-replicate narrowPeak file to run IDR. I check the IDR github and it says "Run idr using an oracle peak list (e.g. peaks called from merged replicates)", which seems the pooled-replicates narrowPeak is merging the bam file first and then use macs2 to call peaks. But I am not sure what is the control file for the merged bam file? Should I merged two replicate control file to get the control file?

Or this pooled replicate is just simply merged the narrowPeak files from replicates.

Thanks a lot!

idr pooled-replicate

• 2,345 views • link

written 7.3 years ago by ben.kunfang •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Best way to combine biological replicate peaks?
written 1.3 years ago by SpartanII •

I have performed ATAC-seq and have three biological replicates and would like to combine their macs2 narrowPeak files to get the common peaks between them, …
Getting differential ChIP-seq peaks between conditions after calling consensus peaks
written 2.2 years ago by Alexandra •

Hello! I have some, pearhaps naive, questions about consensus peaks called using Tao Liu's script on GitHub called merge_then_call_consensus.sh. Say, for example, I have 4 …
ChIP seq peak annotation
written 3.6 years ago by tsomakiank •

Hello! I am trying to annotate H3K9me2 called peaks(macs2) to genes. I am using the homer's annotatePeaks.pl script for the annotation and from the resulting …
Question about ENCODE ATAC-seq bed file, bed narrowPeak vs bed idr_ranked_peak
written 5.4 years ago by yanweng •

I am trying to download all the human ATAC-seq peak from ENCODE. They have a page to allow you filter the samples, and generate a …
IDR on blacklisted region-filtered narrowPeak
written 8.3 years ago by endicytosis •

Hi guys, I am having a problem performing IDR on my blacklist region-filtered file. After I called my narrow peaks using MACS2, I used bedtools …
Does IDR (irreproducible discovery rate) higher mean less 'reproducible' or vice versa?
written 8.5 years ago by sckinta

Applying [idr package][1] to ChIP-seq peaks allows researchers to compare multiple replicates and evaluate each peaks with reproducibility. It produces a output file which mimics …
Biological Replicate Chip-exo
written 9.5 years ago by Sudhir Jadhao

Hello everyone, I have biological replicate, I have called peaks separately and to get final peaks I have used IDR and common peak with at …
ChipSeq: merge bam file before peak calling
written 9.9 years ago by Mike

Hi all, I have chipseq data for 6 samples (treatment/control), with replicate (so total 12 bam files). Now I am looking for peak calling using …
MACS2 Peak Calling : Different number of peaks for each replicate
written 10.5 years ago by AB

Hi All, I am a newbie and I am doing a ChIP-Seq analysis on a small dataset of 6 samples (Wild-type and knockout having 3 …
Apply Idr On List Of Read Counts
written 13.5 years ago by 14134125465346445

<p>Does anybody know how to apply the IDR code on lists of read counts per peak, instead of the usual narrowPeak format:</p> <pre><code>peak-1 20 peak-2 …