goseq: recent update seems to restrict the range of applications
1
0
Entering edit mode
Oleg Moskvin ▴ 60
@oleg-moskvin-4293
Last seen 9.3 years ago
United States

Hi Matt and Nadia,

For almost a year, I was using goseq as a generalized enrichment tool (leveraging the gene length bias correction implemented there, thank you) using my custom gene sets this way:

goseq(myData, genome=NULL, id=NULL, gene2cat = MyVeryOwnGeneCategorization, method="Sampling", repcnt = 100000)

After recent update to version 1.18.0, the pipeline crashed and returned this:

################################
Error in .testForValidKeys(x, keys, keytype) :
 None of the keys entered are valid keys for 'GOID'. Please use the keys method to see a listing of valid arguments.
################################

It looks like the current version became explicitly restricted to use of the pre-defined geneset categories (I see "test.cats" argument that defaults to c("GO:CC", "GO:BP", "GO:MF") in the manual and when I looked at the code itself, it seems to accept "KEGG" as well and return an error in any other case).

This restriction dramatically narrows the applicability of goseq. Could this be changed in the upcoming version?

Thank you,

Oleg

goseq • 1.3k views
ADD COMMENT
0
Entering edit mode

P.S. I happened to still have an older Bioconductor version on one workstation, and exactly the same pipeline runs without any issues there (goseq_1.16.2).

 

ADD REPLY
0
Entering edit mode
@nadia-davidson-5739
Last seen 4.9 years ago
Australia

This response is a bit late, but for those interested, I was not able to recreate this issue.

If you can send me the results of sessionInfo() and a small example with data and commands run, I'll look into it more detail.

 

Cheers,

Nadia.

ADD COMMENT
0
Entering edit mode

I've also been successfully using goseq 1.18.0 for non-GO enrichment analysis.

I have noticed issues when attempting to using `method='Sampling'` in the past, so my first thought would be to try disabling that initially and see if it works.

What does your `MyVeryOwnGeneCategorization` look like?

 

ADD REPLY

Login before adding your answer.

Traffic: 899 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