set.seed(1)
X <- matrix(rnorm(50), nrow = 10)
D <- as.matrix(dist(X))
r <- 3
S <- exp(-D^2)
SNN <- compute.SNN(S, r)
head(SNN, 5)
# inspect average SNN similarity to nearest neighbour
rowMeans(SNN - diag(diag(SNN)))
Run the code above in your browser using DataLab