Question: Issues in generating workable rank file for GSEA
0
gravatar for csijst
12 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 12 months ago by Michael Love24k • written 12 months ago by csijst0
Answer: Issues in generating workable rank file for GSEA
0
gravatar for assaron
12 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 12 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 12 months ago by csijst0
Answer: Issues in generating workable rank file for GSEA
0
gravatar for Michael Love
12 months ago by
Michael Love24k
United States
Michael Love24k wrote:

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

 

 

ADD COMMENTlink written 12 months ago by Michael Love24k
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: 210 users visited in the last hour