Speed up rds creation

jenzopr · Oct 2, 2018 · 5f51977 · 5f51977
1 parent b67a62a
commit 5f51977
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/createSCE_kallisto.R b/src/createSCE_kallisto.R
@@ -41,7 +41,7 @@ expression <- list(umi = umi_counts)
 #
 if (snakemake@config$container$annotate) {
   ec_members %>%
-    lapply(function(m) {unique(tx2gene[m, "gene"])}) -> ec2gene
+    lapply(function(m) {unique(as.character(tx2gene[m, "gene"]))}) -> ec2gene
 
   bm <- biomaRt::getBM(c(snakemake@config$container$biomart$filter, 'external_gene_name'),
                        filters = snakemake@config$container$biomart$filter,
@@ -51,7 +51,7 @@ if (snakemake@config$container$annotate) {
   gene2symbol <- bm[, 'external_gene_name']
   names(gene2symbol) <- bm[, snakemake@config$container$biomart$filter]
 
-  rowData_ <- data.frame(transcript_ids = sapply(ec_members, function(m) paste0(tx2gene[m, "transcript"], collapse = ",")),
+  rowData_ <- data.frame(transcript_ids = sapply(ec_members, function(m) paste0(as.character(tx2gene[m, "transcript"]), collapse = ",")),
                          gene_ids = sapply(ec2gene, function(m) paste0(m, collapse = ",")),
                          symbols = sapply(ec2gene, function(m) { paste0(unname(gene2symbol[m]), collapse = ",") }))
 } else {