0001 function correct_baseline_eeg(p)
0002
0003
0004
0005
0006
0007 disp(mfilename);
0008
0009
0010 twin_base = [-0.5 0];
0011 ch_eeg = 1:63;
0012
0013
0014 file = fullfile(p.proj_root, p.eeg_dirname, p.trial_dirname, ['r_' p.file_list{1} '.eeg.mat']);
0015 [~, ~, time_info] = vb_load_meg_data(file);
0016 [~, from] = min(abs(time_info.time-twin_base(1)));
0017 [~, to] = min(abs(time_info.time-twin_base(2)));
0018
0019 for file = 1:length(p.file_list)
0020
0021 old_file = fullfile(p.proj_root, p.eeg_dirname, p.trial_dirname, ['r_' p.file_list{file} '.eeg.mat']);
0022 load(old_file, 'eeg_data')
0023
0024
0025 for ch = ch_eeg
0026 for tr = 1:size(eeg_data, 3)
0027 eeg_data(ch, :, tr) = eeg_data(ch, :, tr)-mean(eeg_data(ch, from:to, tr), 2);
0028 end
0029 end
0030
0031
0032 new_file = fullfile(p.proj_root, p.eeg_dirname, p.trial_dirname, ['br_' p.file_list{file} '.eeg.mat']);
0033 copyfile(old_file, new_file);
0034 vb_save(new_file, 'eeg_data');
0035 end
0036