% 
% Written by:
% -- 
% John L. Weatherwax                2005-08-12
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

Y1 = [ 6 0 8 0 -3;
       7 -1 12 8 1;
       10 4 6 2 -1;
       -1 6 14 4 -4;
       2 7 10 10 -2;
       2 0 -1 0 1;
       4 7 2.5 7 2;
       2.5 2 2.5 10 2;
       1 5 1  5 2;];
Y2 = [ 4 3 4 4 3;
       -1 1 15 10 -2;
       8 9 8 0 -2 ;
       1 8 2.5 5 -1;
       6 9 10 3 -3;];


M1 = size(Y1,1); M2 = size(Y2,1); 
M = M1+M2; 

% Pt (a):
% 
m1 = mean( Y1 ); var1 = var( Y1 ); 
m2 = mean( Y2 ); var2 = var( Y2 ); 

fisher = (m1 - m2).^2 ./ ( var1 + var2 )

% Plot histograms of each featuer:
% 
for fi=1:size(Y1,2),
  plotTruthFAHistograms( Y1(:,fi), Y2(:,fi), 10 );   
end