R/middlemanPower.R

In O1sims/networkR: A package for analysing social and economic networks

#' Calculate Middleman power for all nodes
#'
#' A function that calculates the middleman power of each node in a network.
#' @param edgeList A dataframe of network data where sources are in the first column and targets are in the second column.
#' @param nodeList A dataframe where all nodes and their respective names are listed.
#' @param adjMatrix The network represented as an adjacency matrix.
#' @keywords brokerage middleman
#' @export
#' @examples
#' middlemanPower()

middlemanPower <- function(edgeList, nodeList, adjMatrix, normalised) {
  if (missing(adjMatrix)) {
    originalAdjMatrix <- adjMatrix <- adjacencyMatrix(edgeList, nodeList)
  } else {
    originalAdjMatrix <- adjMatrix
  }
  if (missing(normalised)) { normalised <- FALSE }
  PS <- predecessorsSuccessors(edgeList = edgeList,
                               nodeList = nodeList,
                               adjMatrix = adjMatrix)
  K <- connectivity(adjMatrix = originalAdjMatrix)
  power <- 0
  for (i in 1:nrow(adjMatrix)) {
    adjMatrix <- originalAdjMatrix
    adjMatrix[i, ] <- adjMatrix[, i] <- 0
    kappa <- connectivity(adjMatrix = adjMatrix)
    power[i] <- K - kappa - PS$noPred[i] - PS$noSucc[i]
  }
  if (normalised == TRUE) {
    potBroker <- potentialBrokerage(edgeList = edgeList,
                                    nodeList = nodeList,
                                    adjMatrix = originalAdjMatrix)
    power <- round(power/as.integer(potBroker),
                   digits = 3)
  }
  return(power)
}