******   Is the model performance critically dependent on the   ******
******   particular parameter values used in the simulations?   ******

%%%%%  June 17, 2003  --- Optimal learning params for cacheO30F10W20
 load('cache\cacheO30F10W20');global PLM_CACHE;PLM_CACHE=cacheO30F10W20;
 load sched;load HUMAN_ZCONGR;

 descr=zeros(9600,5,20);for k=1:10 d=block_from_cache(sched1,PLM_CACHE);descr(:,:,k)=cat(1,d{:});end;
 for k=11:20 d=block_from_cache(sched2,PLM_CACHE);descr(:,:,k)=cat(1,d{:});end

 for n=1:2 for k=1:5 Mdata(n,k).stim=descr(:,:,k+5*n-5);Mdata(n,k).chev=make_chev(9600);Mdata(n,k+5).stim=descr(:,:,k+5*n+5);Mdata(n,k+5).chev=make_chev(9600);end;end;clear n

 Sparams=PLM_search_params =
    lsfun_name: 'PLM_paramsearch'
    sumsq_name: 'PLM_sumsqfit2'
    model_name: 'PLM_Zcongr1'
        params: [1x1 struct]
     p2v_templ: {1x6 cell}
     v2p_templ: {1x6 cell}
        bounds: [2x6 double]
         optns: [1x1 struct]

 Sparams.params =
      rep_size: [7 5 1]
       max_act: 0.5000
      rep_gain: 0.8000
     rep_noise: 0.1000
     out_noise: 0.2000
      W_minmax: [-1 1]
        W_init: 0.1700
       W0_seed: [35x1 double]
    learn_rate: 0.0020
    runav_rate: 0.9800
     criterion: 0.0450
       trigger: 3.3        % set to 0 later
      blk_size: 300

 Sparams.p2v_templ'
    'VAL = PARAMS.learn_rate *100 ;'
    'VAL = PARAMS.out_noise ;'
    'VAL = PARAMS.rep_noise ;'
    'VAL = PARAMS.W_init ;'
    'VAL = PARAMS.rep_gain ;'
    'VAL = PARAMS.criterion *10 ;'

 Sparams.optns.MaxFunEvals=200;

 for k=1:2 disp(['Data set ',int2str(k),':  ',datestr(now)]);[Sres(k).opt_par,Sres(k).ssq,Sres(k).opt_X,Sres(k).det]=paramsearch(Mdata(k,:),Sparams);[Sres(k).opt_X,Sres(k).ssq],save sensitAnal;end
0.0500    0.2344    0.0068    0.3653    0.6027    0.7120   18.6036
0.1661    0.1000    0.1076    0.0954    0.5618    0.2026   10.9855

% Set TRIGGER=0 and re-run both searchs starting from the optimal point of the second
 par=vector2params(Sres(2).opt_X,Sparams.params,Sparams.v2p_templ);par.trigger=0;
 Sparams.optns.MaxFunEvals=120;
 Sparams.params=par;for k=1:2 disp(['Data set ',int2str(k),':  ',datestr(now)]);[Sres(k).opt_par,Sres(k).ssq,Sres(k).opt_X,Sres(k).det]=paramsearch(Mdata(k,:),Sparams);[Sres(k).opt_X,Sres(k).ssq],save sensitAnal;end

0.1804  0.1190  0.0506  0.1093  0.6027  0.2523    9.8411
0.0883  0.1157  0.0264  0.0866  0.8503  0.2594    9.5547


%%%  June 19, 2003  --- The same for O15 F15 W20
 load cacheO15F15W20; global PLM_CACHE;PLM_CACHE=cacheO15F15W20;clear cacheO15F15W20;

 Sparams.params
      rep_size: [7 5 1]
       max_act: 0.5000
      rep_gain: 0.8000
     rep_noise: 0.1000
     out_noise: 0.2000
      W_minmax: [-1 1]
        W_init: 0.1700
       W0_seed: [35x1 double]
    learn_rate: 0.0020
    runav_rate: 0.9800
     criterion: 0.0450
       trigger: 0
      blk_size: 300

 for k=1:2 disp(['Data set ',int2str(k),':  ',datestr(now)]);[Sres(k).opt_par,Sres(k).ssq,Sres(k).opt_X,Sres(k).det]=paramsearch(Mdata(k,:),Sparams);[Sres(k).opt_X,Sres(k).ssq],save sensitAnal;end
0.1897  0.2194  0.0846  0.1630  0.7831  0.4837    9.1338
0.1906  0.2128  0.0679  0.1644  0.6784  0.4293    9.0838

 SresO15F15W20=Sres;clear Sres PLM_CACHE k ans


%%%  June 20, 2003  --- The same for O15 F07 W20
 load cacheO15F07W20; global PLM_CACHE;PLM_CACHE=cacheO15F07W20;clear cacheO15F07W20;
 for k=1:2 disp(['Data set ',int2str(k),':  ',datestr(now)]);[Sres(k).opt_par,Sres(k).ssq,Sres(k).opt_X,Sres(k).det]=paramsearch(Mdata(k,:),Sparams);[Sres(k).opt_X,Sres(k).ssq],save sensitAnal;end
0.1545  0.1845  0.1382  0.1513  0.8300  0.3713    6.8585
0.1851  0.1930  0.0901  0.1715  0.7506  0.4183    6.7889

 SresO15F07W20=Sres;clear Sres PLM_CACHE k ans


%%%  June 22, 2003  --- The same for O20 F10 W20
 load cacheO20F10W20; global PLM_CACHE;PLM_CACHE=cacheO20F10W20;clear cacheO20F10W20;
 for k=1:2 disp(['Data set ',int2str(k),':  ',datestr(now)]);[Sres(k).opt_par,Sres(k).ssq,Sres(k).opt_X,Sres(k).det]=paramsearch(Mdata(k,:),Sparams);[Sres(k).opt_X,Sres(k).ssq],save sensitAnal;end
0.1648  0.1854  0.0765  0.1624  0.7114  0.3891    8.9559
0.1554  0.1711  0.1230  0.1300  0.7369  0.3566    8.7224

 SresO20F10W20=Sres;clear Sres PLM_CACHE k ans


%%%  June 28, 2003  --- The same for O10 F10 W20
 cd cache;load cacheO10F10W20;cd ..; global PLM_CACHE;PLM_CACHE=cacheO10F10W20;clear cacheO10F10W20;
 for k=1:2 disp(['Data set ',int2str(k),':  ',datestr(now)]);[Sres(k).opt_par,Sres(k).ssq,Sres(k).opt_X,Sres(k).det]=paramsearch(Mdata(k,:),Sparams);[Sres(k).opt_X,Sres(k).ssq],save sensitAnal;end
0.1691  0.2050  0.0905  0.1882  0.6798  0.4433    7.5012
0.1874  0.2110  0.1110  0.1786  0.7139  0.4653    7.5019

 SresO10F10W20=Sres;clear Sres PLM_CACHE k ans
