Question: Help with HOMER motif profile plots
0
gravatar for venu
8 months ago by
venu6.3k
Germany
venu6.3k wrote:

Hi,

I'm trying to understand profile plots for motifs. Homer has a function annotatePeaks.pl) to generate motif density profiles that looks like following (lifted from here)

enter image description here

While I understand what these profile plots convey for ChIP-seq data (e.g. those generated by deepTools plotProfile), but I'm not sure what exactly is the y-axis in above plot.

So the question is, is it just the probability of finding given motif at each position (as MEME does here) in a set of given sequences or is it the actual count normalized to something?

I didn't find any explanation on homer website. It would be great if anyone can explain what exactly is going on in above plot.

Thank you.

motifs homer • 491 views
ADD COMMENTlink modified 8 months ago • written 8 months ago by venu6.3k
1

Convert a BED file of motif positions to a bigWig with a 1 where there's a motif and a 0 otherwise and it's the same thing you'd then get with plotProfile I think.

ADD REPLYlink written 8 months ago by Devon Ryan92k

I created a bigiwg ( using bedtools genomecov + kent utils) and plotted with plotProfile. The motif I'm interested in is has very sparse hits, applying loess fit on plotProfile output data is more/less looking closer to the above plot. I think I understood what's happening inside.

ADD REPLYlink written 8 months ago by venu6.3k
3
gravatar for venu
8 months ago by
venu6.3k
Germany
venu6.3k wrote:

For future googlers, here is the response from lead author of HOMER, Chris Benner

The motif plots are normalized by taking the actual motif totals, then normalizing by the bin size used in the histogram (so that the motif densities are comparable across different bin sizes) and the total number of peaks/regions used to create the histogram (so that the motif densities are independent of the number of regions used to search).

ADD COMMENTlink written 8 months ago by venu6.3k
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: 970 users visited in the last hour