Dear Forum,
I’m following this tutorial to compute induced power time course in different frequency bands. I wanted to compute separately for a number of ROIs that I created based on FreeSurfer parcellation (I have verified them in other analyses). As the API indicates, if I pass the label parameter of the following function as a list of multiple labels, I should get the estimated power for all labels (i.e. my ROIs):
stcs = source_band_induced_power(
epochs, inverse_operator1, bands, label=label, n_cycles=2, n_jobs=-1
)
However, this gives me the following error:
ValueError: Vertices must be ordered in increasing order.
What does this error mean?
For the moment I work around the problem by looping over every label and computing the average of all vertices in the label as, for example:
stcs[“theta”].data.mean(axis=0)
However, I’m afraid this is not the most efficient way and would like to understand what the error message means.
A more important question is that I noticed the tutorial code directly used epoch data, i.e. without subtracting evoked response, when calculating the induced power. Is this ideal? Does the function subtract the evoked response automatically (I couldn’t identify from the source code)? And as I understand from an earlier discussion, it is preferred to apply again baseline correction when computing the TFR from already baseline-corrected evoked response. Is this still true for calculating induced power?
Thank you in advance for your input.
Yaqing