Documentation of lagged_corr_pcs

Global Index (short | long) | Local contents | Local Index (short | long)

Cross-Reference Information

Listing of script lagged_corr_pcs

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load RAW_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = pcs;
load LP10_detrend_L1-7_yr101-1000.mat; pcs3 = pcs;
load HP10_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = pcs;
lags = -20:20;
nlag = length(lags);
a = NaN * ones(3, nlag);
for i = 1:nlag;
  a(1,i) = corr(pcs3(:,1), pcs3(:,1), lags(i));
  a(2,i) = corr(pcs3(:,2), pcs3(:,2), lags(i));
  a(3,i) = corr(pcs3(:,1), pcs3(:,2), lags(i));
end
figure(1); 
for i = 1:3
  subplot(3,1,i);
    bar(lags, a(i,:));
    axis([-20 20 -.25 1.1]);
    grid on;
    set(gca, 'XTick', -20:5:20, 'YTick', -1:.2:1);
end
subplot(3,1,1);
  title('Autocorrelation for LP10 PC1');
subplot(3,1,2);
  title('Autocorrelation for LP10 PC2');
subplot(3,1,3);
  title('Lagged Correlation between LP10 PC1 and PC2');
  xlabel('PC1 leads PC2      |     PC2 leads PC1');
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap1/Plots