What parameter values to change in peak calling with MACS2 when calculating IDR afterwards?

To use IDR we need to relax peak calling settings when running MACS2. How do we do that? Is it enough to increase the qvalue threshold? Are there any specific parameter values we need to use?

chip-seq macs2 idr

• 2,315 views • link

written 8.0 years ago by salamandra

For MACS-related questions, you will likely get a response on the group's GitHub page, where the program athour is active: https://github.com/taoliu/MACS

• 0 views • link

written 8.0 years ago by Kevin Blighe

0 answers

No answers yet.

More posts like this

CUT&Tag data processing: peak caller choice and typical peak size profile
written 3.6 years ago by nihilior •

Greetings everyone, I am currently analysing CUT&Tag data (a recent technique, variant of the ChIP-seq), capturing different histone variants. To do so, I tried all …
Low read coverage in Input: Peak calling/Filtering
written 4.8 years ago by Ankit

Hi everyone, I have a datasets where the coverage of input is relatively. Actually many reads were filtered out during deduplication step. When I performed …
Chip-seq peak calling best way usage of parameters
written 5.5 years ago by luckysardar171 •

Hello, currently working on chip-seq analysis calling peaks using MACS2, need suggestions how to set parameters according to my data i have H3K4me1 H3K4me3 and …
how to call weak peak using macs2
written 6.3 years ago by huangweiping.official •

We are interesting in ChIP-seq small peak , we have avidence that the the small or weak peaks are functional sites , how to call …
ChIP-seq: To use IDR, samples have to pass both tests or is it enough to pass the pseudo-replicate …
written 8.0 years ago by salamandra

In order to use IDR threshold to select consistent peaks between two replicates of ChIP-seq, samples have to pass the pseudo-replicate consistency test (Np/Nt < …
Using MACS2 for peak calling
written 8.3 years ago by anais1396 •

Hello everyone !! I am new in dealing with peak calling and I chose macs2 for doing this on ATAC-seq samples. I performed a mapping …
Files to compare ChIP vs Input on genome browser: after alignment or peak calling?
written 8.4 years ago by salamandra

When visualizing ChIP-seq data on genome browser should we use files obtained upon read alignment (bam files) or instead use the files treat_qvalue.bdg and control_lambda.bdg …
If we use MACS2 do we need to remove duplicate sequences with samtools rmdup ?
written 8.5 years ago by salamandra

Hi, have some questions about sequence duplicates in ChIP-seq analysis: a) When running fastQC on my samples, they give the error on 'Sequence duplication levels'. …
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 …
Differential chip-seq peak calling using macs2
written 9.2 years ago by lucas.michel •

I have been using macs2 for peak calling and posterior differential peak calling. I am working with paired-end data. After calling peaks with macs2 callpeak …