Time-frequency beamforming using LCMV

Compute LCMV source power 1 in a grid of time-frequency windows and display results.

References

1

Dalal et al. Five-dimensional neuroimaging: Localization of the time-frequency dynamics of cortical activity. NeuroImage (2008) vol. 40 (4) pp. 1686-1700

# Author: Roman Goj <roman.goj@gmail.com>
#
# License: BSD (3-clause)

import mne
from mne import compute_covariance
from mne.datasets import sample
from mne.event import make_fixed_length_events
from mne.beamformer import tf_lcmv
from mne.viz import plot_source_spectrogram

print(__doc__)

data_path = sample.data_path()
raw_fname = data_path + '/MEG/sample/sample_audvis_raw.fif'
noise_fname = data_path + '/MEG/sample/ernoise_raw.fif'
event_fname = data_path + '/MEG/sample/sample_audvis_raw-eve.fif'
fname_fwd = data_path + '/MEG/sample/sample_audvis-meg-eeg-oct-6-fwd.fif'
subjects_dir = data_path + '/subjects'
label_name = 'Aud-lh'
fname_label = data_path + '/MEG/sample/labels/%s.label' % label_name

Out:


Read raw data, preload to allow filtering

raw = mne.io.read_raw_fif(raw_fname, preload=True)
raw.info['bads'] = ['MEG 2443']  # 1 bad MEG channel

# Pick a selection of magnetometer channels. A subset of all channels was used
# to speed up the example. For a solution based on all MEG channels use
# meg=True, selection=None and add grad=4000e-13 to the reject dictionary.
# We could do this with a "picks" argument to Epochs and the LCMV functions,
# but here we use raw.pick_types() to save memory.
left_temporal_channels = mne.read_selection('Left-temporal')
raw.pick_types(meg='mag', eeg=False, eog=False, stim=False, exclude='bads',
               selection=left_temporal_channels)
reject = dict(mag=4e-12)
# Re-normalize our empty-room projectors, which should be fine after
# subselection
raw.info.normalize_proj()

# Setting time limits for reading epochs. Note that tmin and tmax are set so
# that time-frequency beamforming will be performed for a wider range of time
# points than will later be displayed on the final spectrogram. This ensures
# that all time bins displayed represent an average of an equal number of time
# windows.
tmin, tmax = -0.5, 0.75  # s
tmin_plot, tmax_plot = -0.3, 0.5  # s

# Read epochs. Note that preload is set to False to enable tf_lcmv to read the
# underlying raw object.
# Filtering is then performed on raw data in tf_lcmv and the epochs
# parameters passed here are used to create epochs from filtered data. However,
# reading epochs without preloading means that bad epoch rejection is delayed
# until later. To perform bad epoch rejection based on the reject parameter
# passed here, run epochs.drop_bad(). This is done automatically in
# tf_lcmv to reject bad epochs based on unfiltered data.
event_id = 1
events = mne.read_events(event_fname)
epochs = mne.Epochs(raw, events, event_id, tmin, tmax, proj=True,
                    baseline=None, preload=False, reject=reject)

# Read empty room noise, preload to allow filtering, and pick subselection
raw_noise = mne.io.read_raw_fif(noise_fname, preload=True)
raw_noise.info['bads'] = ['MEG 2443']  # 1 bad MEG channel
raw_noise.pick_types(meg='mag', eeg=False, eog=False, stim=False,
                     exclude='bads', selection=left_temporal_channels)
raw_noise.info.normalize_proj()

# Create artificial events for empty room noise data
events_noise = make_fixed_length_events(raw_noise, event_id, duration=1.)
# Create an epochs object using preload=True to reject bad epochs based on
# unfiltered data
epochs_noise = mne.Epochs(raw_noise, events_noise, event_id, tmin, tmax,
                          proj=True, baseline=None,
                          preload=True, reject=reject)

# Make sure the number of noise epochs is the same as data epochs
epochs_noise = epochs_noise[:len(epochs.events)]

# Read forward operator
forward = mne.read_forward_solution(fname_fwd)

# Read label
label = mne.read_label(fname_label)

Out:

Opening raw data file /home/circleci/mne_data/MNE-sample-data/MEG/sample/sample_audvis_raw.fif...
    Read a total of 3 projection items:
        PCA-v1 (1 x 102)  idle
        PCA-v2 (1 x 102)  idle
        PCA-v3 (1 x 102)  idle
    Range : 25800 ... 192599 =     42.956 ...   320.670 secs
Ready.
Current compensation grade : 0
Reading 0 ... 166799  =      0.000 ...   277.714 secs...
72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Opening raw data file /home/circleci/mne_data/MNE-sample-data/MEG/sample/ernoise_raw.fif...
Isotrak not found
    Read a total of 3 projection items:
        PCA-v1 (1 x 102)  idle
        PCA-v2 (1 x 102)  idle
        PCA-v3 (1 x 102)  idle
    Range : 19800 ... 85867 =     32.966 ...   142.965 secs
Ready.
Current compensation grade : 0
Reading 0 ... 66067  =      0.000 ...   109.999 secs...
110 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 110 events and 751 original time points ...
1 bad epochs dropped
Reading forward solution from /home/circleci/mne_data/MNE-sample-data/MEG/sample/sample_audvis-meg-eeg-oct-6-fwd.fif...
    Reading a source space...
    Computing patch statistics...
    Patch information added...
    Distance information added...
    [done]
    Reading a source space...
    Computing patch statistics...
    Patch information added...
    Distance information added...
    [done]
    2 source spaces read
    Desired named matrix (kind = 3523) not available
    Read MEG forward solution (7498 sources, 306 channels, free orientations)
    Desired named matrix (kind = 3523) not available
    Read EEG forward solution (7498 sources, 60 channels, free orientations)
    MEG and EEG forward solutions combined
    Source spaces transformed to the forward solution coordinate frame

Time-frequency beamforming based on LCMV

# Setting frequency bins as in Dalal et al. 2008 (high gamma was subdivided)
freq_bins = [(4, 12), (12, 30), (30, 55), (65, 299)]  # Hz
win_lengths = [0.3, 0.2, 0.15, 0.1]  # s

# Setting the time step
tstep = 0.05

# Setting the whitened data covariance regularization parameter
data_reg = 0.001

# Subtract evoked response prior to computation?
subtract_evoked = False

# Calculating covariance from empty room noise. To use baseline data as noise
# substitute raw for raw_noise, epochs.events for epochs_noise.events, tmin for
# desired baseline length, and 0 for tmax_plot.
# Note, if using baseline data, the averaged evoked response in the baseline
# period should be flat.
noise_covs = []
for (l_freq, h_freq) in freq_bins:
    raw_band = raw_noise.copy()
    raw_band.filter(l_freq, h_freq, n_jobs=1, fir_design='firwin')
    epochs_band = mne.Epochs(raw_band, epochs_noise.events, event_id,
                             tmin=tmin_plot, tmax=tmax_plot, baseline=None,
                             proj=True)

    noise_cov = compute_covariance(epochs_band, method='shrunk', rank=None)
    noise_covs.append(noise_cov)
    del raw_band  # to save memory

# Computing LCMV solutions for time-frequency windows in a label in source
# space for faster computation, use label=None for full solution
stcs = tf_lcmv(epochs, forward, noise_covs, tmin, tmax, tstep, win_lengths,
               freq_bins=freq_bins, subtract_evoked=subtract_evoked,
               reg=data_reg, label=label, rank=None)

# Plotting source spectrogram for source with maximum activity.
# Note that tmin and tmax are set to display a time range that is smaller than
# the one for which beamforming estimates were calculated. This ensures that
# all time bins shown are a result of smoothing across an identical number of
# time windows.
plot_source_spectrogram(stcs, freq_bins, tmin=tmin_plot, tmax=tmax_plot,
                        source_index=None, colorbar=True)
../../_images/sphx_glr_plot_tf_lcmv_001.png

Out:

Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 4 - 12 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 4.00
- Lower transition bandwidth: 2.00 Hz (-6 dB cutoff frequency: 3.00 Hz)
- Upper passband edge: 12.00 Hz
- Upper transition bandwidth: 3.00 Hz (-6 dB cutoff frequency: 13.50 Hz)
- Filter length: 991 samples (1.650 sec)

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 481 original time points ...
0 bad epochs dropped
Computing data rank from raw with rank=None
    Using tolerance 1.1e-11 (2.2e-16 eps * 13 dim * 4e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using SHRUNK
Done.
Number of samples used : 34632
[done]
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 12 - 30 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 12.00
- Lower transition bandwidth: 3.00 Hz (-6 dB cutoff frequency: 10.50 Hz)
- Upper passband edge: 30.00 Hz
- Upper transition bandwidth: 7.50 Hz (-6 dB cutoff frequency: 33.75 Hz)
- Filter length: 661 samples (1.101 sec)

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 481 original time points ...
0 bad epochs dropped
Computing data rank from raw with rank=None
    Using tolerance 2e-11 (2.2e-16 eps * 13 dim * 7e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using SHRUNK
Done.
Number of samples used : 34632
[done]
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 30 - 55 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 30.00
- Lower transition bandwidth: 7.50 Hz (-6 dB cutoff frequency: 26.25 Hz)
- Upper passband edge: 55.00 Hz
- Upper transition bandwidth: 13.75 Hz (-6 dB cutoff frequency: 61.88 Hz)
- Filter length: 265 samples (0.441 sec)

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 481 original time points ...
0 bad epochs dropped
Computing data rank from raw with rank=None
    Using tolerance 1.9e-11 (2.2e-16 eps * 13 dim * 6.4e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using SHRUNK
Done.
Number of samples used : 34632
[done]
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 65 - 3e+02 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 65.00
- Lower transition bandwidth: 16.25 Hz (-6 dB cutoff frequency: 56.88 Hz)
- Upper passband edge: 299.00 Hz
- Upper transition bandwidth: 1.31 Hz (-6 dB cutoff frequency: 299.65 Hz)
- Filter length: 1517 samples (2.526 sec)

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 481 original time points ...
0 bad epochs dropped
Computing data rank from raw with rank=None
    Using tolerance 2.2e-11 (2.2e-16 eps * 13 dim * 7.5e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using SHRUNK
Done.
Number of samples used : 34632
[done]
Loading data for 72 events and 751 original time points ...
0 bad epochs dropped
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 4 - 12 Hz

IIR filter parameters
---------------------
Butterworth bandpass zero-phase (two-pass forward and reverse) non-causal filter:
- Filter order 16 (effective, after forward-backward)
- Cutoffs at 4.00, 12.00 Hz: -6.02, -6.02 dB

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 751 original time points ...
0 bad epochs dropped
Computing time-frequency LCMV beamformer for time window -500 to -200 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.4e-11 (2.2e-16 eps * 13 dim * 2.6e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -450 to -150 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.3e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -400 to -100 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.2e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -350 to -49 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7e-11 (2.2e-16 eps * 13 dim * 2.4e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -300 to 0 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 6.8e-11 (2.2e-16 eps * 13 dim * 2.4e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -250 to 49 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 6.9e-11 (2.2e-16 eps * 13 dim * 2.4e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -199 to 100 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.1e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -149 to 150 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.5e-11 (2.2e-16 eps * 13 dim * 2.6e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -99 to 200 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.6e-11 (2.2e-16 eps * 13 dim * 2.6e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -49 to 250 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.8e-11 (2.2e-16 eps * 13 dim * 2.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 0 to 300 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.9e-11 (2.2e-16 eps * 13 dim * 2.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 50 to 350 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.8e-11 (2.2e-16 eps * 13 dim * 2.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 100 to 400 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.7e-11 (2.2e-16 eps * 13 dim * 2.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 150 to 450 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.3e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 200 to 500 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.2e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 250 to 550 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.2e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 300 to 600 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.2e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 350 to 650 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.2e-11 (2.2e-16 eps * 13 dim * 2.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 400 to 700 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 7.7e-11 (2.2e-16 eps * 13 dim * 2.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 450 to 750 ms, in frequency range 4 to 12 Hz
Computing data rank from raw with rank=None
    Using tolerance 8.1e-11 (2.2e-16 eps * 13 dim * 2.8e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 13032
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 1.3e-18 (2.2e-16 eps * 13 dim * 0.00046  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 12 - 30 Hz

IIR filter parameters
---------------------
Butterworth bandpass zero-phase (two-pass forward and reverse) non-causal filter:
- Filter order 16 (effective, after forward-backward)
- Cutoffs at 12.00, 30.00 Hz: -6.02, -6.02 dB

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 751 original time points ...
0 bad epochs dropped
Computing time-frequency LCMV beamformer for time window -500 to -300 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.7e-11 (2.2e-16 eps * 13 dim * 2e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -450 to -250 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.5e-11 (2.2e-16 eps * 13 dim * 1.9e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -400 to -200 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.1e-11 (2.2e-16 eps * 13 dim * 1.8e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -350 to -149 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.8e-11 (2.2e-16 eps * 13 dim * 1.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -300 to -99 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.8e-11 (2.2e-16 eps * 13 dim * 1.6e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -250 to -49 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.6e-11 (2.2e-16 eps * 13 dim * 1.6e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -199 to 0 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.2e-11 (2.2e-16 eps * 13 dim * 1.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -149 to 50 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.4e-11 (2.2e-16 eps * 13 dim * 1.5e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -99 to 100 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.7e-11 (2.2e-16 eps * 13 dim * 1.6e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -49 to 150 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.2e-11 (2.2e-16 eps * 13 dim * 1.8e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 0 to 200 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.5e-11 (2.2e-16 eps * 13 dim * 1.9e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 50 to 250 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.5e-11 (2.2e-16 eps * 13 dim * 1.9e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 100 to 300 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.5e-11 (2.2e-16 eps * 13 dim * 1.9e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 150 to 350 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.1e-11 (2.2e-16 eps * 13 dim * 1.8e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 200 to 400 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.1e-11 (2.2e-16 eps * 13 dim * 1.8e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 250 to 450 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.2e-11 (2.2e-16 eps * 13 dim * 1.8e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 300 to 500 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5e-11 (2.2e-16 eps * 13 dim * 1.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 350 to 550 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5.1e-11 (2.2e-16 eps * 13 dim * 1.8e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 400 to 600 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5e-11 (2.2e-16 eps * 13 dim * 1.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 450 to 650 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.5e-11 (2.2e-16 eps * 13 dim * 1.6e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 500 to 700 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 4.8e-11 (2.2e-16 eps * 13 dim * 1.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 550 to 750 ms, in frequency range 12 to 30 Hz
Computing data rank from raw with rank=None
    Using tolerance 5e-11 (2.2e-16 eps * 13 dim * 1.7e+04  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 8712
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.1e-18 (2.2e-16 eps * 13 dim * 0.0014  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 30 - 55 Hz

IIR filter parameters
---------------------
Butterworth bandpass zero-phase (two-pass forward and reverse) non-causal filter:
- Filter order 16 (effective, after forward-backward)
- Cutoffs at 30.00, 55.00 Hz: -6.02, -6.02 dB

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 751 original time points ...
0 bad epochs dropped
Computing time-frequency LCMV beamformer for time window -500 to -350 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.4e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -450 to -300 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.4e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -400 to -250 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -350 to -199 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -300 to -150 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -250 to -100 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -199 to -49 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -149 to 0 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.4e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -99 to 50 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.4e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -49 to 100 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.5e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 0 to 150 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.6e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 50 to 200 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.6e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 100 to 250 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.6e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 150 to 300 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.5e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 200 to 350 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.3e-11 (2.2e-16 eps * 13 dim * 4.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 250 to 400 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 300 to 450 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 350 to 500 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 400 to 550 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.1e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 450 to 600 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.1e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 500 to 650 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.1e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 550 to 700 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 600 to 750 ms, in frequency range 30 to 55 Hz
Computing data rank from raw with rank=None
    Using tolerance 1.2e-11 (2.2e-16 eps * 13 dim * 4.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 6552
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 3.4e-18 (2.2e-16 eps * 13 dim * 0.0012  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 65 - 3e+02 Hz

IIR filter parameters
---------------------
Butterworth bandpass zero-phase (two-pass forward and reverse) non-causal filter:
- Filter order 16 (effective, after forward-backward)
- Cutoffs at 65.00, 299.00 Hz: -6.02, -6.02 dB

72 matching events found
No baseline correction applied
Not setting metadata
Created an SSP operator (subspace dimension = 3)
3 projection items activated
Loading data for 72 events and 751 original time points ...
0 bad epochs dropped
Computing time-frequency LCMV beamformer for time window -500 to -400 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.7e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -450 to -350 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.4e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -400 to -300 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.4e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -350 to -249 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.2e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -300 to -199 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.3e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -250 to -150 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.5e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -199 to -99 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.4e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -149 to -49 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.3e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -99 to 0 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.4e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window -49 to 50 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.4e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 0 to 100 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.4e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 50 to 150 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.3e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 100 to 200 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.3e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 150 to 250 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.5e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 200 to 300 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.3e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 250 to 350 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.2e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 300 to 400 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.2e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 350 to 450 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9e-12 (2.2e-16 eps * 13 dim * 3.1e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 400 to 500 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.1e-12 (2.2e-16 eps * 13 dim * 3.1e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 450 to 550 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.5e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 500 to 600 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.5e-12 (2.2e-16 eps * 13 dim * 3.3e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 550 to 650 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.3e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 600 to 700 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.2e-12 (2.2e-16 eps * 13 dim * 3.2e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]
Computing time-frequency LCMV beamformer for time window 650 to 750 ms, in frequency range 65 to 299 Hz
Computing data rank from raw with rank=None
    Using tolerance 9.1e-12 (2.2e-16 eps * 13 dim * 3.1e+03  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Created an SSP operator (subspace dimension = 3)
    Setting small MAG eigenvalues to zero (without PCA)
Reducing data rank from 13 -> 10
Estimating covariance using EMPIRICAL
Done.
Number of samples used : 4392
[done]
Computing inverse operator with 13 channels.
    13 out of 366 channels remain after picking
Selected 13 channels
Whitening the forward solution.
    Created an SSP operator (subspace dimension = 3)
Computing data rank from covariance with rank=None
    Using tolerance 4.7e-18 (2.2e-16 eps * 13 dim * 0.0016  max singular value)
    Estimated rank (mag): 10
    MAG: rank 10 computed from 13 data channels with 3 projectors
    Setting small MAG eigenvalues to zero (without PCA)
Creating the source covariance matrix
Adjusting source covariance matrix.
[done]

Total running time of the script: ( 0 minutes 22.297 seconds)

Estimated memory usage: 457 MB

Gallery generated by Sphinx-Gallery