LFP recordings and feature extraction
Local field potentials were acquired at a sampling rate of 422 Hz using the Medtronic Activa PC + S systemAll LFP analyses were performed using custom-written scripts in Python (v.3.6) and Matlab (R2018b). LFP recorded within a session was divided into 10-second segments from which spectral power, coherence and phase-amplitude coupling (PAC) were estimated. Spectral power and magnitude-squared coherence were estimated using the Python library Nitime’sPAC was estimated using the PACtools (v.0.3.1) python library