kingjr

import mne
import numpy as np
from scipy.stats import norm
from matplotlib import pyplot as plt

# Load required sample data
data_path = mne.datasets.sample.data_path()
subjects_dir = data_path + '/subjects'
evoked = mne.read_evokeds(data_path + '/MEG/sample/sample_audvis-ave.fif')[0]
fwd = mne.read_forward_solution(data_path + '/MEG/sample/sample_audvis-meg-eeg-oct-6-fwd.fif', surf_ori=False, force_fixed=True)