/tools/expression/expPkg.xml

 1<tool name="Expression CEL file packager" id="expPkg">
 2  <description>can download .cel files from GEO by given GSM IDs and prepare a cel.zip file for expression analysis.</description>
 3  <command interpreter="command">/bin/bash $shscript</command>
 4  <inputs>
 5    <param name="group0name" type="text" label="Control group name"/>
 6    <repeat name="group0gsmid" title="GSM ID of control group">
 7      <param name="gsmid" type="text" label="GEO GSM ID"/>
 8    </repeat>
 9
10    <param name="group1name" type="text" label="Sample group name"/>
11    <repeat name="group1gsmid" title="GSM ID of sample group">
12      <param name="gsmid" type="text" label="GEO GSM ID"/>
13    </repeat>
14  </inputs>
15  <outputs>
16    <data format="cel.zip" name="output" label="Expression CEL file zip" />
17  </outputs>
18  <configfiles>
19    <configfile name="shscript">
20#!/bin/bash
21
22#set $tmp0 = ""
23#for $g0gsmid in $group0gsmid
24#set $tmp0 = $tmp0 + " 0:"+str($group0name).replace(" ","_").replace("\t","_").replace(":","_")+":GSM"+str($g0gsmid.gsmid).upper().lstrip("GSM")+" "
25#end for
26
27#set $tmp1 = ""
28#for $g1gsmid in $group1gsmid
29#set $tmp1 = $tmp1 + " 1:"+str($group1name).replace(" ","_").replace("\t","_").replace(":","_")+":GSM"+str($g1gsmid.gsmid).upper().lstrip("GSM")+" "
30#end for
31
32#set $dollar = chr(36)
33#set $gt = chr(62)
34#set $lt = chr(60)
35#set $ad = chr(38)
36
37expressPkgr.py $tmp0 $tmp1
38mv package.zip $output
39    </configfile>
40  </configfiles>
41  <help>
42This expression CEL file packager can prepare a cel.zip file for Cistrome expression array tools. This tool is only designed to fetch Affymetrix CEL files for expression array from NCBI GEO database. This script was written by Len Taing.
43
44The script can group CEL files into two groups -- a sample group and a control group. Typically, control group can be the wildtype, and sample group can be the samples after certain genes are knocked down/off. When the gene expression index results from the package generated by this script are used in differential expression analysis, the difference represents sample group vs control group, so if a fold change is 1.5, then the gene is highly expressed in sample group.
45
46Possible reasons of Error:
47
481. GSM ID can't be found on GEO site;
492. If you fetch many CEL files or there are many users downloading CEL files from our Cistrome server at the same time, connection to GEO site may be failed;
50
51  </help>
52</tool>