Inverse: Check for bad-channel agreement in noise/data cov and record…

…ings

See: https://neuroimage.usc.edu/forums/t/45191/8

toolbox/process/functions/process_inverse_2018.m

@@ -498,6 +498,12 @@
             errMessage = [errMessage 'The noise covariance contains NaN values. Please re-calculate it after tagging correctly the bad channels in the recordings.' 10];
             break;
         end
+        % Check that bad channels in noise covariance are the same as bad channels in recordings
+        badChNoiseCov_goodChRecs = intersect(find(and(~any(NoiseCovMat.NoiseCov,1), ~any(NoiseCovMat.NoiseCov,2)')), GoodChannel);
+        if ~isempty(badChNoiseCov_goodChRecs)
+            errMessage = [errMessage 'Bad channels in noise covariance are different from bad channels in recordings.' 10 'Please re-calculate it after tagging correctly the bad channels in the recordings.' 10];
+            break;
+        end
 %         % Divide noise covariance by number of trials (DEPRECATED IN THIS VERSION)
 %         if ~isempty(nAvg) && (nAvg > 1)
 %             NoiseCovMat.NoiseCov = NoiseCovMat.NoiseCov ./ nAvg;
@@ -512,6 +518,11 @@
                 errMessage = [errMessage 'The data covariance contains NaN values. Please re-calculate it after tagging correctly the bad channels in the recordings.' 10];
                 break;
             end
+            badChDataCov_goodChRecs = intersect(find(and(~any(DataCovMat.NoiseCov,1), ~any(DataCovMat.NoiseCov,2)')), GoodChannel);
+            if ~isempty(badChDataCov_goodChRecs)
+                errMessage = [errMessage 'Bad channels in data covariance are different from bad channels in recordings.' 10 'Please re-calculate it after tagging correctly the bad channels in the recordings.' 10];
+                break;
+            end
 %             % Divide data covariance by number of trials
 %             if isempty(nAvg) && (nAvg > 1)
 %                 DataCovMat.NoiseCov = DataCovMat.NoiseCov ./ nAvg;