ERP实验设计和数据分析...ppt

ERP实验设计和数据分析,李晓庆中国科学院心理研究所2011-07-01,ERP实验设计中的若干问题如何从ERP数据中得出结论ERP成分(component)与ERP波形(waveform)ERP实验设计中的无关变量及应对策略刺激材料数量和呈现策略数据分析时的基线问题ERP数据分析中的注意事项ERP叠加平均的原理、及伪迹用ICA排除伪迹ERP的PCA分析法及适用范围F.C.DondersCentreforCognitiveNeuroimaging,ERP实验设计中的若干问题,好数据是无可替代的,好数据是无可替代的。相信伪迹提出可以将坏数据变成好数据是愚蠢的。伪迹排除只能偶尔去掉那些偶尔发生的含有伪迹的试次，使好数据变得更好。如果被试在对待特定事件进行反应时总是眨眼，或者有连续的高振幅活动，伪迹排除方法并不能给予补偿，也就是说，那些具有系统噪声和系统伪迹的数据，并不能通过伪迹排除的方法得到更大改善。(Hansen.J.C.,未发表的软件文档),ERP实验设计中的若干问题,(一)如何从ERP数据中得出结论,1.不基于特定ERP成分的推论,(Li,Yang,Hagoort,2008),ERP实验设计中的若干问题,例子，重读(合适vs.不合适)声调(合适vs.不合适),现在正是玫瑰盛开的季节。明天小明去买什么把房间装饰一下？a明天小秦去买花装饰房间（AB）b明天小秦去买花装饰房间（AB）c明天小秦去买画装饰房间（AB）d明天小秦去买画装饰房间（AB）,Lietal.,2008,(d-a)=(c-a)+(b-a):线性相加性说明accent和tone的语义加工存在一定的对应之处,ERP实验设计中的若干问题,(一)如何从ERP数据中得出结论,2.基于特定ERP成分的推论,ERP成分(Component):apartofwaveformwithacircumscribedscalpdistribution(alludingtotheunderlyingneuralconfiguration)andacircumscribedrelationshiptoexperimentalvariables(alludingtothecognitivefunctionservedbytheactivityofthisconfiguration)(Donchinetal.,1978).,N400:central-parietal;semanticintegration,信息结构如何影响语义加工的深度(Wangetal.,2009)焦点(焦点vs.非焦点)语义合适(合适vs.不合适),ERP实验设计中的若干问题,(二)ERP成分(component)与ERP波形(waveform),ERP波峰ERP成分不能从单一的ERP波形估计ERP成分的时间过程和峰潜伏期不要把实验效应与波形的变化进行简单对应波峰振幅的差异成分大小的差异波峰潜伏期的差异成分时程的变化,ERP实验设计中的若干问题,(二)ERP成分(component)与ERP波形(waveform),聚焦于特定的成分使用已成熟的实验操作聚焦于大的成分用差异波来分离成分使用成分独立的实验设计,实验设计减少ERP成分和ERP波形之间的关系不确定性,可数名词,语义相关(plate.cup)可数名词,语义不相关(sock.cup)不可数名词,语义相关(rain.water)不可数名词,语义不相关(garbage.water)(Vogel,Luck,Shapiro,1998),ERP实验设计中的若干问题,(三)ERP实验设计中的无关变量及应对策略,例子,实验目的：刺激可辨性(discriminability)对P3的影响材料：字母表（10%“Q”;90%随即选择的非“Q”;“Q”从不连续出现）任务：字母“Q”时按键条件:容易辨别(亮条件)vs.难于辨别(暗条件)无关变量：刺激的物理属性：“Q”与非“Q”的形状不同;与非“Q”相比,“Q”出现的频率高被试的反应方式：对靶“Q”反应，而对非靶不反应(受与运动有关的ERP的污染)刺激的呈现方式：“Q”从不连续出现(P3后效只影响非“Q”)刺激材料的数量：10%“Q”;90%非“Q”(峰振幅)刺激条件的操作：亮条件vs.暗条件-二次反应(早期成分；加工难度差异),ERP实验设计中的若干问题,(三)ERP实验设计中的无关变量及应对策略,应对策略：刺激的物理属性：使用不同的靶字母，在全体被试间进行平衡；10个字母，1个为靶-每个字母的出现概率都是10%被试的反应方式：对靶“Q”和非靶做不同的反应刺激的呈现方式：分析数据时，去除有靶领先的非靶刺激刺激材料的数量：平均波幅代替峰波幅；随机选择同等数量的非靶刺激刺激条件的操作：增加控制条件；亮刺激与暗刺激在同一区组内随机呈现；额外的数据分析(RTs快vs.RTs慢),ERP实验设计中的若干问题,(三)ERP实验设计中的无关变量及应对策略,1.刺激的物理属性：不同实验条件下，尽量使用物理属性相同的刺激材料；增加控制实验例子，重读(合适vs.不合适)韵律边界(合适vs.不合适),他撞倒谁的车子之后逃跑了？撞倒/肖明的车子（AB）撞倒/肖明的车子（AB）撞倒肖明的/车子（AB）撞倒肖明的/车子（AB）,他撞倒谁的车子之后逃跑了？（AB）他撞倒肖明的什么物品之后逃跑了？（AB）警察把撞倒谁的车子没收了？（AB）警察把撞倒肖明的什么物品没收了？（AB）撞倒/肖明的车子(Lietal.,2011),下划线:表示锁时的关键词;红色表示重读的位置;左侧:物理属性不一样;右侧:物理属性一样,AudiovisualsituationStandard随后(红色)的方形出现了(75%)Deviant1随后(红色)的方形出现了(12.5%)Deviant2随后(红色)的方形出现了(12.5%)AuditorysituationStandard随后红色的方形出现了Deviant随后红色的方形出现了(Lietal.,2009),(三)ERP实验设计中的无关变量及应对策略,1.刺激的物理属性：不同实验条件下，尽量使用物理属性相同的刺激材料；增加控制实验例子，重读语义与视觉背景整合的时间特点,下划线:表示锁时的关键词;红色表示重读的位置;上:正式实验;下:控制实验,重读与视觉背景一致,重读与视觉背景不一致,重读与视觉背景不一致,ERP实验设计中的若干问题,(三)ERP实验设计中的无关变量及应对策略,2.刺激材料的重复3.不同实验条件下，刺激材料的数量不同4.不同实验条件下，被试反应的有无不同；反应的难度不同5.尽量在同一区组内变化实验条件；而不是在不同区组间变化实验条件,ERP实验设计中的若干问题,(四)刺激材料数量和呈现策略,1.刺激材料的数量效应的大小(ERP成分;实验条件操纵)；数据的噪声水平30-60-largecomponent(P3,N400)150-200-medium-sizedcomponent(N2)400-800-smallcomponent(P1)(Luck,2004)2.刺激材料的数量和信噪比,ERP实验设计中的若干问题,(四)刺激材料数量和呈现策略,3.刺激材料-ERP成分重叠4.应对ERP成分重叠的策略a运用ISIjitter(抖动)jitter的范围尽量宽；b不呈现刺激的小概率事件(Pictonetal.,1974)c估计已采集数据的重叠从平均波形中减去估计的重叠ADJAR(adjacentresponse)(Woldorff,1993),ERP实验设计中的若干问题,(五)数据分析时的基线问题-实验设计,1.尽量使锁时点之前的刺激材料在不同实验条件中相同(Li,Chen,andYang.,2011)2.选择恰当的时间段作为基线，并在实验设计时使其在不同实验条件中相同(Nan,Knsche,Friederici,2009)3.利用差异波（把效应进行比较）(Lietal.,minorrevision),ERP实验设计中的若干问题,例子，重读(合适vs.不合适)韵律边界(合适vs.不合适),他撞倒谁的车子之后逃跑了？撞倒/肖明的车子（AB）撞倒/肖明的车子（AB）撞倒肖明的/车子（AB）撞倒肖明的/车子（AB）,他撞倒谁的车子之后逃跑了？（AB）他撞倒肖明的什么物品之后逃跑了？（AB）警察把撞倒谁的车子没收了？（AB）警察把撞倒肖明的什么物品没收了？（AB）撞倒/肖明的车子(Li,Chen,andYang,2011),下划线:表示锁时的关键词;红色表示重读的位置;左侧:基线不一样;右侧:基线一样,ERP实验设计中的若干问题,例，非音乐家对音乐小结边界的感知，及文化差异对其的影响,Baseline:-2000ms-800-300ms(Nan,Knsche,Friederici,2009),Researchdesign,Activesentences(lowworkingmemoryload)highsemanticrelatedlowsemanticrelated萌萌把“新娘”浅色的婚纱弄脏后感到很抱歉。萌萌把“队长”浅色的婚纱弄脏后感到很抱歉。萌萌把“新娘”浅色的婚纱弄脏后感到很抱歉。萌萌把“队长”浅色的婚纱弄脏后感到很抱歉。Passivesentences(highworkingmemoryload)highsemanticrelatedlowsemanticrelated萌萌说“新娘”被弄脏浅色的婚纱后叹了口气。萌萌说“队长”被弄脏浅色的婚纱后叹了口气。萌萌说“新娘”被弄脏浅色的婚纱后叹了口气。萌萌说“队长”被弄脏浅色的婚纱后叹了口气。,MemoryLoad(highvs.low)Accentuation(accentedvs.de-accented)SemanticRelatedness(highvs.low),(Lietal.,minorrevision),下划线:表示锁时的关键词;红色表示重读的位置;右侧减左侧得到差异波,ERP数据分析中的注意事项,（一）ERP叠加的平均的原理、及伪迹,1.ERP的平均的原理,ERP数据分析中的注意事项,（一）ERP叠加平均的原理、及伪迹,2.EEG数据中的伪迹,眨眼眼动(阅读,视觉搜索)波肌电jaw,neck,andfacemusclessquinting,eyebrowmovementssmiling,laughter出声阅读心电(e.g.,参考位于血管上)电极接触不良,剔除含有伪迹的刺激信噪比降低矫正较好的信噪比；剔出不正确反而更歪曲数据,ERP数据分析中的注意事项,（一）ERP叠加平均的原理、及伪迹,2.EEG数据中的伪迹-波,10Hz,100ms一个周期，parietal-occipital很难平均掉：波幅高与刺激材料同步影响数据的解释；特别是始潜时的分析控制策略：ISIjitter，100ms,ERP数据分析中的注意事项,（一）ERP叠加平均的原理、及伪迹,2.EEG数据中的伪迹-慢电位漂移,皮肤电位变化；电极接触不良剔除：计算EEG在所有trials上的斜率，并将那些斜率超出一定阈值的trial排除,ERP数据分析中的注意事项,（一）ERP叠加平均的原理、及伪迹,3.EEG数据中的伪迹-眼动,固定的电压梯度(正极指向眼球正面)；眼动指向的方面，头皮电压变得更正剔除眼动污染的刺激优点：数据更干净缺点：信噪比降低；不可避免会产生眼动；有的眼动难以侦查到；剔除方法错误眼动矫正回归分析；偶极子定位；ICA,ERP数据分析中的注意事项,3.EEG数据中的伪迹-眼动,ERP数据分析中的注意事项,（二）用ICA排除伪迹(Jungetal.,2000,2001),数据驱动data-driving利用数据的相关结构，线形方法和非线性方法相结合；N电极-N个独立的成分与特定伪迹有关的成分具有独特的特点利用逆运算(inversecomputation)-减去有伪迹有关的成分,眼动,肌电,波,坏电极,ERP数据分析中的注意事项,（二）用ICA排除伪迹-具体程序（1）,Neuroscan数据预处理Ocularartifactreduction、re-Reference、filter、epoch、baselinecorrection、Artfactrejection（-75+75）、Deleterejectredsweeps用matlab打开eeglabfile-setpath(包含eeglab程序所在的路径、eeg数据所在的路径)；直接在matlab中输入eeglab把.eeg数据转换为eeglab的数据(.set)file-import-fromNeuroscan.eegfilefile-savecurrentdatasetas(保存为.set数据),ERP数据分析中的注意事项,（二）用ICA排除伪迹-具体程序（2）,RunICAfile-loadexistingdatasetedit-channellocation(一路点击ok)tools-runICAfile-savecurrentdatasetas(保存包含ica信息的数据)用ICA去伪迹tools-rejectdatausingICA-rejectcomponentbymapstools-romovecompomentsfile-savecurrentdatasetas(保存包含去除伪迹后的数据),把.set数据转换为Neuroscan能读取的.avg数据1.files-export-DataandICAactivitytotextfiles对话窗的以下处划勾ExportERPaverageinsteadoftrialsTransposematrix(ele-rows)Exportchannellabes/componentnumbersExporttimevalues2.把导出的数据转换为.xls数据-去除第一行（电极点名称）和第一列（时间点名称）；把.xls数据转换为.dat数据把.dat数据转换为.avg数据.dat数据的打开方式选择Neuroscanchannel62(电极数)acquisitionrate500(采样率)XMimum-100（epoch时的时间最小值）point551（采样点数）,（二）用ICA排除伪迹-具体程序（3）,ERP数据分析中的注意事项,（三）ERP的PCA分析法及适用范围-具体程序,数据驱动data-driving利用数据的相关结构M个连续变量-N个成分，NM很难从ERP波形中确定每个ERP成分的代表性时间窗口,(Lietal.,2010),例子，重读(合适vs.不合适)韵律边界(合适vs.不合适)(Lietal.,2011),ERP数据分析中的注意事项,（三）ERP的PCA分析法及适用范围-基本原理,.avg数据的预处理打开.avg文件，然后另存为.dat文件；把.dat文件的扩展名改为.xls文件，删除基线上的数据、及电极名称；把.xls文件的扩展名改为.txt文件,ERP数据分析中的注意事项,（三）ERP的PCA分析法及适用范围-基本原理,用matlab运行rawdata_to_pca.m来合并数据1.输入是