I’m not exactly sure what you are trying to do, but you can compute correlations with numpy.corrcoef — NumPy v1.26 Manual, and cross-correlation with numpy.correlate — NumPy v1.26 Manual. These functions take two 1D NumPy arrays, which you can generate from your epochs with epochs.get_data(). You have to pick one channel, and you might have to flatten the resulting 2D array.