Hi all, I'm trying to get Gene Ontology and Gene Set Enrichment using topGO. I'm pretty novice to bioinformatics using R, so I need your help!
So from RNA-Seq raw count data (31973 genes), I used DESeq2 to get differentially expressed genes, using results(), and from there I sorted by p-values. I didn't use adjusted p-values, since there are too many genes, when I adjusted p-values, all of them were above 0.9, which did not give differentially expressed genes. So, I sorted by p-values, and then about to look for genes that are less than 0.05, however, since I didn't use adjusted p-values, I reduced the p-value cutoff to 0.01. QUESTION: Is it acceptable to use reduced p-values (0.05 -> 0.01) instead of adjusted p-values?
QUESTION2: What is the best R tools package to use this .csv file (consists of Ensembl gene ID, baseMean, log2FoldChange, lfcSE, stat, pvalue, padj) as an input and give me Gene Ontology analysis and Gene Set Enrichment as outputs?
Thanks!
maybe you can follow this workflow: https://github.com/twbattaglia/RNAseq-workflow