The editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: GOseq output but strange warning message!
gravatar for TB18
6 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 • 174 views
ADD COMMENTlink written 6 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 25 days 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: 219 users visited in the last hour