## Not run:
# \donttest{
if (requireNamespace("phyloseq", quietly = TRUE)) {
data("GlobalPatterns", package = "phyloseq")
## convert phyloseq object to LorMe format
gp_obj <- Trans_from_phylo(GlobalPatterns, outputtax = "standard")
## specify experimental design
gp_plan <- object_config(gp_obj,
treat_location = 6,
rep_location = 7)
## generate community-structure PCoA plot (Genus level)
community_structure <- structure_plot(gp_plan, taxlevel = "Genus")
community_structure$PCoA_Plot
}
# }
Run the code above in your browser using DataLab