#code for my discriminant analysis
hab.lda <- lda(grp ~ ., data=hab_std)

hab.lda.values <- predict(hab.lda, hab_std)
hab.class <- predict(hab.lda)$class

#create a histogram of the discriminant function values
ldahist(data = hab.lda.values$x[,1], g=grp)

#create a scatterplot of the discriminant function values
plot(hab.lda.values$x[,1], type="n", xlim=c(0,30), ylab=c("LDA Axis 1"))
    text(hab.lda.values$x[,1], row.names(hab_std),  col=c(as.numeric(hab.class)+10))
abline(v=0, lty="dotted")
abline(h=0, lty="dotted")