No, it is not possible to orditkplot()
with ordispider()
, it simply doesn't know how to handle arbitrary plotting functions like this.
You don't indicate why you want to use ordispider()
to show your grouped sites in the DCA ordination? You don't need to join them to some centroid or similar just to indicate group membership. Instead you can use plotting symbols to differentiate the groups, e.g.
require("vegan")
data(dune)
data(dune.env)
mod <- decorana(dune)
plot(mod, display = "sites", type = "n")
## colour & shape according to Management
col <- c("red","orange","forestgreen","navy")
pch <- 1:4
## add the points
with(dune.env,
points(mod, display = "sites", col = col[Management],
pch = pch[Management]))
## add a legend
legend("topright",
legend = with(dune.env, levels(Management)),
col = col, pch = pch, title = "Management",
bty = "n")
Alternatively, I suppose you could plot without labels and add them later, perhaps using locator()
to identify clear regions of the plot in which to place the labels, for example:
plot(mod, display = "sites", type = "p")
with(dune.env, ordispider(mod, groups = Management, col = "red"))
## select 4 locations
coords <- locator(with(dune.env, length(levels(Management))))
## now you have to click on the plot where you want the labels
## automagically finishes after you click the 4th label in this case
## draw labels
text(coords, labels = with(dune.env, levels(Management)))