Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

wr_hier.py only return GOterms in namespace 'BP' #163

Closed
NewTearingHeaven opened this issue Apr 26, 2020 · 5 comments
Closed

wr_hier.py only return GOterms in namespace 'BP' #163

NewTearingHeaven opened this issue Apr 26, 2020 · 5 comments

Comments

@NewTearingHeaven
Copy link

Dear author, i was running
$ scripts/wr_hier.py BP MF CC --gene2go=gene2go --taxid=9606 --dash_len=17 --concise -o human_BP_MF_CC.txt
but the results in txt only includes those GOterms in 'BP'. I've checked the gene2go file and it does include MF and CC.
Full results are here:
EXISTS: go-basic.obo
go-basic.obo: fmt(1.2) rel(2020-03-23) 47,232 GO Terms
HMS:0:00:04.080605 330,518 annotations, 20,562 genes, 18,403 GOs, 1 taxids READ: gene2go
18499 IDs in loaded association branch, BP
NoneGoSubDag: 47232 sources in 47232 GOs rcnt(True). 2701 alt GO IDs
NoneGoSubDag: namedtuple fields: NS level depth GO alt GO_name dcnt D1 tcnt tfreq tinfo id
NoneGoSubDag: relationships: set()
NoneGoSubDag: 12378 sources in 15696 GOs rcnt(True). 0 alt GO IDs
NoneGoSubDag: namedtuple fields: NS level depth GO alt GO_name dcnt D1 id
NoneGoSubDag: relationships: set()
WROTE: human_BP_MF_CC_test.txt

@dvklopfenstein
Copy link
Collaborator

dvklopfenstein commented Apr 29, 2020

Thanks for reporting this and thank you for the specific command to replicate the issue.

You are correct.

DETAILS: We changed the underlying infrastructure of the annotation readers so that they all share a common base class and in addressing that and another issue to separate the GO terms more easily into BP, MF, and CC, this command wrongly no longer prints MF and CC.

I'll fix this shortly. Not catching this earlier was due to missing functionality in the tests, which will also be corrected.

Thank so much for being alert and for taking the time to report this issue. And thank you for your interest in GOATOOLS.

@dvklopfenstein
Copy link
Collaborator

@tanghaibao , Can you please create a new version? This is an important fix and should be available to all. Thanks.

@NewTearingHeaven
Copy link
Author

Thanks for your reply

@tanghaibao
Copy link
Owner

@dvklopfenstein @NewTearingHeaven

Updated the PyPi version to v1.0.5.

@NewTearingHeaven
Copy link
Author

NewTearingHeaven commented May 16, 2020 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants