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()

参考文献

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です