%
% 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 the lungB data
% modified from Example 1.1

% First standardize the data such that each row has mean 0 and standard
% deviation 1.
load lungB; 
x = lungB;
[n,p]=size(x);         % [n=675 number of genes; p=number of patients]
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.
imagesc(y);
colormap(gray(256))
colorbar
title('Gene Expression for LungB')
xlabel('Patient')
ylabel('Gene')

saveas( gcf, '../../WriteUp/Graphics/Chapter1/prob_1_3', 'eps' );