Mercurial > repos > bgruening > hicexplorer_hicdetectloops

comparison hicDetectLoops.xml @ 0:2d1f9759fa0e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/hicexplorer commit 3b41d687ff30583540d055f6995de00530cca81d-dirty"
author bgruening
date Mon, 16 Dec 2019 15:49:32 -0500
parents
children ce8d154fc6b5
comparison
equal deleted inserted replaced
-1:000000000000 0:2d1f9759fa0e
1 <tool id="hicexplorer_hicdetectloops" name="@BINARY@" version="@WRAPPER_VERSION@.0">
2 <description>searches for enriched regions</description>
3 <macros>
4 <token name="@BINARY@">hicDetectLoops</token>
5 <import>macros.xml</import>
6 </macros>
7 <expand macro="requirements" />
8 <command detect_errors="exit_code"><![CDATA[
9 ln -s '$matrix_h5_cooler' 'matrix.$matrix_h5_cooler.ext' &&
10 @BINARY@
11
12 --matrix 'matrix.$matrix_h5_cooler.ext'
13
14 #if $peakWidth:
15 --peakWidth $peakWidth
16 #end if
17
18 #if $windowSize:
19 --windowSize $windowSize
20 #end if
21
22 #if $pValuePreselection:
23 --pValuePreselection $pValuePreselection
24 #end if
25
26 #if $peakInteractionsThreshold:
27 --peakInteractionsThreshold $peakInteractionsThreshold
28 #end if
29
30 #if $maximumInteractionPercentageThreshold:
31 --maximumInteractionPercentageThreshold $maximumInteractionPercentageThreshold
32 #end if
33 #if $pValue:
34 --pValue $pValue
35 #end if
36 #if $maxLoopDistance:
37 --maxLoopDistance $maxLoopDistance
38 #end if
39
40 #if $region:
41 --region '$region'
42 #end if
43
44 #if $chromosomes:
45 #set $chromosome = ' '.join([ '\'%s\'' % $chrom for $chrom in str($chromosomes).split(' ') ])
46 --chromosomes $chromosome
47 #end if
48
49 --statisticalTest $statisticalTest_selector
50
51 --outFileName output_loop.bedgraph
52
53 --threads @THREADS@
54 ]]>
55 </command>
56 <inputs>
57 <expand macro="matrix_h5_cooler_macro" />
58 <param argument="--peakWidth" type="integer" optional='true' label="Peak width" help= "The width of the peak region in bins. The square around the peak will include (2 * peakWidth)^2 bins." />
59 <param argument="--windowSize" type="integer" optional='true' label="Window size" help= "The window size for the neighborhood region the peak is located in. All values from this region (exclude the values from the peak
60 region) are tested against the peak region for significant difference. The square will have the size of (2 * windowSize)^2 bins" />
61 <param argument="--pValuePreselection" type="float" label="P-value preselection" help= "Only candidates with p-values less the given threshold will be considered as candidates.
62 For each genomic distance a negative binomial distribution is fitted and for each pixel a p-value given by the cumulative density function is given.
63 This does NOT influence the p-value for the neighborhood testing." value='0.05'/>
64 <param argument="--peakInteractionsThreshold" type="integer" label="Minimum interaction number" help= "The minimum number of interactions a detected peaks needs to have to be considered." value='5' />
65 <param argument="--maximumInteractionPercentageThreshold" type="float" value='0.1' label="Maximum interaction share" help= "For each genomic distance the maximum value is considered and all candidates need to have at least \'max_value * maximumInteractionPercentageThreshold\' interactions." />
66 <param argument="--pValue" type="float" label="P-value" help= "Rejection level for the statistical test for H0. H0 is peak region and background have the same distribution." value='0.05'/>
67 <param argument="--maxLoopDistance" optional='true' type="integer" label="Maximal loop distance" help= "Maximum genomic distance of a loop, usually loops are within a distance of ~2MB." value='2000000'/>
68 <param argument="--minLoopDistance" optional='true' type="integer" label="Maximal loop distance" help= "Minimum genomic distance of a loop to be considered." value='100000'/>
69 <param argument="--chromosomes" optional='true' type="text" label="Chromosomes to include" help= "Chromosomes to include in the analysis. If not set, all chromosomes are included." />
70 <param argument="--region" optional='true' type="text" label="Chromosomes to include" help= "The format is chr:start-end." />
71 <param name="statisticalTest_selector" type="select" label="Stistical test">
72 <option value="wilcoxon-rank-sum" selected="True">Wilcoxon rank-sum'</option>
73 <option value="anderson-darling">Anderson-Darling</option>
74 </param>
75 </inputs>
76 <outputs>
77 <data name='output_loops' from_work_dir='output_loop.bedgraph' format='bedgraph' label='Computed loops'/>
78 </outputs>
79 <tests>
80 <test>
81 <param name="matrix_h5_cooler" value="small_test_matrix.cool"/>
82 <param name="maxLoopDistance" value="30000000"/>
83 <param name="windowSize" value="5"/>
84 <param name="peakWidth" value="2"/>
85 <param name="pValue" value="0.5"/>
86 <param name="pValuePreselection" value="0.55"/>
87 <!-- <param name="chromosomes" value="X"/> -->
88 <output name="output_loops" file="hicDetectLoops/loops.bedgraph" ftype="bedgraph" compare="sim_size"/>
89 </test>
90 </tests>
91 <help><![CDATA[
92
93 Loop detection
94 ==============
95
96 Computes enriched regions (peaks) or long range contacts on the given contact matrix.
97
98 For more information about HiCExplorer please consider our documentation on readthedocs.io_.
99
100 .. _readthedocs.io: http://hicexplorer.readthedocs.io/en/latest/index.html
101
102 ]]></help>
103 <expand macro="citations" />
104 </tool>