linearRule(location, D, dataElem, aprioriProb, cvCoef)

function  [rule, constCoef, linCoef] = linearRule(location, D, dataElem, aprioriProb, cvCoef)

% Linear Discriminant Rule

% Johnson, R.A. & Wichern, D.W. Applied Multivariate Statistical Analysis,
% Prentice -Hall

constCoef = -0.5 * location * D * location';
linCoef = location * D;

% cvCoeff: Cross Validation Coefficient
rule = cvCoef * (linCoef * dataElem' + constCoef) + log(aprioriProb);

constCoef = constCoef + log(aprioriProb);