title: using ordiellipse in NMDS in R on axes 2 and 3
body: <p>I am trying to apply NMDS by vegan package on a dataset by the following functions:</p>
  
  <pre><code>library(vegan)
  library(MASS)
  library(MVA)
  library(ggplot2)
  all_d &lt;- vegdist(landmetr, method= "hor", binary=FALSE, diag=TRUE,
               upper=FALSE, na.rm = FALSE)
  allmds3 &lt;- isoMDS(all_d, k=3)
  gr0 &lt;- factor(france100_7c$LAND_CODE)
  plot(allmds3$points[,2:3], col=gr0, bg="grey", cex=0.2,pch=20) 
  # regress environmental variables on configuration of site points
  provaenv3 &lt;- envfit(allmds3$points[,2:3], landmetr,999, na.rm = FALSE, france100_7c$LAND_CODE)
  # add the environmental regressions to the area plot
  plot(provaenv3,add=T, cex=0.5)
  #plot0
  ordiellipse (allmds3, gr0, kind = "sd", conf=0.95, label = T)
  </code></pre>
  
  <p>However, despite a this biplot is constructed on axes  2 and 3, the ordiellipse function write the ellipses which identify clusters on the axes 1 and 2 biplot.
  Could anyone explain me why this happens and how could I plot the right ellipses on axes 2 and 3 biplot?</p>