计算时代的平均值
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)
我正在使用 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)