Pythonで脳を研究する 10
引き続きPython MNEで色々と試している.
今回はEEGのevoked応答プロットのサンプルコードを出してみた.
import mne
import matplotlib.pyplot as plt
from mne.datasets import sample
from pathlib import Path
# サンプルデータのパスを取得
data_path = Path(mne.datasets.sample.data_path())
evoked_file = data_path / "MEG" / "sample" / "sample_audvis-ave.fif"
# evokedオブジェクトを読み込み(baseline補正およびプロジェクター適用済み)
evokeds = mne.read_evokeds(evoked_file, baseline=(None, 0), proj=True, verbose=False)
# 条件名に合わせた辞書を作成
conds = ("auditory_left", "auditory_right", "visual_left", "visual_right")
evks = dict(zip(conds, evokeds))
# EEGチャネルのみでトポロジカルサブプロットを作成
mne.viz.plot_compare_evokeds(
evks,
picks="eeg", # EEGチャネルに限定
colors=dict(auditory_left=0, auditory_right=0, visual_left=1, visual_right=1),
linestyles=dict(auditory_left="solid", auditory_right="solid",
visual_left="dashed", visual_right="dashed"),
axes="topo", # トポロジカルな配置でプロット
styles=dict(
auditory_left=dict(linewidth=1),
auditory_right=dict(linewidth=1),
visual_left=dict(linewidth=1),
visual_right=dict(linewidth=1),
),
)
plt.show()

参考文献