In allanecology/BetaDivMultifun: Testing effects of multitrophic beta-diversity and land-use intensity (LUI) on ecosystem multifunctionality

#TODO : maybe delete this script?

Depends on spdiv which is loaded in calc_betadiversities.Rmd

Exploring dataset

trlevels <- sort(unique(spdiv$Trophic_level))
overview <- matrix(NA, 1, 2, dimnames = list(c(), c("Trophic_level", "years")))
for(t in trlevels){
  y <- unique(spdiv[Trophic_level == t, Year])
  del <- cbind("Trophic_level" = rep(t, length(y)), "years" = y)
  overview <- rbind(overview, del)
}
overview <- overview[-1,]
overview <- data.table::data.table(overview)
overview[ years == "2011_2012", years := "2011.5"]
overview[years == "2014_2015", years := "2014.5"]
overview[, years := as.numeric(years)]

Years

library(cowplot)
ggplot(overview, aes(years, Trophic_level)) +
        geom_point() +
  theme(axis.text = element_text(angle = 45, hjust = 1, size=10))

unique(spdiv[Trophic_level == "ant.omnivore", .(DataID, Year)])

Number of species

trlevels <- sort(unique(spdiv$Trophic_level))
overview <- matrix(NA, 1, 2, dimnames = list(c(), c("Trophic_level", "Species")))
for(t in trlevels){
  y <- length(unique(spdiv[Trophic_level == t, Species]))
  del <- cbind("Trophic_level" = t, "Species" = y)
  overview <- rbind(overview, del)
}
overview <- overview[-1,]
overview <- data.table::data.table(overview)

allanecology/BetaDivMultifun documentation built on Nov. 9, 2023, 8:47 p.m.