reject_trial_meg

 Reject trials in which blinks occurred

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

0001 function reject_trial_meg(p)
0002 % Reject trials in which blinks occurred
0003 %
0004 % Copyright (C) 2011, ATR All Rights Reserved.
0005 % License : New BSD License(see VBMEG_LICENSE.txt)
0006 
0007 disp(mfilename);
0008 
0009 % Set input file
0010 data_file = fullfile(p.proj_root, p.meg_dirname, p.trial_dirname, 'ri_all.info.mat');
0011 
0012 % Set output file
0013 save_file = fullfile(p.proj_root, p.meg_dirname, p.trial_dirname, 'rri_all.info.mat');
0014 
0015 % Load data
0016 load(data_file, 'fileinfo', 'Measurement')
0017 
0018 % Reject trials in which blinks occurred
0019 blink_trial = find(fileinfo.blink == 1);
0020 fileinfo.ActiveTrial(blink_trial) = 0;
0021 disp([num2str(length(blink_trial)) '/' num2str(fileinfo.Ntotal) ' trials are rejected.'])
0022 
0023 % Save result
0024 save(save_file, 'fileinfo', 'Measurement')
0025