%
%   read the strain data and compute the correlation matrix
%
load becker.dat
load bird.dat
load bormann_hammond.dat
load hackl.dat
load holt.dat
load kreemer.dat
load loveless_meade.dat
load mccaffrey.dat
load parsons.dat
load gpsgridder.dat
load shen.dat
load smith_konter.dat
load tape.dat
load tong.dat
load zeng.dat
load SEISM.dat
%
%   compute the variance of each
%
std(tape)		% 74 01
std(becker)		% 82 02
std(shen)		% 91 03
std(parsons)		% 92 04
std(mccaffrey)		% 93 05
std(bormann_hammond)	% 106 06
std(gpsgridder)		% 109 07
std(holt)		% 126 08
std(hackl)		% 137 09
std(zeng)		% 139 10
std(loveless_meade)	% 149 11
std(kreemer)		% 152 12
std(smith_konter)	% 175 13
std(tong)		% 178 14
std(bird)		% 234 15
std(SEISM)		% 86 16

%
%   load the A-matrix and compute the correlations
%
A=[tape,becker,shen,parsons,mccaffrey,bormann_hammond,gpsgridder,holt,hackl,zeng,loveless_meade,kreemer,smith_konter,tong,bird,SEISM];
R=corrcoef(A);
%
figure(1)
imagesc(R,[0.4,0.9]);colorbar;colormap(gray);axis('square')
%
%  load just the rough models
%
A2=[shen,mccaffrey,gpsgridder,holt,zeng,loveless_meade,kreemer,smith_konter,tong,bird,SEISM];
R2=corrcoef(A2);
figure(2)
imagesc(R2,[0.4,0.9]);colorbar;colormap(gray);axis('square')

