if (requireNamespace("zetadiv")) {
data(otutab, package = "pcutils")
zeta_decay_result <- z_diversity_decay(otutab, metadata[, c("lat", "long")],
metadata["Group"],
zetadiv_params = list(sam = 10)
)
plot(zeta_decay_result)
}
Run the code above in your browser using DataLab