Question: GOseq output but strange warning message!
gravatar for TB18
10 months ago by
TB180 wrote:


Dear all,

I am running GOseq on a custom annotated dataset. I am getting a nicely looking output: list of enriched and depleted GOterms for my different groups, just as I want. But after calculating pwf I get the following warning message:

> pwf=nullp(cat_genes_vec,
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​

But pwf is calculated:

> head(pwf)
  DEgenes         pwf
1       0       243 0.001000788
2       0       560 0.001000780
3       0       262 0.001000788
4       0      1176 0.001000628
5       0       248 0.001000788
6       0       291 0.001000788

Further, the plot from nullp is just a flat line on 0 (or more likely 0.001) hard to tell. I cannot clarify this. I think the warning message and plot result are related; I do get an output when I continue. 

Can someone help me understanding this?

Thanks in advance!

go goseq pwf • 236 views
ADD COMMENTlink written 10 months ago by TB180

Hi there. I am getting the same warning message while executing broadenrich() for GO and KEGG pathway enrichment.

ADD REPLYlink written 4 months ago by Shaurya Jauhari0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 281 users visited in the last hour