clusterProfiler - GeneRatio & nb of genes always positively correlated in dotplot
3
1
Entering edit mode
@jane-merlevede-5019
Last seen 6.1 years ago

Hello,

Sorry for asking a lot of questions about clusterProfiler. I try to fully understand the results and graphs I want to use.

When using dotplot() on the result of enrichGO() and enrichDO(), I always observe that the most significant terms are the one with the biggest number of genes and the highest gene ratio. Whereas, we can imagine to get a very significant term with very high gene ratio but small (at least not the biggest) number of genes...

I observed that also on all the examples I saw on the vignette and on the web.

What is the rational behind this? This way, the small gene sets can never appear significant...

Jane

clusterProfiler • 4.9k views
ADD COMMENT
0
Entering edit mode
Guangchuang Yu ★ 1.2k
@guangchuang-yu-5419
Last seen 23 days ago
China/Guangzhou/Southern Medical Univer…

pvalue is not correlated with gene ratio, see the color of the plot here.

ADD COMMENT
0
Entering edit mode
@jane-merlevede-5019
Last seen 6.1 years ago

Indeed, pvalue is not correlated with gene ratio, in my examples as well. But I did not mention this relashionship. I speak about the positive correlation between the number of genes and the gene ratio. Can you answer this please?

Your example illustrates perfectly my point: one more case with no high gene ratio and small number of genes, as we can expect.

ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode
@jane-merlevede-5019
Last seen 6.1 years ago

Thank you for your answer. I thought n was the number of genes in a category and not the size of the input list of genes of interest. That is why I did not expect correlation between the number of genes and the gene ratio.

It would be great to have the possibility to plot this information as well (ratio of genes of interest over number of genes in a category), that is adding a third possible value to "x".

ADD COMMENT
0
Entering edit mode

you mean k/M?

ADD REPLY
0
Entering edit mode

yes, exactly

ADD REPLY
2
Entering edit mode

I will not change gene ratio (k/n, data from input gene list) and background ratio (M/N, data from background annotation), as they exist quite a long time in the package.

Maybe I can introduce another column, namely odd ratio by (gene ratio)/(background ratio), in future release.

 

 

ADD REPLY

Login before adding your answer.

Traffic: 795 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6