# NOT RUN {
data(shoshveg) # returns vegetation data set
data(shoshsite) # returns site data
elev.clust <- as.numeric(factor(cut(shoshsite$elevation,5)))
# 5 elevation bands
elev.const <- const(shoshveg,elev.clust)
compare(elev.const,1,2) # identify diagnostic species
# }
Run the code above in your browser using DataLab