Warning messages: gseaKEGG
Entering edit mode
adR ▴ 10
Last seen 2 days ago
Germany, M√ľnchen

Hello Dear friends,

I had a question and was already posted in the Biostars platform but recommended by a friend to bring it here for further feedback. The questions are here bellow in the link. https://www.biostars.org/p/454938/#454949 Could you please give me some more tips on how to handle these? The input data has only geneIDs(enterZ) and log2fold change values of the genes. Thank you so much!

clusterProfiler r Tutorial • 660 views
Entering edit mode
Last seen 1 hour ago
Republic of Ireland

The warning messages say more about your upstream analysis, as opposed to this being an issue with clusterProfiler. You may want to figure out why ties may exist in your [assumed ranked] input data. Ties of this nature may exist in studies with small sample n, imbalanced comparison groups, or, generally studies of poor design and/or poor model assumptions.

If you want to show a reproducible example, that would help us very much.


Entering edit mode

Thanks, Kevin, I would like to share with you the top 50 in the following picture. In the picture, you see the named gene with the log2FC values. The code and link for the picture are below I used:

z <- gseKEGG(geneList = y, 
    pvalueCutoff = 0.5, # pvalue cutoff 
    organism = "hsa", # supported organisms listed below
    keyType = "kegg",
    exponent = 1,
    eps = 0,
    pAdjustMethod = "none",
    minGSSize = 0,
    by = "fgsea",
    use_internal_data = FALSE)

top 50 genes with their LFC values


Login before adding your answer.

Traffic: 408 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6