# NOT RUN {
iris.dist <- dist(iris[, 1:4], method="manhattan")
iris.dist[iris.dist == 0] <- abs(jitter(0))
library(MASS)
iris.m <- isoMDS(iris.dist)
cor(iris[, 1:4], iris.m$points) # MDS loadings surrogate
MDSv(iris.m$points) # MDS explained variance surrogate
# }
Run the code above in your browser using DataLab