How do I analyze phase-amplitude coupling (PAC)?

Share icon

Fri, 06/30/2023 - 07:37 am

Hi all.

I'm currently using Igor Pro v9.02 (64-bit) with the NeuroMatic v3.0a package. 

I want to analyze phase-amplitude coupling (PAC) with my raw wave data file. I'm expecting to visualize the results in a spectrogram (or sonogram), with each vertical and horizontal axis showing amplitude frequency (Hz) and phase frequency (Hz). I would also like to quantify the magnitude of PAC in arbitrary units. As far as I know, I am supposed to perform Hilbert transform with my data, but I am not quite aware of this signal processing method nor how to execute this function in Igor Pro. 

My raw wave is 1.5 second-long, where the specific interval of interest is 0.5 – 1.0 seconds (500 milliseconds).

I attached my current working file (.pxp), where my target for analysis is Graph4, and Table0 shows this wave's time points (in seconds). The raw wave trace recording is in the attached (.csv) file.

Could you please explain the GUI-oriented step-by-step procedures or provide the codes for this analysis to achieve the result as described above?

(For example, how to analyze PAC between 8 Hz and 40 Hz.)

Many thanks,
Sung-Jun.

 working file (target of analysis: Graph4) Raw wave trace (1st column: time points, 2nd column: recorded signal)

Hello Sung-Jun,

I'm not aware of an Igor package that computes PAC but that does not mean that one does not exist.  If you want to implement this yourself in Igor you would have to write procedures or execute instructions from the command line; a "GUI-oriented" approach is not likely here. 

If you would like to write your own procedures to compute PAC I recommend starting by choosing a paper that presents the theory of the method you would like to follow.  Once you break the method into computational steps you can get help on how to implement them in Igor.

AG

June 30, 2023 at 04:52 pm - Permalink