# Draw a plot for a non-vegan ordination (cmdscale). data(dune) dune.dis <- vegdist(wisconsin(dune)) dune.mds <- cmdscale(dune.dis, eig = TRUE) dune.mds$species <- wascores(dune.mds$points, dune, expand = TRUE) fig <- ordiplot(dune.mds, type = "none") points(fig, "sites", pch=21, col="red", bg="yellow") text(fig, "species", col="blue", cex=0.9) # Default plot of the previous using identify to label selected points ## Not run: # fig <- ordiplot(dune.mds) # identify(fig, "spec")## End(Not run)
Run the code above in your browser using DataCamp Workspace