Selecting discriminant function models for predicting the expected richness of aquatic macroinvertebrates