Mercurial > repos > xuebing > sharplabtool

comparison tools/rgenetics/rgHaploView.xml @ 0:9071e359b9a3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author xuebing
date Fri, 09 Mar 2012 19:37:19 -0500
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:9071e359b9a3
1 <tool id="rgHaploView1" name="LD plots:" version="0.3">
2
3 <description>and comparisons with HapMap data</description>
4
5 <command interpreter="python">
6 rgHaploView.py "$ucsc_region" "$rslist" "$title" "$out_file1"
7 "$lhistIn.extra_files_path" "$lhistIn.metadata.base_name"
8 "$minmaf" "$maxdist" "$ldtype" "$hires" "$memsize" "$out_file1.files_path"
9 "$infoTrack" "$tagr2" "$hmpanel" ${GALAXY_DATA_INDEX_DIR}/shared/jars/haploview.jar
10 </command>
11
12 <inputs>
13
14 <param name="lhistIn" type="data" format="lped"
15 label="Current history lPed format data"
16 size="80" help="Linkage Ped format data from your current history" />
17
18 <param name="ucsc_region" type="text" label="Optional subset region (blank=ALL. WARNING: doing this will fail if >1 chromosome in input data!)"
19 size="80" optional="true"
20 help="Region eg: chr9:119,506,000-119,518,000 Leave blank for all or to extract the rs list supplied below."/>
21
22 <param name="rslist" type="text" area='true' size='5x20' label="rs list" optional="true"
23 help="List of rs numbers to select - cut and paste or type, use space delimiters. Leave blank to extract region supplied above." />
24
25 <param name="title" type="text" size="80" label="Title for output files" optional="true"
26 help="Descriptive title for new genotype/map files" value="LD Plots" />
27
28 <param name="ldtype" type="select" label="Type of LD measure to estimate and plot"
29 size="80" help="" >
30 <option value="RSQ" selected="True">rsquared (default)</option>
31 <option value="DEFAULT">D prime</option>
32 <option value="DPALT">D prime alternative</option>
33 <option value="GAB">Gabriel</option>
34 <option value="GAM">4 Gamete test</option>
35 </param>
36
37 <param name="minmaf" type="float" label = "Minimum minor allele frequency to use" value="0.05"
38 help="If &gt; 0.0, markers below this MAF will be ignored for calculations"/>
39
40 <param name="maxdist" type="integer" label = "Maximum distance (kbp) between markers for LD estimate"
41 value="200" help="If &lt; &gt; 0, only marker pairs at or below this distance will have LD calculated"/>
42
43 <param name="hmpanel" type="select" multiple="true" label="Hapmap panels to compare"
44 size="40" help="HapMap data LD plots will also be produced for each selected population panel" >
45 <option value='CEU' selected="True">CEPH (European) (default)</option>
46 <option value='YRI'>Yoruba (African)</option>
47 <option value='CHB+JPT'>Chinese + Japanese</option>
48 <option value="">(None - no comparison)</option>
49 </param>
50 <param name="tagr2" type="float" label = "rsquared threshold for tagging outputs" value="0.8"
51 help="Tagging output will use this value as the minimum rsquared threshold"/>
52
53 <param name="infoTrack" type="select" label="Add Hapmap information track to image"
54 help="Refseq genes and snp density can be added to the plot if desired for orientation" >
55 <option value="info">Add Information track (DISABLED! Awaiting bug fix from Haploview authors since reported in October 2009)</option>
56 <option value="noinfo" selected="True">No Information track</option>
57 </param>
58
59 <param name="hires" type="select" label="High resolution plots"
60 help="A high resolution plot file may be possible but only for small regions - not reliable &gt;100's of snps">
61 <option value="hi">High resolution - only a few (hundreds of) markers</option>
62 <option value="lo" selected="True">Low resolution - large number of markers</option>
63 </param>
64
65 <param name="memsize" type="select" label="System RAM to allocate"
66 size="80" help="Very large files will need extra memory (java is a bit of a pig)" >
67 <option value="1024">1GB</option>
68 <option value="2048" selected="True">2GB (default)</option>
69 <option value="4096">4GB</option>
70 <option value="6144">6GB</option>
71 <option value="8192">8GB</option>
72 </param>
73
74 </inputs>
75
76 <outputs>
77 <data format="html" name="out_file1" label="${title}.html" />
78 </outputs>
79
80 <!-- python $TOOLPATH/$TOOL.py "" "rs2283802Xrs2267000Xrs16997606Xrs4820537Xrs3788347Xrs756632Xrs4820539Xrs2283804Xrs2267006Xrs4822363X" \
81 "$NPRE" $OUTPATH/${NPRE}.html "test" "" "$INPATH" "tinywga" 0.0 200000 "RSQ" "lo" "2048" "$OUTPATH" "hg18" "noinfo" "0.8" \
82 "['CEU','YRI','CHB+JPT']" $BINPATH/haploview.jar -->
83 <tests>
84 <test>
85 <param name='lhistIn' value='tinywga' ftype='lped' >
86 <metadata name='base_name' value='tinywga' />
87 <composite_data value='tinywga.ped' />
88 <composite_data value='tinywga.map' />
89 <edit_attributes type='name' value='tinywga' />
90 </param>
91 <param name='ucsc_region' value='' />
92 <param name='title' value='rgHaploViewtest1' />
93 <param name='rslist' value="rs2283802 rs2267000 rs16997606 rs4820537 rs3788347 rs756632Xrs4820539 rs2283804 rs2267006 rs4822363" />
94 <param name='ldtype' value='RSQ' />
95 <param name='minmaf' value='0.0' />
96 <param name='maxdist' value='200000' />
97 <param name='tagr2' value='0.8' />
98 <param name='hmpanel' value="YRI" />
99 <param name='infoTrack' value='noinfo' />
100 <param name='hires' value='lo' />
101 <param name='memsize' value='2048' />
102 <output name='out_file1' file='rgtestouts/rgHaploView/rgHaploViewtest1.html' ftype='html' lines_diff="60">
103 <extra_files type="file" name='alljoin.pdf' value="rgtestouts/rgHaploView/alljoin.pdf" compare="sim_size" delta="50000"/>
104 <extra_files type="file" name='allnup.pdf' value="rgtestouts/rgHaploView/allnup.pdf" compare="sim_size" delta="50000" />
105 <extra_files type="file" name='Log_rgHaploViewtest1.txt' value="rgtestouts/rgHaploView/Log_rgHaploViewtest1.txt" compare="diff" lines_diff="50"/>
106 <extra_files type="file" name='rgHaploViewtest1.ped.TESTS' value="rgtestouts/rgHaploView/rgHaploViewtest1.ped.TESTS" compare="diff"
107 lines_diff="20"/>
108 <extra_files type="file" name='rgHaploViewtest1.ped.TAGS' value="rgtestouts/rgHaploView/rgHaploViewtest1.ped.TAGS" compare="diff"
109 lines_diff="20" />
110 </output>
111 </test>
112 </tests>
113
114 <help>
115
116 .. class:: infomark
117
118 **Note**
119
120 The input file must be in linkage ped format. A suitable file can be chosen from the system library,
121 or from the files already imported into your current history. Use either one of the selection boxes to
122 make your choice.
123
124 -----
125
126 **Syntax**
127
128 - **Library Linkage Ped** is a linkage format pedigree file chosen from the system file Library
129 - **History Linkage Ped** is a linkage format pedigree file chosen from your current Galaxy History
130 - **Region** is the genomic region cut and paste from a UCSC browser location window
131 - **Genome Build** is the version of the genome your markers are from - use hg18 for CAMP illumina data
132
133 -----
134
135 **Summary**
136
137 This tool is a special purpose tool to estimate and plot linkage disequilibrium estimated
138 from genotype data in linkage pedigree format (separate map file). All markers in the input file
139 are used as the default. To limit the calculations to a subset of the input data, supply
140 a specified genomic region in UCSC browser location format or a list of specific marker IDs.
141
142 Note that you can choose either a file of the correct type (linkage pedigree - lped) from
143 your current history **or** from the system library
144
145 This tool currently calls Haploview for estimation and plots. For full attribution, source code and documentation, see
146 http://www.broad.mit.edu/mpg/haploview/index.php
147
148 Copyright, Ross Lazarus, April 2008 for the Rgenetics project
149 Released under the LGPL. See http://www.gnu.org/licenses/lgpl.html for license terms.
150
151 </help>
152 </tool>