For group analysis, save evoked data in BIDS-derivatives:
t_obs, clusters, p_values, H0 = cluster_stats Use a configuration file (YAML) # config.yaml subjects: ['001', '002', '003'] task: 'visual' preprocessing: l_freq: 0.1 h_freq: 40 notch: 50 epochs: tmin: -0.2 tmax: 0.8 baseline: [-0.2, 0] Python script with argparse import yaml, argparse from mne_bids import BIDSPath, read_raw_bids def main(subject, config): # load config # run pipeline for one subject pass mne bids pipeline
Run in parallel:
# 4. Set average reference (EEG) if 'eeg' in raw: raw.set_eeg_reference('average', projection=False) For group analysis, save evoked data in BIDS-derivatives:
pip install mne mne-bids pybv from pathlib import Path import mne from mne_bids import BIDSPath, write_raw_bids, make_dataset_description Define your project root bids_root = Path('/path/to/your/bids_dataset') bids_root.mkdir(exist_ok=True) Create a dataset description (required for BIDS) make_dataset_description( path=bids_root, name="My MEG/EEG Study", authors=["Your Name", "Collaborator"], dataset_doi="", funding="Grant #", ) Define a subject and session subject_id = '001' session_id = '01' # optional task = 'visual' Convert a single raw file (e.g., BrainVision .vhdr) raw_path = Path('/raw_data/sub-001/session_1/eeg.vhdr') bids_path = BIDSPath( subject=subject_id, session=session_id, task=task, suffix='eeg', root=bids_root, ) Write to BIDS (copies and anonymizes) raw = mne.io.read_raw_brainvision(raw_path, preload=False) write_raw_bids( raw, bids_path, overwrite=False, verbose=True, ) For group analysis
# Read events from BIDS events, event_id = mne.events_from_annotations(raw_clean) selected_events = ['stimulus/face', 'stimulus/car'] event_id_sel = k: v for k, v in event_id.items() if k in selected_events Create epochs epochs = mne.Epochs( raw_clean, events, event_id=event_id_sel, tmin=-0.2, # 200 ms pre-stimulus tmax=0.8, # 800 ms post-stimulus baseline=(-0.2, 0), reject=dict(eeg=100e-6), # reject epochs with >100 µV peak-to-peak preload=True, ) Drop bad epochs automatically epochs.drop_bad() print(f"Retained len(epochs) / len(epochs.events) epochs") Step 5: Sensor-Level Analysis – Evoked Responses # Compute evoked for each condition evoked_face = epochs['stimulus/face'].average() evoked_car = epochs['stimulus/car'].average() Plot butterfly plot evoked_face.plot_joint(title='Face condition') Difference wave evoked_diff = mne.combine_evoked([evoked_face, evoked_car], weights=[1, -1]) evoked_diff.plot_joint(title='Face - Car')