diff --git a/brainiak/funcalign/srm.py b/brainiak/funcalign/srm.py
index 023e71872..e436ad7f3 100644
--- a/brainiak/funcalign/srm.py
+++ b/brainiak/funcalign/srm.py
@@ -261,7 +261,8 @@ def transform(self, X, y=None):
         s = [None] * len(X)
         for subject in range(len(X)):
             if X[subject] is not None:
-                s[subject] = self.w_[subject].T.dot(X[subject])
+                s[subject] = self.w_[subject].T.dot(
+                    X[subject] - self.mu_[subject][:, np.newaxis])
 
         return s