diff --git a/dge_workflow/dge_utils.sh b/dge_workflow/dge_utils.sh index 38e4285a28276ad20af4ea11dee1736942ab0c60..8d43570eb0a0fc56a4a9f1aa77087509b1b4279a 100755 --- a/dge_workflow/dge_utils.sh +++ b/dge_workflow/dge_utils.sh @@ -459,14 +459,14 @@ devtools::source_url("https://raw.githubusercontent.com/holgerbrandl/datautils/v exprCounts <- list.files(".", "ReadsPerGene.out.tab") %>% ldply(function(countFile){ read.delim(countFile, header=F) %>% select(V1, V2) %>% - set_names("gene_id", "num_alignments") %>% - filter(!str_detect(gene_id, "^N_")) %>% + set_names("ensembl_gene_id", "num_alignments") %>% + filter(!str_detect(ensembl_gene_id, "^N_")) %>% mutate(sample=trim_ext(countFile, ".ReadsPerGene.out.tab")) }, .progress="text") countMatrix <- spread(exprCounts, sample, num_alignments) -write.delim(countMatrix, "star_count_matrix.txt") +write.delim(countMatrix, "star_counts_matrix.txt") ' | R --vanilla -q }