Can I return a result for every GO term tested? I'd hoped increasing the p value to 1 would mean no results would be filtered out but it doesn't seem to have worked that way.
geneList below is a vector of fold changes with ensembl IDs as names
> library(clusterProfiler)
> library(org.Hs.eg.db)
> data(geneList, package = "DOSE")
> z <- gseGO(geneList, "BP", org.Hs.eg.db, pvalueCutoff = 1)
using 'fgsea' for GSEA analysis, please cite Korotkevich et al (2019).
preparing geneSet collections...
GSEA analysis...
leading edge analysis...
done...
Warning message:
In fgseaMultilevel(pathways = pathways, stats = stats, minSize = minSize, :
For some pathways, in reality P-values are less than 1e-10. You can set the `eps` argument to zero for better estimation.
> zz <- gseGO(geneList, "BP", org.Hs.eg.db)
using 'fgsea' for GSEA analysis, please cite Korotkevich et al (2019).
preparing geneSet collections...
GSEA analysis...
leading edge analysis...
done...
Warning message:
In fgseaMultilevel(pathways = pathways, stats = stats, minSize = minSize, :
For some pathways, in reality P-values are less than 1e-10. You can set the `eps` argument to zero for better estimation.
> summary(z$pvalue)
Min. 1st Qu. Median Mean
0.00000 0.04104 0.24027 0.33392
3rd Qu. Max.
0.58711 1.00000
> summary(zz$pvalue)
Min. 1st Qu. Median
0.000e+00 9.487e-06 5.611e-04
Mean 3rd Qu. Max.
1.450e-03 2.417e-03 6.513e-03
> length(z$ID)
[1] 5777
> length(zz$ID)
[1] 753