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

rainfallTransform() for circos.genomicRainfall() does not work for single locus per chromosome #390

Open
evan-john opened this issue Jun 3, 2024 · 0 comments

Comments

@evan-john
Copy link

Hi, I am trying to plot a genomicRainfall track.
Some data have only 1 locus per chromosome (I have some important data points to present that are low frequency occurrence).
After some investigating, it seems the rainfallTransform() function to obtain the plotting dist returns 'NA' for those single-locus-per-chromosome. See the example below for Chr1.

head(rainfallTransform(x))
chr start end dist
1 Chr1 2474424 2476277 NA
2 Chr2 52994 55607 86197.0
3 Chr2 141804 143777 86197.0
4 Chr3 76845 79978 4656.0
5 Chr3 84634 86360 2871.5
6 Chr3 87447 88624 810.0

The issue is, when I plot, I can not observe in Chr1, but Ch2/Chr3 data points are fine. I couldn't figure out a workaround. Is there a fix you can suggest?

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

1 participant