I am using morlet wavelets for my data. I am trying to find a way to adjust the code in such a way that I can visualize the power per frequency band [12-17Hz], [17-23Hz], [23-29 Hz]. Right now it seems that power is computed for every single frequency from 12 up to 29.
I am following this tutorial:
This is a piece from my code I am trying to adjust:
###determining the frequencies### freqs = np.arange(12, 29, 1) # full range freqsbeta1 = np.arange(12, 17, 1) # beta1 (12-17Hz) freqsbeta2 = np.arange(17, 23, 1) # beta2 (17-23Hz) freqsbeta3 = np.arange(23, 29, 1) # beta3 (23-29Hz) #different number of cycle per frequency n_cycles = freqs / 2. n_cycles1 = freqsbeta1 / 2. n_cycles2 = freqsbeta2 / 2. n_cycles3 = freqsbeta3 / 2. baseline = [-6.5, -5.5] # baseline interval (in s) baseline2 = [-7, -6] baseline3 = [-0.1, -0.01] powerprep1allbands, itc1 = tfr_morlet(concatpreparationb1, freqs=freqs, n_cycles=n_cycles, use_fft=True, return_itc=True, n_jobs=1)
And here is how I visualize it:
powerprep1allbands.plot(['C3'], tmin=-0.5, tmax=1.5, baseline=baseline, mode='percent', title=powerprep1allbands.ch_names)
Can someone tell me how I can get it seperated on the basis of the three frequency bands?