comparison phyloseq_plot_richness.R @ 0:c0101c72b8af draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 5ec9f9e81bb9a42dec5c331dd23215ca0b027b2b
author iuc
date Sat, 16 Mar 2024 07:56:17 +0000
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:c0101c72b8af
1 #!/usr/bin/env Rscript
2
3 suppressPackageStartupMessages(library("optparse"))
4 suppressPackageStartupMessages(library("phyloseq"))
5
6 option_list <- list(
7 make_option(c("--input"), action = "store", dest = "input", help = "Input RDS file containing a phyloseq object"),
8 make_option(c("--output"), action = "store", dest = "output", help = "Output PDF")
9 )
10
11 parser <- OptionParser(usage = "%prog [options] file", option_list = option_list)
12 args <- parse_args(parser, positional_arguments = TRUE)
13 opt <- args$options
14 phyloseq_obj <- readRDS(opt$input)
15 # Start PDF device driver and generate the plot.
16 dev.new()
17 pdf(file = opt$output)
18 plot_richness(phyloseq_obj, x = "samples", color = "samples")
19 dev.off()