extract_indices_change_log(ind)

Log-ratio change in functional indices across two conditions.

statisticsfunctional diversity

Args:ind — named list of index matrices

extract_indices_change_log <- 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 = log(ind[2, i] / ind[1, i]),
      rand = log(ind[-c(1:2), i] / ind[1, i]),
      nreps = length(ind[, 1]) - 2,
      probs = c(0.025, 0.975),
      rnd = 2
    )
  }
  rownames(indSES) <- gsub("\\.", " ", rownames(indSES))
  return(indSES)
}