Documentation of plot_regressions

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

Cross-Reference Information

Listing of script plot_regressions

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs2/data
load heat_pcs.mat
ct = getct;
ct = detrend(ct);
lims = [90 310 -50 50];
[heat1, lat, lon, depth, middepth] = getheat(lims, 1:3, tim, 'temp');
heat2 = getheat(lims, 4:7, tim, 'temp');
[b1, a1] = butter(9, 2/18); [b2, a2] = butter(9, 2/40);
tseries = filtfilt(b1, a1, ct)-filtfilt(b2, a2, ct);
lags = [-10:1:1];
lags = [-1:1:10];
lags = -12:3:12;
[b1, a1] = butter(9, 2/20); [b2, a2] = butter(9, 2/30);
tseries = filtfilt(b1, a1, ct)-filtfilt(b2, a2, ct);
lags = [-10:1:1];
lags = [-1:1:10];
lags = -12:3:12;
tseries = bpcs1540(:,1);
lags = [-12:3:12];
tseries = lpcs15(:,1);
lags = [-12:3:12];
tseries = hpcs7(:,1);
lags = [-3:1:3];
tseries = bpcs412(:,1);
lags = [-4:1:4];
[reg1, c1] = regress_eof(heat1, tseries, lags);
[reg2, c2] = regress_eof(heat2, tseries, lags);
[tem, tem, tem, lm] = getll('temp', lims);
figure(1); fo(1); clf;
default_global; FRAME = [100 300 -45 45];
figure(1); fo(1); clf;
for i = 1:12;
  subplot(6,2,i);
    gcont(2e-7*reg1(i,:,:), 1.5, 0);
    dc2(lm, 0.7, 1000);
    title(['lag ' num2str(lags(i))]);
end
figure(2); fo(2); clf;
for i = 1:12;
  subplot(6,2,i);
    gcont(2e-7*reg2(i,:,:), 2.5, 0);
    dc2(lm, 0.7, 1000);
    title(['lag ' num2str(lags(i))]);
end
figure(1); fo(1); clf;
for i = 1:12;
  subplot(6,2,i);
    gcont(c1(i,:,:), 0.1, 0);
    dc2(lm, 0.7, 1000);
    title(['lag ' num2str(lags(i))]);
end
figure(2); fo(2); clf;
for i = 1:12;
  subplot(6,2,i);
    gcont(c2(i,:,:), 0.1, 0);
    dc2(lm, 0.7, 1000);
    title(['lag ' num2str(lags(i))]);
end