Hi,
I'm generating a density plot using the ggplot2 library:
ggplot(df_chr20, aes(x=start1)) + geom_density() + theme_bw()
ggplot(df_chr21, aes(x=start1)) + geom_density() + theme_bw()
The data frame for df_chr20 has 1573 rows, and I wanted to plot the probability density of genomic regions as x = start1. The same for df_chr21, with 1271 rows. What I see in the plot is what I expect, but is there a reason why the y-axis probablilty densities are so low? The area under the curve is surely not equal to 1.
I expected the high localised high density in chr 21, and both the plots have similar scales in terms of y-axis values. But since they are so low, how can I interpret this plot?
I found some information here, but my values are much lower than the examples.
Any help would b every kind!
If the plot isn't visible, please click here or the link: https://freeimage.host/i/JhowWN
0 answers
No answers yet.
Log in to answer this question.
Why do you think they are low? Your chromosome has millions of bases, but you only have values for about a thousand of them. In other words, most values are 0, so the average should be very low.
Please see How to add images to a Biostars post to add your images properly. You need the direct link to the image, not the link to the webpage that has the image embedded (which is what you have used here)
Hello rishi!
It appears that your post has been cross-posted to another site: https://bioinformatics.stackexchange.com/questions/13381/very-low-probability-densities
This is typically not recommended as it runs the risk of annoying people in both communities.
Hello rishi!
We believe that this post does not fit the main topic of this site.
Asked and answered here: https://bioinformatics.stackexchange.com/questions/13381/very-low-probability-densities. OP, please don't do this, it's immensely disrespectful and annoying.
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!