#' Calculate species richness
#'
#' @description Calculates species richness, by year, site, and transect, for a specific location.
#'
#' @param data A dataframe that contains at least columns of Ano, Zona (e.g. "Zona de pesca" or "Zona de no pesca"), Sitio, Transecto, GeneroEspecie.
#' @param location A quoted string that indicates the location.
#'
#'@return S A dataframe with columns for Ano, Zona, Sitio, Transecto, and Species richness (S; number of species).
#'
#'@export
richness <- function(data, location){
library(dplyr) #Load dplyr
library(tidyr) #Load tidyr
S <- data %>% #Set S equal to data
filter(Comunidad == location) %>% #Filter by location
filter(Abundancia > 0) %>%
group_by(Ano,
Zona,
Sitio,
Transecto,
Temperatura,
Visibilidad,
ProfundidadInicial,
GeneroEspecie) %>%
summarize(N = n()) %>%
group_by(Ano,
Zona,
Sitio,
Transecto,
Temperatura,
Visibilidad,
ProfundidadInicial) %>%
summarize(S = n()) %>%
ungroup() %>%
select(Ano, Zona, Sitio, Transecto, Temperatura, Visibilidad, Profundidad = ProfundidadInicial, Indicador = S)
return(as.data.frame(S))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.