library(mva)
library(MASS)
data(varespec)
vare.dist <- vegdist(varespec, "bray")
mds.null <- isoMDS(vare.dist)
## See if you can get any better
## (This is a small data set, and an easy case: probably you can't
## improve much).
stress <- rep(NA,10)
for (i in 1:10)
stress[i] <- isoMDS(vare.dist, initMDS(vare.dist), maxit=200, trace=F)
Run the code above in your browser using DataLab