bayesdemo2.m

BAYESDEMO1 demo how to display discriminat function for Bayes classifier.

% BAYESDEMO2 demo on Bayes classifiation and error estimation.

% Statistical Pattern Recognition Toolbox, Vojtech Franc, Vaclav Hlavac
% (c) Czech Technical University Prague, http://cmp.felk.cvut.cz
% Modifications:
% 9-June-2001, V.Franc 

echo on;

% load model (see unsudemo).
pause;  % press anykey
model = load('gaussc2b');

% generate randomly testing set (200 patterns) accroding to the model.
pause;  % press anykey
[Xtst,Itst]=nmix(model.MI, model.SIGMA, model.Pk, 200);

% display testing set and Bayes discriminat function
pause;  % press anykey
figure; hold on;
ppoints( Xtst, Itst);
pbayescln( model.MI, model.SIGMA, model.Pk, 0, 'k--');

% classify testing points
pause;  % press anykey
[Ipred]=bayescln(Xtst,model.MI,model.SIGMA,model.Pk);

% compute classification error on testing points
pause;  % press anykey
100*cerror(Ipred,Itst)

% estimate upper limit on mean classification error 
pause;  % press anykey
mi1 = model.MI(:,1); % 1st mean value vector
mi2 = model.MI(:,2); % 2nd mean value vector
c1 = acov( model.SIGMA, 1 ); % 1st covariance matrix
c2 = acov( model.SIGMA, 2 ); % 2nd covariance matrix
100*bhattach(mi1,mi2,c1,c2,model.Pk(1),model.Pk(2))

echo off;