计算时代的平均值

Computing a mean of epochs

我正在使用 MNE Toolbox 在 Python 上进行 EEG 数据分析。我是 Python 新手,我想知道是否有办法计算纪元平均值? "epoch mean",我的意思是获取每个时期,并找到它给出的平均曲线。 (英语不是我的第一语言,所以我希望它是清楚的) 感谢您的帮助!

假设每个epoch总共有100个数据点,而你有这样的20个epoch。然后您可以将此数据重塑为 (20,100):20 行和 100 列。你想要的是每个时期的平均值。我假设您不需要滚动平均值(移动平均线:MA);如果您需要MA,请在评论区留言。

让我们制作一些虚拟数据并将此逻辑应用于它。

import numpy as np
import matplotlib.pyplot as plt

%matplotlib inline
%config InlineBackend.figure_format = 'svg' # 'svg', 'retina'
plt.style.use('seaborn-white')

x = np.random.randn(20*100) # assume this is your original data
x = x.reshape((20,100)) # now we reshape it: here each row is an epoch
epoch_average = x.mean(axis=1)

# Plot figure to show results
show_figure = True
if show_figure:
    fig, axs = plt.subplots(nrows=5, ncols=4, figsize=(12,15), sharey='row')
    for i, (ax, x_epoch) in enumerate(zip(axs.flatten(), x)):
        plt.sca(ax)
        plt.plot(np.arange(x.shape[1]), x_epoch, 'k-', label='epoch-{}'.format(i))
        plt.axhline(epoch_average[i], label='epoch-average', color='red', alpha=0.8, lw=2.0, ls='--')
        plt.legend()
        plt.title('Epoch-{} Average: {:.3f}'.format(str(i).zfill(2), epoch_average[i]))

    plt.tight_layout()
    plt.show()
    fig.savefig('output.png', dpi=300)