gravatar for karlaarz

2 hours ago by

Hello,

I have the following human gene list:

CCT3 ABCA4 HEG1 PPARGC1B ADD1 CEP85 SLC1A4 DUSP10 PLAGL2 UBE2G2 NTRK2 PPIP5K1 DDB1 PRPH2 OAZ2 PEA15 ICMT KDM4A NCOA6 ZNF609 AKAP1 SYNE3 CAMSAP1 POLE4 ZDHHC5 ANGEL1 KCNJ14 NDUFA8 SIPA1L2 BTD CCT7 ANO2

I did the enrichment analysis using clusterProfiler

gobp <-enrichGO(glist,OrgDb=org.Hs.eg.db, ont = "BP",  pAdjustMethod = "fdr", keyType = 'SYMBOL', pvalueCutoff = 0.05)

head(gobp)
                   ID                                                              Description GeneRatio   BgRatio       pvalue   p.adjust    qvalue                              geneID Count
GO:1903405 GO:1903405                                     protein localization to nuclear body      2/36  10/18670 0.0001611010 0.02351193 0.0197561                           CCT3/CCT7     2
GO:1904851 GO:1904851 positive regulation of establishment of protein localization to telomere      2/36  10/18670 0.0001611010 0.02351193 0.0197561                           CCT3/CCT7     2
GO:1904867 GO:1904867                                       protein localization to Cajal body      2/36  10/18670 0.0001611010 0.02351193 0.0197561                           CCT3/CCT7     2
GO:0060249 GO:0060249                                         anatomical structure homeostasis      6/36 439/18670 0.0001748712 0.02351193 0.0197561 ABCA4/CCT3/CCT7/POLE4/ADD1/PPARGC1B     6
GO:0070203 GO:0070203          regulation of establishment of protein localization to telomere      2/36  11/18670 0.0001966624 0.02351193 0.0197561                           CCT3/CCT7     2
GO:0070202 GO:0070202        regulation of establishment of protein localization to chromosome      2/36  12/18670 0.0002357086 0.02351193 0.0197561                           CCT3/CCT7     2

However, if I do the same analysis using Panther or David, I don't find statistically relevant results as all FDR values are equal to 1. I know that cluster profiler calculates FDR values following the Storey, 2002 paper, but I wouldn't expect to see a big difference.

1) Why does clusterProfiler show different FDR values from Panther and David?

Thanks!



Source link