This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Selecting DEGs for gene ontology using clusterprofiler and DESeq2

Hi all,

I am using DESeq2 to generate a DEG list and then feeding significant genes into enrichGO function of clusterprofiler (using below code).

My question is if there are any standard definitions of what constitutes a "significant" gene for GO analysis on the basis of L2FC, padj, basemean, etc.

Thanks much!

res<-results(dds)
sigs<-na.omit(res)
sigs<-sigs[sigs$padj < 0.05 & sigs$baseMean >50,]
sigspos<-sigs[sigs$log2FoldChange>0,]
sigsposname<-sigspos$symbol
sigsneg<-sigs[sigs$log2FoldChange< 0,]
sigsnegname<-sigsneg$symbol

SZup<- enrichGO(gene = sigsposname, OrgDb = "org.Hs.eg.db", keyType ="SYMBOL", ont = "BP")
r deseq2 rnaseq

depends on your results. typically padj <=0.05, abs(logFC)>=1

0 answers

No answers yet.

Log in to answer this question.