Mercurial > repos > iuc > phyloseq_plot_ordination

annotate phyloseq_plot_ordination.R @ 2:dfe800a3faaf draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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:05 +0000
parents 92e77800ef2c
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
11d43fa12aab "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
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
2
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
3 suppressPackageStartupMessages(library("optparse"))
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
4 suppressPackageStartupMessages(library("phyloseq"))
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
5
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
6 option_list <- list(
11d43fa12aab "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 file containing a phyloseq object"),
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
8 make_option(c("--method"), action = "store", dest = "method", help = "Ordination method"),
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
9 make_option(c("--distance"), action = "store", dest = "distance", help = "Distance method"),
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
10 make_option(c("--type"), action = "store", dest = "type", help = "Plot type"),
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
11 make_option(c("--output"), action = "store", dest = "output", help = "Output")
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
12 )
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
13
1
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
14 parser <- OptionParser(usage = "%prog [options] file", option_list = option_list)
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
15 args <- parse_args(parser, positional_arguments = TRUE)
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
16 opt <- args$options
0
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
17 # Construct a phyloseq object.
1
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
18 phyloseq_obj <- readRDS(opt$input)
0
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
19 # Transform data to proportions as appropriate for
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
20 # Bray-Curtis distances.
1
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
21 proportions_obj <- transform_sample_counts(phyloseq_obj, function(otu) otu / sum(otu))
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
22 ordination_obj <- ordinate(proportions_obj, method = opt$method, distance = opt$distance)
0
11d43fa12aab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit d1004c06207be773c278e12745aada276b63172e"
iuc
parents:
diff changeset
23 # Start PDF device driver and generate the plot.
1
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
24 dev.new()
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
25 pdf(file = opt$output)
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
26 plot_ordination(proportions_obj, ordination_obj, type = opt$type)
92e77800ef2c planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/phyloseq commit 7df921baa7aa8680421b9440a1cd6eaab1a15ce2
iuc
parents: 0
diff changeset
27 dev.off()