diff --git a/mne_nirs/statistics/_statsmodels.py b/mne_nirs/statistics/_statsmodels.py
index 826c7cdb0..19854f59f 100644
--- a/mne_nirs/statistics/_statsmodels.py
+++ b/mne_nirs/statistics/_statsmodels.py
@@ -112,7 +112,7 @@ def statsmodels_to_results(model, order=None):
             for j in range(i + 1):
                 sdf[jj, 0] = np.asarray(model.cov_re)[i, j]
                 sdf[jj, 1] = np.sqrt(model.scale) * \
-                             model.bse.iloc[model.k_fe + jj]
+                    model.bse.iloc[model.k_fe + jj]
                 jj += 1
 
         # Variance components