Displaying 1 result from an estimated 1 matches for "comm4".
Did you mean:
comms
2011 Aug 10
1
Plotting Ellipses and Points of Matching Colors in an Ordination
...ctional code is below:
library(vegan)
comvars <- read.csv("commvars9trim2.csv")
com.mds <- metaMDS(comvars, trace=FALSE)
firefactors <- read.csv("commvars9factors.csv")
plot(com.mds, type = "n")
points(com.mds, dis="si", pch = as.numeric(firefactors$comm4), col =
as.numeric(firefactors$comm4))
This code will generate the ordination and plot it, using the factor levels
from firefactors$comm4 to determine the colors and characters to use in the
plot.
What I would like to do next is plot the ellipses denoting the 95%
confidence intervals of the group...