Mercurial > repos > iuc > phyloseq_plot_ordination
changeset 1:92e77800ef2c draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
| author | iuc |
|---|---|
| date | Fri, 09 Feb 2024 21:42:09 +0000 |
| parents | 11d43fa12aab |
| children | dfe800a3faaf |
| files | macros.xml phyloseq_from_dada2.R phyloseq_plot_ordination.R phyloseq_plot_ordination.xml phyloseq_plot_richness.R |
| diffstat | 5 files changed, 32 insertions(+), 36 deletions(-) [+] |
line wrap: on
line diff
--- a/macros.xml Thu Mar 03 13:28:30 2022 +0000 +++ b/macros.xml Fri Feb 09 21:42:09 2024 +0000 @@ -2,6 +2,11 @@ <token name="@TOOL_VERSION@">1.38.0</token> <token name="@VERSION_SUFFIX@">0</token> <token name="@PROFILE@">21.01</token> + <xml name="bio_tools"> + <xrefs> + <xref type="bio.tools">phyloseq</xref> + </xrefs> + </xml> <xml name="requirements"> <requirements> <requirement type="package" version="@TOOL_VERSION@">bioconductor-phyloseq</requirement>
--- a/phyloseq_from_dada2.R Thu Mar 03 13:28:30 2022 +0000 +++ b/phyloseq_from_dada2.R Fri Feb 09 21:42:09 2024 +0000 @@ -10,30 +10,25 @@ make_option(c("--output"), action = "store", dest = "output", help = "RDS output") ) -parser <- OptionParser(usage = "%prog [options] file", option_list = option_list); -args <- parse_args(parser, positional_arguments = TRUE); -opt <- args$options; - +parser <- OptionParser(usage = "%prog [options] file", option_list = option_list) +args <- parse_args(parser, positional_arguments = TRUE) +opt <- args$options # The input sequence_table is an integer matrix # stored as tabular (rows = samples, columns = ASVs). -seq_table_numeric_matrix <- data.matrix(read.table(opt$sequence_table, sep = "\t")); - +seq_table_numeric_matrix <- data.matrix(read.table(opt$sequence_table, sep = "\t")) # The input taxonomy_table is a table containing # the assigned taxonomies exceeding the minBoot # level of bootstrapping confidence. Rows correspond # to sequences, columns to taxonomic levels. NA # indicates that the sequence was not consistently # classified at that level at the minBoot threshold. -tax_table_matrix <- as.matrix(read.table(opt$taxonomy_table, header = FALSE, sep = "\t")); - +tax_table_matrix <- as.matrix(read.table(opt$taxonomy_table, header = FALSE, sep = "\t")) # Construct a tax_table object. The rownames of # tax_tab must match the OTU names (taxa_names) # of the otu_table defined below. -tax_tab <- tax_table(tax_table_matrix); - +tax_tab <- tax_table(tax_table_matrix) # Construct an otu_table object. -otu_tab <- otu_table(seq_table_numeric_matrix, taxa_are_rows = TRUE); - +otu_tab <- otu_table(seq_table_numeric_matrix, taxa_are_rows = TRUE) # Construct a phyloseq object. -phyloseq_obj <- phyloseq(otu_tab, tax_tab); -saveRDS(phyloseq_obj, file = opt$output, compress = TRUE); +phyloseq_obj <- phyloseq(otu_tab, tax_tab) +saveRDS(phyloseq_obj, file = opt$output, compress = TRUE)
--- a/phyloseq_plot_ordination.R Thu Mar 03 13:28:30 2022 +0000 +++ b/phyloseq_plot_ordination.R Fri Feb 09 21:42:09 2024 +0000 @@ -11,20 +11,17 @@ make_option(c("--output"), action = "store", dest = "output", help = "Output") ) -parser <- OptionParser(usage = "%prog [options] file", option_list = option_list); -args <- parse_args(parser, positional_arguments = TRUE); -opt <- args$options; - +parser <- OptionParser(usage = "%prog [options] file", option_list = option_list) +args <- parse_args(parser, positional_arguments = TRUE) +opt <- args$options # Construct a phyloseq object. -phyloseq_obj <- readRDS(opt$input); - +phyloseq_obj <- readRDS(opt$input) # Transform data to proportions as appropriate for # Bray-Curtis distances. -proportions_obj <- transform_sample_counts(phyloseq_obj, function(otu) otu / sum(otu)); -ordination_obj <- ordinate(proportions_obj, method = opt$method, distance = opt$distance); - +proportions_obj <- transform_sample_counts(phyloseq_obj, function(otu) otu / sum(otu)) +ordination_obj <- ordinate(proportions_obj, method = opt$method, distance = opt$distance) # Start PDF device driver and generate the plot. -dev.new(); -pdf(file = opt$output); -plot_ordination(proportions_obj, ordination_obj, type = opt$type); -dev.off(); +dev.new() +pdf(file = opt$output) +plot_ordination(proportions_obj, ordination_obj, type = opt$type) +dev.off()
--- a/phyloseq_plot_ordination.xml Thu Mar 03 13:28:30 2022 +0000 +++ b/phyloseq_plot_ordination.xml Fri Feb 09 21:42:09 2024 +0000 @@ -3,6 +3,7 @@ <macros> <import>macros.xml</import> </macros> + <expand macro="bio_tools"/> <expand macro="requirements"/> <command detect_errors="exit_code"><![CDATA[ Rscript '${__tool_directory__}/phyloseq_plot_ordination.R'
--- a/phyloseq_plot_richness.R Thu Mar 03 13:28:30 2022 +0000 +++ b/phyloseq_plot_richness.R Fri Feb 09 21:42:09 2024 +0000 @@ -8,14 +8,12 @@ make_option(c("--output"), action = "store", dest = "output", help = "Output PDF") ) -parser <- OptionParser(usage = "%prog [options] file", option_list = option_list); -args <- parse_args(parser, positional_arguments = TRUE); -opt <- args$options; - -phyloseq_obj <- readRDS(opt$input); - +parser <- OptionParser(usage = "%prog [options] file", option_list = option_list) +args <- parse_args(parser, positional_arguments = TRUE) +opt <- args$options +phyloseq_obj <- readRDS(opt$input) # Start PDF device driver and generate the plot. -dev.new(); -pdf(file = opt$output); -plot_richness(phyloseq_obj, x = "samples", color = "samples"); +dev.new() +pdf(file = opt$output) +plot_richness(phyloseq_obj, x = "samples", color = "samples") dev.off()