NMDS (Non-metric multidimensional scaling)
So this works for my data! I had started with the commands under the section NMDS from http://userweb.eng.gla.ac.uk/umer.ijaz/bioinformatics/ecological.html the sample and data samples were prepared according to the above. and proceed with the following: library("vegan", lib.loc="~/R/win-library/3.2") library("lattice", lib.loc="C:/Program Files/R/R-3.2.5/library") set.seed(2) abund_table<-read.csv("Sample.csv",row.names=1,check.names=FALSE) #Transpose the data to have sample names on rows abund_table<-t(abund_table) abund_table<-t(abund_table) example_NMDS=metaMDS(abund_table, k=2) plot(example_NMDS) ordiplot(example_NMDS,type="n") orditorp(example_NMDS,display="species",col="red",air=0.01) orditorp(example_NMDS,display="sites",cex=1.25,air=0.01) And I changed "Treatment 1" and "Treatment 2" to: treat=c(rep("pH",5),rep("Cond