Warning messages: gseaKEGG
1
0
Entering edit mode
adR ▴ 30
@do-it-23093
Last seen 7 months 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 • 1.2k views
ADD COMMENT
0
Entering edit mode
Kevin Blighe ★ 3.9k
@kevin
Last seen 1 day 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.

Kevin

ADD COMMENT
0
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

ADD REPLY

Login before adding your answer.

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