apply_ica_meg

 Apply ICA to MEG data

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

0001 function apply_ica_meg(p)
0002 % Apply ICA to 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 % Load MEG data
0010 data_file = fullfile(p.proj_root, p.meg_dirname, p.trial_dirname, 'all.info.mat');
0011 data = vb_load_meg_data(data_file);
0012 
0013 [Nch, Nsample, Ntr] = size(data);
0014 data2 = reshape(data, [Nch, Nsample*Ntr]);
0015 
0016 % Apply ICA to MEG data
0017 [w, s, v] = runica(data2);
0018 unmix = w*s;% IC = unmix x MEG
0019 mix = inv(unmix);% MEG = mix x IC
0020 
0021 % Save ICA result
0022 ica_dir = fullfile(p.proj_root, p.meg_dirname, p.ica_dirname);
0023 if exist(ica_dir, 'dir') ~= 7
0024     vb_mkdir(ica_dir);
0025 end
0026 output_file = fullfile(ica_dir, [mfilename '.mat']);
0027 save(output_file, 'w', 's', 'v', 'unmix', 'mix')
0028 
0029