Mercurial > repos > iuc > raceid_clustering
diff scripts/cluster.R @ 9:0bff0ee0683a draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/raceid3 commit 82a18e57158136e265a26f27feb40f8dc13437bf
author | iuc |
---|---|
date | Wed, 24 Aug 2022 18:09:06 +0000 |
parents | a4b734cd253b |
children |
line wrap: on
line diff
--- a/scripts/cluster.R Mon Dec 20 10:12:47 2021 +0000 +++ b/scripts/cluster.R Wed Aug 24 18:09:06 2022 +0000 @@ -1,7 +1,7 @@ #!/usr/bin/env R VERSION <- "0.5" # nolint -args <- commandArgs(trailingOnly = T) +args <- commandArgs(trailingOnly = TRUE) if (length(args) != 1) { message(paste("VERSION:", VERSION)) @@ -52,8 +52,8 @@ if (filt.use.ccorrect) { par(mfrow = c(2, 2)) sc <- do.call(CCcorrect, c(sc, filt.ccc)) - print(plotdimsat(sc, change = T)) - print(plotdimsat(sc, change = F)) + print(plotdimsat(sc, change = TRUE)) + print(plotdimsat(sc, change = FALSE)) } return(sc) } @@ -62,8 +62,8 @@ sc <- do.call(compdist, c(sc, clust.compdist)) sc <- do.call(clustexp, c(sc, clust.clustexp)) if (clust.clustexp$sat) { - print(plotsaturation(sc, disp = F)) - print(plotsaturation(sc, disp = T)) + print(plotsaturation(sc, disp = FALSE)) + print(plotsaturation(sc, disp = TRUE)) } print(plotjaccard(sc)) return(sc) @@ -125,7 +125,7 @@ print(do.call(mtext, c(paste(buffer, "(fc > ", genelist.foldchange, ")"), test))) }) - write.table(df, file = out.genelist, sep = "\t", quote = F) + write.table(df, file = out.genelist, sep = "\t", quote = FALSE) } @@ -138,7 +138,7 @@ is.outlier = names(dat) %in% sc@out$out) write.table(tab, file = out.assignments, sep = "\t", - quote = F, row.names = F) + quote = FALSE, row.names = FALSE) } @@ -158,7 +158,7 @@ getfdata(sc))) / ncol(sc@expdata)), "% of cells remain")) write.table(as.matrix(sc@ndata), file = out.table, col.names = NA, - row.names = T, sep = "\t", quote = F) + row.names = TRUE, sep = "\t", quote = FALSE) } if (use.cluster) {