Warning when pspm_load1 ignores NaNs in condition statistics

Author: 4gwe

src/pspm_load1.m

@@ -269,6 +269,18 @@
     else
       for iCond = 1:numel(indata.(mdltype).condnames)
         condindx = strcmpi(indata.(mdltype).condnames{iCond}, indata.(mdltype).trlnames);
+
+        cond_stats = indata.(mdltype).stats(condindx, :);
+        nan_rows = find(any(isnan(cond_stats), 2));
+
+        if ~isempty(nan_rows)
+          warning('ID:invalid_condition_values', ['Condition ''%s'' in first-level model ''%s'' contains %d row(s) with NaN values. ', 'NaN values were omitted when computing condition statistics.'], ...
+            indata.(mdltype).condnames{iCond}, ...
+            indata.(mdltype).modelfile, ...
+            numel(nan_rows));
+        end
+
+
         data.stats(iCond, :) = mean(indata.(mdltype).stats(condindx, :), 1, 'omitnan');
       end
       data.names = indata.(mdltype).names;

test/pspm_load1_test.m

@@ -267,6 +267,19 @@ function test_action_con(this)
         end
       end
     end
+    function test_action_cond_invalid_values(this)
+        model.modelfile = 'test.mat';
+        model.modeltype = 'dcm';
+        model.modality = 'scr';
+        model.stats = [1 2; NaN 4; 3 4];
+        model.names = {'parameter 1', 'parameter 2'};
+        model.trlnames = {'CS+', 'CS+', 'CS-'};
+    
+        [sts, data] = this.verifyWarning( @() pspm_load1(model, 'cond'), 'ID:invalid_condition_values');
+    
+        this.verifyEqual(sts, 1);
+        this.verifyEqual(data.stats, [1 3; 3 4]);
+    end
     function test_action_all(this)
       for i=1:numel(this.modelfiles)
         f = this.modelfiles{i};