ica.decompose_epochs(epochs)
# plot components for one epoch of interest
# A distinct cardiac component should be visible
ica.plot_sources_epochs(epochs, epoch_idx=13, n_components=25)
###############################################################################
# Automatically find the ECG component using correlation with ECG signal
# As we don't have an ECG channel we use one that correlates a lot with heart
# beats: 'MEG 1531'. We can directly pass the name to the find_sources method.
# In our example, the find_sources method returns and array of correlation
# scores for each ICA source.
ecg_scores = ica.find_sources_epochs(epochs,
target='MEG 1531',
score_func='pearsonr')
# get maximum correlation index for ECG
ecg_source_idx = np.abs(ecg_scores).argmax()
print '#%i -- ICA component resembling the ECG' % ecg_source_idx
###############################################################################
# Automatically find the EOG component using correlation with EOG signal
# As we have an EOG channel, we can use it to detect the source.
eog_scores = ica.find_sources_epochs(epochs,
target='EOG 061',
score_func='pearsonr')
# fit sources from epochs or from raw (both works for epochs)
ica.decompose_epochs(epochs)
# plot components for one epoch of interest
# A distinct cardiac component should be visible
ica.plot_sources_epochs(epochs, epoch_idx=13, n_components=25)
###############################################################################
# Automatically find the ECG component using correlation with ECG signal
# As we don't have an ECG channel we use one that correlates a lot with heart
# beats: 'MEG 1531'. We can directly pass the name to the find_sources method.
# In our example, the find_sources method returns and array of correlation
# scores for each ICA source.
ecg_scores = ica.find_sources_epochs(epochs, target='MEG 1531',
score_func='pearsonr')
# get maximum correlation index for ECG
ecg_source_idx = np.abs(ecg_scores).argmax()
print '#%i -- ICA component resembling the ECG' % ecg_source_idx
###############################################################################
# Automatically find the EOG component using correlation with EOG signal
# As we have an EOG channel, we can use it to detect the source.
eog_scores = ica.find_sources_epochs(epochs, target='EOG 061',
score_func='pearsonr')
# get maximum correlation index for EOG
