计算一个年代的平均值

Question

我正在使用MNE工具箱对Python进行脑电数据分析。我是个Python新手，我想知道是否有一种方法来计算一个时代的平均值？所谓“时代”，我的意思是把每一个时代都取下来，然后找出它能给出的平均曲线。(英语不是我的第一语言，所以我希望它是明确的)谢谢你的帮助！

Answer 1

假设每个时代总共有100个数据点，而你有这样的20个数据点。然后，您可以将这些数据重组为(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)    

​

​