一种基于配对样本T检验的事件相关电位分析方法技术

一种基于配对样本T检验的事件相关电位分析方法：设计含两种刺激的脑电诱发实验，并利用脑电采集设备记录多个导联的头皮脑电信号，进行初步的预处理；提取两种刺激下的ERP信号；对两种刺激下的ERP信号进行配对样本T检验，确定具有显著差异的时段；计算两种刺激下的ERP信号在显著差异时段内的差异面积，并绘制脑地形图，确定差异脑区。本发明专利技术确定了两种刺激下的ERP显著差异时段，并绘制了基于ERP波形差异面积的脑地形图，进而得到在显著差异时段内的差异脑区分布。本发明专利技术对于信噪比较差、且单个成分并不明显的ERP研究具有重要意义，并且为ERP信号和自发脑电的剥离提供了新的思路。

【技术实现步骤摘要】
一种基于配对样本T检验的事件相关电位分析方法
本专利技术涉及一种事件相关电位分析方法。特别是涉及一种包含多种外界刺激、且单种刺激重复次数较少的基于配对样本T检验的事件相关电位分析方法。
技术介绍
脑电信号是通过电极记录下来的脑细胞群的自发性、节律性电活动，根据是否含有外部刺激，可分为自发脑电(Electroencephalo-graph,EEG)和事件相关电位(Event-relatedPotentials,ERP)两种。事件相关电位是人们对特定刺激事件进行感知加工或执行某种认知任务时诱发出来的一种脑电信号，常用于反映刺激发生前后脑电位的变化情况，与大脑注意资源分配、客体记忆、思维决策、认知加工等相关。由于ERP信号具有毫秒级的时间分辨率、良好的非侵入性，并且采集设备操作较简单，该信号在脑功能研究和脑疾病预诊方面都有颇多应用。多年来，ERP研究的一个重大难点就是与自发脑电的剥离。研究显示，大脑无时无刻不在运转，即使在不给任何外界刺激的情况下，中枢神经系统亦存在着节律性、自发性放电现象，而外部事件诱发的ERP信号幅值远小于自发脑电，且通常被淹没在自发脑电中。由于自发脑电具有很大的个体差异性和随机性，因此不能形成一个固定的自发脑电模板，使得ERP信号便于剥离。实际过程中，通常采用多次重复施加外部刺激、再求平均的方式提高ERP信号的幅值和纯度，进而将其与自发脑电剥离。为了得到信噪比较好的ERP信号，通常需要几十甚至几百次的重复外界刺激，一方面，多次重复刺激必然会引起感官系统的疲劳，且难以保持完全一致的重复性；另一方面，大量刺激材料的准备并不容易，尤其是对于特定含义的图片、声音等较复杂的刺激。另外，以往的ERP分析多集中于某个或者某几个ERP成分(如P1，N1，P3等)的分析，然而对于重复刺激次数较少、信噪比欠佳的ERP信号，具有明确物理意义的ERP成分往往难以识别，也造成不同刺激下ERP特征提取的困难。基于配对样本T检验的事件相关电位分析方法从显著差异的角度分析ERP信号，能够提取出有显著意义的差异特征，避开了单个ERP成分提取的困难，是ERP对比分析的新思路。另外，由于两种不同外界刺激下的自发脑电虽不完全一致，但也不具有显著的差异性，若对两种刺激下的ERP信号进行配对样本T检验，得到的显著差异时段必为两真实ERP信号具有显著差异的时段，可以实现ERP信号与自发脑电的间接剥离。
技术实现思路
本专利技术所要解决的技术问题是，提供一种可用于信噪比较差、且单个成分并不明显的ERP研究的基于配对样本T检验的事件相关电位分析方法。本专利技术所采用的技术方案是：一种基于配对样本T检验的事件相关电位分析方法，包括如下步骤：1)设计含两种刺激的脑电诱发实验，并利用脑电采集设备记录多个导联的头皮脑电信号，进行初步的预处理；2)提取两种刺激下的ERP信号；3)对两种刺激下的ERP信号进行配对样本T检验，确定具有显著差异的时段；4)计算两种刺激下的ERP信号在显著差异时段内的差异面积，并绘制脑地形图，确定差异脑区。步骤1)所述的初步的预处理，是为去除头皮脑电信号记录过程中的低频漂移、高频干扰以及眼动电生理信号的干扰，对原始脑电信号进行变平均参考、0.5-10Hz带通滤波以及独立成分分析去眼电的预处理操作。步骤2)所述的提取两种以上刺激下的ERP信号，包括如下过程：(1)对初步预处理之后的整段头皮脑电信号进行分割，得到20个时长为4s的静息脑电片段和20个时长为1s的诱发脑电片段，其中，两种刺激对应的诱发脑电片段各10个；(2)选取刺激呈现之前的200ms，即静息期的后200ms为基准脑电，并计算基准脑电的平均幅值，将每个诱发脑电片段减去基准脑电平均幅值，实现去基线操作；(3)分别对去基线后的两种诱发脑电片段进行叠加平均，得到每位被试者的每个导联在两种刺激下的ERP信号，分别表示为X＝{Xijk}和Y＝{Yijk}，其中，i＝1,2,……,N1，N1＝15，N1是被试者数目；j＝1,2,……,N2，N2＝32，N2是导联数目，k＝1,2,……,N3，N3＝1024，N3是数据点数。步骤3)中所述的对两种刺激下的ERP信号进行配对样本T检验，是对每个导联中的每个数据点对应的ERP序列进行配对样本T检验，对于第j个导联第k个数据点，首先建立一个新变量Z＝{Zijk}，Zijk＝Xijk-Yijk，i＝1,2,……,N1，计算新变量的均值和方差构造统计量检验tjk是否服从自由度为N1-1的T分布，并计算出对应的显著程度Pjk，若Pjk<0.05，则序列和序列具有显著性差异，即两种刺激下的第j个导联第k个数据点具有显著性差异，否则，两种刺激下的第j个导联第k个数据点不具有显著性差异。步骤3)中所述的确定具有显著差异的时段，是在已确定两种刺激下的第j个导联第k个数据点是否具有显著性差异的基础上进行的，包含如下过程：(1)对于第j个导联的N3个数据点，若存在10个以上的连续数据点k，使得序列和序列具有显著性差异，那么这些连续的数据点k对应的时段就是第j个导联在两种刺激下的ERP信号具有显著差异的时段，若不存在10个以上的连续数据点k，使得序列和序列具有显著性差异，则第j个导联在两种刺激下的ERP信号不具有显著差异时段；(2)根据所有导联的显著差异时段分布，选择一个以上相对大的时段，使得尽可能包含多数导联的显著差异时段，所选的显著差异时段对应的数据点集标记为其中，r＝1,2,…,m，m为所选的显著差异时段数，可选范围为{1,2,…,100}，nr分别是每个显著差异时段所对应的数据点数，取值均大于10。步骤4)中所述的计算两种刺激下的ERP信号在显著差异时段内的差异面积，是对每位被试者的每个导联均计算两种刺激下的ERP信号在所选显著差异时段内的差异面积，其中，在显著差异时段内的两种刺激下的ERP信号的差异面积为并对所有被试者的数据进行叠加平均，得到第j个导联在所选显著差异时段内两种刺激下的ERP信号差异面积步骤4)中所述的绘制脑地形图、确定差异脑区，是对于每一个显著差异时段，均根据所有导联的两种刺激下的ERP信号差异面积绘制脑地形图，进而分析两种刺激所诱发ERP信号差异的空间分布状况，得到在显著差异时段内的主要激活脑区分布。本专利技术的一种基于配对样本T检验的事件相关电位分析方法，从配对样本T检验出发，确定了两种刺激下的ERP显著差异时段，并绘制了基于ERP波形差异面积的脑地形图，进而得到在显著差异时段内的差异脑区分布。本专利技术主要针对包含多种外界刺激、且单种刺激重复次数较少的诱发脑电研究，对于信噪比较差、且单个成分并不明显的ERP研究具有重要意义，并且为ERP信号和自发脑电的剥离提供了新的思路。附图说明图1是一种基于配对样本T检验的事件相关电位分析方法的流程图。具体实施方式下面结合实施例和附图对本专利技术的一种基于配对样本T检验的事件相关电位分析方法做出详细说明。本专利技术的一种基于配对样本T检验的事件相关电位(Event-relatedPotentials,ERP)分析方法，首先利用脑电采集设备记录两种刺激下的多导头皮脑电信号，并进行初步的预处理；其次提取两种刺激下的ERP信号；再次分别对每一导联在两种外界刺激下的ERP信号进行配对样本T检验，得到每本文档来自技高网...

【技术保护点】
一种基于配对样本T检验的事件相关电位分析方法，其特征在于，包括如下步骤：1)设计含两种刺激的脑电诱发实验，并利用脑电采集设备记录多个导联的头皮脑电信号，进行初步的预处理；2)提取两种刺激下的ERP信号；3)对两种刺激下的ERP信号进行配对样本T检验，确定具有显著差异的时段；4)计算两种刺激下的ERP信号在显著差异时段内的差异面积，并绘制脑地形图，确定差异脑区。

【技术特征摘要】
1.一种基于配对样本T检验的事件相关电位分析方法，其特征在于，包括如下步骤：1)设计含两种刺激的脑电诱发实验，并利用脑电采集设备记录多个导联的头皮脑电信号，进行初步的预处理；2)提取两种刺激下的ERP信号；所述的提取两种以上刺激下的ERP信号，包括如下过程：(1)对初步预处理之后的整段头皮脑电信号进行分割，得到20个时长为4s的静息脑电片段和20个时长为1s的诱发脑电片段，其中，两种刺激对应的诱发脑电片段各10个；(2)选取刺激呈现之前的200ms，即静息期的后200ms为基准脑电，并计算基准脑电的平均幅值，将每个诱发脑电片段减去基准脑电平均幅值，实现去基线操作；(3)分别对去基线后的两种诱发脑电片段进行叠加平均，得到每位被试者的每个导联在两种刺激下的ERP信号，分别表示为X＝{Xijk}和Y＝{Yijk}，其中，i＝1,2,……,N1，N1＝15，N1是被试者数目；j＝1,2,……,N2，N2＝32，N2是导联数目，k＝1,2,……,N3，N3＝1024，N3是数据点数；3)对两种刺激下的ERP信号进行配对样本T检验，确定具有显著差异的时段；所述的对两种刺激下的ERP信号进行配对样本T检验，是对每个导联中的每个数据点对应的ERP序列进行配对样本T检验，对于第j个导联第k个数据点，首先建立一个新变量Z＝{Zijk}，Zijk＝Xijk-Yijk，i＝1,2,……,N1，计算新变量的均值和方差构造统计量检验tjk是否服从自由度为N1-1的T分布，并计算出对应的显著程度Pjk，若Pjk<0.05，则序列和序列具有显著性差异，即两种刺激下的第j个导联第k个数据点具有显著性差异，否则，两种刺激下的第j个导联第k个数据点不具有显著性差异；所述的确定具有显著差异的时段，是在已确定两种刺激下的第j个导联第k个数据点是否具有显著性差异的基础上进行的，包含...

【专利技术属性】
技术研发人员：殷涛，王欣，刘志朋，
申请(专利权)人：中国医学科学院生物医学工程研究所，
类型：发明
国别省市：天津;12