R/RichnessStatistics.R

In mikemahoney218/BeaverForaging:

#' A function to calculate species richness across all plots.
#'
#' @import dplyr
#'

RichnessStatistics <- function(){
  SiteRichness <- InventoryData %>%
    group_by(Site) %>%
    summarise(n = n_distinct(Species))

  RichnessCorrelation <- InventoryData %>%
    group_by(Site, Transect, PlotCode) %>%
    summarise(n = n_distinct(Species))

  RichnessCorrelation <- cor.test(RichnessCorrelation$PlotCode, RichnessCorrelation$n, method = "spearman")
  Richness <<- list()
  Richness$SiteRichness <<- SiteRichness
  Richness$Correlation <<- RichnessCorrelation
  return("Type Richness$SiteRichness to retrieve a table of species richness by site, \nor Richness$Correlation for the output of the Pearson correlation test.")
}