Mercurial > repos > testtool > find_dmp

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

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author testtool
date Sun, 11 Jun 2017 17:29:26 -0400
parents
children
line wrap: on
line source

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)