# Iris data
# Loading the numeric variables of iris data
iris <- as.matrix(iris[,-5])
# standardizing the data
iris <- scale(iris)
# double k-means with 3 unit-clusters and 2 components for the variables
p1 <- redkm(iris, K = 3, Q = 2, Rndstart = 10)
p2 <- doublekm(iris, K=3, Q=2, Rndstart = 10)
mri <- mrand(t(p1$U)%*%p2$U)
Run the code above in your browser using DataLab