calc_power_itc_PAS.py

from my_settings import *
import mne
import sys
import numpy as np

from mne.minimum_norm import read_inverse_operator, source_induced_power

subject = sys.argv[1]

epochs = mne.read_epochs(epochs_folder + "%s_trial_start-epo.fif"
                         % subject)

inv = read_inverse_operator(mne_folder + "%s-inv.fif" % subject)

labels = mne.read_labels_from_annot(subject, parc='PALS_B12_Brodmann',
                                    regexp="Bro",
                                    subjects_dir=subjects_dir)
labels_selc = labels[6], labels[7]

frequencies = np.arange(8, 13, 1)  # define frequencies of interest
n_cycles = 4 # frequencies / 3.  # different number of cycle per frequency
method = "dSPM"

sides = ["left", "right"]
conditions = ["ctl", "ent"]

for label in labels_selc:
    for cond in conditions:
        for j, side in enumerate(sides):
            power, itc = source_induced_power(epochs[cond + "/" + side],
                                              inv,
                                              frequencies,
                                               label=label,
                                               method=method,
                                               pick_ori=None,
                                               use_fft=True,
                                               baseline=(-0.2, 0),
                                               baseline_mode='zscore',
                                               n_cycles=n_cycles,
                                               pca=True,
                                               n_jobs=1)
            np.save(tf_folder + "%s_pow_%s_%s_%s_%s_bs.npy" % (subject,
                                                                  cond,
                                                                  side,
                                                                  method,
                                                                  label.name),
                power)
            np.save(tf_folder + "%s_itc_%s_%s_%s_%s_bs.npy" % (subject,
                                                                  cond,
                                                                  side,
                                                                  method,
                                                                  label.name),
                itc)