diff findDMP/findDMP.R @ 0:8cf72c018b42 draft default tip

Uploaded

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/findDMP/findDMP.R	Sun Jun 11 17:29:26 2017 -0400
@@ -0,0 +1,23 @@
+require("minfi", quietly = TRUE)
+require("data.table", quietly = TRUE)
+
+options(warn = -1)
+options("download.file.method"="wget")
+
+args <- commandArgs(trailingOnly = TRUE)
+
+input1 = args[1]
+input2 = args[2]
+output = args[3]
+
+GRset <- get(load(input1))
+
+pheno <- read.csv(input2)
+
+beta <- getBeta(GRset)
+
+pData  <- pData(GRset)
+
+dmp <- dmpFinder(beta, pheno$Phenotype, type = "continuous")
+
+write.table(dmp, output)