In diazrenata/neonbecs: What the Package Does (One Line, Title Case)

library(replicatebecs)
communities <- load_paper_data()
communities_energy <- lapply(communities, add_energy_sizeclass)
isds <- lapply(communities_energy, make_isd)

Exploring the multimodality of rodent ISDs using Gaussian mixture models (after Thibault et al 2011).

isd_plots <- plot_paper_dists(isds, "isd")
invisible(isd_plots)

gmms <- lapply(isds, fit_gmm)
names(gmms) <- names(communities)

for(i in 1:9) {
  print(names(gmms)[i])
  plot_gmm(gmms[[i]])
}

nmodes <- vapply(gmms, FUN = get_nmodes, FUN.VALUE = 3)

nmodes

diazrenata/neonbecs documentation built on July 6, 2019, 9:25 a.m.