clusterProfiler GeneRatio minGSSize
0
0
Entering edit mode
@zhangjianhai-12955
Last seen 4.8 years ago

Hello,

I am using "enrichKEGG" on rice genes. The entriz ids are entr <- c("4337100", "4336350", "4324731", "4336393", "4347034", "4345789", "4336561", "107276072", "4324047", "4352785", "4344600", "4327916", "4328844", "4336915", "4338072", "4339251", "4325336", "4329887", "107276640", "9266093", "9270198", "4346096", "4343818"). kk <- enrichKEGG(gene=entr, organism="osa", minGSSize=5, pvalueCutoff=0.05, pAdjustMethod = "BH").

The result is data.frame(kk):

               ID                        Description GeneRatio BgRatio
osa00430 osa00430 Taurine and hypotaurine metabolism       1/4 21/4948
osa00941 osa00941             Flavonoid biosynthesis       1/4 25/4948
osa03020 osa03020                     RNA polymerase       1/4 46/4948

             pvalue   p.adjust     qvalue  geneID Count
osa00430 0.01687387 0.04012714 0.01055977 4344600     1
osa00941 0.02006357 0.04012714 0.01055977 4337100     1
osa03020 0.03668235 0.04890979 0.01287100 4325336     1

My questions are: 1. The minGSSize=5, but why only 1 gene is enriched, as shown in the "GeneRatio" "1/4"? 2. The pAdjustMethod is "BH", then it produces "qvalue", so how does "p.adjust" come out?

Regards

enrichKEGG GeneRatio • 894 views
ADD COMMENT

Login before adding your answer.

Traffic: 667 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