Nothing
wdnc <- function(vnetwork){ # vnetwork is a weighted adjacency matrix
if (!is.matrix(vnetwork)){
stop("Input data must be an adjacency matrix.")
}
if (is.matrix(vnetwork)&(dim(vnetwork)[1]!=dim(vnetwork)[2])){
stop("Input data must be an adjacency matrix (square).")
}
network <- vnetwork
network[network>0] <- 1
n <- dim(network)[1]
W <- network # number of walks for distance i
Sigma <- network # last line values for ditance i-1
wd <- colSums(vnetwork, na.rm=TRUE)
i <- 2
while(min(Sigma)==0){ # finishes when the last shortest path is taken into account
Y <- vnetwork %*% W # last line values for distance i
X <- network %*% W # number of walks between two vertices
TM <- X # temp matrix to extract the shortest paths with length i
TM[TM>0] <- 1
TM <- TM-Sigma # substract the no. of shortest paths till distance i-1
TM[TM<0] <- 0
X <- TM * X
diag(Y) <- 0
Z <- (Y/X)/i # calculate the ewc for distance i
Z[Z == "Inf" | Z == "-Inf"] <- 0
wd <- cbind(wd,colSums(Z, na.rm=TRUE))
Sigma <- Sigma + TM # update no. of shortest paths till distance i for (the next cycle)
W <- W %*% network
i <- i + 1
}
wd <- wd/(max(vnetwork)*(n-1)) # standardization of centralities
ewcmatrix <- t(wd)
ewcmatrix[ewcmatrix=="Inf"] <- 0
cluster <- apply(ewcmatrix, 2, which.max)
names(cluster) <- NULL
cluster
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.