Mercurial > repos > iuc > phyloseq_plot_richness
annotate phyloseq_plot_richness.R @ 10:7ed99cb39c9d draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit fa72d860839082a926004d8a97a03a3e27701333
| author | iuc | 
|---|---|
| date | Fri, 04 Apr 2025 10:15:49 +0000 | 
| parents | 1ff178d1757e | 
| children | 
| 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 | 
| 
1
 
1ff178d1757e
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
 
iuc 
parents: 
0 
diff
changeset
 | 
11 parser <- OptionParser(usage = "%prog [options] file", option_list = option_list) | 
| 
 
1ff178d1757e
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
 
iuc 
parents: 
0 
diff
changeset
 | 
12 args <- parse_args(parser, positional_arguments = TRUE) | 
| 
 
1ff178d1757e
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
 
iuc 
parents: 
0 
diff
changeset
 | 
13 opt <- args$options | 
| 
 
1ff178d1757e
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
 
iuc 
parents: 
0 
diff
changeset
 | 
14 phyloseq_obj <- readRDS(opt$input) | 
| 
0
 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 
iuc 
parents:  
diff
changeset
 | 
15 # Start PDF device driver and generate the plot. | 
| 
1
 
1ff178d1757e
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
 
iuc 
parents: 
0 
diff
changeset
 | 
16 dev.new() | 
| 
 
1ff178d1757e
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
 
iuc 
parents: 
0 
diff
changeset
 | 
17 pdf(file = opt$output) | 
| 
 
1ff178d1757e
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
 
iuc 
parents: 
0 
diff
changeset
 | 
18 plot_richness(phyloseq_obj, x = "samples", color = "samples") | 
| 
0
 
92b82deaaed1
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
 
iuc 
parents:  
diff
changeset
 | 
19 dev.off() | 
