Mercurial > repos > testtool > block_finder

view findBlocks/findBlocks.R @ 0:fea57210de6c draft

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

require("IlluminaHumanMethylation450kanno.ilmn12.hg19", quietly = TRUE)
require("minfi", quietly = TRUE)

options(warn = -1)
options("download.file.method"="wget")

args <- commandArgs(trailingOnly = TRUE)

input = args[1]
output = args[2]

GRset <- get(load(input))

blocks <- blockFinder(GRset, design, coef = 2, what = c("Beta", "M"),
            cluster = NULL, cutoff = NULL,
            pickCutoff = FALSE, pickCutoffQ = 0.99,
            nullMethod = c("permutation","bootstrap"),
            smooth = TRUE, smoothFunction = locfitByCluster,
            B = ncol(permutations), permutations = NULL,
            verbose = TRUE)

write.table(blocks, output)