I am trying to run a regression model on my time-frequency EEG data
(AverageTFR or EpochsTFR). I am wondering if there is any way to have
"metadata" for time-frequency data to not lose the behavioral information
like reaction time for each trial.
However, I am using ITC (inter-trial-coherence) and not power data,
therefore I have to set "return_itc=True, average=True" in the
*tfr_morlet() *function and cannot be able to have my output as
"mne.time_frequency.EpochsTFR". For each subject/participant, I got
averaged over all trials/epochs. My final output is a numpy.array with a
shape "n_subjects X n_freqs X n_time X n_channels". I can save this
result as EpochsTFR type, however, I will lose all metadata information.