Question on statistical comparison between PSD of 5 subjects

💬 Support & Discussions
1

I want to compare the PSD in alpha band of 5 subjects. For that I need alpha band power for each subject to be a distribution. I am doing the following steps:

for sub, val in epochs_audi.items():
    psds, freqs = mne.time_frequency.psd_multitaper(val, fmin=0.0, fmax=4.0)
    psds = pd.DataFrame(np.log10(psds.mean(0)))

where epochs_audi is dictionary of epoched data with subject as key. I do not know how to proceed from there.

Alternatively, I also found the average PSD in the frequency band alpha (result - no of channels x 1) . Now I do not know what to do. Please suggest.

Thanks

2

Hello! Looks like you’re on a reasonable track, why can’t you just apply statistics (like ANOVA) from there?

1 Like