hello ....it would be helpful if someone share code for feature extraction using dpli in time frequency domain, how can we get connectivity measures for each band?

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])