# \donttest{
if (requireNamespace("zetadiv")) {
data(otutab, package = "pcutils")
zeta_result <- z_diversity(otutab, metadata["Group"], zetadiv_params = list(sam = 10))
plot(zeta_result, lm_model = "exp", text = TRUE)
}
# }
Run the code above in your browser using DataLab