[BioC] ChIPpeakAnno, callEnrichedRegions results have "holes"

Susan Wilson smwilson at hpc.unm.edu
Wed Sep 5 23:18:57 CEST 2012


I have some ChIP-seq data which returns a fairly uniform set of scores 
with the computeMATScore function. When I export (rtracklayer) that data 
as a .wig file, IGV visually reinforces that fact. However, when I run 
callEnrichedRegions, no matter how I try and "relax" the thresholds, 
using either method = score, pValue, or FDR, the resulting (bedgraph) 
data has single contiguous "holes" in regions of yeast chromosomes IV, 
X, XII, XIII, XV, and XVI (note, the chromosomes are not completely 
blank). The wiggle plot shows no obvious qualitative differences in 
those regions. Has anyone seen this behavior?


More information about the Bioconductor mailing list