annotate minfi_dmp.xml @ 20:5aaf04dbc8bc draft

Deleted selected files
author kpbioteam
date Fri, 22 Feb 2019 08:24:52 -0500
parents cffa2871fe8b
children 3d083550977f
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
17
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
1 <tool id="minfi_dmp" name="Minfi DMP" version="@MINFI_VERSION@">
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
2 <description>to find differentially methylated positions</description>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
3 <macros>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
4 <import>macros.xml</import>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
5 </macros>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
6 <expand macro="requirements" />
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
7 <command detect_errors="exit_code"><![CDATA[
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
8 Rscript '$minfi_dmp_script'
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
9 ]]></command>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
10 <configfiles>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
11 <configfile name="minfi_dmp_script"><![CDATA[
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
12 require("minfi", quietly = TRUE)
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
13
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
14 set <- get(load('$grset'))
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
15
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
16 genomeranges <- as.data.frame(ranges(set))
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
17
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
18 beta <- getBeta(set)
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
19
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
20 pheno <- read.table('$phenotype_table',skip=1)
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
21
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
22 type <- '$phenotype'
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
23
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
24 qCutoff <- as.numeric('$q_cutoff')
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
25
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
26 shrinkVar <- '$variance_shrinkage'
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
27
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
28 tab <- read.table('$ucsc_genome')
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
29
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
30 tab <- tab[,-(11:14),drop=FALSE]
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
31
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
32 colnames(tab) <- c("seqname","source","feature","start","end","score","strand", "frame","attributes", "names")
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
33
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
34 tab[,"source"] <- NULL
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
35
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
36 tab[,"frame"] <- NULL
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
37
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
38 tab[,"attributes"] <- NULL
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
39
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
40 dmp <- dmpFinder(beta, pheno[,"V2"], type = type, qCutoff = qCutoff, shrinkVar = shrinkVar)
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
41
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
42 dmp[,"names"] <- rownames(dmp)
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
43
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
44 data <- merge(dmp, tab, by="names",sort = TRUE)
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
45
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
46 data <- data[c("seqname","start","end","names","score","strand", "feature","intercept", "f", "pval","qval")]
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
47
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
48 write.table(data, file= '$dmp', quote = FALSE,col.names = FALSE, row.names = FALSE, sep = "\t")
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
49 ]]>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
50 </configfile>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
51 </configfiles>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
52 <inputs>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
53 <param type="data" name="grset" format="rdata" label="Input Set" help="MethylSet, RatioSet or GenomicRatioSet" />
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
54 <param type="data" name="phenotype_table" format="tabular" label="Phenotype Table"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
55 <param name="phenotype" type="select" label="Phenotype Type">
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
56 <option value="categorical">categorical</option>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
57 <option value="continuous">continuous</option>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
58 </param>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
59 <param name="q_cutoff" type="integer" value="1" label="qCutoff Size" help="DMPs with an FDR q-value greater than this will not be returned."/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
60 <param name="variance_shrinkage" type="boolean" truevalue="TRUE" falsevalue="FALSE" label="Variance Shrinkage"
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
61 help="Enable variance shrinkage is recommended when sample sizes are small"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
62 <param type="data" name="ucsc_genome" format="gtf" label="Genome Table" help="UCSC genome data Methyl450"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
63 </inputs>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
64 <outputs>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
65 <data name="dmp" format="interval" label="Differentially Methylated Positions"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
66 </outputs>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
67 <tests>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
68 <test>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
69 <param name="grset" value="GRSet_without_SNPs.rdata"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
70 <param name="phenotype_table" value="phenotypeTable.txt"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
71 <param name="phenotype" value="categorical"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
72 <param name="q_cutoff" value="1"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
73 <param name="variance_shrinkage" value="FALSE"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
74 <param name="ucsc_genome" value="ucsc.gtf"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
75 <output name="dmp" file="Differentially_Methylated_Positions.interval"/>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
76 </test>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
77 </tests>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
78 <help><![CDATA[
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
79 This tool creates differentially methylated positions (DMPs) with respect to a phenotype covariate. The phenotype may be categorical (e.g. cancer vs. normal) or continuous (e.g. blood pressure).
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
80
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
81 The output is an interval file with the following columns:
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
82
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
83 - seqnames
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
84 - start
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
85 - end
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
86 - cpg
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
87 - names
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
88 - score
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
89 - strand
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
90 - feature
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
91 - intercept
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
92 - f pval
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
93 - qval
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
94
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
95 ]]></help>
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
96 <expand macro="citations" />
cffa2871fe8b Uploaded
kpbioteam
parents:
diff changeset
97 </tool>