extract_indices_change(ind)

Absolute difference in functional indices across two conditions.

statisticsfunctional diversity

Args:ind — named list of index matrices (before / after)

extract_indices_change <- function(ind) {
  indSES <- matrix(
    NA, nr = ncol(ind), nc = 7,
    dimnames = list(
      colnames(ind),
      c("Observed", "SES", "MeanRd", "CI025Rd", "CI975Rd", "Pval", "Nreps")
    )
  )
  for (i in 1:nrow(indSES)) {
    indSES[i, ] <- sesandpvalues_bis(
      obs = ((ind[2, i] - ind[1, i]) / ind[1, i]) * 100,
      rand = ((ind[-c(1:2), i] - ind[1, i]) / ind[1, i]) * 100,
      nreps = length(ind[, 1]) - 2,
      probs = c(0.025, 0.975),
      rnd = 2
    )
  }
  rownames(indSES) <- gsub("\\.", " ", rownames(indSES))
  return(indSES)
}