# HG changeset patch
# User iuc
# Date 1540940222 14400
# Node ID 5f589c91566e2700b0d1bee23c5e6e8aee199df0
# Parent e41ec5af7472be072410c13cf014696499b38b5d
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/heinz commit a4b0969b33a68a0ea9ba12291f6694aec24f13ed
diff -r e41ec5af7472 -r 5f589c91566e heinz.xml
--- a/heinz.xml Thu Aug 02 11:57:44 2018 -0400
+++ b/heinz.xml Tue Oct 30 18:57:02 2018 -0400
@@ -4,13 +4,13 @@
heinz
'$subnetwork'
+heinz -m \${GALAXY_SLOTS:-2} -n '$score' -e '$edge' > '$subnetwork'
]]>
-
-
+
+
@@ -43,16 +43,15 @@
the second, Heinz score.
========= ===================
-BRCA2 -6.991782933819368
+BRCA2 -6.991782933819368
--------- -------------------
-BRCA1 -5.206139799106934
+BRCA1 -5.206139799106934
--------- -------------------
-AACS -0.9299868303078357
+AACS -0.9299868303078357
--------- -------------------
ABCC11 -5.845009850430119
========= ===================
-
Edge file: the background network Heinz uses in the form of a list of edges; each line is made up of
two node identifiers (e.g., genes, KEGG ORTHOLOGY (KO)) delimited by a tab. In practice, we could
prepare this file using different pathway databases, such as Reactome, STRING and KEGG. Which database
@@ -67,8 +66,6 @@
--------- ---------
NPTX1 CIAO1
========= =========
-
-
]]>
10.1093/bioinformatics/btn161
diff -r e41ec5af7472 -r 5f589c91566e heinz_scoring.py
--- a/heinz_scoring.py Thu Aug 02 11:57:44 2018 -0400
+++ b/heinz_scoring.py Tue Oct 30 18:57:02 2018 -0400
@@ -85,10 +85,7 @@
# Calculate the scores
input_pvalues = pd.read_csv(args.nodes, sep='\t', names=['node', 'pvalue'])
input_pvalues.loc[:, 'score'] = input_pvalues.pvalue.apply(lambda x:
- (alpha - 1) *
- (np.log(x) -
- np.log(
- p_threshold)))
+ (alpha - 1) * (np.log(x) - np.log(p_threshold)))
# print(input_pvalues.loc[:, ['node', 'score']])
input_pvalues.loc[:, ['node', 'score']].to_csv(args.output, sep='\t',
index=False, header=False)