//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// Reading the raw data:
raw = mne.io.read_raw_fif(raw_fname)
raw.info["bads"] = ["MEG 2443"] // 1 bad MEG channel
// Set picks
picks = mne.pick_types(raw.info, meg="grad", eeg=False, eog=False,
After Change
//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// We are interested in the beta band. Define a range of frequencies, using a
// log scale, from 12 to 30 Hz.
freqs = np.logspace(np.log10(12), np.log10(30), 9)
//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// Computing the cross-spectral density matrix for the beta frequency band, for
// different time intervals. We use a decim value of 20 to speed up the
