目录
一、为什么要对ERP进行时域分析
时域分析:时间序列
脑电数据
- 用EEG代表静息态的数据
- 用ERP代表任务态的数据
时域分析只针对ERP(event-related poteneials)的数据
时域分析分析什么:事件诱发的电位是随着时间如何变化的?
预处理之后的数据的组成:
- 信号(ERP):锁时锁相
- 噪声:随机
因此,可以对预处理后数据进行叠加平均,以便抵消噪声、凸显信号
二、ERP进行时域分析:叠加平均
- 数据形状:通道(channel)x 时间(time)x 分段(epoch)
- 叠加平均
- 叠加平均是相同通道下,不同epoch(一个trial一个epoch)在同一时刻的波形进行相加,取平均值。
- 叠加平均的质量取决于:
- 试次(trial)的数量:校正优于拒绝接受
- 试次的质量(由‘数据采集+预处理’控制)
- 叠加平均的水平:
- 个体水平:抵消试次间的差异
- 提值去做统计
- Peak
- Latency
- Mean(平均幅值)
- Aplitute(幅值)
- 提值去做统计
- 组水平:抵消个体差异
- 画图
- 波形图:展示时间信息。选通道,画某个通道随时间变化的ERP波形
- 地形图:展示空间信息。选时间,画在这个时间段内,某些通道的地形图。某些通道在该时间段内各自幅值的均值
- 波形图:展示时间信息。选通道,画某个通道随时间变化的ERP波形
- 画完波形图之后要找成分,定时间窗,去画地形图。
- 画图
- 个体水平:抵消试次间的差异
- 叠加平均是相同通道下,不同epoch(一个trial一个epoch)在同一时刻的波形进行相加,取平均值。
具体的数据形状间的变化
dir函数在筛选数据时,会按照升序排列,注意0的影响
一些看形状的函数
length()
size()
part1:计算组平均 ERP
Part2:绘制组平均ERP
Part3:绘制一定时间范围内的地形图
Part4:绘制一系列地形图
导出并保存生成的图片
组水平叠加平均 AND 个体水平的统计分析
个体水平的统计:t检验、 ANOVA
三、ERP分析结果的报告
1、组水平叠加平均的一般流程:
- 画图
- 找感兴趣的成分、定时间窗(用于个体水平的比较统计)
2、组水平:画出不同条件下的波形图
运行结果:Cz通道下,组件叠加平均后,四种条件的不同波形随时间变化比较
3、根据波形图确定感兴趣的时间段范围、画地形图
4、如何自动求出地形图的maplimits
5、个体水平的叠加平均,提取指标后用来做统计
一般假设驱动更好,因为有理论支撑,
6、个体水平的指标统计的一些知识
1)、ANOVA
方差分析(ANOVA)又称“变异数分析”或“F检验”,用于两个及两个以上样本均数差别的显著性检验。
2)、T检验:
在matlab里面查看文档说明
输出结果参数
一个例子说明T检验
t检验公式:
ps:
- P越小,不是说明实际差别越大,而是说越有理由拒绝H0 ,越有理由说明两者有差异,差别有无统计学意义和有无专业上的实际意义并不完全相同
3)t检验和ANOVA的代码
代码:
ANOVA