Mercurial > repos > kpbioteam > minfi_qc

annotate minfi_qc.R @ 0:7bd0d3c7a8d3 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload
author kpbioteam
date Tue, 06 Feb 2018 05:40:33 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
1 require("minfi", quietly = TRUE)
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
2
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
3 args <- commandArgs(trailingOnly = TRUE)
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
4
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
5 input = args[1]
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
6 output1 = args[2]
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
7 output2 = args[3]
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
8
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
9 MSet <- get(load(input))
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
10
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
11 qc <- getQC(MSet)
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
12
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
13 write.table(qc, output1)
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
14
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
15 png(output2)
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
16 plotQC(qc)
7bd0d3c7a8d3 planemo upload
kpbioteam
parents:
diff changeset
17 dev.off()