# HG changeset patch # User artbio # Date 1670115661 0 # Node ID fb2f1b8b00134c9cb618d71fb924b4fdf0dc74ff # Parent 252eded618481a609f8c9049c5ce124e1132efaf planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_scran_normalize commit e0357f07fdabee1ec6614aca6f7b51095111e0d5 diff -r 252eded61848 -r fb2f1b8b0013 scran-normalize.R --- a/scran-normalize.R Thu Sep 26 10:50:55 2019 -0400 +++ b/scran-normalize.R Sun Dec 04 01:01:01 2022 +0000 @@ -1,59 +1,65 @@ # load packages that are provided in the conda env -options( show.error.messages=F, - error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } ) +options(show.error.messages = FALSE, + error = function() { + cat(geterrmessage(), file = stderr()) + q("no", 1, FALSE) + } +) loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8") warnings() library(optparse) library(scran) -#Arguments -option_list = list( +# Arguments +option_list <- list( make_option( c("-d", "--data"), default = NA, - type = 'character', + type = "character", help = "Input file that contains count values to transform" ), make_option( c("-s", "--sep"), - default = '\t', - type = 'character', + default = "\t", + type = "character", help = "File separator [default : '%default' ]" ), make_option( "--cluster", - default=FALSE, - action="store_true", - type = 'logical', + default = FALSE, + action = "store_true", + type = "logical", help = "Whether to calculate the size factor per cluster or on all cell" ), make_option( c("-m", "--method"), - default = 'hclust', - type = 'character', + default = "hclust", + type = "character", help = "The clustering method to use for grouping cells into cluster : hclust or igraph [default : '%default' ]" ), make_option( "--size", default = 100, - type = 'integer', + type = "integer", help = "Minimal number of cells in each cluster : hclust or igraph [default : '%default' ]" ), make_option( c("-o", "--out"), default = "res.tab", - type = 'character', + type = "character", help = "Output name [default : '%default' ]" ) ) -opt = parse_args(OptionParser(option_list = option_list), +opt <- parse_args(OptionParser(option_list = option_list), args = commandArgs(trailingOnly = TRUE)) -if (opt$sep == "tab") {opt$sep = "\t"} +if (opt$sep == "tab") { + opt$sep <- "\t" + } -data = read.table( +data <- read.table( opt$data, check.names = FALSE, header = TRUE, @@ -62,10 +68,9 @@ ) ## Import data as a SingleCellExperiment object -sce <- SingleCellExperiment(list(counts=as.matrix(data))) +sce <- SingleCellExperiment(list(counts = as.matrix(data))) - -if(opt$cluster){ +if (opt$cluster) { clusters <- quickCluster(sce, min.size = opt$size, method = opt$method) ## Compute sum factors @@ -78,14 +83,14 @@ sce <- normalize(sce) -logcounts <- data.frame(genes = rownames(sce), round(logcounts(sce), digits=5), check.names = F) +logcounts <- data.frame(genes = rownames(sce), round(logcounts(sce), digits = 5), check.names = FALSE) write.table( logcounts, opt$out, - col.names = T, - row.names = F, - quote = F, + col.names = TRUE, + row.names = FALSE, + quote = FALSE, sep = "\t" ) diff -r 252eded61848 -r fb2f1b8b0013 scran_normalize.xml --- a/scran_normalize.xml Thu Sep 26 10:50:55 2019 -0400 +++ b/scran_normalize.xml Sun Dec 04 01:01:01 2022 +0000 @@ -1,4 +1,4 @@ - + Normalize raw counts expression values using deconvolution size factors r-optparse @@ -16,7 +16,7 @@ --method '$metacell.method' --size '$metacell.size' #end if - -o ${output} + -o '${output}' ]]> @@ -92,16 +92,6 @@ - - @Article{, - author = {Aaron T. L. Lun and Davis J. McCarthy and John C. Marioni}, - title = {A step-by-step workflow for low-level analysis of single-cell RNA-seq data with Bioconductor}, - journal = {F1000Res.}, - year = {2016}, - volume = {5}, - pages = {2122}, - doi = {10.12688/f1000research.9501.2}, - } - + 10.12688/f1000research.9501.2