annotate findPeaks.xml @ 16:687df269e597 draft

Uploaded
author kevyin
date Wed, 19 Dec 2012 17:28:55 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
16
687df269e597 Uploaded
kevyin
parents:
diff changeset
1 <tool id="homer_findPeaks" name="homer_findPeaks" version="0.1.2">
687df269e597 Uploaded
kevyin
parents:
diff changeset
2 <requirements>
687df269e597 Uploaded
kevyin
parents:
diff changeset
3 <requirement type="package" version="4.1">homer</requirement>
687df269e597 Uploaded
kevyin
parents:
diff changeset
4 </requirements>
687df269e597 Uploaded
kevyin
parents:
diff changeset
5 <description>Homer's peakcaller. Requires tag directories (see makeTagDirectory)</description>
687df269e597 Uploaded
kevyin
parents:
diff changeset
6 <!--<version_command></version_command>-->
687df269e597 Uploaded
kevyin
parents:
diff changeset
7 <command>
687df269e597 Uploaded
kevyin
parents:
diff changeset
8 findPeaks $tagDir.extra_files_path $options -o $outputPeakFile
687df269e597 Uploaded
kevyin
parents:
diff changeset
9
687df269e597 Uploaded
kevyin
parents:
diff changeset
10 #if $control_tagDir:
687df269e597 Uploaded
kevyin
parents:
diff changeset
11 -i $control_tagDir.extra_files_path
687df269e597 Uploaded
kevyin
parents:
diff changeset
12 #end if
687df269e597 Uploaded
kevyin
parents:
diff changeset
13
687df269e597 Uploaded
kevyin
parents:
diff changeset
14 2&gt; $out_log || echo "Error running findPeaks." >&amp;2
687df269e597 Uploaded
kevyin
parents:
diff changeset
15 </command>
687df269e597 Uploaded
kevyin
parents:
diff changeset
16 <inputs>
687df269e597 Uploaded
kevyin
parents:
diff changeset
17 <param format="homerTagDirectory" name="tagDir" type="data" label="tag directory" help="Must be made with homer_makeTagDirectory" />
687df269e597 Uploaded
kevyin
parents:
diff changeset
18 <param format="homerTagDirectory" name="control_tagDir" type="data" optional="True" label="Control tag directory" help="Must be made with homer_makeTagDirectory" />
687df269e597 Uploaded
kevyin
parents:
diff changeset
19 <param type="text" name="options" label="Extra options" value="" help="See link below for more options">
687df269e597 Uploaded
kevyin
parents:
diff changeset
20 <sanitizer>
687df269e597 Uploaded
kevyin
parents:
diff changeset
21 <valid initial="string.printable">
687df269e597 Uploaded
kevyin
parents:
diff changeset
22 <remove value="&apos;"/>
687df269e597 Uploaded
kevyin
parents:
diff changeset
23 <remove value="/"/>
687df269e597 Uploaded
kevyin
parents:
diff changeset
24 </valid>
687df269e597 Uploaded
kevyin
parents:
diff changeset
25 <mapping initial="none">
687df269e597 Uploaded
kevyin
parents:
diff changeset
26 <add source="&apos;" target="__sq__"/>
687df269e597 Uploaded
kevyin
parents:
diff changeset
27 </mapping>
687df269e597 Uploaded
kevyin
parents:
diff changeset
28 </sanitizer>
687df269e597 Uploaded
kevyin
parents:
diff changeset
29 </param>
687df269e597 Uploaded
kevyin
parents:
diff changeset
30 </inputs>
687df269e597 Uploaded
kevyin
parents:
diff changeset
31 <outputs>
687df269e597 Uploaded
kevyin
parents:
diff changeset
32 <!--<data format="html" name="html_outfile" label="index" />-->
687df269e597 Uploaded
kevyin
parents:
diff changeset
33 <!--<data format="html" hidden="True" name="html_outfile" label="index.html" />-->
687df269e597 Uploaded
kevyin
parents:
diff changeset
34 <data format="txt" name="outputPeakFile" label="${tool.name} on #echo os.path.splitext(str($tagDir.name))[0]#.txt" />
687df269e597 Uploaded
kevyin
parents:
diff changeset
35 <data format="txt" name="out_log" label="${tool.name} on #echo os.path.splitext(str($tagDir.name))[0]#.log" />
687df269e597 Uploaded
kevyin
parents:
diff changeset
36 </outputs>
687df269e597 Uploaded
kevyin
parents:
diff changeset
37 <tests>
687df269e597 Uploaded
kevyin
parents:
diff changeset
38 <test>
687df269e597 Uploaded
kevyin
parents:
diff changeset
39 <!--<param name="input_file" value="extract_genomic_dna.fa" />-->
687df269e597 Uploaded
kevyin
parents:
diff changeset
40 <!--<output name="html_file" file="sample_output.html" ftype="html" />-->
687df269e597 Uploaded
kevyin
parents:
diff changeset
41 </test>
687df269e597 Uploaded
kevyin
parents:
diff changeset
42 </tests>
687df269e597 Uploaded
kevyin
parents:
diff changeset
43
687df269e597 Uploaded
kevyin
parents:
diff changeset
44 <help>
687df269e597 Uploaded
kevyin
parents:
diff changeset
45
687df269e597 Uploaded
kevyin
parents:
diff changeset
46 .. class:: infomark
687df269e597 Uploaded
kevyin
parents:
diff changeset
47
687df269e597 Uploaded
kevyin
parents:
diff changeset
48 **Homer findPeaks**
687df269e597 Uploaded
kevyin
parents:
diff changeset
49
687df269e597 Uploaded
kevyin
parents:
diff changeset
50 For more options, look under: "Command line options for findPeaks"
687df269e597 Uploaded
kevyin
parents:
diff changeset
51
687df269e597 Uploaded
kevyin
parents:
diff changeset
52 http://biowhat.ucsd.edu/homer/ngs/peaks.html
687df269e597 Uploaded
kevyin
parents:
diff changeset
53
687df269e597 Uploaded
kevyin
parents:
diff changeset
54 TIP: use homer_bed2pos and homer_pos2bed to convert between the homer peak positions and the BED format.
687df269e597 Uploaded
kevyin
parents:
diff changeset
55
687df269e597 Uploaded
kevyin
parents:
diff changeset
56 **Parameter list**
687df269e597 Uploaded
kevyin
parents:
diff changeset
57
687df269e597 Uploaded
kevyin
parents:
diff changeset
58 Command line options (not all of them are supported)::
687df269e597 Uploaded
kevyin
parents:
diff changeset
59
687df269e597 Uploaded
kevyin
parents:
diff changeset
60 Usage: findPeaks &lt;tag directory&gt; [options]
687df269e597 Uploaded
kevyin
parents:
diff changeset
61
687df269e597 Uploaded
kevyin
parents:
diff changeset
62 Finds peaks in the provided tag directory. By default, peak list printed to stdout
687df269e597 Uploaded
kevyin
parents:
diff changeset
63
687df269e597 Uploaded
kevyin
parents:
diff changeset
64 General analysis options:
687df269e597 Uploaded
kevyin
parents:
diff changeset
65 -o &lt;filename|auto&gt; (file name for to output peaks, default: stdout)
687df269e597 Uploaded
kevyin
parents:
diff changeset
66 &quot;-o auto&quot; will send output to &quot;&lt;tag directory&gt;/peaks.txt&quot;, &quot;.../regions.txt&quot;,
687df269e597 Uploaded
kevyin
parents:
diff changeset
67 or &quot;.../transcripts.txt&quot; depending on the &quot;-style&quot; option
687df269e597 Uploaded
kevyin
parents:
diff changeset
68 -style &lt;option&gt; (Specialized options for specific analysis strategies)
687df269e597 Uploaded
kevyin
parents:
diff changeset
69 factor (transcription factor ChIP-Seq, uses -center, output: peaks.txt, default)
687df269e597 Uploaded
kevyin
parents:
diff changeset
70 histone (histone modification ChIP-Seq, region based, uses -region -size 500 -L 0, regions.txt)
687df269e597 Uploaded
kevyin
parents:
diff changeset
71 groseq (de novo transcript identification from GroSeq data, transcripts.txt)
687df269e597 Uploaded
kevyin
parents:
diff changeset
72 tss (TSS identification from 5&apos; RNA sequencing, tss.txt)
687df269e597 Uploaded
kevyin
parents:
diff changeset
73 dnase (Hypersensitivity [crawford style (nicking)], peaks.txt)
687df269e597 Uploaded
kevyin
parents:
diff changeset
74
687df269e597 Uploaded
kevyin
parents:
diff changeset
75 chipseq/histone options:
687df269e597 Uploaded
kevyin
parents:
diff changeset
76 -i &lt;input tag directory&gt; (Experiment to use as IgG/Input/Control)
687df269e597 Uploaded
kevyin
parents:
diff changeset
77 -size &lt;#&gt; (Peak size, default: auto)
687df269e597 Uploaded
kevyin
parents:
diff changeset
78 -minDist &lt;#&gt; (minimum distance between peaks, default: peak size x2)
687df269e597 Uploaded
kevyin
parents:
diff changeset
79 -gsize &lt;#&gt; (Set effective mappable genome size, default: 2e9)
687df269e597 Uploaded
kevyin
parents:
diff changeset
80 -fragLength &lt;#|auto&gt; (Approximate fragment length, default: auto)
687df269e597 Uploaded
kevyin
parents:
diff changeset
81 -inputFragLength &lt;#|auto&gt; (Approximate fragment length of input tags, default: auto)
687df269e597 Uploaded
kevyin
parents:
diff changeset
82 -tbp &lt;#&gt; (Maximum tags per bp to count, 0 = no limit, default: auto)
687df269e597 Uploaded
kevyin
parents:
diff changeset
83 -inputtbp &lt;#&gt; (Maximum tags per bp to count in input, 0 = no limit, default: auto)
687df269e597 Uploaded
kevyin
parents:
diff changeset
84 -strand &lt;both|separate&gt; (find peaks using tags on both strands or separate, default:both)
687df269e597 Uploaded
kevyin
parents:
diff changeset
85 -norm # (Tag count to normalize to, default 10000000)
687df269e597 Uploaded
kevyin
parents:
diff changeset
86 -region (extends start/stop coordinates to cover full region considered &quot;enriched&quot;)
687df269e597 Uploaded
kevyin
parents:
diff changeset
87 -center (Centers peaks on maximum tag overlap and calculates focus ratios)
687df269e597 Uploaded
kevyin
parents:
diff changeset
88 -nfr (Centers peaks on most likely nucleosome free region [works best with mnase data])
687df269e597 Uploaded
kevyin
parents:
diff changeset
89 (-center and -nfr can be performed later with &quot;getPeakTags&quot;
687df269e597 Uploaded
kevyin
parents:
diff changeset
90
687df269e597 Uploaded
kevyin
parents:
diff changeset
91 Peak Filtering options: (set -F/-L/-C to 0 to skip)
687df269e597 Uploaded
kevyin
parents:
diff changeset
92 -F &lt;#&gt; (fold enrichment over input tag count, default: 4.0)
687df269e597 Uploaded
kevyin
parents:
diff changeset
93 -P &lt;#&gt; (poisson p-value threshold relative to input tag count, default: 0.0001)
687df269e597 Uploaded
kevyin
parents:
diff changeset
94 -L &lt;#&gt; (fold enrichment over local tag count, default: 4.0)
687df269e597 Uploaded
kevyin
parents:
diff changeset
95 -LP &lt;#&gt; (poisson p-value threshold relative to local tag count, default: 0.0001)
687df269e597 Uploaded
kevyin
parents:
diff changeset
96 -C &lt;#&gt; (fold enrichment limit of expected unique tag positions, default: 2.0)
687df269e597 Uploaded
kevyin
parents:
diff changeset
97 -localSize &lt;#&gt; (region to check for local tag enrichment, default: 10000)
687df269e597 Uploaded
kevyin
parents:
diff changeset
98 -inputSize &lt;#&gt; (Size of region to search for control tags, default: 2x peak size)
687df269e597 Uploaded
kevyin
parents:
diff changeset
99 -fdr &lt;#&gt; (False discovery rate, default = 0.001)
687df269e597 Uploaded
kevyin
parents:
diff changeset
100 -poisson &lt;#&gt; (Set poisson p-value cutoff, default: uses fdr)
687df269e597 Uploaded
kevyin
parents:
diff changeset
101 -tagThreshold &lt;#&gt; (Set # of tags to define a peak, default: 25)
687df269e597 Uploaded
kevyin
parents:
diff changeset
102 -ntagThreshold &lt;#&gt; (Set # of normalized tags to define a peak, by default uses 1e7 for norm)
687df269e597 Uploaded
kevyin
parents:
diff changeset
103 -minTagThreshold &lt;#&gt; (Absolute minimum tags per peak, default: expected tags per peak)
687df269e597 Uploaded
kevyin
parents:
diff changeset
104
687df269e597 Uploaded
kevyin
parents:
diff changeset
105 GroSeq Options: (Need to specify &quot;-style groseq&quot;):
687df269e597 Uploaded
kevyin
parents:
diff changeset
106 -tssSize &lt;#&gt; (size of region for initiation detection/artifact size, default: 250)
687df269e597 Uploaded
kevyin
parents:
diff changeset
107 -minBodySize &lt;#&gt; (size of regoin for transcript body detection, default: 1000)
687df269e597 Uploaded
kevyin
parents:
diff changeset
108 -maxBodySize &lt;#&gt; (size of regoin for transcript body detection, default: 10000)
687df269e597 Uploaded
kevyin
parents:
diff changeset
109 -tssFold &lt;#&gt; (fold enrichment for new initiation dectection, default: 4.0)
687df269e597 Uploaded
kevyin
parents:
diff changeset
110 -bodyFold &lt;#&gt; (fold enrichment for new transcript dectection, default: 4.0)
687df269e597 Uploaded
kevyin
parents:
diff changeset
111 -endFold &lt;#&gt; (end transcript when levels are this much less than the start, default: 10.0)
687df269e597 Uploaded
kevyin
parents:
diff changeset
112 -fragLength &lt;#&gt; (Approximate fragment length, default: 150)
687df269e597 Uploaded
kevyin
parents:
diff changeset
113 -uniqmap &lt;directory&gt; (directory of binary files specifying uniquely mappable locations)
687df269e597 Uploaded
kevyin
parents:
diff changeset
114 Download from http://biowhat.ucsd.edu/homer/groseq/
687df269e597 Uploaded
kevyin
parents:
diff changeset
115 -confPvalue &lt;#&gt; (confidence p-value: 1.00e-05)
687df269e597 Uploaded
kevyin
parents:
diff changeset
116 -minReadDepth &lt;#&gt; (Minimum initial read depth for transcripts, default: auto)
687df269e597 Uploaded
kevyin
parents:
diff changeset
117 -pseudoCount &lt;#&gt; (Pseudo tag count, default: 2.0)
687df269e597 Uploaded
kevyin
parents:
diff changeset
118 -gtf &lt;filename&gt; (Output de novo transcripts in GTF format)
687df269e597 Uploaded
kevyin
parents:
diff changeset
119 &quot;-o auto&quot; will produce &lt;dir&gt;/transcripts.txt and &lt;dir&gt;/transcripts.gtf
687df269e597 Uploaded
kevyin
parents:
diff changeset
120 </help>
687df269e597 Uploaded
kevyin
parents:
diff changeset
121 </tool>
687df269e597 Uploaded
kevyin
parents:
diff changeset
122