Selecting DEGs for gene ontology using clusterprofiler and DESeq2
0
2
Entering edit mode
5 weeks ago

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 • 464 views
ADD COMMENT
2
Entering edit mode

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

ADD REPLY
0
Entering edit mode

Thank you!

ADD REPLY
1
Entering edit mode

Would greatly appreciate any help!!

ADD REPLY

Login before adding your answer.

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