diff --git a/mne_connectivity/effective.py b/mne_connectivity/effective.py
index 6edd36cf..a7a9eff2 100644
--- a/mne_connectivity/effective.py
+++ b/mne_connectivity/effective.py
@@ -383,7 +383,7 @@ def phase_slope_index_time(data,
     bands = list(zip(np.asarray((fmin,)).ravel(), np.asarray((fmax,)).ravel()))
     n_bands = len(bands)
 
-    freq_dim = -2
+    freq_dim = -1 
 
     # allocate space for output
     out_shape = list(cohy.shape)