% Correct baseline for each channel:
    for i = 1:size(eeg,1)
        baseline = mean(eeg(i,:),2);
        data_shifted(i,:) = eeg(i,:) - baseline;
    end
    
    % Rereference to average:
    data_avg = mean(data_shifted,1);
    
    for i = 1:size(eeg,1)
        data_rerefed(i,:) = data_shifted(i,:) - data_avg;
    end