Mercurial > repos > gaelcge > r_signac_galaxy
comparison signac-find_clusters.R @ 0:6e0b320d8b6a draft default tip
"planemo upload commit dc808171975d0012e25bd7b32adc7a5a5c56a145-dirty"
| author | gaelcge |
|---|---|
| date | Tue, 02 Aug 2022 19:11:27 +0000 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:6e0b320d8b6a |
|---|---|
| 1 #!/usr/bin/env Rscript | |
| 2 | |
| 3 # Load optparse we need to check inputs | |
| 4 | |
| 5 suppressPackageStartupMessages(require(optparse)) | |
| 6 | |
| 7 # Load common functions | |
| 8 | |
| 9 suppressPackageStartupMessages(require(workflowscriptscommon)) | |
| 10 | |
| 11 # parse options | |
| 12 | |
| 13 option_list = list( | |
| 14 make_option( | |
| 15 c("--signac-object"), | |
| 16 action = "store", | |
| 17 default = NA, | |
| 18 type = 'character', | |
| 19 help = "" | |
| 20 ), | |
| 21 make_option( | |
| 22 c("--algorithm"), | |
| 23 action = "store", | |
| 24 default = NA, | |
| 25 type = 'character', | |
| 26 help = "." | |
| 27 ), | |
| 28 make_option( | |
| 29 c("--output-object-file"), | |
| 30 action = "store", | |
| 31 default = NA, | |
| 32 type = 'character', | |
| 33 help = "File name in which to store serialized R matrix object." | |
| 34 ) | |
| 35 ) | |
| 36 | |
| 37 opt <- wsc_parse_args(option_list) | |
| 38 | |
| 39 suppressPackageStartupMessages(require(Seurat)) | |
| 40 suppressPackageStartupMessages(require(Signac)) | |
| 41 | |
| 42 set.seed(1234) | |
| 43 | |
| 44 # extract gene annotations from EnsDb | |
| 45 signac_object <- readRDS(file = opt$signac_object) | |
| 46 | |
| 47 signac_object <- FindClusters(object = signac_object, verbose = FALSE, algorithm = as.numeric(opt$algorithm)) | |
| 48 | |
| 49 # Output to a serialized R object | |
| 50 saveRDS(signac_object, file = opt$output_object_file) |