# 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 @@
-
+
-
+
-
+
-
-
+
+
-
+
-
+
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