Dear all,
I am currently trying to run a GO enrichment analysis using GOseq but I am slightly puzzled if I am doing things correctly given my (zero) output. My endgoal is to get results I can use in GOgadget.
I have the following files: - Factor labels (either differently expressed or not, so two groups) - Gene lenghts file - GO annotations linked to the gene IDs
However calculating pwf results in the following warning:
Warning message: In newton(lsp = lsp, X = G$X, y = G$y, Eb = G$Eb, UrS = G$UrS, L = G$L, : Fitting terminated with step failure - check results carefully
I created a plot showing the "proportion of DE", like on page 8 of the GOseq manual. I get a flat line on 0. Meaning that there is no significant difference?
Does anyone know what potentially could be the problem here?
0 answers
No answers yet.
Log in to answer this question.
-
Unable to fetch fasta files of contigs using Edirect in Linux
written by Mohamed Samir •
-
cellranger count error
written by connorjfausto •
-
Goseq - Pwf: Error In If (Hi <= Low) { : Missing Value Where True/False Needed
written by ryme •
-
gogadget use with non-model organism to visualize goseq results
written by arroyo.sr •
-
GOseq gene2cat format errors
written by et395 •
-
decreased genotyping rate after the conversion of lgen to ped and map
written by mingchen720 •
-
GOseq analysis with non-native organism
written by Bin •
-
topGO is performing enrichment test on GO terms that are not present in my dataset
written by Fabio Marroni
-
How to do gene length bias corrected KEGG pathway analysis for non-native species in goseq?
written by CandiceChuDVM
-
Goseq: Pwf: Error In If (Hi <= Low) { : Missing Value Where True/False Needed
written by liupfskygre
Can you provide the code, and plot? I think the failure is the reason why you have a flat line, do you have RNAseq data?
Yes, it's RNAseq data. I'm currently home, but will provide the code when I'll be in the office tomorrow. Thanks so far.