Mercurial > repos > iuc > ruvseq

comparison get_deseq_dataset.R @ 4:3a083c78896e draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/ruvseq commit 30117fce22f3771c9c0c13e70c3ad14b694de6e2
author iuc
date Fri, 21 Apr 2023 14:09:17 +0000
parents fed9d0350d72
children
comparison
equal deleted inserted replaced
3:d1f7fa5bb3cb 4:3a083c78896e
12 gff_file <- NULL 12 gff_file <- NULL
13 tx2gene <- read.table(tx2gene, header = has_header) 13 tx2gene <- read.table(tx2gene, header = has_header)
14 } 14 }
15 } 15 }
16 16
17 if (!use_txi & has_header) { 17 if (!use_txi && has_header) {
18 countfiles <- lapply(as.character(sample_table$filename), read.delim, row.names = 1) 18 countfiles <- lapply(as.character(sample_table$filename), read.delim, row.names = 1)
19 tbl <- do.call("cbind", countfiles) 19 tbl <- do.call("cbind", countfiles)
20 colnames(tbl) <- rownames(sample_table) # take sample ids from header 20 colnames(tbl) <- rownames(sample_table) # take sample ids from header
21 21
22 # check for htseq report lines (from DESeqDataSetFromHTSeqCount function) 22 # check for htseq report lines (from DESeqDataSetFromHTSeqCount function)
33 dds <- DESeqDataSetFromMatrix( 33 dds <- DESeqDataSetFromMatrix(
34 countData = tbl, 34 countData = tbl,
35 colData = subset(sample_table, select = -filename), 35 colData = subset(sample_table, select = -filename),
36 design = design_formula 36 design = design_formula
37 ) 37 )
38 } else if (!use_txi & !has_header) { 38 } else if (!use_txi && !has_header) {
39 39
40 # construct the object from HTSeq files 40 # construct the object from HTSeq files
41 dds <- DESeqDataSetFromHTSeqCount( 41 dds <- DESeqDataSetFromHTSeqCount(
42 sampleTable = sample_table, 42 sampleTable = sample_table,
43 directory = dir, 43 directory = dir,