Mercurial > repos > testtool > find_dmr

view findDMR/findDMR.R @ 9:bae8d90853b7 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author testtool
date Tue, 25 Apr 2017 13:35:53 -0400
parents a0e3f4efb03f
children 7df2b7d79391
line wrap: on
line source

require("IlluminaHumanMethylation450kanno.ilmn12.hg19", quietly = TRUE)
require("data.table", quietly = TRUE)
require("minfi", 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 <- fread(input2)

designMatrix <- model.matrix(~ pheno$Phenotype)

dmrs <- bumphunter(GRset, design = designMatrix,
                   cutoff = 0.2, B=0, type="Beta")
write.table(dmrs, output)