fmin, fmax = 0.5, 30
sfreq = raw.info[‘sfreq’] # the sampling frequency
tmin = 0.0 # exclude the baseline period
con = spectral_connectivity_epochs(epochs1, method=‘dpli’, mode=‘multitaper’, sfreq=sfreq, fmin=fmin, fmax=fmax, faverage=True, tmin=tmin, mt_adaptive=False, n_jobs=1)
Now, visualize the connectivity in 3D:
plot_sensors_connectivity(epochs.info,con.get_data(output=‘dense’)[:,:,0])