diff --git a/scripts/counts_matrix_from_SALMON.R b/scripts/counts_matrix_from_SALMON.R index 46dadd3..f9d9f1d 100644 --- a/scripts/counts_matrix_from_SALMON.R +++ b/scripts/counts_matrix_from_SALMON.R @@ -37,7 +37,11 @@ writeCounts <- function(colDataFile, salmon_output_folder, counts_dir, type) { names(files) <- rownames(colData) txi <- tximport::tximport(files, type = "salmon", txOut = TRUE) - dds <- DESeq2::DESeqDataSetFromTximport(txi, colData, ~group) + if(length(unique(colData$sample_type)) < 2) { + dds <- DESeq2::DESeqDataSetFromTximport(txi, colData, ~1) + } else { + dds <- DESeq2::DESeqDataSetFromTximport(txi, colData, ~group) + } #save raw counts write.table(x = DESeq2::counts(dds),