ChIP-seq heatmap: Which MACS2 output file should we give as input into deeptools computeMatrix?
Entering edit mode
3.1 years ago
salamandra ▴ 360

Hi, I want to generate an heatmap with peak intensity values for genomic regions.

I'm using the 'computeMatrix' followed by 'plotHeatmap' from 'deeptools'.

MACS2 outputs 3 bedGraph files: 'treat_pileup.bdg', 'control_lambda.bdg' and 'qvalue.bdg', all them i converted to bigWig. Which of those bigWig files should we use in computeMatrix? Which of those shows peak intensity as difference from ChIP and Input?

computeMatrix reference-point -S -R regions.bed -out output.gz --referencePoint center --outFileSortedRegions sortedreg.bed -b 4000 -a 4000 --missingDataAsZero

plotHeatmap -m output.gz -out figure.eps --dpi 300 --colorList '#FFFFFF,#FDF0E3,#FCE1C7,#FBD3AB,#FAC48F,#F8B573,#F7A757,#F6983B,#F58A20' --dpi 300 --whatToShow 'heatmap only' --refPointLabel 'Peak center'
deeptools computeMatrix RNA-Seq ChIP-Seq MACS2 • 1.7k views

Login before adding your answer.

Traffic: 1916 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6