correct_baseline_meg

 Correct baseline of MEG data

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

0001 function correct_baseline_meg(p)
0002 % Correct baseline of MEG data
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 parameters
0010 twin_base = [-0.5 0];% time window of baseline
0011 
0012 % Load time information
0013 file = fullfile(p.proj_root, p.meg_dirname, p.trial_dirname, ['r_' p.file_list{1} '.meg.mat']);
0014 [~, ~, time_info] = vb_load_meg_data(file);
0015 [~, from] = min(abs(time_info.time-twin_base(1)));
0016 [~, to] = min(abs(time_info.time-twin_base(2)));
0017 
0018 for file = 1:length(p.file_list)
0019     % Load data
0020     old_file = fullfile(p.proj_root, p.meg_dirname, p.trial_dirname, ['r_' p.file_list{file} '.meg.mat']);
0021     load(old_file, 'bexp')
0022     
0023     % Make average of pre-stimulus  MEG to 0
0024     for ch = 1:size(bexp, 1)
0025         for tr = 1:size(bexp, 3)
0026             bexp(ch, :, tr) = bexp(ch, :, tr)-mean(bexp(ch, from:to, tr), 2);
0027         end
0028     end
0029     
0030     % Save baseline-corrected MEG data
0031     new_file = fullfile(p.proj_root, p.meg_dirname, p.trial_dirname, ['br_' p.file_list{file} '.meg.mat']);
0032     copyfile(old_file, new_file);
0033     vb_save(new_file, 'bexp');
0034 end
0035