The support.bioconductor.org editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: Issues in generating workable rank file for GSEA
0
gravatar for csijst
8 months ago by
csijst0
Singapore/National University of Singapore
csijst0 wrote:

Hi,

I'm trying to generate a ranking file for my differential expressed genes that was generated via DESeq2. I saw some websites like Genome Spot that taught me how to use the p-values and log2foldchange to rank the genes. But in DESeq2 output, some of them appear as "NA". And this will reappear in the ranking file. Apparently, pre-ranked GSEA either java Desktop application (from Broad Institute) or GenePattern are not able to process the NA entries.

As such, what should I do for genes that display NA? As in what is the usual procedure in handling such information? Remove them?

Thank you.

ADD COMMENTlink modified 8 months ago by Michael Love22k • written 8 months ago by csijst0
Answer: Issues in generating workable rank file for GSEA
0
gravatar for assaron
8 months ago by
assaron150
assaron150 wrote:

Hi,

In DESeq2 you can use t-statistic for ranking, there are no NA-s. If you still want to use p-values, set DESeq2 option cooksCutoff=FALSE for the results function.

ADD COMMENTlink written 8 months ago by assaron150

Hi,

I tried the command,


dDif_res <- results(dLRT, contrast=c("treatment", "shRNA", "ctrl"), alpha=0.05, altHypothesis="greaterAbs", cooksCutoff=FALSE)

head(dDif_res)


But the output still shows NA for some of the genes. These genes do not have any baseMean as well.

ADD REPLYlink written 8 months ago by csijst0
Answer: Issues in generating workable rank file for GSEA
0
gravatar for Michael Love
8 months ago by
Michael Love22k
United States
Michael Love22k wrote:

You can also set NA to 1. There’s example code for this in the vignette FAQ

 

 

ADD COMMENTlink written 8 months ago by Michael Love22k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 289 users visited in the last hour