fix for identifying if one or multiple durations are employed to describe the condition

Author: MartinPerez

nipype/algorithms/modelgen.py

@@ -462,6 +462,10 @@ def _concatenate_info(self, infolist):
         # names, onsets, durations, amplitudes, pmod, tmod, regressor_names,
         # regressors
         infoout = infolist[0]
+        for j, val in enumerate(infolist[0].durations):
+            if len(infolist[0].onsets[j]) > 1 and len(val) == 1:
+                infoout.durations[j] = (infolist[0].durations[j] *
+                                        len(infolist[0].onsets[j]))
         for i, info in enumerate(infolist[1:]):
             #info.[conditions, tmod] remain the same
             if info.onsets:
@@ -476,8 +480,15 @@ def _concatenate_info(self, infolist):
                                  sum(nscans[0:(i + 1)])
                         infoout.onsets[j].extend(onsets.tolist())
                 for j, val in enumerate(info.durations):
-                    if len(val):
+                    if len(info.onsets[j]) > 1 and len(val) == 1:
+                        infoout.durations[j].extend(info.durations[j] *
+                                                    len(info.onsets[j]))
+                    elif len(info.onsets[j]) == len(val):
                         infoout.durations[j].extend(info.durations[j])
+                    else:
+                        raise ValueError('Mismatch in number of onsets and \
+                                         durations for run {0}, condition \
+                                         {1}'.format(i+2, j+1))
                 if hasattr(info, 'amplitudes') and info.amplitudes:
                     for j, val in enumerate(info.amplitudes):
                         infoout.amplitudes[j].extend(info.amplitudes[j])

nipype/algorithms/tests/test_modelgen.py

@@ -79,6 +79,7 @@ def test_modelgen_spm_concat():
     yield assert_equal, np.sum(res.outputs.session_info[0]['regress'][0]['val']), 30
     yield assert_equal, len(res.outputs.session_info[0]['cond']), 1
     yield assert_almost_equal, np.array(res.outputs.session_info[0]['cond'][0]['onset']), np.array([2.0, 50.0, 100.0, 170.0, 210.0, 220.0, 280.0, 330.0])
+    yield assert_almost_equal, np.array(res.outputs.session_info[0]['cond'][0]['duration']), np.array([1., 1., 1., 1., 1., 1., 1., 1.])
     # Test case of scans as output units instead of seconds
     setattr(s.inputs, 'output_units', 'scans')
     yield assert_equal, s.inputs.output_units, 'scans'