so far you start with already extracted PSD spectra + corresponding frequency bins. all functions work with 3d (trials x channels x frequencies), 2d (trials or channels x frequencies) and 1d data ...