nem <- read_nem2(tab = nemtab, tax = nemtax, meta = nemmeta)
nem_index <- nem |>
calc_alpha() |>
calc_nemindex() |>
calc_lm(group = Treatments,
x = Chao1,
y = TotalBiomass)
nem_index
Run the code above in your browser using DataLab