Skip to content

Commit

Permalink
revert previous change
Browse files Browse the repository at this point in the history
  • Loading branch information
tanghaibao committed Sep 27, 2023
1 parent 8e86e44 commit 13d7cc4
Show file tree
Hide file tree
Showing 2 changed files with 27 additions and 25 deletions.
50 changes: 26 additions & 24 deletions goatools/anno/init/reader_idtogos.py
Original file line number Diff line number Diff line change
Expand Up @@ -64,31 +64,33 @@ def _init_w_godag(self, prt=stdout):
for itemid, gos in self.id2gos.items():
to_add = set()
for goid in gos:
if goid in s_godag:
goobj = s_godag[goid]
if goobj.is_obsolete:
if self.obsolete == "keep":
prt.write(f"**WARNING: {goid} obsolete in DAG, kept\n")
to_add.add(goid)
elif self.obsolete == "replace":
to_replace = set()
if "replaced_by" in goobj.__dict__ and goobj.replaced_by:
to_replace |= set(goobj.replaced_by.split(","))
if "consider" in goobj.__dict__ and goobj.consider:
to_replace |= goobj.consider
if to_replace:
prt.write(
f"**WARNING: {goid} obsolete in DAG, replaced by {to_replace}\n"
)
else:
prt.write(
f"**WARNING: {goid} obsolete in DAG, no replacement\n"
)
to_add |= to_replace
elif self.obsolete == "skip":
prt.write(f"**WARNING: {goid} obsolete in DAG, skipped\n")
else:
if goid not in s_godag:
prt.write(f"**WARNING: {goid} NOT FOUND IN DAG\n")
continue
goobj = s_godag[goid]
if goobj.is_obsolete:
if self.obsolete == "keep":
prt.write(f"**WARNING: {goid} obsolete in DAG, kept\n")
to_add.add(goid)
elif self.obsolete == "replace":
to_replace = set()
if "replaced_by" in goobj.__dict__ and goobj.replaced_by:
to_replace |= set(goobj.replaced_by.split(","))
if "consider" in goobj.__dict__ and goobj.consider:
to_replace |= goobj.consider
if to_replace:
prt.write(
f"**WARNING: {goid} obsolete in DAG, replaced by {to_replace}\n"
)
else:
prt.write(
f"**WARNING: {goid} obsolete in DAG, no replacement\n"
)
to_add |= to_replace
elif self.obsolete == "skip":
prt.write(f"**WARNING: {goid} obsolete in DAG, skipped\n")
else:
to_add.add(goid)
for goid in to_add:
goobj = s_godag[goid]
namespace = goobj.namespace
Expand Down
2 changes: 1 addition & 1 deletion goatools/cli/find_enrichment.py
Original file line number Diff line number Diff line change
Expand Up @@ -354,7 +354,7 @@ def __init__(self, args):
self.godag = GODag(
obo_file=self.args.obo,
optional_attrs=godag_optional_attrs,
load_obsolete=self.args.obsolete in ("keep", "replace"),
load_obsolete=True,
)
# GET: Gene2GoReader, GafReader, GpadReader, or IdToGosReader
self.objanno = self._get_objanno(self.args.filenames[2])
Expand Down

0 comments on commit 13d7cc4

Please sign in to comment.