Mercurial > repos > iuc > ampvis2_boxplot
view test-data/generate.R @ 8:5e4e2b612a4f draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/ampvis2 commit 9fe04d29ea604a152144908dbd20c0754695a025
author | iuc |
---|---|
date | Sat, 16 Nov 2024 19:57:11 +0000 (5 weeks ago) |
parents | 0c069fc2471c |
children |
line wrap: on
line source
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") # construct data for merge_ampvis2 t <- amp_load(otutable = "AalborgWWTPs.otu.csv", metadata = "AalborgWWTPs.tsv", taxonomy = "AalborgWWTPs.tax.tsv", fasta = "AalborgWWTPs.fa", tree = "AalborgWWTPs.nwk") d_2010 <- amp_subset_samples(t, Period %in% "2010") # 60 samples and 109 OTUs have been filtered # Before: 67 samples and 200 OTUs # After: 7 samples and 91 OTUs d_2011 <- amp_subset_samples(t, Period %in% "2011") # 61 samples and 94 OTUs have been filtered # Before: 67 samples and 200 OTUs # After: 6 samples and 106 OTUs saveRDS(d_2010, "AalborgWWTPs.2010.rds") saveRDS(d_2011, "AalborgWWTPs.2011.rds") # => merging should give 13 samples