EEGLAB数据分析：预处理与后续处理

来源：http://blog.sina.com.cn/s/blog_13171a73d0102v4zx.html

数据预处理主要包括数据导入、电极定位、电极返回、滤波、去除伪迹、重建参考、分段、叠加平均等步骤。只有通过预处理的数据，才能进行特征值提取以进一步进行方差分析等操做。EEGLAB对数据处理的优势无需赘述。不少状况下，后续被试的预处理与第一个被试的处理过程是同样的。在这种状况下就能够采用批处理的方式进行，以便减小处理数据的时间并保持数据间参数的一致性，有利于后续数据的统计分析。在此主要有两种方法。第一，最简单的方法就是采用EEGLAB里面的插件Batch。该插件无需编程知识，在第一个数据处理后，其余数据均可以根据对第一个数据的处理步骤依次进行，方便快捷。这里须要注意的是，第一个被试的处理必须经过GUI（界面）实现，而不能经过Script实现。第二，就是经过编写简单的程序进行批处理。两种方法都是数据处理的好方法。

方法一：经过Batch插件进行批处理

• 预处理

在批处理以前，咱们必须对第一个数据进行完整且适当的分析。由于后续全部被试的分析都是和第一个被试的分析同样的，所以第一个数据的分析是很是关键的，它决定了后续被试的数据处理是否得当以及可否保持数据间参数的一致性。这里提醒一点，就是ERP实验前，必须打好分段标记！！这是后续数据处理中一个基本且重要的方面。

1 在第一个被试预处理所有作完后，在Batch中保存，默认选中第一个后缀名为set。

2 将刚才处理的过程保存为一个历史过程文件，路径为存储当前数据的文件夹。

3 将刚才处理的数据在EEGLAB里面清除（切记）。

4 打开Batch选项，选择历史文件，就是刚才处理第一个数据是保存的文件，第二打开要处理的其余文件，直接点击就进行批处理。另外旁边有一个选项是并行处理，选择这个选项能够加快处理速度，可是对内存要求很大。

Batch批处理，第一个数据

1 数据导入；

2 将采样率降为250 Hz: (非必须步骤)；

3 导入电极位置；

4 声明原始参考点：能够在导入电极位置以后，在最后一个电极点后增长参考点，电极自动寻找，就生成原来的参考电极的位置FCz。同时，须要制定这个电极点为参考点，其余全部的电极点都以这点为参考(1:63)；

5 在重建参考，进行平均参考（加入FCz，去除两个眼电电极点）；而后进行双耳参考（TP9， TP10），两步以后，就生成了FCz的数据值且是以TP9和TP10为参考点的矩阵；

6 自动检测电极点，坏电极点采用用插值法进行差值处理。

7 对连续数据删除正式试次开始和中间间隔的数据，减小数据处理对内存的要求。

8 采用ASR方法对上述连续数据进行伪迹矫正，这里注意不使用删除电极功能（off便可），在ASR以前最好用自动检测电极点程序。

9 滤波：高通、低通分别进行；(注意高通滤波)

10 分段：选择须要分析的全部marker；

11 基线矫正(自动步骤)；

12 删除错误试次和波幅过大试次；这步骤能够在下面操做；

上述就是采用Batch进行批处理的第一步；而后依次进行第二到第四步，处理剩余的数据。所有被试的数据进行批处理完毕后，对每一个被试的数据进行以下操做：

13 删除特定试次：这里的主要是指波幅范围标准和极端值分布标准；

请注意有些步骤不能用批处理进行，好比坏电极的处理，这是属于个体差别

14 检查是否有坏导：

① 检查是否有坏导; 若是有坏导，进行第 ②步操做；

② 替换坏电极：采用命令的方法或者采用差值方法进行替换，如下是命令方法：

EEG.data(62,:,:)=(EEG.data(4,:,:)+EEG.data(9,:,:)+EEG.data(13,:,:)+EEG.data(40,:,:))/4  

 EEG.data(63,:) = mean(EEG.data([4,9,13,41],:),1);

• 后处理

插值法可用工具菜单的interploration命令进行。通过上述处理的步骤，只是初步进行了处理，主要是滤波，重建参考，替换坏点，删除变异试次等操做。尚未去除眼电(上述盲源分析方法有效，若是已经处理过，就略过这步)，所以还需进一步进行去除除眼电的操做ICA。方法同上述操做步骤同样，过程以下：

1 先进行第一个数据的操做，而后进行批处理。

2 最后对每一个被试的数据进行去除眼电成分，再次进行基线矫正，获得干净的数据。

3 从每一个被试的数据中提取各个条件的数据；创建dataset。

进过上述两步的操做后，全部被试的数据都已经预处理完毕，就能够对这些数据进行后期的特征值提取，主要是指数据统计分析等工做。

方法二：经过编程进行批处理

EEGLAB最方便之处在于，通过每步处理后，均可以经过脚本的历史命令查询出相应的脚本，这样的话就能够在这些脚本的基础上进行编程，从而减小数据处理时间。能够经过EEG.history命令也能够经过菜单将处理过的数据的脚本进行保存。编程的批处理能够分两步进行：

1 所有数据进行预处理（无ICA）；

2 把第一步进行处理过的数据进行ICA，目的是去除眼电、肌电和其余噪音，获得干净的数据。通过第1步的预处理后，统一对已经生成的文件进行ICA。经过ICA，能够分别对生成的文件去除眼电、肌电和其余噪音ICA，这里笔者根据以往的经验，只去除眼电的噪音，并且对数据依次处理。不推荐使用独立成分自动去除的程序（如AJUST等自动软件）。

方法三：经过编程进行批处理

第二种方法中，咱们首先对全部数据进行了批处理，而后经过ICA的方法去除伪迹。实际上，有不少种去除伪迹的方法，一种比较直接的方法就是对原始的数据进行去除。这里推荐使用ASR方法，即伪迹空间重建的方法。这种方法只能用在原始数据即没有分段前的数据进行处理。步骤包括导入数据、电极定位、返回参考（可选择）、重建参考、滤波、ASR自动校订、分段、去除Epoch以及叠加平均过程。若是原始数据较大，则能够在预处理前分别对全部的原始数据进行处理，删除在刺激前和实验中途休息的是大段数据，达到减小数据量的目的。

---

EEGLAB数据分析

• 时间/波形分析

EEGLAB中，有个独特的模块就是study模块，能够方便的进行预处理后的数据分析。可是目前来看，经过这种方法虽然可行，可是较为麻烦。所以，更简便的方法就是不利用study编写脚本程序的方法进行后续的数据分析。所以，能够充分利用matlab的函数进行波幅和潜伏期分析。须要注意的是，在提取波幅和潜伏期以前，必须创建每种条件下的全部被试的总平均文件，利用总平均文件才能够进行后续的波幅和潜伏期特征值的提取。利用总平均文件，能够方便的提取波幅和潜伏期，进行方差分析；另外，总平均文件的创建还有利于之后论文投稿中的波形图和地形图的准备，一石二鸟。

画出总平均波形图/画出总平均地形图

总平均波形图 要画出总平均波形图，必须首先创建总平均的文件。这里经过GUI操做的方式经过Grand-average选项，对每一种条件下全部被试的.set数据进行总平均，创建一个N个不一样条件总平均文件。须要注意的是，总平均后须要再一次进行电极点定位，这样才能进行后续操做，不然没法画出总平均地形图。而后经过脚本命令方式的进行画出总平均图。该命令能够根据输出的电极点的数量，以及投稿时的要求进行适当的修改，以输出符合要求的图片(分辨率和大小)。利用上述语句能够在九个子图里分别画出两种条件以及差别波的总平均图，优势是方便快捷，无需对图片进行额外的准备，方便后续的投稿和图表准备。

图片制做分为三步：

(1) 完成画图及相关设置(字体大小、线宽、图例大小也是正常尺寸)，

(2) 此时WindowStyle is 'docked'，要改成normal，有两种操做：在Figure properties——more properties中找到Windowstyle，而后用鼠标改成normal；或者直接用命令：set (gcf,'windowstyle','normal')

(3) 将图片保存成TIFF格式图片。

(4) 将图片拷贝到word里，并生成PDF文件，用Photoshop软件打开刚才生成的PDF文件，而后对之进行操做，设置最终的图片大小和分辨率。

总平均地形图 前提是首先利用Grand-average插件进行不一样条件的总平均，获得两个总平均文件。再利用相关的EEGLAB函数进行画图。如须要改变电极点，在上述命令中修改须要的点便可，方便快捷。利用上一步生成的总平均文件和总平均地形图的命令，能够画出所须要的总平均地形图。能够根据某些特定要求简单的改变某些参数的设置便可。该命令能够很轻易的调整要输出的总平均地形图的数量，地形图的时间点和是否画出总平均差别地形图等。利用这个命令能够容易的画出总平均地形图，与上述总平均波形图同样。

当总平均文件创建之后，若是不利用脚本提取波幅和潜伏期的值，则能够采用第二种方法。全部数据导入letswave中，很方便进行提取，而后导入SPSS，进一步进行方差分析（参见相关文献）。该软件的优势之一是，在进行传统的波幅和潜伏期分析后，还能够进一步进行时间-频率分析；并且这个软件是线性的，无需过多的编程知识，有利于初学者学习。

• 时间/频率分析

利用每一个被试的处理事后的文件，能够进行时间-频率分析。这里主要采用的是EEGLAB自带的newtimef()函数，该函数是进行时频分析的主要函数，这里咱们主要采用的是FFT函数，能够对感兴趣的较低的频率进行分析。

总平均时间频率图/ROI频率的时间频率图 在这里，须要辨别两个图，一个是总平均时间-频率图，第二个是感兴趣的频率的时间能量波形图。前者是指根据newtimef()进行参数设置后，计算每一个被试每种条件下的时频图，而后进行平均得出的。后者是根据本身研究的兴趣点，只计算和分析本身感兴趣的频带的能量分布图。要画出上述两个重要的图，首先必须创建总平均时间频率文件，这个文件是进行后续数据处理的基础，从中能够提取出本身感兴趣的频带的ERSP的值，进行进一步的ANOVA等操做。