%
% Example epage: 28
% Problem epage: 40 
% 
% Written by:
% -- 
% John L. Weatherwax                2008-02-20
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

clear all; 
close all; drawnow; 
clc; 

addpath('../../Code/eda_toolbox');
addpath('../../Code/eda_data');

% This shows the normalized genes for ALL and AML. The data are standardized.
% similar to Example 1.1

% First standardize the data such that each row has mean 0 and standard
% deviation 1.
load leukemia
%x = leukemia(:,1:38);   % <- take the training set of patients ... 
x = leukemia(:,39:end);  % <- take the testing set of patients ... 
x = leukemia;            % <- take all of the patients ... 
[n,p]=size(x);
y = zeros(n,p);
for i = 1:n
  sig = std(x(i,:));
  mu = mean(x(i,:));
  y(i,:)= (x(i,:)-mu)/sig;
end

% Now do the image of the data.
pcolor(y)
colormap(gray(256))
colorbar
title('Gene Expression for Leukemia')
xlabel('Leukema Patient')
ylabel('Gene')

if( 0 ) saveas( gcf, '../../WriteUp/Graphics/Chapter1/prob_1_2', 'eps' ); end