Is this code for calculating the absolute and relative power of each electrode correct

💬 Support & Discussions
1 
# Calculate the absolute and relative power of each electrode
psds, freqs = mne.time_frequency.psd_welch(epochs, fmin=0, fmax=30, tmin=-1, tmax=60,
                                           n_fft=1280, n_overlap=250, n_per_seg=1280, picks=range(64), proj=False,
                                           n_jobs=1,  average='mean', window='hamming',
                                           verbose=None)
psds_mean = psds.mean(0).mean(0)
total_power = simps(psds_mean, dx=freq_res)
# Define frequency bands
freq_bands = {'delta': [0.5, 4],
              'theta': [4, 8],
              'alpha': [8, 13],
              'beta': [13, 30],
              }

# Initialize dictionaries to store results for each channel
abs_power = {}
rel_power = {}
# Loop through each channel
for i, ch_name in enumerate(raw.ch_names[:64]):

    # Compute absolute power for all frequency bands
    abs_power[ch_name] = {}
    for band, freq_range in freq_bands.items():
        idx_band = np.logical_and(freqs >= freq_range[0], freqs <= freq_range[1])
        band_power = simps(psds[:, i, idx_band].mean(axis=1), dx=freq_res)
        abs_power[ch_name][band] = '{:.20f} '.format(band_power*1e12)

    # Compute relative power for all frequency bands
    rel_power[ch_name] = {}
    for band, freq_range in freq_bands.items():
        idx_band = np.logical_and(freqs >= freq_range[0], freqs <= freq_range[1])
        band_power = simps(psds[:, i, idx_band].mean(axis=1), dx=freq_res)
        rel_power[ch_name][band] = '{:.20f} '.format(band_power / total_power)
2

I don’t know whether the code for calculating the absolute power and relative power of each electrode is correct, because the calculated value is very small, such as 0.00000000034427678985, and I don’t know what the unit of calculated data is, can someone teach me, thank you very much