Mercurial > repos > iuc > ampvis2_export_fasta

diff test-data/generate.R @ 0:3c3c359d9b85 draft

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/ampvis2 commit 7c0ecbffdb5e993f5af7e3b52c424c2761fb91d3"
author iuc
date Mon, 04 Apr 2022 10:22:24 +0000
parents
children e10b5eebd415
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/generate.R	Mon Apr 04 10:22:24 2022 +0000
@@ -0,0 +1,10 @@
+library(ampvis2)
+
+# subset taxa using 200 random OTUs
+aalborgwwtps <- amp_subset_taxa(AalborgWWTPs, tax_vector = sample(AalborgWWTPs$tax$OTU, 200))
+
+ape::write.tree(aalborgwwtps$tree, "AalborgWWTPs.nwk")
+amp_export_fasta(aalborgwwtps, "AalborgWWTPs.fa", tax = F)
+amp_export_otutable(aalborgwwtps, "AalborgWWTPs.otu")
+write.table(aalborgwwtps$tax, file = "AalborgWWTPs.tax", quote = F, sep = "\t", row.names = F)
+write.table(aalborgwwtps$metadata, file = "AalborgWWTPs.tsv", quote = F, sep = "\t")