dat <- rbind(matrix(rnorm(100, mean = 0, sd = 0.3), ncol = 2),
matrix(rnorm(100, mean = 2, sd = 0.3), ncol = 2))
sim <- computeGaussianSimilarity(dat, 1)
res <- bipartitionShi(sim)
plot(dat[,1], dat[,2], type = "p", xlab = "x", ylab = "y",
col = res$label, main = "Initial features space")
plot(res$eigenvector, type = "p", xlab = "Indices", ylab = "1st eigenvector",
col = res$label, main = "Spectral embedding")
Run the code above in your browser using DataLab