vegan - using ordiellipse in NMDS in R on axes 2 and 3 -

i trying apply nmds vegan package on dataset following functions:

library(vegan) library(mass) library(mva) library(ggplot2) all_d <- vegdist(landmetr, method= "hor", binary=false, diag=true,              upper=false, na.rm = false) allmds3 <- isomds(all_d, k=3) gr0 <- 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 <- envfit(allmds3$points[,2:3], landmetr,999, na.rm = false, france100_7c$land_code) # add environmental regressions area plot plot(provaenv3,add=t, cex=0.5) #plot0 ordiellipse (allmds3, gr0, kind = "sd", conf=0.95, label = t) 

however, despite biplot constructed on axes 2 , 3, ordiellipse function write ellipses identify clusters on axes 1 , 2 biplot. explain me why happens , how plot right ellipses on axes 2 , 3 biplot?

you use alternative plot function ordiplot():

require(vegan) data(dune) data(dune.env) nmds <- metamds(dune, k = 3) pl<-ordiplot(nmds, choices = c(2, 3), display = 'sites', type = 'n') points(pl, = 'sites', col = dune.env$management, pch = 16) ordiellipse(pl, dune.env$management)