view findBlocks/findBlocks.R @ 1:7ce7fb0f8d05 draft default tip

Uploaded
author testtool
date Sun, 11 Jun 2017 08:59:12 -0400
parents fea57210de6c
children
line wrap: on
line source

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)