文章来源于微信公众号(茗创科技),欢迎有兴趣的朋友搜索关注。

1、安装与预设

图片

实战数据

图片

图片

图片

图片

2、背景概念

DataSet:被试数据导入EEGLAB称之为DataSet,各个预处理步骤会更新DataSet

ERPsets:对各个Bin对应的分段进行叠加平均,生成的平均ERP波形对应于一个ERPset。ERPset可以存储在磁盘上,包含了各种条件下的ERP波形

Bins:一种条件(或条件的组合)的刺激对应于一个Bin

3、ERPLAB处理数据的步骤

(1)导入数据

eeglab

(2)定位电极

eeglab

(3)去除无用电极

eeglab

(4)滤波

erplab

(5)降采样率

eeglab

(6)创建EVENTLIST

erplab

(7)创建bin

erplab

(8)分段

erplab

(9)去坏段、替换坏导联

eeglab

(10)去伪迹成分

erplab

(11)去除伪迹较大的分段

erplab

(12)重参考

erplab

(13)计算平均ERPs

erplab

(14) Plot波形图、地形图

erplab

(15)测量ERP与导出数据

erplab

(1)导入数据(厂商的数据格式;EEGLAB自带格式;文本文档;.mat文件)

图片

2)定位电极(标准位置模板;自定义文件;手工查找位置)

图片

图片

图片

(3)去除无用电极

BP:HEO, VEO

Neuroscan: HEO, VEO, CB1, CB2, Trigger

图片

图片

 (4)滤波

图片

图片

图片

图片

(5)降采样率

图片

(6)创建EVENTLIST(是为了后面创建bin

图片

图片

图片

图片

必须注意的是:在创建EVENTLIST时,实验过程中的所有MARK都要写入,以防后面建BIN出现问题。

图片

图片

有两种方法可以查看刚刚建立的list(可以用于检查,不是必须要操作的步骤)

第一种:在MATLAB命令窗口输入“EEG.EVENTLIST”,会出现如下内容

图片

第二种:MATLAB的文本编辑器中直接打开创建好的文本文件(list.txt)。如下所示

图片

(7)创建bin

图片

打开的空白编辑器是这样的

图片

但是建bin需要注意的地方有:

1、必须是在英文字符状态下编写。

2、bin编号必须从1开始,并且编号是连续的,如从bin1开始,按bin1,bin2,bin3的顺序往下写,不要间断。

3、写bin分为3步:第1行写bin的编号;第二行写一句话,用来描述该刺激条件;第3行写条件锁时和反应限制,并且“.”代表以哪个Mark为零点。

4、一般会有3种锁时情况:限刺激锁时;限刺激锁时的反应时间;限反应锁时。接下来一一列举。

限刺激锁时:

图片

限刺激锁时的反应时间:

图片

限反应锁时:

图片

本手册中,示例数据的bin是这样写的

图片

写完之后保存为bin.txt文件或者bin.m格式的文件

图片

图片

图片

(8)分段和基线校正(基于bin进行分段)

图片

图片

图片

(9)去坏段替换坏导联

图片

a) 去坏段

图片

b) 插值坏导

图片

图片

图片

(10去伪迹成分(运行ICA;去伪迹成分) 

a) 运行ICA

图片

图片

图片

图片

否则ICA运行会终止!

图片

运行完ICA对数据进行保存

图片

b) 去伪迹成分(水平眼电、垂直眼电)

图片

图片

图片

图片

 把要删除的成分(水平眼电、垂直眼电等)打开,点击ACCEPT变成REJECT,再点OK。

图片

11去除伪迹较大的分段(峰-峰检验;平均波幅检验;阶跃伪迹)

a) 以峰-峰检验为例

图片

图片

图片

图片

图片

导出伪迹拒绝率

图片

b) 以平均波幅为例

图片

图片

c) 阶跃伪迹

用峰-峰检验标记伪迹后,再做阶跃伪迹检验,在这个实验例子中,检测HEOG通道(通道16)中的水平眼动。

图片

图片

可以在EventList中查看伪迹信息,先把EventList导出

图片

直接打开导出的文本文档,或者用Matlab文本编辑器打开文本文件,显示如下:

图片

伪迹拒绝(a_flags)这一列,每个被拒绝的试次在这列最右侧会看到标识1。峰-峰检验中有些events设置了标识4。其他events在阶跃函数中设置了标识3,还有些events同时设置了标识3和标识4,表明这两种类型的伪迹都被检测到了。这一信息对于行为分析中排除特定伪迹的试次是有用的。

(12重参考(双侧颞乳突参考;全脑平均参考)

a) 双侧颞乳突参考

图片

b) 全脑平均参考

图片

(13计算平均ERPs

a) 单被试平均

图片

图片

图片

图片

图片

b) 所有被试总平均

图片

图片

图片

接下来把总平均的数据重命名,进行保存即可

(14plot 波形图、地形图

a) plot ERP波形图

图片

图片

图片

图片

b) plot ERP地形图

图片

图片

(15测量ERP指标

从ERPLAB中选择Load existing ERPset打开ERPset数据文件;如果已经加载了,就在ERPset菜单栏中选择它。

图片

加载之后,在ERPLAB中选择ERP Measurement Tool

图片

图片

点击Viewer  

图片

图片

Create plot(导出的图形显示如下)

图片

整个ERP Measurement Tool 设置完成后,点击RUN,对数据进行保存

图片

Logo

开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!

更多推荐