# 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)