Irreproducible Discovery Rate
1
2
Entering edit mode
10.2 years ago
Dataminer ★ 2.7k

Hi!

I am computing irreproducible discovery rate (see here https://sites.google.com/site/anshulkundaje/projects/idr ) for my two ChIP seq profiles for Pol II on same cell line.

My peaks were called using MACS2 and peaks for both ChIP seq profiles were of width 500-520 bases.

I am using following command to compute IDR:

Rscript batch-consistency-analysis2.r ./encodepeaks_1 ./encodepeaks_2 -1 ./Analysis 0 F p.value

If you have used this IDR package you must be knowing that you get a file which has genomic co-ordinates from two ChIP seq profiles which are overlapping with a IDR value in last column (lower the IDR, better the peak).

My problem is few of the genomic co-ordinate in the resulting file are much bigger, in the range of 20041463 bps. Which is very abnormal because none of my peaks are that broad in any of the profile original profile (given as input, i.e., broadpeak1 & broadpeak2).

Have you ever faced such a problem while using IDR? if yes, how you had solved it? I cannot discard these peaks because they fall on my target gene, but they are just too broad.

Please give your advice or suggestion or solution.

Thank you

chip-seq genomics bioinformatics r • 8.6k views
ADD COMMENT
0
Entering edit mode

The old version of IDR (1.x) sometimes mixed and matched starts/ends from input MACS2 peaks, resulting in these multi-Mb monsters. In addition, it regularly paired peaks that were on different chromosomes!!

Be sure to use the new IDR (2.x), available on GitHub here. None of the problems of the old, except the output fold-change value must be divided by 2 (they add replicate fold-changes, instead of averaging??).

ADD REPLY
0
Entering edit mode

This version seems far easier to use, but the old idr is not restricted to only two replications

ADD REPLY
1
Entering edit mode
10.2 years ago
Dataminer ★ 2.7k

Edited the Command

Rscript batch-consistency-analysis2.r ./broadpeaks_1 ./broadpeaks_2 -1 ./Analysis 0 F p.value

This is was one of the things which needed to be fixed, and be careful of the human genome table you use.

ADD COMMENT
0
Entering edit mode

Be sure to mark you question as answered--answering your own question is allowed and encouraged.

ADD REPLY
0
Entering edit mode

The site won't allow me to mark my answer as accepted :(

ADD REPLY
0
Entering edit mode

I guess that's also falls under the "voting on own post" prohibition since it also nets reputation points.

ADD REPLY
0
Entering edit mode

Oops--sorry about the misinformation. I think this could be an exception?

ADD REPLY

Login before adding your answer.

Traffic: 2521 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6