Hi, all!
I'm trying to plot the DNA methylation profiles using deepTools and having annoying peaks on the gene body region as you can see the picture below.
The code I used;
$ computeMatrix scale-regions -S Sample1.bw Sample2.bw Sample3.bw -R genes.bed -a 5000 -b 5000 -o Samples.matrix --regionBodyLength 15000 --skipZeros
$ plotProfile -m Samples.matrix -out Samples.profile.png --perGroup
How can I remove them?
Thank you!
•
written
by
mbk0asis
1 answer
Give plotHeatmap a try, I expect you have some outlier bin in one gene (just remove that gene, though in the future you'll want to blacklist that region when making bigWig files) that will be apparent in the heatmap.
I am stuck at the first step which is ID the outliers. I was going to export the matrix using "–outFileNameMatrix" option and edit the matrix. Then, how can I make the edited file back to the matrix file 'plotHeamap' recognizes? Should I just use other tools like R to draw the heatmap?
If you want to actually edit the heatmap then don't bother with --outFileNameMatrix, just gunzip the output from computeMatrix, remove the problematic line, change the group_boundaries so they limits are correct, and then gzip the results back up.
In general, it's probably simpler to just figure out where the problem region is and remake the bigWig files with that region blacklisted.
Log in to answer this question.
-
reference-point setting from computeMatrix
written by QX •
-
DiffBind issue in identification of differentially enriched peaks between two conditions
written by accibio •
-
Are gene order maintained across heatmaps in the perGroup option in plotHeatmap deeptools?
written by glyptonotus •
-
Several labels in refPointLabel plotHeatmap Deeptools
written by glyptonotus •
-
Problem with the output of Deeptools PlotProfile: a strange pattern of repetitive summits
written by Alireza •
-
Calculate the ratio of H3K4me3 to H3K4me1 using output of deeptools computeMatrix
written by jia •
-
deeptools plotProfile plots inconsistently
written by Ekalavya •
-
plot a smoth TSS profile in R
written by Lila M
Hello.
Not sure if my problem aligns perfectly to this thread, but I want to share some weird plots from the computeMatrix-plotProfile workflow. Not sure about "flatness" of the plot.
The corresponding plotProfile script is as under.
The plot is very boxy and abrupt. I have checked the bedGraph values and they aren't uniform as such.
Any suggestions on what could be erroneous here.
Use the directions: How to add images to a Biostars post
Create a new thread if this is not directly related to the original question.