%
% Problem Epage 200 
% Example 5.3 Epage 133
% 
% 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; 
clc; 

addpath('../../Code/CSTool');
load tibetan;   % <- size(tibetan) = [32 skulls,5 features each]
%tibetan = randn(1,100); 

% sort our data (ascending order): 
tibetan_s = sort(tibetan,1); 

% get the equivalent data for the standard normal: 
ntib        = size(tibetan,1);
norm_data   = randn(1,ntib); 
norm_data_s = sort(norm_data); 

% plot the q-q plot for both sets of data: 
figure; plot( norm_data_s, tibetan_s, 'x' ); grid on; 
title( 'q-q plots for tibetan skull measurements' ); 
xlabel( 'Standard Normal Quantiles' ); ylabel( 'Skull Quantiles' ); 
saveas( gcf, 'prob_5_3_plt.eps', 'epsc' );