R/NetworkOverlap.R

In stnava/itkImageR: Examples using Rcpp to interface R and ITK

# NetworkOverlap.R - Essentially the DICE Metric, but for network edge overlap n1 - network matrix n2 -
# network matrix
networkOverlap <- function(n1, n2) {
  
  intersect <- sum(n1[upper.tri(n1)] * n2[upper.tri(n2)])
  union <- sum(n1[upper.tri(n1)]) + sum(n2[upper.tri(n2)])
  
  overlap <- 0
  if (union > 0) 
    overlap <- 2 * intersect/union
  
  return(overlap)
}