Goseq problem with "missing value"
1
0
Entering edit mode
@mgierlinski-7369
Last seen 8 weeks ago
United Kingdom

I'm running goseq on a sequence of bootstraps where a DE test is performed on a random selection of replicates. THis way, I have a constant list of background genes and a selection of significant genes which is slightly changing from bootstrap to bootstrap. goseq works nicely on most of these bootstraps, but every now and then it crashes. Here is an example:

> library(goseq)
> library(org.Sc.sgd.db)
> all.genes = toupper(scan("test_all.txt", what="character"))
Read 7126 items
> sel.genes = toupper(scan("test_sel.txt", what="character"))
Read 630 items
> gene.vector = as.integer(all.genes%in%sel.genes)
> names(gene.vector) = all.genes
> pwf=nullp(gene.vector,  "sacCer2", "ensGene")
Loading sacCer2 length data...
Error in if (min(fv) < lower_bound) fv = fv - min(fv) + lower_bound : 
  missing value where TRUE/FALSE needed

This happens only to a few bootstraps, must be some peculiar selection of genes causing it.

Here are the text files (gene lists) I used in this example (they crash goseq):

test_all.txt   test_sel.txt

 

 

 

goseq • 1.3k views
ADD COMMENT
1
Entering edit mode
@nadia-davidson-5739
Last seen 4.9 years ago
Australia

Hi,

Which versions of goseq and mgcv are you using? There was an update made to goseq from 1.16.2, which may fix this issue, so consider updating it if it's older. You'll may also need to update mgcv to the latest available version. If this still doesn't work, can you post the result of "sessionInfo()". I haven't be able to reproduce the error with the versions of packages I have.

 

Cheers,

Nadia.

ADD COMMENT
0
Entering edit mode

OK, this helped! I updated goseq to version 1.18.0 and the error disappeared.

Thanks a lot,

Marek

ADD REPLY

Login before adding your answer.

Traffic: 760 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6