R/Vegsoup-VisualisationMethods.R

In vegsoup: Classes and Methods for Phytosociology

plotPCO <- function (x, display = "sites", ...) {
	ord <- capscale(as.matrix(x) ~ 1, data = sites(x))	
	fig <- ordiplot(ord, display = display, type = "n") # , ...
	ol <- outlier(x, ...)
	points(fig, select = !ol, col = 1)
	points(fig, select = ol, col = 2, pch = "+")
	#orditorp(fig, select = ol, col = 2)	
}