# NOT RUN {
iris.dist <- dist(iris[, -5], method="manhattan")
library(MASS)
iris.m <- isoMDS(iris.dist + 1e-9) # isoMDS() does not like zero distances
(vv <- MDSv(iris.m$points)) # MDS explained variance surrogate
xlab <- paste0("Dim 1 (", round(vv[1], 2), "%)")
ylab <- paste0("Dim 2 (", round(vv[2], 2), "%)")
plot(iris.m$points, col=as.numeric(iris$Species), xlab=xlab, ylab=ylab)
iris.cmd <- cmdscale(iris.dist)
MDSv(iris.cmd)
iris.p <- prcomp(iris[, -5], scale=TRUE)
MDSv(iris.p$x)
100*summary(iris.p)$importance[2, ] # compare with MDSv() results
# }
Run the code above in your browser using DataLab