R/nodal.local.R

In Blaunet: Calculate and Analyze Blau Statuses for Measuring Social Distance

nodal.local <-
function(blauObj, focal.niche=NULL, dev.range, ecologies.off=FALSE){
  
  if (is.null(focal.niche)){
    return("Primary Membership needed for nodal.local")
  }

  blauObj$primMemCol <- correctFormat(focal.niche, blauObj$memberships)

  if (ecologies.off == TRUE){
    blauObj <- niches(blauObj, dev.range, ecologies.off)
  }

  uniqueEcologies <- unique(blauObj$ids[,2])
  if(length(uniqueEcologies) == 1 || ecologies.off == TRUE){ #if there's only one ecology and, we don't have isInNiche

    if(is.null(blauObj$isInNiche)){
      blauObj <- niches(blauObj, dev.range)
    }

    blauObj <- calc.nodal(blauObj, mode = "local") #has isInNiches now; does a bunch of stuff if a primaryMembership is specified
  }

  else if(length(uniqueEcologies) > 1) { #if there's more than one ecology, we need to split up primary membership, weights, dimensions and memberships
    if(is.null(blauObj$isInNiche)){
      blauObj <- niches(blauObj, dev.range)
    }    
    blauObj <- calc.nodal.ecology(blauObj, uniqueEcologies, mode = "local")
  }

  colnames(blauObj$nodalLocal) <- c("FocNicher", "Nicher", "MemNotNiche")
  
  return(blauObj)

}