Applying idr package to ChIP-seq peaks allows researchers to compare multiple replicates and evaluate each peaks with reproducibility. It produces a output file which mimics the input narrowPeak file with additional two columns stating localIDR and globalIDR. Both localIDR and globalIDR represents -log10 IDR values. global IDR is more informative since it is more like multiple test corrected value of local IDR, according to this post (https://groups.google.com/forum/#!topic/idr-discuss/FY2K5VKx8AQ)
So if I want to find the most reproducible peaks, should I select the peaks with higher globalIDR or lower globalIDR? In other words, Does IDR (irreproducible discovery rate) higher mean less 'reproducible' or vice versa?
Another question I have is about pvalue/qvalue in 12-col report. Is the pvalue/qvalue for the peak prediction or for IDR test? In the original individual input narrowPeak file, you have pvalue representing the peak calling confidence for each peak in each file. When you combine in idr, how do you get pvalue/qvalue for it? average the original pvalue from all inputs?