Methods are provided for determining discriminant functions of minimum risk linear classification systems, wherein a discriminant function is represented by a geometric locus of a principal eigenaxis of a linear decision boundary. A geometric locus of a principal eigenaxis is determined by solving a system of fundamental locus equations of binary classification, subject to geometric and statistical conditions for a minimum risk linear classification system in statistical equilibrium. Feature vectors and machine learning algorithms are used to determine discriminant functions and ensembles of discriminant functions of minimum risk linear classification systems, wherein distributions of the feature vectors have similar covariance matrices, and wherein a discriminant function of a minimum risk linear classification system exhibits the minimum probability of error for classifying given collections of feature vectors and unknown feature vectors related to the collections.