coefficent case study(liton)

clc;
clear all;
data= load('sediments.txt');
for i=1:10
sample(i,:)=['sample',sprintf('%02.of',i)];
end
clear i

mineral=['amp';'pyr';'pla';'ksp';'qtz';'cla';'flu';'sph';'gsl']
corrmatrix= corrcoef(data);
corrmatrix= flipud(corrmatrix);
imagesc(corrmatrix),colormap(hot)
title('correlation matrix');
axis square,colorbar,hold
set(gca,'xticklabel',mineral,'yTicklabel',flipud(mineral))

[pcs,newdata,variances]=princomp(data)
pcs(:,1:5)

subplot(2,2,1),plot(1:9,pcs(:,1),'0'),axis([1 9 -1 1])
text((1:9)+0.2,pcs(:,1),minerals,'fontsize',8),hold
plot(1:9,zeros(9,1),'r'),title('pc 1')

subplot(2,2,2),plot(1:9,pcs(:,2),'0'),axis([1 9 -1 1])
text((1:9)+0.2,pcs(:,2),minerals,'fontsize',8),hold
plot(1:9,zeros(9,1),'r'),title('pc 2')

subplot(2,2,3),plot(1:9,pcs(:,3),'0'),axis([1 9 -1 1])
text((1:9)+0.2,pcs(:,3),minerals,'fontsize',8),hold
plot(1:9,zeros(9,1),'r'),title('pc 3')

subplot(2,2,4),plot(1:9,pcs(:,4),'0'),axis([1 9 -1 1])
text((1:9)+0.2,pcs(:,4),minerals,'fontsize',8),hold
plot(1:9,zeros(9,1),'r'),title('pc 4')

plot(pcs(:,1),pcs(:,2),'o')
text((:,1)+0.2,pcs(:,2),minerals,'fontsize';14),hold
x=get(gca,'xlim');y=get(gca,'ylim');
plot(x,zeros(size((x)),'r')
plot(zeros(size((y)),y,'r')

plot(newdata(:,1),newdata(:,2),'+')
text(newdata(:,1)+0.01,newdata(:,2),sample),hold
x=get(gca,'xlim');y=get(gca,'ylim');
plot(x,zeros(size((x)),'r')
plot(x,zeros(size((y)),y,'r')