annotate profile_cl/getOntology.R @ 0:8547aedf1350 draft

Uploaded
author immport-devteam
date Mon, 27 Feb 2017 13:03:19 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
1 # Cell Ontology Module for Galaxy
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
2 # FlowCL
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
3 ######################################################################
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
4 # Copyright (c) 2016 Northrop Grumman.
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
5 # All rights reserved.
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
6 ######################################################################
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
7 #
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
8 # Version 1
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
9 # Cristel Thomas
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
10 #
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
11 #
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
12
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
13 library(flowCL)
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
14 library(base)
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
15
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
16 getOntology <- function(output_file, markers) {
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
17 res <- flowCL(markers, ResetArch = TRUE)
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
18 if (length(res) == 6) {
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
19 report <- capture.output(res$Table)
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
20 sink(output_file)
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
21 cat(report, sep = "\n")
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
22 sink()
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
23 }
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
24 }
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
25
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
26 args <- commandArgs(trailingOnly = TRUE)
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
27 markers <- paste(args[3:length(args)], collapse="")
8547aedf1350 Uploaded
immport-devteam
parents:
diff changeset
28 getOntology(args[2], markers)