Skip to content
Snippets Groups Projects
Commit 63b01794 authored by Arthur Le Bars's avatar Arthur Le Bars
Browse files

Added Homo sapiens deletion from chado database

parent 0b4c80b7
No related branches found
No related tags found
1 merge request!1Release 1.0
#!/usr/bin/env bash
\ No newline at end of file
...@@ -9,10 +9,7 @@ import json ...@@ -9,10 +9,7 @@ import json
import urllib3 as ul import urllib3 as ul
from chado import ChadoInstance from chado import ChadoInstance
from workflow import Workflow from workflow import Workflow
from filetransfer import FileTransfer
from toolrunner import ToolRunner from toolrunner import ToolRunner
from webscrap import WebScrap
""" """
TODO: script description TODO: script description
python3 ~/PycharmProjects/ggauto/gga_load_data/main.py ~/PycharmProjects/ggauto/gga_load_data/dataloader.json python3 ~/PycharmProjects/ggauto/gga_load_data/main.py ~/PycharmProjects/ggauto/gga_load_data/dataloader.json
...@@ -207,6 +204,19 @@ def main(): ...@@ -207,6 +204,19 @@ def main():
gi.histories.upload_dataset_from_library(history_id=current_hi_id, lib_dataset_id=datasets["transcripts_file"]) gi.histories.upload_dataset_from_library(history_id=current_hi_id, lib_dataset_id=datasets["transcripts_file"])
gi.histories.upload_dataset_from_library(history_id=current_hi_id, lib_dataset_id=datasets["proteins_file"]) gi.histories.upload_dataset_from_library(history_id=current_hi_id, lib_dataset_id=datasets["proteins_file"])
# Delete Homo sapiens from Chado database
toolrunner = ToolRunner(parameters_dict=sp_dict, instance=gi, history=current_hi_id)
sapiens_id = None
sapiens = toolrunner.get_sapiens_id()
sapiens_job_out = sapiens["outputs"][0]["id"]
sapiens_json_output = gi.datasets.download_dataset(dataset_id=sapiens_job_out)
try:
sapiens_output = json.loads(sapiens_json_output)[0]
sapiens_id = str(sapiens_output["organism_id"]) # needs to be str to be recognized by the chado tool
toolrunner.delete_sapiens(hs_id=sapiens_id)
except bb.ConnectionError:
print("Homo sapiens isn't in the database")
# Workflow generation # Workflow generation
workflow = Workflow(parameters_dict=sp_dict, instance=gi, history_id = current_hi_id) workflow = Workflow(parameters_dict=sp_dict, instance=gi, history_id = current_hi_id)
wf_dict_json = workflow.generate(working_directory=wd, main_directory=main_dir, workflow_name="preset_workflow") wf_dict_json = workflow.generate(working_directory=wd, main_directory=main_dir, workflow_name="preset_workflow")
......
...@@ -45,6 +45,7 @@ class ToolRunner: ...@@ -45,6 +45,7 @@ class ToolRunner:
self.instance.tools.run_tool(tool_id="toolshed.g2.bx.psu.edu/repos/gga/chado_organism_delete_organisms/organism_delete_organisms/2.3.2", self.instance.tools.run_tool(tool_id="toolshed.g2.bx.psu.edu/repos/gga/chado_organism_delete_organisms/organism_delete_organisms/2.3.2",
history_id=self.history, history_id=self.history,
tool_inputs={"name": self.organism}) tool_inputs={"name": self.organism})
def purge_analyses(self): def purge_analyses(self):
return None return None
...@@ -101,6 +102,7 @@ class ToolRunner: ...@@ -101,6 +102,7 @@ class ToolRunner:
"sourcename": "Genoscope", "sourcename": "Genoscope",
"date_executed": self.date "date_executed": self.date
}) })
def add_genome(self): def add_genome(self):
""" """
Tool invocation must be as return (or else it executes but returns nothing when called) Tool invocation must be as return (or else it executes but returns nothing when called)
...@@ -144,16 +146,7 @@ class ToolRunner: ...@@ -144,16 +146,7 @@ class ToolRunner:
""" """
return self.instance.tools.run_tool(tool_id="toolshed.g2.bx.psu.edu/repos/gga/chado_organism_delete_organisms/organism_delete_organisms/2.3.2", return self.instance.tools.run_tool(tool_id="toolshed.g2.bx.psu.edu/repos/gga/chado_organism_delete_organisms/organism_delete_organisms/2.3.2",
history_id=self.history, history_id=self.history,
tool_inputs={"organism": hs_id}) tool_inputs={"organism": str(hs_id)})
# def load_fasta(self, inputs_dict):
# """
#
# :return:
# """
# return self.instance.tools.run_tool(tool_id="toolshed.g2.bx.psu.edu/repos/gga/chado_feature_load_fasta/feature_load_fasta/2.3.2",
# history_id=self.history,
# tool_inputs=)
def show_pannel(self): def show_pannel(self):
print(self.instance.tools.get_tool_panel()) print(self.instance.tools.get_tool_panel())
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment