forgot to regress pre-computing std

Author: Shoshana Berleant

nipype/algorithms/compcor.py

@@ -118,13 +118,16 @@ class TCompCor(CompCor):
 
     def _run_interface(self, runtime):
         imgseries = nb.load(self.inputs.realigned_file).get_data()
-        time_voxels = imgseries.T
-        num_voxels = np.prod(time_voxels.shape[1:])
 
         # From the paper:
         # "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 = imgseries - np.mean(imgseries, axis=1)[:, np.newaxis]
+
+        time_voxels = imgseries.T
+        num_voxels = np.prod(time_voxels.shape[1:])
 
         # "To construct the tSTD noise ROI, we sorted the voxels by their
         # temporal standard deviation ..."

nipype/algorithms/tests/test_compcor.py

@@ -61,9 +61,9 @@ def test_compcor_with_extra_regressors(self, mock_add_extras):
 
     def test_tcompcor(self):
         ccinterface = TCompCor(realigned_file=self.realigned_file)
-        self.run_cc(ccinterface, [['-0.2846272268'], ['0.7115680670'],
-                                  ['-0.6048328569'], ['0.2134704201'],
-                                  ['-0.0355784033']])
+        self.run_cc(ccinterface, [['-0.2500000000'], ['0.2500000000'],
+                                  ['-0.2500000000'], ['0.7500000000'],
+                                  ['-0.5000000000']])
 
     def test_tcompcor_with_percentile(self):
         ccinterface = TCompCor(realigned_file=self.realigned_file, percentile_threshold=0.2)