'Compute power spectrum density (PSD) using a multi-taper method.
Parameters
epochs_data : array, shape (n_epochs, n_channels, n_times)
The data
y : None | array, shape (n_epochs,)
The label for each epoch.
If None not used. Defaults to None.
Returns
psd : array, shape (n_signals, len(freqs)) or (len(freqs),)
The computed PSD.'

def transform(self, epochs_data, y=None):

    if (not isinstance(epochs_data, np.ndarray)):
        raise ValueError(('epochs_data    should    be    of    type    ndarray    (got    %s).' % type(epochs_data)))
    (psd, _) = psd_array_multitaper(epochs_data, sfreq=self.sfreq, fmin=self.fmin, fmax=self.fmax, bandwidth=self.bandwidth, adaptive=self.adaptive, low_bias=self.low_bias, normalization=self.normalization, n_jobs=self.n_jobs)
    return psd
