Compute the power spectral density (PSD) using multitapers.
Calculates spectral density for orthogonal tapers, then averages them together for each channel/epoch. See [1] for a description of the tapers and [2] for the general method.
Epochs
or Raw
or Evoked
The data for PSD calculation.
float
Min frequency of interest.
float
Max frequency of interest.
float
| None
Min time of interest.
float
| None
Max time of interest.
float
The bandwidth of the multi taper windowing function in Hz. The default value is a window half-bandwidth of 4.
Use adaptive weights to combine the tapered spectra into PSD (slow, use n_jobs >> 1 to speed up computation).
Only use tapers with more than 90% spectral concentration within bandwidth.
Normalization strategy. If “full”, the PSD will be normalized by the
sampling rate as well as the length of the signal (as in
Nitime). Default is 'length'
.
str
| list
| slice
| None
Channels to include. Slices and lists of integers will be interpreted as
channel indices. In lists, channel type strings (e.g., ['meg',
'eeg']
) will pick channels of those types, channel name strings (e.g.,
['MEG0111', 'MEG2623']
will pick the given channels. Can also be the
string values “all” to pick all channels, or “data” to pick data
channels. None (default) will pick good data channels (excluding reference
MEG channels). Note that channels in info['bads']
will be included if
their names or indices are explicitly provided.
Apply SSP projection vectors. If inst is ndarray this is not used.
int
| None
The number of jobs to run in parallel. If -1
, it is set
to the number of CPU cores. Requires the joblib
package.
None
(default) is a marker for ‘unset’ that will be interpreted
as n_jobs=1
(sequential execution) unless the call is performed under
a joblib.parallel_backend()
context manager that sets another
value for n_jobs
.
Whether to omit bad segments from the data before fitting. If True
(default), annotated segments whose description begins with 'bad'
are
omitted. If False
, no rejection based on annotations is performed.
Has no effect if inst
is not a mne.io.Raw
object.
str
| int
| None
Control verbosity of the logging output. If None
, use the default
verbosity level. See the logging documentation and
mne.verbose()
for details. Should only be passed as a keyword
argument.
Notes
New in version 0.12.0.
References
mne.time_frequency.psd_multitaper
#Frequency and time-frequency sensor analysis