# HG changeset patch # User iuc # Date 1580328458 18000 # Node ID cdf4443d56250a58bd3cddf58400a8e97e9fb658 # Parent 1797cbe967d12c927a8e3780c8dc484a17d7b9c4 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/dropletutils/ commit ea9a67b5711b53fb911a0c50b6b49eb9db2366a7" diff -r 1797cbe967d1 -r cdf4443d5625 dropletutils.xml --- a/dropletutils.xml Thu Dec 05 01:32:29 2019 -0500 +++ b/dropletutils.xml Wed Jan 29 15:07:38 2020 -0500 @@ -3,7 +3,7 @@ Utilities for handling droplet-based single-cell RNA-seq data 1.2.1 - 5 + galaxy6 tenx.input tenx.output @@ -80,8 +80,8 @@ out.type='$operation.outformat' #if str($operation.outformat) == 'directory': files\$out='@TXOUT@' - #else if str($operation.outformat) == 'h5ad': -files\$out='$fileout_h5ad' + #else if str($operation.outformat) == 'h5': +files\$out='$fileout_h5' #else if str($operation.outformat) == 'tsv': files\$out='$fileout_tsv' #end if @@ -99,7 +99,7 @@ - + @@ -107,8 +107,8 @@ - - + + @@ -137,7 +137,7 @@ - + @@ -151,8 +151,8 @@ operation['use']=='filter' and operation['outformat'] == 'tsv' - - operation['use']=='filter' and operation['outformat'] == 'h5ad' + + operation['use']=='filter' and operation['outformat'] == 'h5' operation['use']=='filter' and operation['outformat'] == 'directory' @@ -177,9 +177,9 @@ - + - + @@ -210,14 +210,14 @@ - + - + @@ -229,22 +229,22 @@ - + - - + + - + - + diff -r 1797cbe967d1 -r cdf4443d5625 scripts/dropletutils.Rscript --- a/scripts/dropletutils.Rscript Thu Dec 05 01:32:29 2019 -0500 +++ b/scripts/dropletutils.Rscript Wed Jan 29 15:07:38 2020 -0500 @@ -34,7 +34,7 @@ if (typeout == "tsv"){ writeTSV(fileout, getCounts(object)) } - else if (typeout == "h5ad"){ + else if (typeout == "h5"){ write10xCounts(fileout, getCounts(object), type="HDF5", gene.symbol=determineGeneIDs(object), @@ -55,7 +55,7 @@ ## - solution is to use the readSparseCounts function from scater sce <- SingleCellExperiment(assays = list(counts = readSparseCounts(filein))) } - else if (typein == "h5ad"){ + else if (typein == "h5"){ sce <- read10xCounts(filein, col.names=T, type="HDF5") # use barcodes.tsv as column names } else if (typein == "directory"){ @@ -69,7 +69,7 @@ ## Methods -doEmptyDrops <- function(files, eparams, in.type="directory", out.type="h5ad", fdr_threshold = 0.01){ +doEmptyDrops <- function(files, eparams, in.type="directory", out.type="h5", fdr_threshold = 0.01){ sce <- read10xFiles(files$infile, in.type) eparams$... <- NULL ## hack @@ -115,7 +115,7 @@ } -doDefaultDrops <- function(files, dparams, in.type="directory", out.type="h5ad"){ +doDefaultDrops <- function(files, dparams, in.type="directory", out.type="h5"){ sce <- read10xFiles(files$infile, in.type) dparams$m = counts(sce) diff -r 1797cbe967d1 -r cdf4443d5625 test-data/defs_defaultdrops.h5 Binary file test-data/defs_defaultdrops.h5 has changed diff -r 1797cbe967d1 -r cdf4443d5625 test-data/defs_defaultdrops.h5ad Binary file test-data/defs_defaultdrops.h5ad has changed diff -r 1797cbe967d1 -r cdf4443d5625 test-data/defs_emptydrops_150_0002.h5 Binary file test-data/defs_emptydrops_150_0002.h5 has changed diff -r 1797cbe967d1 -r cdf4443d5625 test-data/defs_emptydrops_150_0002.h5ad Binary file test-data/defs_emptydrops_150_0002.h5ad has changed diff -r 1797cbe967d1 -r cdf4443d5625 test-data/defs_emptydrops_150_0002a.h5 Binary file test-data/defs_emptydrops_150_0002a.h5 has changed diff -r 1797cbe967d1 -r cdf4443d5625 test-data/defs_emptydrops_150_0002a.h5ad Binary file test-data/defs_emptydrops_150_0002a.h5ad has changed diff -r 1797cbe967d1 -r cdf4443d5625 test-data/in_matrix.h5 Binary file test-data/in_matrix.h5 has changed diff -r 1797cbe967d1 -r cdf4443d5625 test-data/in_matrix.h5ad Binary file test-data/in_matrix.h5ad has changed