R/networkOverlap.R

In neuroconductor-devel/ANTsR: ANTs in R: Quantification Tools for Biomedical Images

# 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)
}