# HG changeset patch # User ebi-gxa # Date 1695777092 0 # Node ID 046d8ff974ff487e5b69328f7f5a4c2b1eff2ca7 # Parent 59a7f3f83aec21fee119799579d36dffd2166716 planemo upload for repository https://github.com/ebi-gene-expression-group/container-galaxy-sc-tertiary/ commit 45145f380b27c3092e1fa2249adc36d7d6fdf5fe diff -r 59a7f3f83aec -r 046d8ff974ff decoupler_pseudobulk.py --- a/decoupler_pseudobulk.py Sun Sep 24 08:44:24 2023 +0000 +++ b/decoupler_pseudobulk.py Wed Sep 27 01:11:32 2023 +0000 @@ -58,15 +58,15 @@ # avoid dash that is read as point on R colnames. obs_for_deseq.index = obs_for_deseq.index.str.replace("-", "_") obs_for_deseq.index = obs_for_deseq.index.str.replace(" ", "_") - col_metadata_file = f"{output_dir}col_metadata.csv" + col_metadata_file = f"{output_dir}col_metadata.tsv" # write obs to a col_metadata file if factor_fields: # only output the index plus the columns in factor_fields in that order - obs_for_deseq[factor_fields].to_csv(col_metadata_file, sep=",", index=True) + obs_for_deseq[factor_fields].to_csv(col_metadata_file, sep="\t", index=True) else: - obs_for_deseq.to_csv(col_metadata_file, sep=",", index=True) + obs_for_deseq.to_csv(col_metadata_file, sep="\t", index=True) # write var to a gene_metadata file - pdata.var.to_csv(f"{output_dir}gene_metadata.csv", sep=",", index=True) + pdata.var.to_csv(f"{output_dir}gene_metadata.tsv", sep="\t", index=True) # write the counts matrix of a specified layer to file if layer is None: # write the X numpy matrix transposed to file @@ -75,7 +75,7 @@ df = pd.DataFrame( pdata.layers[layer].T, index=pdata.var.index, columns=obs_for_deseq.index ) - df.to_csv(f"{output_dir}counts_matrix.csv", sep=",", index_label="") + df.to_csv(f"{output_dir}counts_matrix.tsv", sep="\t", index_label="") def plot_pseudobulk_samples( diff -r 59a7f3f83aec -r 046d8ff974ff decoupler_pseudobulk.xml --- a/decoupler_pseudobulk.xml Sun Sep 24 08:44:24 2023 +0000 +++ b/decoupler_pseudobulk.xml Wed Sep 27 01:11:32 2023 +0000 @@ -1,4 +1,4 @@ - + aggregates single cell RNA-seq data for running bulk RNA-seq methods decoupler @@ -74,9 +74,9 @@ produce_anndata - - - + + + produce_plots @@ -105,19 +105,22 @@ - + + - + + - + + @@ -153,9 +156,11 @@ The tool will output the filtered AnnData, count matrix, samples metadata, genes metadata (in DESeq2 format), and the pseudobulk plot and filter by expression plot (if enabled). + You can obtain more information about Decoupler pseudobulk at the developers documentation `here `_ . + ]]> - doi.org/10.1093/bioadv/vbac016 + 10.1093/bioadv/vbac016