# HG changeset patch # User iuc # Date 1707514936 0 # Node ID 1ff178d1757e3f4bdb96646d70c6fd26d7f95db7 # Parent 92b82deaaed10ec0d027415ed4173a64227c4d66 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2 diff -r 92b82deaaed1 -r 1ff178d1757e macros.xml --- a/macros.xml Thu Mar 03 13:28:55 2022 +0000 +++ b/macros.xml Fri Feb 09 21:42:16 2024 +0000 @@ -2,6 +2,11 @@ 1.38.0 0 21.01 + + + phyloseq + + bioconductor-phyloseq diff -r 92b82deaaed1 -r 1ff178d1757e phyloseq_from_dada2.R --- a/phyloseq_from_dada2.R Thu Mar 03 13:28:55 2022 +0000 +++ b/phyloseq_from_dada2.R Fri Feb 09 21:42:16 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) diff -r 92b82deaaed1 -r 1ff178d1757e phyloseq_plot_ordination.R --- a/phyloseq_plot_ordination.R Thu Mar 03 13:28:55 2022 +0000 +++ b/phyloseq_plot_ordination.R Fri Feb 09 21:42:16 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() diff -r 92b82deaaed1 -r 1ff178d1757e phyloseq_plot_richness.R --- a/phyloseq_plot_richness.R Thu Mar 03 13:28:55 2022 +0000 +++ b/phyloseq_plot_richness.R Fri Feb 09 21:42:16 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() diff -r 92b82deaaed1 -r 1ff178d1757e phyloseq_plot_richness.xml --- a/phyloseq_plot_richness.xml Thu Mar 03 13:28:55 2022 +0000 +++ b/phyloseq_plot_richness.xml Fri Feb 09 21:42:16 2024 +0000 @@ -3,6 +3,7 @@ macros.xml +