Entering edit mode
5.3 years ago
Jingyue
▴
70
Hi,
I'm wondering how to plot the heatmap around the peak center in ChIPseeker? Can I plot 3 samples in 1 plot? Do I need a merged peak file?
Here is the code I use to plot around TSS, but don't know how to do it for peak center. Thank you!
promoter <- getPromoters(TxDb=txdb, upstream=1000, downstream=1000)
tagMatrixList <- lapply(files, getTagMatrix, windows=promoter)
hi Jingyue,
In your example you plot the profile of ChIP peaks binding to TSS region (distance of peaks from TSS).
What is exactly the other plot, you are talking about?
Using the center of the peak as reference and plot the signal of bam or bigwig file?
p.s. sorry editors, I accidentally used answer button
I would like to plot the ChIP peaks binding enrichment around its binding center.
I plotted this using deeptools, but would like to see how it look like in ChIP-seeker.
ChIP-seq inputs are peak.bed files, so I don't think is over bam/bigwig files
Thank you!