Mercurial > repos > proteore > proteore_reactome
changeset 11:19195d1a4063 draft default tip
"planemo upload commit a2de8dd8f2468046e787f83812d201bb191edab8-dirty"
| author | proteore |
|---|---|
| date | Mon, 10 May 2021 15:30:34 +0000 |
| parents | ef31b5ac28d7 |
| children | |
| files | reactome_analysis.py reactome_analysis.xml test-data/Reactome_Analysis_Tools.html test-data/Reactome_Analysis_Tools_invalid_ID.txt |
| diffstat | 4 files changed, 35 insertions(+), 26 deletions(-) [+] |
line wrap: on
line diff
--- a/reactome_analysis.py Fri Jan 24 10:45:56 2020 -0500 +++ b/reactome_analysis.py Mon May 10 15:30:34 2021 +0000 @@ -1,7 +1,12 @@ -import os, re, json, argparse, csv +import argparse +import csv +import json +import os +import re CURRENT_DIR = os.path.dirname(os.path.abspath(__file__)) + def id_valid(identifiers): """ Validate IDs if they contain special characters @@ -15,7 +20,8 @@ else: remove.append(id) return res, remove - + + def isnumber(format, n): """ Check if an variable is numeric @@ -32,6 +38,7 @@ else: return False + def data_json(identifiers): """ Submit IDs list to Reactome and return results in json format @@ -40,67 +47,70 @@ trash = [] if identifiers[1] == "list": ids = identifiers[0].split() - ids = [x.split(";") for x in ids] - ids = [item.strip() for sublist in ids for item in sublist if item != ''] - json_string = os.popen("curl -H \"Content-Type: text/plain\" -d \"$(printf '%s')\" -X POST --url www.reactome.org/AnalysisService/identifiers/\?pageSize\=1\&page\=1" % ids).read() + ids = [x.split(";") for x in ids] + ids = [item.strip() for sublist in ids for item in sublist if item != ''] # noqa 501 + json_string = os.popen("curl -H \"Content-Type: text/plain\" -d \"$(printf '%s')\" -X POST --url www.reactome.org/AnalysisService/identifiers/\?pageSize\=1\&page\=1" % ids).read() # noqa 501 if len(id_valid(identifiers[0].split())[1]) > 0: trash = id_valid(identifiers[0].split())[1] elif identifiers[1] == "file": header = identifiers[2] - with open(identifiers[0],"r") as mq : - file_content = csv.reader(mq,delimiter="\t") - file_content = list(file_content) #csv object to list + with open(identifiers[0], "r") as mq: + file_content = csv.reader(mq, delimiter="\t") + file_content = list(file_content) # csv object to list ncol = identifiers[3] if isnumber("int", ncol.replace("c", "")): if header == "true": - idens = [x for x in [line[int(ncol.replace("c", ""))-1].split(";") for line in file_content[1:]]] + idens = [x for x in [line[int(ncol.replace("c", ""))-1].split(";") for line in file_content[1:]]] # noqa 501 else: - idens = [x for x in [line[int(ncol.replace("c", ""))-1].split(";") for line in file_content]] - - idens = [item.strip() for sublist in idens for item in sublist if item != ''] #flat list of list of lists, remove empty items + idens = [x for x in [line[int(ncol.replace("c", ""))-1].split(";") for line in file_content]] # noqa 501 + # flat list of list of lists, remove empty items + idens = [item.strip() for sublist in idens for item in sublist if item != ''] # noqa 501 ids = "\n".join(id_valid(idens)[0]) - json_string = os.popen("curl -H \"Content-Type: text/plain\" -d \"$(printf '%s')\" -X POST --url www.reactome.org/AnalysisService/identifiers/\?pageSize\=1\&page\=1 2> stderr" % ids).read() + json_string = os.popen("curl -H \"Content-Type: text/plain\" -d \"$(printf '%s')\" -X POST --url www.reactome.org/AnalysisService/identifiers/\?pageSize\=1\&page\=1 2> stderr" % ids).read() # noqa 501 if len(id_valid(idens)[1]) > 0: trash = id_valid(idens)[1] - #print(json_string) + # print(json_string) j = json.loads(json_string) - print ("Identifiers not found: " + str(j["identifiersNotFound"])) - print ("Pathways found: " + str(j["pathwaysFound"])) + print("Identifiers not found: " + str(j["identifiersNotFound"])) + print("Pathways found: " + str(j["pathwaysFound"])) return json_string, trash + def write_output(filename, json_string, species, trash_file, trash): """ Replace json result in template and print to output """ template = open(os.path.join(CURRENT_DIR, "template.html")) output = open(filename, "w") - try: + try: for line in template: if "{token}" in line: line = line.replace("{species}", species) - line = line.replace("{token}", json.loads(json_string)["summary"]["token"]) + line = line.replace("{token}", json.loads(json_string)["summary"]["token"]) # noqa 501 output.write(line) except ValueError: - output.write("An error occurred due to unavailability of Reactome web service. Please return later.") + output.write("An error occurred due to unavailability of Reactome web service. Please return later.") # noqa 501 template.close() output.close() - + if trash: - #print(trash) + # print(trash) trash_out = open(trash_file, "w") trash_out.write("\n".join(trash)) trash_out.close() + def options(): parser = argparse.ArgumentParser() argument = parser.add_argument("--json", nargs="+", required=True) argument = parser.add_argument("--output", default="output.html") argument = parser.add_argument("--trash", default="trash.txt") - argument = parser.add_argument("--species", default="48887") + argument = parser.add_argument("--species", default="48887") # noqa 841 args = parser.parse_args() filename = args.output json_string, trash = data_json(args.json) write_output(filename, json_string, args.species, args.trash, trash) + if __name__ == "__main__": options()
--- a/reactome_analysis.xml Fri Jan 24 10:45:56 2020 -0500 +++ b/reactome_analysis.xml Mon May 10 15:30:34 2021 +0000 @@ -1,4 +1,4 @@ -<tool id="reactome_analysis" name="Pathway enrichment analysis" version="2020.01.23"> +<tool id="reactome_analysis" name="Pathway enrichment analysis" version="2021.04.20"> <description>[Reactome]</description> <requirements> </requirements> @@ -65,7 +65,7 @@ <param name="ncol" value="c1" /> </conditional> <output name="output" file="Reactome_Analysis_Tools.html" ftype="html" /> - <!-- <output name="trash" file="Reactome_Analysis_Tools_invalid_ID.txt" ftype="tab" hidden="true" /> --> + <output name="trash" file="Reactome_Analysis_Tools_invalid_ID.txt" ftype="tabular" /> </test> </tests> <help><![CDATA[
--- a/test-data/Reactome_Analysis_Tools.html Fri Jan 24 10:45:56 2020 -0500 +++ b/test-data/Reactome_Analysis_Tools.html Mon May 10 15:30:34 2021 +0000 @@ -35,7 +35,7 @@ <p>Please click the button to execute the analysis:</p> - <form action="http://www.reactome.org/PathwayBrowser/#/DTAB=AN&ANALYSIS=MjAxODAxMTYxMTQ5NDZfNjczMg%3D%3D" target="_blank"> + <form action="http://www.reactome.org/PathwayBrowser/#/SPECIES=48887&DTAB=AN&ANALYSIS=MjAyMTA0MjAwODQyNDJfMjc1NTM%3D" target="_blank"> <button>Analyse</button>