Question: ChIP-seq heatmap: Which MACS2 output file should we give as input into deeptools computeMatrix?
0
gravatar for salamandra
9 months ago by
salamandra170
salamandra170 wrote:

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 intensity.bw -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'
ADD COMMENTlink written 9 months ago by salamandra170
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1159 users visited in the last hour