Mercurial > repos > kpbioteam > chipeakanno_annopeaks

comparison chipeakanno_annopeaks.xml @ 0:95ad07b229df draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author kpbioteam
date Sun, 28 Apr 2019 03:43:31 -0400
parents
children 6f52c0dc56d4
comparison
equal deleted inserted replaced
-1:000000000000 0:95ad07b229df
1 <tool id="chipeakanno_annopeaks" name="ChIPpeakAnno annoPeaks" version="0.1.0">
2 <description>annotate peaks by annoGR object in the given range</description>
3 <requirements>
4 <requirement type="package" version="3.16.0-0">bioconductor-chippeakanno</requirement>
5 <requirement type="package" version="2.99.0-3">bioconductor-ensdb.hsapiens.v75</requirement>
6 <requirement type="package" version="1.42.1-1">bioconductor-rtracklayer</requirement>
7 </requirements>
8 <command detect_errors="exit_code"><![CDATA[
9 Rscript '$chipeakanno_annopeaks_script'
10 ]]></command>
11 <configfiles>
12 <configfile name="chipeakanno_annopeaks_script"><![CDATA[
13 require("ChIPpeakAnno", quietly = TRUE)
14 require("EnsDb.Hsapiens.v75", quietly = TRUE)
15 require("rtracklayer", quietly = TRUE)
16
17 options(warn = -1)
18
19 data <- read.table('$dmp')
20
21 peaks <- GRanges(seqnames = data[, 1],
22 ranges = IRanges
23 (start = data[, 2], end = data[, 3]),score=data[, as.numeric('$score')])
24
25 db <- toGRanges(EnsDb.Hsapiens.v75)
26 seqlevelsStyle(peaks) <- seqlevelsStyle(db)
27
28 start <- as.numeric('$start')
29 end <- as.numeric('$end')
30
31 anno <- annoPeaks(peaks, annoData=db, bindingType =c('$type'), bindingRegion=c(start, end))
32
33 write.table(anno, file= '$tab', quote = FALSE, row.names = FALSE, sep = "\t")
34 ]]>
35 </configfile>
36 </configfiles>
37 <inputs>
38 <param type="data" name="dmp" format="interval" label="Differentially Methylated Data"/>
39 <param name="type" type="select" label="bindingType" help="Specifying the criteria to associate peaks with annotation e.g. startSite start position of the feature strand">
40 <option value="startSite">startSite</option>
41 <option value="endSite">endSite</option>
42 <option value="fullRange">fullRange</option>
43 <option value="nearestBiDirectionalPromoters">nearestBiDirectionalPromoters</option>
44 </param>
45 <param name="start" type="text" value="-5000" label="bindingRegionStart" help="e.g. to obtain peaks within 5kb upstream and up to 3kb downstream of TSS within the gene body, set bindingType = startSite and bindingRegion = c(-5000, 3000)"/>
46 <param name="end" type="text" value="3000" label="bindingRegionEnd" help="e.g. to obtain peaks with nearest bi-directional promoters within 5kb upstream and 3kb downstream of TSS, set bindingType = nearestBiDirectionalPromoters and bindingRegion = c(-5000, 3000)"/>
47
48 <param name="score" type="text" value="5" label="Additional Column of Score" help="Position of Additional Column of the Score" />
49 </inputs>
50 <outputs>
51 <data name="tab" format="tabular" label="Table of Annotated Peaks"/>
52 </outputs>
53 <tests>
54 <test>
55 <param name="dmp" value="Differentially_Methylated_Regions.interval"/>
56 <param name="type" value="startSite"/>
57 <param name="start" value="-5000"/>
58 <param name="end" value="5000"/>
59 <param name="score" value="5"/>
60 <output name="tab" file="annoTab.tabular"/>
61 </test>
62 </tests>
63 <help><![CDATA[
64 This tool annotate peaks from bed file.
65 bindingType
66 Specifying the criteria to associate peaks with annotation. Here is how to use it together with the parameter bindingRegion
67 To obtain peaks within 5kb upstream and up to 3kb downstream of TSS within the gene body, set bindingType = "startSite" and bindingRegion = c(-5000, 3000)
68 To obtain peaks up to 5kb upstream within the gene body and 3kb downstream of gene/Exon End, set bindingType = "endSite" and bindingRegion = c(-5000, 3000)
69 To obtain peaks from 5kb upstream to 3kb downstream of genes/Exons , set bindingType = "fullRange" and bindingRegion = c(-5000, 3000)
70 To obtain peaks with nearest bi-directional promoters within 5kb upstream and 3kb downstream of TSS, set bindingType = "nearestBiDirectionalPromoters" and bindingRegion = c(-5000, 3000)
71 startSite
72 start position of the feature (strand is considered)
73 endSite
74 end position of the feature (strand is considered)
75 fullRange
76 whole range of the feature
77 nearestBiDirectionalPromoters
78 nearest promoters from both direction of the peaks (strand is considered). It will report bidirectional promoters if there are promoters in both directions in the given region (defined by bindingRegion). Otherwise, it will report the closest promoter in one direction.
79 bindingRegion
80 Annotation range used together with bindingType, which is a vector with two integer values, default to c (-5000, 5000). The first one must be no bigger than 0. And the sec ond one must be no less than 1. For details, see bindingType.
81 ignore.peak.strand
82 ignore the peaks strand or not.
83 ]]></help>
84 <citations>
85 <citation type="doi">10.18129/B9.bioc.ChIPpeakAnno</citation>
86 </citations>
87 </tool>