title: Labelling NMDS plots in R
body: <p>[enter image description here][1]
  I am having an issue with metaMDS and ordiplot from the vegan package. I am new to this so I can't come up with an ultra flashy example. </p>
  
  <p>Using ordiplot I would like to label plots according to a type (the groups A, B, C, D, E). I have a data frame "df" with a column of plot names and a column of the type. I NO LONGER care whether the type is illustrated by color or with different types of points. 
  my very basic example:</p>
  
  <pre><code>data(dune)
  data(dune.env)
  dist.dune&lt;-distShell(dune, distSimRatio)
  ord&lt;-metaMDS(dist.dune)
  metaMDS(comm=dune)
  plot(ord, type="n")
  colfactor&lt;-factor(dune.env$Management)
  points(ord, display="sites", cex=2, pch=21, col=c(1,2,3,4,5)[colfactor])
  text(ord, display="sites", cex=2, col=c(1,2,3,4,5)[colfactor])
  </code></pre>
  
  <p>UPDATE: Checked through a friends script and found a slightly different solution. This works on the dune data and is easily verfiable.
  The same solution does NOT work on my data however. I can pretty easily check my NMDS results (the attribute I am trying to color by is part of my point names) and their coloring makes no sense, it is however correct in the sense that the correct amount of plots is colored in one color. It is also always the same points that are grouped into one color, seemingly by random. </p>
  
  <p>UPDATE: warning message has disappeared. Now I only get the warning that the species scores are not available, which is clear, since I used a distance matrix for the calculation</p>
  
  <p>UPDATE: [Look at the image below, the letter after the underscore is the plot type, same as the type I used in colfactor. Clearly the results are off :/][2]</p>
  
  <p>UPDATE: [the last letter of the labels specifies the type like in my table, so strR should be where all the triangles are, since both of the labelling commands are identical][1]</p>