In [1]:

    

cd ..
clc; clear all; close all;
addpath(genpath(pwd));
zones_Prec = get_zones_Prec_weekly();
OISST      = dlmread('data/OISST_19811101-20161116.dat');
uwnd       = dlmread('data/uwnd_WIND_CCMP_merge_OISST_weekly.dat');
vwnd       = dlmread('data/vwnd_WIND_CCMP_merge_OISST_weekly.dat');

% Intercept SST to make them time consistent
OISST = OISST(1:size(zones_Prec, 1), :);
uwnd  = uwnd(1:size(zones_Prec, 1), :);
vwnd  = vwnd(1:size(zones_Prec, 1), :);
x_train = [OISST uwnd vwnd];
y_train = zones_Prec;


    

In [2]:

    

whos


    

    




Name               Size                   Bytes  Class     Attributes

  OISST           1688x20800            280883200  double              
  uwnd            1688x20800            280883200  double              
  vwnd            1688x20800            280883200  double              
  x_train         1688x62400            842649600  double              
  y_train         1688x9                   121536  double              
  zones_Prec      1688x9                   121536  double

In [ ]:

    

% idx   hiddenLayerSize delay_weeks trainPerformance    valPerformance  testPerformance performance trainR  valR    testR   R   region1 region2 region3 region4 region5 region6 region7 region8 region9
% result = tdnn_train(X, Y, delay_times, hiddenLayerSize)
result = [];
for delay_times = 1:30
    for repeat = 1:4
        result_one = tdnn_train(x_train, y_train, delay_times, 30);
        result = [result; result_one];
    end
end
dlmwrite('result/20161218_1_result_repeat_4_delay_1-30_hid_30.dat', result, 'delimiter', '\t');


    

In [7]:

    

whos


    

    




Name                Size                   Bytes  Class     Attributes

  OISST            1688x20800            280883200  double              
  ans                 1x1                        8  double              
  delay_times         1x1                        8  double              
  repeat              1x1                        8  double              
  result             80x20                   12800  double              
  result_one          1x20                     160  double              
  uwnd             1688x20800            280883200  double              
  vwnd             1688x20800            280883200  double              
  x_train          1688x62400            842649600  double              
  y_train          1688x9                   121536  double              
  zones_Prec       1688x9                   121536  double

In [ ]:

    

% idx   hiddenLayerSize delay_weeks trainPerformance    valPerformance  testPerformance performance trainR  valR    testR   R   region1 region2 region3 region4 region5 region6 region7 region8 region9
% result = tdnn_train(X, Y, delay_times, hiddenLayerSize)
result = [];
for delay_times = 21:30
    for repeat = 1:10
        result_one = tdnn_train(x_train, y_train, delay_times, 30);
        result = [result; result_one];
    end
end
dlmwrite('result/20161220_1_result_repeat_10_delay_21-30_hid_30.dat', result, 'delimiter', '\t');


    

In [9]:

    

whos


    

    




Name                Size                   Bytes  Class     Attributes

  OISST            1688x20800            280883200  double              
  ans                 1x1                        8  double              
  delay_times         1x1                        8  double              
  repeat              1x1                        8  double              
  result              0x0                        0  double              
  result_one          1x20                     160  double              
  uwnd             1688x20800            280883200  double              
  vwnd             1688x20800            280883200  double              
  x_train          1688x62400            842649600  double              
  y_train          1688x9                   121536  double              
  zones_Prec       1688x9                   121536  double

In [10]:

    

dlmwrite('result/20161220_1_result_repeat_10_delay_21-30_hid_30.dat', result, 'delimiter', '\t');


    

In [14]:

    

result_one(2:3)


    

    




ans =

    30    20

In [ ]: