# HG changeset patch # User galaxyp # Date 1565870872 14400 # Node ID 33c911e66d8e3f56cee09160a7e45edbf1ba6362 # Parent 461ba35a313fd5f31165dc546ce811b9dbf8cce3 "planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/maxquant commit 74b5aa29e49deaaebe19ce2355a70d4f570f4951" diff -r 461ba35a313f -r 33c911e66d8e init.py --- a/init.py Wed Aug 07 10:49:12 2019 -0400 +++ b/init.py Thu Aug 15 08:07:52 2019 -0400 @@ -45,8 +45,7 @@ standard_mods = [] label_mods = [] for m in mods: - if (m.findtext('type') == 'Standard' - or m.findtext('type') == 'AaSubstitution'): + if (m.findtext('type') == 'Standard' or m.findtext('type') == 'AaSubstitution'): standard_mods.append(m.get('title')) elif m.findtext('type') == 'Label': label_mods.append(m.get('title')) diff -r 461ba35a313f -r 33c911e66d8e maxquant.xml --- a/maxquant.xml Wed Aug 07 10:49:12 2019 -0400 +++ b/maxquant.xml Thu Aug 15 08:07:52 2019 -0400 @@ -1,4 +1,4 @@ - + macros.xml @@ -310,17 +310,29 @@ - Thermo raw file or mzXML file - The datatype has to be 'thermo.raw' or 'mzXML'. Make sure to specify the correct datatype either during upload to Galaxy or afterwards (edit attributes --> datatypes) - Optional files: - - Tabular file with experimental design template: - - Currently four columns are needed: Name, Fraction, Experiment and PTM. The headers must have this exact naming. Name and Fraction are abitrary strings, Experiment is an integer, PTM is either True or False. + - Tabular file with experimental design template: + - Currently four columns are needed: Name, Fraction, Experiment and PTM. The headers must have this exact naming. Name and Experiment are abitrary strings, Fraction is an integer or emtpy, PTM is either 'True', 'False' or empty. Consider you uploaded files named File1.mzxml, ..., File5.mzxml. This is a (syntactically) correct experimental design template: :: Name Fraction Experiment PTM File1 1 E1 False File2 2 E1 False + ghost 234 none File3 3 E1 False - ... - ... + File4 E2 true + File5 E1 + + - This is the counter-example with one error per line: + + :: + + Name Fraction Experiment PTM + File1 1 E1 no (wrong PTM value) + File2.mzxml 1 E2 (filename with extension) + File3 f3 E1 (fraction not an integer) + File4 1 (missing experiment) + (File5 misssing) **Parameter Options** diff -r 461ba35a313f -r 33c911e66d8e mqparam.py --- a/mqparam.py Wed Aug 07 10:49:12 2019 -0400 +++ b/mqparam.py Thu Aug 15 08:07:52 2019 -0400 @@ -65,6 +65,27 @@ child = ET.SubElement(el, name, attrib=attrib if attrib else {}) child.text = str(text) + @staticmethod + def _check_validity(design, len_infiles): + "Perform some checks on the exp. design template" + design_len = len(design['Name']) + match = len(list(filter(lambda x: bool(x), design['Name']))) + if match < len_infiles: + raise Exception("Error parsing experimental design template: " + + "Found only {} matching entries ".format(design_len) + + "for {} input files".format(len_infiles)) + for i in range(0, design_len): + msg = "Error in line " + str(i + 2) + " of experimental design: " + if not (design['Name'][i] and design['Experiment'][i]): + raise Exception(msg + " Name or Experiment is empty.") + if design['PTM'][i].lower() not in ('true', 'false'): + raise Exception(msg + "Defines invalid PTM value, " + + "should be 'True' or 'False'.") + try: + int(design['Fraction'][i]) + except ValueError as e: + raise Exception(msg + str(e)) + def _make_exp_design(self, infiles): """Create a dict representing an experimental design from an experimental design template and a list of input files. @@ -80,6 +101,7 @@ >>> design['Fraction'] ['1', '2'] """ + design = {s: [] for s in ("Name", "PTM", "Fraction", "Experiment")} if not self.exp_design: design["Name"] = infiles @@ -94,11 +116,16 @@ if i in design: index.append(i) else: - raise Exception("Invalid comlumn index in experimental" + raise Exception("Invalid column index in experimental" + " design template: {}".format(i)) + for line in design_file: row = line.strip().split('\t') for e, i in zip_longest(row, index): + if i == "Fraction" and e == '': + e = 32767 + elif i == "PTM" and not e: + e = 'False' design[i].append(e) # map infiles to names in exp. design template @@ -117,6 +144,7 @@ else None) # replace orig. file names with matching links to galaxy datasets design['Name'] = names + MQParam._check_validity(design, len(infiles)) return design @@ -223,6 +251,8 @@ for index in range(len(files)): filepath = '' + files[index] + identifier = identifier.replace('<', '<') + description = description.replace('<', '<') fasta = self.fasta_template.replace('', filepath) fasta = fasta.replace('', '' + identifier) diff -r 461ba35a313f -r 33c911e66d8e mqwrapper.py --- a/mqwrapper.py Wed Aug 07 10:49:12 2019 -0400 +++ b/mqwrapper.py Thu Aug 15 08:07:52 2019 -0400 @@ -35,7 +35,7 @@ 'summary') # arguments for mqparam -## global +# global global_flags = ('calc_peak_properties', 'write_mztab', 'ibaq', @@ -51,7 +51,7 @@ 'min_peptide_len', 'max_peptide_mass') -## parameter group specific +# parameter group specific param_group_flags = ('lfq_skip_norm',) param_group_simple_args = ('missed_cleavages',