discrimRuleCoef(m_model, m_dap, classifyRuleCoefList)

function [discrimRuleCoefList, groupToGroupDistances, numGroupToGroupDistances] = discrimRuleCoef(m_model, m_dap, classifyRuleCoefList)

% Computes the Discrimination Rule Coefficients, that is, the Discriminant
% Regions Sample Space Partition

groupToGroup = 2;
groupToGroupDistances = combnk(1:m_dap.constants.numGrp, groupToGroup);

numGroupToGroupDistances = size(groupToGroupDistances,1);
for i = 1:numGroupToGroupDistances
    discrimRuleCoefList{i}{1} = classifyRuleCoefList{groupToGroupDistances(i, 1)}{1} - classifyRuleCoefList{groupToGroupDistances(i, 2)}{1};
	discrimRuleCoefList{i}{2} = classifyRuleCoefList{groupToGroupDistances(i, 1)}{2} - classifyRuleCoefList{groupToGroupDistances(i, 2)}{2};
	if ~strcmp(m_model, 'linear')
        discrimRuleCoefList{i}{3} = classifyRuleCoefList{groupToGroupDistances(i, 1)}{3} - classifyRuleCoefList{groupToGroupDistances(i, 2)}{3};
	end
end