mne.time_frequency.psd_array_multitaper¶

mne.time_frequency.psd_array_multitaper(x, sfreq, fmin=0, fmax=inf, bandwidth=None, adaptive=False, low_bias=True, normalization=’length’, n_jobs=1, verbose=None)[source]¶

Compute power spectrum density (PSD) using a multi-taper method.

Parameters:

Parameters:	x : array, shape=(…, n_times) The data to compute PSD from. sfreq : float The sampling frequency. fmin : float The lower frequency of interest. fmax : float The upper frequency of interest. bandwidth : float The bandwidth of the multi taper windowing function in Hz. adaptive : bool Use adaptive weights to combine the tapered spectra into PSD (slow, use n_jobs >> 1 to speed up computation). low_bias : bool Only use tapers with more than 90% spectral concentration within bandwidth. normalization : str Either “full” or “length” (default). If “full”, the PSD will be normalized by the sampling rate as well as the length of the signal (as in nitime). n_jobs : int Number of parallel jobs to use (only used if adaptive=True). verbose : bool, str, int, or None If not None, override default verbose level (see `mne.verbose()` and Logging documentation for more).
Returns:	psds : ndarray, shape (…, n_freqs) or The power spectral densities. All dimensions up to the last will be the same as input. freqs : array The frequency points in Hz of the PSD.

x : array, shape=(…, n_times)

The data to compute PSD from.

sfreq : float

The sampling frequency.

fmin : float

The lower frequency of interest.

fmax : float

The upper frequency of interest.

bandwidth : float

The bandwidth of the multi taper windowing function in Hz.

adaptive : bool

Use adaptive weights to combine the tapered spectra into PSD (slow, use n_jobs >> 1 to speed up computation).

low_bias : bool

Only use tapers with more than 90% spectral concentration within bandwidth.

normalization : str

Either “full” or “length” (default). If “full”, the PSD will be normalized by the sampling rate as well as the length of the signal (as in nitime).

n_jobs : int

Number of parallel jobs to use (only used if adaptive=True).

verbose : bool, str, int, or None

If not None, override default verbose level (see mne.verbose() and Logging documentation for more).

Returns:

psds : ndarray, shape (…, n_freqs) or

The power spectral densities. All dimensions up to the last will be the same as input.

freqs : array

The frequency points in Hz of the PSD.