A module to calculate phase-amplitude coupling in Python.
Pacpy is no longer actively supported, but it is in a stable state. For an actively maintained package, we recommend pactools. Note that these two packages may give different results when using the default filter parameters.
A Binder demo, complete with simulated data, can be found here.
pip install pacpy
Tested on Linux (Ubuntu 4.10), OS X (10.10.4), and Windows 9.
- numpy
- scipy
- pytest (optional)
That is , we assume Anaconda is installed.
The wrapper for MATLAB can be found at, https://github.com/voytekresearch/pacmat
An example of calculating PAC from two simulated voltage signals using the phase-locking value (PLV) method:
import numpy as np
from scipy.signal import hilbert
from pacpy.pac import plv
t = np.arange(0, 10, .001) # Define time array
lo = np.sin(t * 2 * np.pi * 6) # Create low frequency carrier
hi = np.sin(t * 2 * np.pi * 100) # Create modulated oscillation
hi[np.angle(hilbert(lo)) > -np.pi*.5] = 0 # Clip to 1/4 of cycle
plv(lo, hi, (4,8), (80,150)) # Calculate PAC0.99863308613553081