Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

plot showing broken peak, not continuous #51

Open
wangxu1990 opened this issue Jan 31, 2025 · 1 comment
Open

plot showing broken peak, not continuous #51

wangxu1990 opened this issue Jan 31, 2025 · 1 comment

Comments

@wangxu1990
Copy link

Hi, thank you for this wonderful package that helped me a lot to visualize my chip-seq data.
Now I have no idea why the plot shows the broken peak as shown in the figure. My bigwig file was produced by deeptools bamcoverage with binsize 10 and smooth 30. Does it result from the big binsize? I guess not, because previously, I tried with another peak visualizing package, BioSeqUtils, https://github.com/junjunlab/BioSeqUtils, which gave me the non-broken plot.
Thank you!
Image

@m-jahn
Copy link
Collaborator

m-jahn commented Feb 3, 2025

Dear @wangxu1990, this is a known bug/phenomenon of plotting many small "bars" (essentially lines) next to each other. The geom_coverage functions do not plot a continuous horizon like graph, but each bar separately. If you have many tiny bars, the png plotting device tries to interpolate the observed thousands of bars to fewer dominant ones, and the interlacing (merging points together) produces these stripes as an artifact.

Needs to be handled somehow in nex releases. You can try to bin your data und reduce the number of vertical bars until a fix is out.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants