R/intraguild_in_strength.R

In RadicalCommEcol/MultitrophicFun: Different functions for working with multitrophic communities

#' species-level metrics
#'
#' different functions with the same structure, metric is in function name
#'
#' @param A interaction matrix
#' @param sp.names which sp to calculate, must be in names of A
#'
#' @return dataframe
#' @export
intraguild_in_strength <- function(A, sp.names){

  result <- data.frame(sp = rownames(A),
                       value = NA_real_,
                       stringsAsFactors = FALSE)

  plant.pos <- which(rownames(A) %in% sp.names[[1]]$plants)
  fv.pos <- which(rownames(A) %in% sp.names[[1]]$floral.visitors)
  herb.pos <- which(rownames(A) %in% sp.names[[1]]$herbivores)

  for(i.sp in 1:nrow(A)){
    if(result$sp[i.sp] %in% sp.names[[1]]$plants){
      result$value[i.sp] <- sum(A[i.sp,plant.pos]) - A[i.sp,i.sp]
    }else if(result$sp[i.sp] %in% sp.names[[1]]$herbivores){
      result$value[i.sp] <- sum(A[i.sp,herb.pos]) - A[i.sp,i.sp]
    }else if(result$sp[i.sp] %in% sp.names[[1]]$floral.visitors){
      result$value[i.sp] <- sum(A[i.sp,fv.pos]) - A[i.sp,i.sp]
    }

  }# for i.sp

  return(result)

}