%
% Problem Epage 201 
% Example 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 forearm;   % <- size(forearm) = [140,1]

% produce a q-q plot of this data to study/investigate if it came from a normal distribution 
n  = length(forearm); 
fs = sort(forearm);

norm_s = sort(randn(1,n));
figure; plot( norm_s, fs, '-x' ); grid on; 

% sphere this data and replot: 
fs = fs - mean(fs); 
fs = fs/std(fs); 
figure; plot( norm_s, fs, '-x' ); grid on; 
title( 'q-q plot of sphered forearm data' ); 
xlabel('quantiles of the standard normal'); ylabel('quantiles of the sphered forearm data'); 
saveas( gcf, 'forearm_q_q', 'epsc' );