Mercurial > repos > kpbioteam > ewastools
annotate minfi_dmr.xml @ 75:9c6fbb7d5a2a draft
planemo upload for repository https://github.com/kpbioteam/ewas_galaxy commit 9ce5c86f66f4a0cb64a1b8f48a2a937268b62064-dirty
| author | kpbioteam |
|---|---|
| date | Mon, 20 May 2019 07:14:26 -0400 |
| parents | 1334312dc198 |
| children |
| rev | line source |
|---|---|
| 16 | 1 <tool id="minfi_dmr" name="Minfi DMR" version="@MINFI_VERSION@"> |
| 2 <description>to find differentially methylated regions</description> | |
| 3 <macros> | |
| 4 <import>macros.xml</import> | |
| 5 </macros> | |
| 6 <expand macro="requirements" /> | |
| 7 <command detect_errors="exit_code"> | |
| 8 <![CDATA[ | |
| 9 Rscript '$minfi_dmr_script' | |
| 10 ]]> | |
| 11 </command> | |
| 12 <configfiles> | |
| 13 <configfile name="minfi_dmr_script"> | |
| 14 <![CDATA[ | |
| 15 require("minfi", quietly = TRUE) | |
| 16 options(warn = -1) | |
| 17 | |
| 18 GRSet <- get(load('$grset')) | |
| 19 | |
| 20 pheno <- read.table('$phenotype_table',skip = 1) | |
| 21 | |
| 22 group <- pheno\$V2 | |
| 23 | |
| 24 pair <- factor(pheno\$V3) | |
| 25 | |
| 26 design.matrix <- model.matrix(~ group + pair) | |
| 27 | |
| 28 maxGap <- as.numeric('$maxgap_size') | |
| 29 | |
| 30 if(is.null(GRSet\$cluster)){ | |
| 31 cluster = NULL | |
| 32 maxGap = maxGap | |
| 33 } else { | |
| 34 cluster = GRSet\$cluster | |
| 35 maxGap = NULL | |
| 36 } | |
| 37 | |
| 38 cutoff <- as.numeric('$cutoff_size') | |
| 39 B <- as.numeric('$number_of_resamples') | |
| 40 nullMethod <- '$null_method' | |
| 41 coef <- 2 | |
| 42 | |
| 43 dmrs <- bumphunter(GRSet, | |
| 44 design = design.matrix, | |
| 45 cluster = cluster, | |
| 46 maxGap = maxGap, | |
| 47 cutoff = cutoff, | |
| 48 nullMethod = nullMethod, | |
| 49 B = B) | |
| 50 | |
| 63 | 51 dmrGR <- dmrs\$table[,c(1,2,3)] |
| 52 colnames(dmrGR) <- c("chr","start","end") | |
| 16 | 53 |
| 63 | 54 write.table(dmrGR, file= '$dmr', quote = FALSE,col.names = TRUE, row.names = FALSE, sep = "\t") |
| 16 | 55 ]]> |
| 56 </configfile> | |
| 57 </configfiles> | |
| 58 <inputs> | |
|
75
9c6fbb7d5a2a
planemo upload for repository https://github.com/kpbioteam/ewas_galaxy commit 9ce5c86f66f4a0cb64a1b8f48a2a937268b62064-dirty
kpbioteam
parents:
71
diff
changeset
|
59 <param type="data" name="grset" format="rdata" label="Data Mapped To The Genome" |
|
9c6fbb7d5a2a
planemo upload for repository https://github.com/kpbioteam/ewas_galaxy commit 9ce5c86f66f4a0cb64a1b8f48a2a937268b62064-dirty
kpbioteam
parents:
71
diff
changeset
|
60 help="This class (GenomicRatioSet) holds preprocessed data for Illumina methylation microarrays, mapped to a genomic location." /> |
| 16 | 61 <param type="data" name="phenotype_table" format="tabular" label="Phenotype Table" |
| 62 help="Phenotype Table must include the following information: sampleID, phenotype and paird or unpaired samples column"/> | |
| 63 <param name="maxgap_size" type="integer" value="250" label="maxGap Size" | |
| 64 help="If cluster is not provided this maximum location gap will be used to define cluster"/> | |
| 65 <param name="cutoff_size" type="float" value="0.1" label="Cutoff Size" | |
| 66 help="A numeric value. Values of the estimate of the genomic profile above the cutoff or below the negative of the cutoff will be used as candidate regions. It is possible to give two separate values (upper and lower bounds). If one value is given, the lower bound is minus the value."/> | |
| 67 <param name="number_of_resamples" type="integer" value="0" label="Number of Resamples" | |
| 68 help="An integer denoting the number of resamples to use when computing null distributions. This defaults to 0. If permutations is supplied that defines the number of permutations/bootstraps and B is ignored."/> | |
| 69 <param name="null_method" type="select" label="null Method" help="Method used to generate null candidate regions."> | |
| 70 <option value="permutation" selected="True">permutation</option> | |
| 71 <option value="bootstrap">bootstrap</option> | |
| 72 </param> | |
| 73 </inputs> | |
| 74 <outputs> | |
| 75 <data name="dmr" format="interval" label="Differentially Methylated Regions"/> | |
| 76 </outputs> | |
| 77 <tests> | |
| 78 <test> | |
| 79 <param name="grset" value="GRSet_without_SNPs.rdata"/> | |
| 80 <param name="phenotype_table" value="phenotypeTable.txt"/> | |
| 81 <param name="maxgap_size" value="250"/> | |
| 82 <param name="cutoff_size" value="0.1"/> | |
| 83 <param name="number_of_resamples" value="0"/> | |
| 84 <param name="null_method" value="permutation"/> | |
| 85 <output name="dmr" file="Differentially_Methylated_Regions.interval"/> | |
| 86 </test> | |
| 87 </tests> | |
| 88 <help><![CDATA[ | |
| 89 Genomic Ratio Set should contain 450k arrays reads for each sample listed in Phenotype Table: sampleID phenotype paired e.g. GSM1588708 sensitive 1 | |
| 90 | |
| 91 This tool is a version of the bump hunting algorithm (Andrew E Jaffe et al. 2012) adapted to the 450k array finding find differentially methylated regions DMR. | |
| 92 | |
| 93 The output is an interval file with the following columns: | |
| 94 | |
| 63 | 95 - chr |
| 16 | 96 - start |
| 97 - end | |
| 98 | |
| 99 ]]></help> | |
| 100 <expand macro="citations" /> | |
| 101 </tool> |