root/galaxy-central/tools/rgenetics/rgHaploView.xml

リビジョン 2, 7.4 KB (コミッタ: hatakeyama, 14 年 前)

import galaxy-central

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