Mercurial > repos > iuc > edger

comparison edger.R @ 13:0cb907a2a810 draft

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/edger commit 95a36f5b0bd139f161d131b4c2c9c36940e749c6
author iuc
date Wed, 27 Sep 2023 19:53:04 +0000
parents a8305cf0c951
children 070900306913
comparison
equal deleted inserted replaced
12:a8305cf0c951 13:0cb907a2a810
281 if (!any(factordata[, 1] %in% colnames(counts))) { 281 if (!any(factordata[, 1] %in% colnames(counts))) {
282 stop("Sample IDs in factors file and count matrix don't match") 282 stop("Sample IDs in factors file and count matrix don't match")
283 } 283 }
284 # order samples as in counts matrix 284 # order samples as in counts matrix
285 factordata <- factordata[match(colnames(counts), factordata[, 1]), ] 285 factordata <- factordata[match(colnames(counts), factordata[, 1]), ]
286 factors <- factordata[, -1, drop = FALSE] 286 factors <- data.frame(sapply(factordata[, -1, drop = FALSE], make.names))
287 } else { 287 } else {
288 factors <- unlist(strsplit(opt$factInput, "|", fixed = TRUE)) 288 factors <- unlist(strsplit(opt$factInput, "|", fixed = TRUE))
289 factordata <- list() 289 factordata <- list()
290 for (fact in factors) { 290 for (fact in factors) {
291 newfact <- unlist(strsplit(fact, split = "::")) 291 newfact <- unlist(strsplit(fact, split = "::"))