[FIX] do not compute subject level contrast when running dataset level (#1102)

* do not compute subject level contrast when running dataset level

* update change log

Author: Remi-Gau

CHANGELOG.md

@@ -24,7 +24,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ### Added
 
-* [DOC] support label of activations with all atlases [1100](https://github.com/cpp-lln-lab/bidspm/pull/1100) by [Remi-Gau](https://github.com/Remi-Gau)
+* [ENH] support label of activations with all atlases [1100](https://github.com/cpp-lln-lab/bidspm/pull/1100) by [Remi-Gau](https://github.com/Remi-Gau)
 
 ### Changed
 
@@ -38,7 +38,8 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ### Fixed
 
-* [INFRA] copy `RepetitionTime` in sidecar JSON after running smoothing in #1099 by @Remi-Gau
+* [FIX] do not compute subject level contrast when running dataset level #1102 by @Remi-Gau
+* [FIX] copy `RepetitionTime` in sidecar JSON after running smoothing in #1099 by @Remi-Gau
 
 ### Security
 

demos/openneuro/ds000001_aroma_run.m

@@ -47,3 +47,9 @@
        'participant_label', participant_label, ...
        'preproc_dir', preproc_dir, ...
        'model_file', model_file);
+
+%% dataset level
+bidspm(bids_dir, output_dir, 'dataset', ...
+       'action', 'results', ...
+       'preproc_dir', preproc_dir, ...
+       'model_file', model_file);

demos/openneuro/models/model-balloonanalogrisktaskAroma_smdl.json

@@ -85,6 +85,36 @@
       "DummyContrasts": {
         "Test": "t"
       }
+    },
+    {
+      "Level": "Dataset",
+      "Name": "dataset_level",
+      "GroupBy": [
+        "contrast"
+      ],
+      "Model": {
+        "X": [
+          1
+        ],
+        "Type": "glm",
+        "Software": {
+          "bidspm": {
+            "Results": [
+              {
+                "name": [
+                  "cash_demean",
+                  "control_pumps_demean",
+                  "explode_demean",
+                  "pumps_demean"
+                ]
+              }
+            ]
+          }
+        }
+      },
+      "DummyContrasts": {
+        "Test": "t"
+      }
     }
   ]
 }

src/cli/cliStats.m

@@ -88,7 +88,8 @@ function cliStats(varargin)
     end
 
     if results
-      bidsResults(opt, 'nodeName', nodeName);
+      bidsResults(opt, 'nodeName', nodeName, ...
+                  'analysisLevel', analysisLevel);
     end
 
   end

src/workflows/stats/bidsResults.m

@@ -143,12 +143,14 @@
 
   args.addRequired('opt', @isstruct);
   args.addParameter('nodeName', '');
+  args.addParameter('analysisLevel', '', @ischar);
 
   args.parse(varargin{:});
 
   opt =  args.Results.opt;
 
   nodeName =  args.Results.nodeName;
+  analysisLevel =  args.Results.analysisLevel;
 
   %%
   currentDirectory = pwd;
@@ -177,10 +179,15 @@
     keep = ismember(listNodeNames, nodeName);
     opt.results = opt.results(keep);
   end
+  listNodeLevels = returnlistNodeLevels(opt);
+  if ~isempty(analysisLevel)
+    keep = ismember(listNodeLevels, analysisLevel);
+    opt.results = opt.results(keep);
+    listNodeLevels = listNodeLevels(keep);
+  end
 
   % skip data indexing if we are only at the group level
   indexData = true;
-  listNodeLevels = returnlistNodeLevels(opt);
   if all(ismember(listNodeLevels, 'dataset'))
     indexData = false;
   end