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?