neurofractal / pacmeg Goto Github PK

I want to extend these scripts, perhaps with Fieldtrip integration... However this will require lots of fundamental changes to the code. To keep this repo "future-proof" for the Frontiers article, should I create a new repo.. or integrate with the current repo?

See:

• https://bitbucket.org/ramonmc/pop_pac/src/master/
• https://www.sciencedirect.com/science/article/pii/S1053811918320056#fn1

Dear sir
Thank your for your great PAC codes for analyzing phase amplitude coupling.
I have a question:
we observe the following error when we consider the frequency lower than 6 Hz in the Calc_MI function.
Could you please help us?

Error using filter_with_correction (line 51)
Calculated filter coefficients have poles on or outside the unit circle and will not be stable. Try a higher cutoff
frequency or a different type/order of filter.

Error in filter_with_correction (line 51)
error('Calculated filter coefficients have poles on or outside the unit circle and will not be stable. Try a higher
cutoff frequency or a different type/order of filter.');

Error in ft_preproc_bandpassfilter (line 286)
filt = filter_with_correction(B,A,dat,dir,usefftfilt);

Error in preproc (line 324)
if strcmp(cfg.bpfilter, 'yes'), dat = ft_preproc_bandpassfilter(dat, fsample, cfg.bpfreq, cfg.bpfiltord,
cfg.bpfilttype, cfg.bpfiltdir, cfg.bpinstabilityfix, cfg.bpfiltdf, cfg.bpfiltwintype, cfg.bpfiltdev, cfg.plotfiltresp,
cfg.usefftfilt); end

Error in ft_preprocessing (line 358)
[dataout.trial{i}, dataout.label, dataout.time{i}, cfg] = preproc(data.trial{i}, data.label, data.time{i}, cfg,
begpadding, endpadding);

Error in calc_MI (line 70)
[virtsens_phase] = ft_preprocessing(cfg, virtsens);

Error in simulated_PAC_analysis1 (line 90)
tort_PAC = calc_MI(VE_PAC,[0.001 10],[5 20],[30 70],'no','no','tort');

All the best
Reza

1. According to the calc_MI function (virtsens = MEG data (1 channel)), can I use continuous EEG data as input? and would you please explain this input (toi = times of interest in seconds)?

See: https://doi.org/10.1016/j.neuroimage.2017.07.051

Also see: https://github.com/brainstorm-tools/brainstorm3/blob/bc6ccdc8e93be58bece50394234c3712e7e4c0cd/toolbox/process/functions/process_pac_dynamic.m

Where a portion of the comodulogram would not allow PAC to be detected, instead of throwing up an error make these values equal to NaN.

Implement approaches for surrogate analysis

For trial data, the function should be able to accept the indices of time-points between which PAC should be estimated