fixed TCompCor (bug found while refactoring)

Author: Shoshana Berleant

nipype/algorithms/confounds.py

@@ -407,8 +407,7 @@ def _run_interface(self, runtime):
         # "For each voxel time series, the temporal standard deviation is
         # defined as the standard deviation of the time series after the removal
         # of low-frequency nuisance terms (e.g., linear and quadratic drift)."
-        imgseries = regress_poly(2, imgseries)
-        imgseries = regress_poly(0, imgseries, remove_mean=True, axis=1) #this is super weird
+        imgseries = regress_poly(2, imgseries, remove_mean=True)
 
         time_voxels = imgseries.T
         num_voxels = np.prod(time_voxels.shape[1:])

nipype/algorithms/tests/test_compcor.py

@@ -52,11 +52,11 @@ def test_compcor(self):
 
     def test_tcompcor(self):
         ccinterface = TCompCor(realigned_file=self.realigned_file, percentile_threshold=0.75)
-        self.run_cc(ccinterface, [['-0.1535587949', '-0.4318584065'],
-                                  ['0.4286265207', '0.7162580102'],
-                                  ['-0.6288093202', '0.0465452630'],
-                                  ['0.5859742580', '-0.5144309306'],
-                                  ['-0.2322326636', '0.1834860639']])
+        self.run_cc(ccinterface, [['-0.1114536190', '-0.4632908609'],
+                                  ['0.4566907310', '0.6983205193'],
+                                  ['-0.7132557407', '0.1340170559'],
+                                  ['0.5022537643', '-0.5098322262'],
+                                  ['-0.1342351356', '0.1407855119']])
 
     def test_tcompcor_with_default_percentile(self):
         ccinterface = TCompCor(realigned_file=self.realigned_file)