Mercurial > repos > iuc > phyloseq_plot_richness
annotate phyloseq_plot_richness.R @ 0:92b82deaaed1 draft
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
| author | iuc | 
|---|---|
| date | Thu, 03 Mar 2022 13:28:55 +0000 | 
| parents | |
| children | 1ff178d1757e | 
| rev | line source | 
|---|---|
| 0 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 1 #!/usr/bin/env Rscript | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 2 | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 3 suppressPackageStartupMessages(library("optparse")) | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 4 suppressPackageStartupMessages(library("phyloseq")) | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 5 | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 6 option_list <- list( | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 7 make_option(c("--input"), action = "store", dest = "input", help = "Input RDS file containing a phyloseq object"), | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 8 make_option(c("--output"), action = "store", dest = "output", help = "Output PDF") | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 9 ) | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 10 | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 11 parser <- OptionParser(usage = "%prog [options] file", option_list = option_list); | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 12 args <- parse_args(parser, positional_arguments = TRUE); | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 13 opt <- args$options; | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 14 | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 15 phyloseq_obj <- readRDS(opt$input); | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 16 | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 17 # Start PDF device driver and generate the plot. | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 18 dev.new(); | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 19 pdf(file = opt$output); | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 20 plot_richness(phyloseq_obj, x = "samples", color = "samples"); | 
| 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 iuc parents: diff
changeset | 21 dev.off() | 
