job_plot_meg_search_peak

 Search peak


 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function job_plot_meg_search_peak(inst_id)
0002 %
0003 % Search peak
0004 %
0005 %
0006 % Copyright (C) 2011, ATR All Rights Reserved.
0007 % License : New BSD License(see VBMEG_LICENSE.txt)
0008 
0009 global vbmeg_inst;
0010 define_job_plot_meg_const;
0011 
0012 % Global variables
0013 n_trial = vbmeg_inst.plotmeg{inst_id}.n_trial;
0014 T = vbmeg_inst.plotmeg{inst_id}.T;
0015 ed4 = vbmeg_inst.plotmeg{inst_id}.ed4;
0016 twin_meg = vbmeg_inst.plotmeg{inst_id}.twin_meg;
0017 
0018 % Time window size
0019 Tsize = str2num(get(ed4,'String'));
0020 if Tsize<0, 
0021   Tsize=10; 
0022   set(ed4,'Value','10');
0023 end
0024 
0025 % Search peak
0026 ix = twin_meg(1):twin_meg(2);
0027 tmp = sum((vbmeg_inst.plotmeg{inst_id}.bexp(:,ix,n_trial)).^2,1);
0028 [tmp,ix] = max(tmp);
0029 ix = ix+twin_meg(1)-1;
0030 
0031 % Update time window
0032 vbmeg_inst.plotmeg{inst_id}.t(1) = T(ix)-Tsize;
0033 vbmeg_inst.plotmeg{inst_id}.t(2) = T(ix)+Tsize;
0034 
0035 % Update spatiotemporal pattern
0036 job_plot_meg_update_timecourse(inst_id,FORCE_UPDATE);
0037 job_plot_meg_update_spatialpattern(inst_id);