Mercurial > repos > testtool > find_snp

diff findSNP/findSNPs.R @ 0:803decff7ac8 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 10:00:48 -0400
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/findSNP/findSNPs.R	Sun Jun 11 10:00:48 2017 -0400
@@ -0,0 +1,17 @@
+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))
+
+snps <- getSnpInfo(GRset)
+
+write.table(snps, output)
+
+