#NMDS
library(vegan)
otu=read.table("clipboard", header=T,row.names=1)
otu=t(otu)
dist=vegdist(otu,method="bray")
nmds=monoMDS(dist)
nmds
scores(nmds, choices = c(1,2))