data(varespec)
vare.dist <- vegdist(wisconsin(varespec), "bray")
mds.null <- isoMDS(vare.dist)
## See if you can get any better.
stress <- rep(NA,10)
for (i in 1:10)
stress[i] <- isoMDS(vare.dist, initMDS(vare.dist), maxit=200, trace=F)$stress
rbind(stress, stress < mds.null$stress)
Run the code above in your browser using DataLab