Repository 'find_dmr'
hg clone https://toolshed.g2.bx.psu.edu/repos/testtool/find_dmr

Changeset 12:083895dbc289 (2017-06-12)
Previous changeset 11:24ac6f93cc3e (2017-05-11) Next changeset 13:3735f4bd0a2d (2017-06-12)
Commit message:
Uploaded
modified:
findDMR/.Rhistory
findDMR/findDMR.R
findDMR/findDMR.xml
removed:
findDMR/test-data/._input.csv
findDMR/test-data/DMR.bed
findDMR/test-data/MetaTable.csv
findDMR/test-data/input.csv
findDMR/tool_dependencies.xml
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/.Rhistory
--- a/findDMR/.Rhistory Thu May 11 11:26:02 2017 -0400
+++ b/findDMR/.Rhistory Mon Jun 12 11:05:04 2017 -0400
[
@@ -1,200 +1,1 @@
-clusterSize=2
-class(clusterSize)
-type(clusterSize)
-require("minfi", quietly = TRUE)
-??getGenomicRatioSetFromGEO
-getGenomicRatioSetFromGEO
-getwd()
-GEODataTable <- getGEO(GSE51547)
-require("BiocGenerics", quietly = TRUE)
-require("data.table", quietly = TRUE)
-require("GEOquery", quietly = TRUE)
-require("rtracklayer", quietly = TRUE)
-require("FDb.InfiniumMethylation.hg19", quietly = TRUE)
-GEODataTable <- getGEO(GSE51547)
-GEODataTable <- getGEO("GSE51547")
-IlmnIDTable <- Table(GEODataTable)
-MetaData <- data.frame(Meta(GEODataTable))
-IlmnIDTable <- Table(GEODataTable)
-GEODataTable
-GSE<-"GSE51547"
-function (GSE = NULL, path = NULL, array = "IlluminaHumanMethylation450k",
-annotation = .default.450k.annotation, what = c("Beta", "M"),
-mergeManifest = FALSE, i = 1)
-{
-what <- match.arg(what)
-if (is.null(GSE) && is.null(path))
-stop("Either GSE or path must be supplied.")
-if (!is.null(GSE))
-gset <- GEOquery::getGEO(GSE)
-else gset <- GEOquery::getGEO(filename = file.path(path,
-list.files(path, pattern = ".soft")))
-if (length(gset) == 0)
-stop("Empty list retrieved from GEO.")
-if (length(gset) > 1) {
-warning("More than one ExpressionSet found:\n", names(gset),
-"\nUsing entry ", i)
-gset <- gset[[i]]
-}
-else gset <- gset[[1]]
-platform <- annotation(gset)
-if (platform != "GPL13534")
-warning(sprintf("%s is not the platform ID associated with IlluminaHumanMethylation450k. Should be GPL13534.",
-platform))
-if (what == "Beta" & (min(exprs(gset)[, 1], na.rm = TRUE) <
-0 | max(exprs(gset)[, 1], na.rm = TRUE) > 1))
-warning("Values outside [0,1] detected. 'what' argument should not be Beta.")
-ann <- .getAnnotationString(c(array = array, annotation = annotation))
-if (!require(ann, character.only = TRUE))
-stop(sprintf("cannot load annotation package %s", ann))
-object <- get(ann)
-gr <- getLocations(object, mergeManifest = mergeManifest,
-orderByLocation = TRUE)
-locusNames <- names(gr)
-sampleNames(gset) <- gset$title
-common <- intersect(locusNames, fData(gset)$Name)
-if (length(common) == 0)
-stop("No rowname matches. 'rownames' need to match IlluminaHumanMethylation450k probe names.")
-ind1 <- match(common, fData(gset)$Name)
-ind2 <- match(common, locusNames)
-preprocessing <- c(rg.norm = paste0("See GEO ", GSE, " for details"))
-if (what == "Beta") {
-out <- GenomicRatioSet(gr = gr[ind2, ], Beta = exprs(gset)[ind1,
-, drop = FALSE], M = NULL, CN = NULL, pData = pData(gset),
-annotation = c(array = array, annotation = annotation),
-preprocessMethod = preprocessing)
-}
-else {
-out <- GenomicRatioSet(gr = gr[ind2, ], Beta = NULL,
-M = exprs(gset)[ind1, , drop = FALSE], CN = NULL,
-pData = pData(gset), annotation = c(array = array,
-annotation = annotation), preprocessMethod = preprocessing)
-}
-return(out)
-}
-<environment: namespace:minfi>
-function (GSE = NULL, path = NULL, array = "IlluminaHumanMethylation450k",
-annotation = .default.450k.annotation, what = c("Beta", "M"),
-mergeManifest = FALSE, i = 1)
-{
-what <- match.arg(what)
-if (is.null(GSE) && is.null(path))
-stop("Either GSE or path must be supplied.")
-if (!is.null(GSE))
-gset <- GEOquery::getGEO(GSE)
-else gset <- GEOquery::getGEO(filename = file.path(path,
-list.files(path, pattern = ".soft")))
-if (length(gset) == 0)
-stop("Empty list retrieved from GEO.")
-if (length(gset) > 1) {
-warning("More than one ExpressionSet found:\n", names(gset),
-"\nUsing entry ", i)
-gset <- gset[[i]]
-}
-else gset <- gset[[1]]
-platform <- annotation(gset)
-if (platform != "GPL13534")
-warning(sprintf("%s is not the platform ID associated with IlluminaHumanMethylation450k. Should be GPL13534.",
-platform))
-if (what == "Beta" & (min(exprs(gset)[, 1], na.rm = TRUE) <
-0 | max(exprs(gset)[, 1], na.rm = TRUE) > 1))
-warning("Values outside [0,1] detected. 'what' argument should not be Beta.")
-ann <- .getAnnotationString(c(array = array, annotation = annotation))
-if (!require(ann, character.only = TRUE))
-stop(sprintf("cannot load annotation package %s", ann))
-object <- get(ann)
-gr <- getLocations(object, mergeManifest = mergeManifest,
-orderByLocation = TRUE)
-locusNames <- names(gr)
-sampleNames(gset) <- gset$title
-common <- intersect(locusNames, fData(gset)$Name)
-if (length(common) == 0)
-stop("No rowname matches. 'rownames' need to match IlluminaHumanMethylation450k probe names.")
-ind1 <- match(common, fData(gset)$Name)
-ind2 <- match(common, locusNames)
-preprocessing <- c(rg.norm = paste0("See GEO ", GSE, " for details"))
-if (what == "Beta") {
-out <- GenomicRatioSet(gr = gr[ind2, ], Beta = exprs(gset)[ind1,
-, drop = FALSE], M = NULL, CN = NULL, pData = pData(gset),
-annotation = c(array = array, annotation = annotation),
-preprocessMethod = preprocessing)
-}
-else {
-out <- GenomicRatioSet(gr = gr[ind2, ], Beta = NULL,
-M = exprs(gset)[ind1, , drop = FALSE], CN = NULL,
-pData = pData(gset), annotation = c(array = array,
-annotation = annotation), preprocessMethod = preprocessing)
-}
-return(out)
-}
-GSE<-"GSE51547"
-function (GSE = NULL, path = NULL, array = "IlluminaHumanMethylation450k",
-annotation = .default.450k.annotation, what = c("Beta", "M"),
-mergeManifest = FALSE, i = 1)
-{
-what <- match.arg(what)
-if (is.null(GSE) && is.null(path))
-stop("Either GSE or path must be supplied.")
-if (!is.null(GSE))
-gset <- GEOquery::getGEO(GSE)
-else gset <- GEOquery::getGEO(filename = file.path(path,
-list.files(path, pattern = ".soft")))
-if (length(gset) == 0)
-stop("Empty list retrieved from GEO.")
-if (length(gset) > 1) {
-warning("More than one ExpressionSet found:\n", names(gset),
-"\nUsing entry ", i)
-gset <- gset[[i]]
-}
-else gset <- gset[[1]]
-platform <- annotation(gset)
-if (platform != "GPL13534")
-warning(sprintf("%s is not the platform ID associated with IlluminaHumanMethylation450k. Should be GPL13534.",
-platform))
-if (what == "Beta" & (min(exprs(gset)[, 1], na.rm = TRUE) <
-0 | max(exprs(gset)[, 1], na.rm = TRUE) > 1))
-warning("Values outside [0,1] detected. 'what' argument should not be Beta.")
-ann <- .getAnnotationString(c(array = array, annotation = annotation))
-if (!require(ann, character.only = TRUE))
-stop(sprintf("cannot load annotation package %s", ann))
-object <- get(ann)
-gr <- getLocations(object, mergeManifest = mergeManifest,
-orderByLocation = TRUE)
-locusNames <- names(gr)
-sampleNames(gset) <- gset$title
-common <- intersect(locusNames, fData(gset)$Name)
-if (length(common) == 0)
-stop("No rowname matches. 'rownames' need to match IlluminaHumanMethylation450k probe names.")
-ind1 <- match(common, fData(gset)$Name)
-ind2 <- match(common, locusNames)
-preprocessing <- c(rg.norm = paste0("See GEO ", GSE, " for details"))
-if (what == "Beta") {
-out <- GenomicRatioSet(gr = gr[ind2, ], Beta = exprs(gset)[ind1,
-, drop = FALSE], M = NULL, CN = NULL, pData = pData(gset),
-annotation = c(array = array, annotation = annotation),
-preprocessMethod = preprocessing)
-}
-else {
-out <- GenomicRatioSet(gr = gr[ind2, ], Beta = NULL,
-M = exprs(gset)[ind1, , drop = FALSE], CN = NULL,
-pData = pData(gset), annotation = c(array = array,
-annotation = annotation), preprocessMethod = preprocessing)
-}
-return(out)
-}
-out
-GSE<-"GSE51547"
-gset <- getGEO(GSE)
-ann <- .getAnnotationString(c(array = array, annotation = annotation))
-platform <- annotation(gset)
-hm450.hg19 <- getPlatform()
-IlmnInfo <-
-data.table(
-IlmnID = names(hm450.hg19),
-CHR = as.data.frame(hm450.hg19@seqnames)$value,
-BP = as.numeric(hm450.hg19@elementMetadata$probeStart)
-)
-IlmnIDTable <- Table(gset)
-??getGenomicRatioSetFromGEO()
-getGenomicRatioSetFromGEO()
-getGenomicRatioSetFromGEO
+??fread
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/findDMR.R
--- a/findDMR/findDMR.R Thu May 11 11:26:02 2017 -0400
+++ b/findDMR/findDMR.R Mon Jun 12 11:05:04 2017 -0400
[
@@ -1,42 +1,36 @@
-require("BiocGenerics", quietly = TRUE)
-require("data.table", quietly = TRUE)
-require("bumphunter", quietly = TRUE)
+require("minfi", quietly = TRUE)
+
+options(warn = -1)
+options("download.file.method"="wget")
 
 args <- commandArgs(trailingOnly = TRUE)
-GSMTable = args[1]
-platform = args[2]
-Data_Table = args[3]
-cutoff = as.numeric(args[4])
-clusterSize = as.numeric(args[5])
-DMR = args[6]
+input1 = args[1]
+input2 = args[2]
+cutoff = as.numeric(args[3])
+B = as.numeric(args[4]) 
+pickCutoffQ = as.numeric(args[5])
+output1 = args[6]
+output2 = args[7]
 
-TAB = fread(GSMTable)
-
-IlmnInfo = fread(platform)
+GRset <- get(load(input1))
 
-gmSet = fread(Data_Table)
+pheno <- fread(input2)
 
-# bumphunter Run with processed data
-designMatrix <- model.matrix( ~ TAB$Phenotype)
+designMatrix <- model.matrix(~ pheno$Phenotype)
 
-bumps <- bumphunter(
-  as.matrix(gmSet),
-  design = designMatrix,
-  pos = IlmnInfo$BP,
-  cutoff = cutoff,
-  chr = IlmnInfo$CHR
-)
+dmrs <- bumphunter(GRset, design = designMatrix,
+                   cutoff =cutoff, B=B, type="Beta",  pickCutoff=TRUE,
+                   pickCutoffQ=pickCutoffQ)
+
+DMRTable <- dmrs$table
+
+write.table(DMRTable, output1)
 
-# choose DMR's of a certain length threshold
-DMRTable <- bumps$table[which(bumps$table$L >= clusterSize), ]
-DMRInfo <- data.table(DMRTable$chr, DMRTable$start, DMRTable$end)
-
+sign=sign(DMRTable$value)
 
-write.table(
-  DMRInfo,
-  DMR,
-  quote = F,
-  sep = "\t",
-  row.names = F,
-  col.names = F
-)
+sign[sign==-1]="-"
+sign[sign==1]="+"
+
+dmr_track<-cbind(as.character(DMRTable$chr),DMRTable$start,DMRTable$end,sign)
+
+write.table(dmr_track,output2,quote = FALSE, sep = "\t\t\t\t",row.names = FALSE, col.names = FALSE,append=TRUE)
\ No newline at end of file
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/findDMR.xml
--- a/findDMR/findDMR.xml Thu May 11 11:26:02 2017 -0400
+++ b/findDMR/findDMR.xml Mon Jun 12 11:05:04 2017 -0400
[
@@ -1,57 +1,34 @@
-<tool id="DMR" name="findDMR" version="1.16.2">
-  <description>
-from series of samples
-  </description>
+<tool id="findDMR" name="findDMR" version="1.16.2">
   <requirements>
-    <requirement type= "package" version="0.0.1">450kdependency</requirement>
-     <requirement type="package" version="3.2.1">R</requirement>
-     </requirements>
+   <requirement name="package_r_3_2_1" type="package" version="3.2.1">R</requirement>
+</requirements>
 <stdio>
    <exit_code range="1:" />
 </stdio>
-  <command> Rscript $__tool_directory__/findDMR.R "$GSMTable" "$platform" "$Data_Table" "$cutoff" "$clusterSize" "$DMR"</command>
+  <command> Rscript $__tool_directory__/findDMR.R "$input1" "$input2" "$output" </command>
   <inputs>
-    <param optional="false" format="txt" name="GSMTable" type="data" value="" help="GSM ID and phenotype table." label="[required] GSMTable">
-      <validator type="empty_field" message="This field is required."/>
-    </param>
-     <param optional="false" format="txt" name="platform" type="data" value="" help="CG, CHR, BP table." label="[required] platform">
-      <validator type="empty_field" message="This field is required."/>
-    </param>
-     <param optional="false" format="txt" name="Data_Table" type="data" value="" help="Matrix of data." label="[required] Data_Table">
-      <validator type="empty_field" message="This field is required."/>
-    </param>
-    <param name="cutoff" type="float" value=""
-        label="Enter cutoff size (number)"
-        help="e.g. '0.2'" >
-      </param>
-    <param name="clusterSize" type="float" value=""
-        label="Enter cluster size (number)"
-        help="e.g. '2'" >
-      </param>
+    <param name="input1" format="Rdata" label="GenomicRatioSet.Rdata" help="e.g. 'GSE51547'"/>
+    <param name="input2" format="txt" label="PhenoTab.txt" help="ID and Phenotype table"/>
   </inputs>
   <outputs>
-    <data format="bed" name="DMR" label="DMR.bed"/>
+    <data format="txt" name="output" label="findDMR.txt"/>
   </outputs>
   <tests>
     <test>
       <param name="test">
-      <element name="test-data">
-          <collection type="data">
-              <element name="GSMTable" value="test-data/input.txt"/>
-              <element name="platform" value="test-data/platform.txt"/>
-              <element name="Data_Table" value="test-data/Data_Table.txt"/>
-              <element name="cutoff" value="0.2"/>
-              <element name="clusterSize" value="2"/>
+      <element name="test-data"/>
+           <collection type="data">
+                <element format="Rdata" name="input1" label="test-data/GenomicRatioSet.Rdata"/>
+                <element format="txt" name="input2" label="test-data/PhenoTab.txt"/>
           </collection>
-        </element>
         </param>
-        <output format="bed" name="DMR" label="test-data/DMR.bed"/>
+        <output format="txt" name="output" label="test-data/findDMR.txt"/>
         </test>
     </tests>
   <help>
-**Description**
+Identify CpGs where methylation is associated with a continuous or categorical phenotype.
 </help>
 <citations>
-DMR
+Aryee, Martin J., et al. "Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays." Bioinformatics 30.10 (2014): 1363-1369.
 </citations>
 </tool>
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/test-data/._input.csv
b
Binary file findDMR/test-data/._input.csv has changed
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/test-data/DMR.bed
--- a/findDMR/test-data/DMR.bed Thu May 11 11:26:02 2017 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
b'@@ -1,9881 +0,0 @@\n-chrX\t106959609\t106960186\n-chr1\t54203844\t54204501\n-chr5\t112073385\t112073454\n-chr18\t56338460\t56338659\n-chr6\t29943460\t29944045\n-chr11\t64405452\t64405993\n-chr4\t1363998\t1364344\n-chr16\t78079753\t78080193\n-chr3\t38070886\t38071301\n-chr2\t175351922\t175352168\n-chr4\t1806633\t1807332\n-chr6\t28793329\t28793608\n-chr4\t1205447\t1206004\n-chr3\t113775302\t113775388\n-chr16\t3085934\t3086940\n-chr2\t27485967\t27486122\n-chr1\t2918058\t2918526\n-chr17\t422325\t422568\n-chr5\t141031071\t141031094\n-chr6\t17987573\t17987659\n-chr17\t2952307\t2952767\n-chr19\t16653151\t16653581\n-chr3\t193921489\t193922037\n-chr2\t70418309\t70418827\n-chr16\t1095625\t1096754\n-chr11\t2161355\t2161411\n-chr12\t51664503\t51664607\n-chr16\t1705991\t1706153\n-chr14\t23067033\t23067188\n-chr3\t173113110\t173113519\n-chr7\t50444039\t50444317\n-chr3\t100053070\t100053140\n-chr17\t6918792\t6918952\n-chr2\t175350573\t175351649\n-chr6\t33216768\t33217142\n-chr16\t85646762\t85646832\n-chr16\t11439254\t11439785\n-chr2\t201981786\t201981873\n-chr10\t128993810\t128994432\n-chr7\t116593215\t116593288\n-chr16\t29873565\t29874506\n-chr22\t40742327\t40742374\n-chr21\t46031707\t46032038\n-chr6\t31940192\t31940302\n-chr2\t128784344\t128784923\n-chr19\t45927001\t45927018\n-chr16\t24550834\t24552079\n-chr6\t26284890\t26285184\n-chr22\t47158324\t47158521\n-chr3\t63850035\t63850278\n-chr16\t672235\t672882\n-chr14\t57274763\t57275234\n-chr16\t3766844\t3767679\n-chr3\t44036098\t44036453\n-chr15\t32907294\t32907402\n-chr17\t48637584\t48637666\n-chr5\t137800779\t137801033\n-chr22\t20004369\t20004611\n-chr19\t35531102\t35531591\n-chr18\t32397976\t32398246\n-chr22\t30162987\t30162991\n-chr17\t78032296\t78032492\n-chr10\t131927297\t131927449\n-chr19\t1947674\t1947912\n-chr1\t36106988\t36107155\n-chr9\t101568927\t101569648\n-chr18\t48346660\t48346795\n-chr1\t211652228\t211652741\n-chr4\t1294432\t1294566\n-chr1\t12079449\t12079505\n-chr14\t103740037\t103740476\n-chr12\t108237572\t108237937\n-chr13\t113097259\t113097633\n-chr18\t63417557\t63417735\n-chr3\t51976003\t51976171\n-chr12\t123237034\t123237126\n-chr9\t88555654\t88555863\n-chr3\t53289290\t53289728\n-chr16\t31214346\t31214378\n-chr12\t8185029\t8185344\n-chr10\t38146893\t38147096\n-chr11\t62341426\t62341521\n-chr14\t65239309\t65239653\n-chr20\t42219469\t42219655\n-chr12\t10364612\t10365082\n-chr3\t48487981\t48488082\n-chr7\t45614242\t45614720\n-chr2\t208632784\t208633283\n-chr9\t75136064\t75136701\n-chr12\t12966076\t12966271\n-chr1\t16563734\t16563906\n-chr16\t2255332\t2255883\n-chr9\t77642922\t77643323\n-chr14\t56047764\t56047916\n-chr19\t10138244\t10138615\n-chr17\t73901016\t73901339\n-chr16\t14397688\t14397766\n-chr16\t30457360\t30457524\n-chr1\t155243076\t155243300\n-chr16\t32936895\t32937009\n-chr7\t18067530\t18067535\n-chr11\t1304252\t1304348\n-chr16\t56677510\t56677585\n-chr20\t2821434\t2821485\n-chr1\t153914312\t153914453\n-chr16\t2569925\t2570212\n-chr15\t68870979\t68871357\n-chr17\t7517617\t7517649\n-chr16\t28608846\t28609531\n-chr7\t2563136\t2563622\n-chr1\t75139583\t75139958\n-chr7\t44104935\t44105118\n-chr16\t4852372\t4852777\n-chr6\t101328458\t101329304\n-chr20\t35504198\t35504371\n-chr6\t33656083\t33656500\n-chr12\t120739972\t120740007\n-chr1\t247579560\t247580517\n-chr8\t37824481\t37824496\n-chr12\t108079012\t108079503\n-chr19\t10735958\t10736001\n-chr7\t154706045\t154706916\n-chr16\t8962538\t8962970\n-chr10\t45914477\t45914640\n-chr16\t12211050\t12211322\n-chr16\t31366118\t31366418\n-chr14\t58893467\t58893740\n-chr18\t23806180\t23806263\n-chr7\t82541\t82911\n-chr15\t77287290\t77287608\n-chr15\t69706335\t69706469\n-chr3\t195622522\t195622598\n-chr12\t24102533\t24103021\n-chr10\t31607330\t31607651\n-chr3\t53164962\t53165031\n-chr16\t19566218\t19566823\n-chr19\t18451444\t18451525\n-chr7\t100493811\t100494332\n-chr16\t3550032\t3550792\n-chr16\t20359396\t20359923\n-chr2\t27994492\t27994636\n-chr16\t8715546\t8715882\n-chr16\t30379343\t30380183\n-chr16\t29296186\t29296693\n-chr15\t45879564\t45879610\n-chr11\t70289501\t70289510\n-chr20\t56935373\t56935483\n-chr10\t48355005\t48355070\n-chr14\t63568147\t63568558\n-chr1\t224518178\t224518192\n-chr10\t99447118\t99447212\n-chr2\t130763379\t130763641\n-chr1\t27901757\t27902687\n-chr16\t19124971\t19125137\n-chr17\t74236897\t74237210\n-chr1\t228782809\t228783335\n-chr15\t85144175\t85144257\n-chr5\t122181105\t122181680\n-chr2\t143886278\t143886375\n-chr20\t31349704\t31349891\n-chr19\t7964567\t7964703\n-chr16\t338531'..b'33247885\t33247984\n-chr2\t113398727\t113399144\n-chr6\t139695804\t139695808\n-chr6\t73935125\t73935200\n-chr1\t1228942\t1229105\n-chr5\t156773546\t156773616\n-chr3\t18480249\t18480341\n-chr2\t242566905\t242566940\n-chr2\t240239755\t240239767\n-chr17\t48858588\t48858603\n-chr16\t89342351\t89342629\n-chr13\t28544814\t28545214\n-chr1\t26233404\t26233435\n-chr13\t114062198\t114062507\n-chr1\t33235571\t33235717\n-chr3\t48594424\t48594454\n-chrX\t42637449\t42637506\n-chr5\t131347720\t131347725\n-chr15\t93257477\t93257579\n-chr1\t9380708\t9380800\n-chr12\t110888158\t110888232\n-chr13\t53776211\t53776396\n-chr9\t139840384\t139840630\n-chr3\t96495326\t96495603\n-chr4\t7287567\t7287737\n-chr11\t65819726\t65819742\n-chr20\t25371652\t25371671\n-chr19\t5567758\t5568093\n-chr6\t33272178\t33272186\n-chr7\t156801539\t156801818\n-chr17\t48698713\t48698751\n-chr12\t131258228\t131258361\n-chr9\t5831674\t5831951\n-chrX\t154842736\t154842793\n-chr15\t43513504\t43513563\n-chr8\t30890469\t30890553\n-chr8\t145019491\t145019563\n-chr11\t65196125\t65196174\n-chr21\t46019376\t46019570\n-chr22\t41777755\t41777779\n-chr1\t43147939\t43148010\n-chr2\t178257542\t178257632\n-chr11\t1446386\t1446802\n-chr6\t32820153\t32820214\n-chr19\t15343347\t15343377\n-chr9\t4679736\t4679986\n-chr10\t90343387\t90343641\n-chr8\t54852115\t54852232\n-chr9\t95947098\t95947348\n-chr16\t88761048\t88761284\n-chr20\t20257813\t20258001\n-chr13\t30948114\t30948122\n-chr7\t95026538\t95026672\n-chr7\t44530455\t44530521\n-chr5\t2752553\t2752834\n-chr16\t30798770\t30798776\n-chr19\t46713398\t46713416\n-chr4\t1352452\t1352787\n-chr8\t61566034\t61566391\n-chr6\t170065843\t170065938\n-chr16\t21294899\t21295207\n-chr20\t23028842\t23028977\n-chr6\t133562475\t133562479\n-chr2\t62732892\t62733114\n-chr11\t7951168\t7951173\n-chr17\t79792936\t79792949\n-chr13\t96705664\t96705738\n-chr5\t132386262\t132386298\n-chr12\t7126208\t7126344\n-chr3\t119421553\t119421604\n-chr8\t72754421\t72754905\n-chr6\t30134908\t30135065\n-chr22\t20130575\t20130698\n-chr1\t53542934\t53543090\n-chr8\t1765679\t1765772\n-chr5\t180615482\t180615543\n-chr7\t158857749\t158858057\n-chr16\t89557908\t89557923\n-chr10\t11783974\t11784058\n-chr19\t10444117\t10444479\n-chr11\t102669042\t102669291\n-chr4\t1604438\t1604576\n-chr10\t50747266\t50747328\n-chr6\t33265322\t33265386\n-chr3\t49142215\t49142226\n-chr8\t143193033\t143193131\n-chr16\t23193800\t23193885\n-chr17\t40332910\t40332961\n-chr3\t54157647\t54157677\n-chr6\t30881579\t30881595\n-chr6\t27100067\t27100370\n-chr16\t31383055\t31383088\n-chr5\t175223598\t175223982\n-chr11\t102745886\t102746236\n-chr17\t56296115\t56296444\n-chr16\t4511869\t4511939\n-chr4\t103790056\t103790299\n-chr15\t31284413\t31284569\n-chr11\t104034282\t104034571\n-chr20\t18569003\t18569196\n-chr19\t22019044\t22019059\n-chr8\t61429582\t61429779\n-chr11\t66335975\t66336077\n-chr10\t133758671\t133758839\n-chr6\t31856959\t31857034\n-chr3\t51428372\t51428392\n-chr14\t104647027\t104647046\n-chr20\t60737957\t60738265\n-chr2\t119606698\t119607049\n-chr11\t2813073\t2813298\n-chr11\t2020129\t2020231\n-chr14\t104175847\t104176333\n-chr11\t1785547\t1785552\n-chr1\t161009415\t161009644\n-chr7\t39170539\t39170715\n-chr17\t45855505\t45855691\n-chr20\t21686205\t21686264\n-chr5\t141693817\t141694185\n-chr17\t840927\t841205\n-chr2\t3750517\t3750952\n-chrX\t118708041\t118708065\n-chr17\t43238138\t43238230\n-chr6\t169226282\t169226337\n-chr3\t160283421\t160283489\n-chr7\t1959941\t1960073\n-chr1\t231762664\t231762689\n-chr10\t7860354\t7860463\n-chr7\t157367159\t157367462\n-chr17\t76498943\t76498964\n-chr14\t102027124\t102027514\n-chr17\t26941944\t26942145\n-chr8\t96036568\t96036755\n-chr4\t15964902\t15965016\n-chrX\t48329927\t48329937\n-chr11\t10472101\t10472108\n-chr4\t39460330\t39460342\n-chr6\t32549313\t32549496\n-chr7\t104654852\t104654870\n-chr6\t159240766\t159240774\n-chr19\t1079986\t1080466\n-chr1\t39339553\t39339683\n-chr4\t74486409\t74486510\n-chr6\t29999567\t29999708\n-chr1\t84543108\t84543369\n-chr1\t8877774\t8877813\n-chr12\t39298591\t39298776\n-chr6\t42946146\t42946178\n-chr14\t77279079\t77279247\n-chr1\t2565265\t2565729\n-chr17\t72368116\t72368179\n-chr5\t138612986\t138613030\n-chr6\t26502058\t26502192\n-chr7\t2598460\t2598473\n-chr19\t8649796\t8649942\n-chr3\t127006239\t127006414\n-chr5\t132072786\t132072925\n-chr4\t870161\t870443\n-chr11\t1365592\t1365668\n-chr16\t1575153\t1575219\n-chr1\t76189753\t76189897\n-chr2\t130691926\t130691932\n-chr11\t1413842\t1413990\n-chr7\t1881390\t1881839\n'
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/test-data/MetaTable.csv
--- a/findDMR/test-data/MetaTable.csv Thu May 11 11:26:02 2017 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
@@ -1,2 +0,0 @@
-"channel_count","characteristics_ch1","contact_address","contact_city","contact_country","contact_department","contact_institute","contact_name","contact_zip.postal_code","data_processing","data_row_count","description","extract_protocol_ch1","geo_accession","hyb_protocol","label_ch1","label_protocol_ch1","last_update_date","molecule_ch1","organism_ch1","platform_id","scan_protocol","series_id","source_name_ch1","status","submission_date","supplementary_file","taxid_ch1","title","type"
-"1","sample type: melanoma cell line","Barngatan 2B","Lund","Sweden","Dept of Oncology","Lund University","Martin,,Lauss","22185","Raw intensities for methylated (M) and unmethylated (U) signal were extracted from Illumina’s GenomeStudio. Beta-values were calculated as M/(M+U). A total of 496 missing values (melanomas) were imputed using k-nearest neighbor imputation (k=10). For each sample we performed a peak-based correction of Illumina I and II chemical assays. For both assays we smoothed the beta values (Epanechnikov smoothing kernel) to estimate unmethylated and methylated peaks, respectively; and the unmethylated peak was moved to 0 and the methylated peak to 1 using linear scaling, with beta-values in between stretched accordingly. Beta-values below 0 were set back to 0 and values above 1 were set to 1.","485577","melanoma cell line","Genomic DNA was extracted from the biopsies using QIAamp DNA Mini Kit (Qiagen). A total of 500 ng of DNA were used for bisulfite treatment, using the EZ DNA Methylation Kit (Zymo). We hybridized 200 ng in 4 μl to the Infinium Human Methylation450 BeadChip array.","GSM1247787","Bisulphite converted DNA was amplified, fragmented and hybridised to Illumina Infinium Human Methylation450 Beadchip using standard Illumina protocol","Cy5 and Cy3","Standard Illumina Protocol","May 17 2015","genomic DNA","Homo sapiens","GPL13534","Arrays were imaged using BeadArray Reader using standard recommended Illumina scanner setting","GSE51547","SKMEL3","Public on May 17 2015","Oct 22 2013","NONE","9606","genomic DNA from Sample SKMEL3","genomic"
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/test-data/input.csv
--- a/findDMR/test-data/input.csv Thu May 11 11:26:02 2017 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
@@ -1,4 +0,0 @@
-ID,Phenotype
-GSM1247787,melanoma
-GSM1247784,melanoma
-GSM1247733,healthy
b
diff -r 24ac6f93cc3e -r 083895dbc289 findDMR/tool_dependencies.xml
--- a/findDMR/tool_dependencies.xml Thu May 11 11:26:02 2017 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
@@ -1,10 +0,0 @@
-<?xml version="1.0"?>
-<tool_dependency>
-    <package name="450kdependency" version="0.0.1">
-     # TODO:
-        <repository changeset_revision="6ab5f9bf0d30" name="450kdependency" owner="testtool" toolshed="https://toolshed.g2.bx.psu.edu" />
-    </package>
-    <package name="R" version="3.2.1">
-        <repository changeset_revision="d0bf97420fb5" name="package_r_3_2_1" owner="iuc" toolshed="https://toolshed.g2.bx.psu.edu" />
-    </package>
-</tool_dependency>