0001 function trialdata = reference_eeg_for_meg(trialdata,eeg_channels,mask,EEG_reference); 0002 0003 if strcmp(EEG_reference,'mastoids'); 0004 0005 temp = trialdata(:,eeg_channels) - 0.5*trialdata(:,eeg_channels(4))*ones(1,size(eeg_channels,2)); 0006 trialdata(:,eeg_channels) = temp; 0007 elseif strcmp(EEG_reference,'bipolar1'); 0008 % 0009 % bipolars are along one dimension in anterior-posterior directions 0010 % there will be other bipolars for other montages 0011 % 0012 temp = zeros(size(trialdata,1),size(eeg_channels,2)); 0013 temp(:,1) = trialdata(eeg_channels); 0014 % 0015 % to be finished later 0016 % 0017 0018 end; 0019 0020 0021