if (requireNamespace("NST") && requireNamespace("geosphere")) {
library(ggplot2)
data(otutab, package = "pcutils")
metadata[, c("lat", "long")] -> geo
geo_sim(otutab, geo) -> geo_res
pcutils::my_lm(geo_res[4], "dis.geo", geo_res)
}
Run the code above in your browser using DataLab