使用sLORETA进行源定位

1、、溯源数据生成溯源数据生成的步骤如下：1.对数据进行处理之后，如下图所示，对平均过后的ERP数据导出。-i l ctZlh 越 Raw Data& T1 CT21O:右睜Raw Dn匕 SQ Edt CMannete 自 莎 New Reference6 j$ Ocuhr Corroction 1CA自 妙 Ew Data InspectorEJ ： FitersO 总X yahi自 Baseline Correction 自*妙 Aitifact Rejectton;:O Average-；PPRk DGtQEonr. t2n ltt212t2L3 t ：214Ct215 tt216L ,

2、tt217 Q Ct218 -I Ct219少 Generic Data Export 占妙 buyt2hi ；謂 Baseline Con-ecton g Artifact Rejecton6 ：坊 Avera oe2Peak Detect io n2 Generic Data ExportGeneric Data Export2在Export里选择Generic,如下图所示，将方框里的文件格式改成txt,其它选项保持默 认状态点击nextoGeneric Data ExportEnotSdededElxi:VWrteHeadrfite_ WrtcMitaRc印xxtFteSbpccrott

3、r F:細EPlacedJi=Htor/rfeeBasePe:-湎血-3:2-7:ea电血亦 Nartcerrfeftrrrfi Ti FonTd l;/Mr. .mk)C MfcnafMxrrt:；Rj 刪 HcMcrFkPgHarker ft fetwi!213 a mt;Nod Caned3在Da File Fomat里面选择输出用于溯源数据的格式，切记一定要选择Text Fomat这 个选项，若选择Binary Fomat,那么生成的数据将是一堆不能用于分析的乱码。然后点击n extoZq Fie Fomiat0閔約Forrnm: Toi FomNO VedonrodGeneric D

4、ata ExportLre MnenjComnon E oei-tcro for Fzr. nke and w： foma* ds!a: PC Forrn：FHun / Se Feed)O UNIX Fcm (U Fed)O b*c Ferret Caago Rrtum) FvexasCancel4出现下面的窗口，保持各个选择的默认状态，点击nextoGeneric Dat Export口 SAJ Qianrei 2s胪e 9 So Oau F4e I1 OgrwrtioWCm屋 Symbol C )WKOadmal Sy; Pd 7&g* 10 Coum wan Seoarace Value

5、 by SoaoeU S*pjrat Vais, by Cutfom Sapaolor匚Ua Costom Preoanrrwoaon35出现下面的窗口，选择你需要的电极，眼电就不用选了，之后点击FinishoGeneric Data Export匚Cham 曲A/ciabeChanrefc-Sslecied Channes:HEOGVKGSebi rfenoveC3C4CP2QBRL2FP1FP2FP2RO1O2OZP3P4叶FirishCcTd下面这些容就是我们将ERP数据导出这个过程中，后台记录的操作步骤。生成溯源数据之后, 查看你的数据格式是否为ASCII,这个很重要。* Generi

6、c Data Export *File name parameter： $h SnFile extensiom txtWrite header file： yesWrite marker file： noFormat： ASCIIOrientation： MULTIPLEXEDLine Delimiter： CRLF (PC style)Add channel names： noOverwrite default decimal symbol: noExport al 1 channels： noExport the following channels：C3 C4 CPz Cz F3 F4

7、FCz Fpl Fp2 Fpz Fz 01 02 Oz P3 P4P07 P08 POz PzExported to： F：溯源Ett210 Average, txt Header file： F：溯源Ett210 Average, vhdrCreated Using Component Version 2.1.0.12577.下面我们就在输出路径的文件夹中找到我们“tt210”格式的两个文件。一个是txt格式， 打开之后如图1所示。另一个是vhdr格式，里面包含了导出数据的一些信息，包括txt格 式中每个数据对应的电极等，如图2。例如：若选用的电极数为20个，采样率为lOOOHZ/s, 所分

8、析的时间轴（算上基线）为550ms,那么txt文档中打开的数据就是有20列，550行（每 一个电极下面所采集的550个点，若分析的时间轴是1000ms的话，那么每一个电极下面就 会有1000个点，以此类推。），图2显示的就是20个电极的电极顺序。电极的排列顺序一定 不能乱。一0 17SO一0 1854一0 2043-0. 2372一0 2842-0. 3426-0. 4076-0. 4726-0. 5305-0. 5744-0. 5988一0 6004一0. 5785-0. 5353-0. 4741-0. 3993-0. 3154-0. 2266-0. 1368-0. 0497O. 5318O

9、. 53650. 54210. 54810. 55370 55800 56050. 56120. 56060 5592O. 5572O. 55420. 54930. 54150. 52980. 51310 49060 46170 42610 3839O. 2478O. 22690. 21240. 20400. 20050. 20140 20630. 21540 22900 2465O. 2665O. 28690. 30520. 31910. 32670. 32710. 31980. 30470 28250 25370. 29700. 28860. 28500. 28610. 29110 299

10、20 30950. 32110 33350. 34610. 35830. 36920. 37810. 38430. 38730 38680. 38270. 37510 36430 35060. 13990. 12670. 11690. 11050. 10780. 10970. 11770. 13330. 15810. 19210. 23440. 28210. 33160. 37910. 42060. 45330. 47510. 48510. 48340. 47050. 06470. 05070. 04450. 04620. 05570. 07200 09380 11990. 14910 179

11、90. 21080. 23970. 26450 28330. 29450. 29760 29260. 28040. 26270. 2411一O. 051 亠 一0 065 一0 073 一0. 075 一0. 072 一0 066 一0 057 一0 047 一0 035 一0 022 -O. 009 O. 0034 0. 0152 0. 0256 0. 0340 0 0402 0 0442 0. 0463 0 0472 0 04750. 0312O.3355O. 21950. 33460.44740. 21770.04860. 10300.28220. 18160. 31660. 41580

12、. 1945O. 05210. 16400.22570. 14230. 29740. 37750. 17220. 0590v;, , Futuro oxtcnainj .; Fl ol da orc dolinixod by coxnaa, soxiQ fields; rdght bo ril ttod .p ty.; Conasa In cbamG nas ara coded as 】Chl=Fpl,阳Ch2=Fz,册Ch3=F3,哦Ch4=Cz, , , .VCh5=C3 .Ch6=Pz. . .Ch7=P3. . .ChB=PO7.盼Ch9=Oz .Chi0=01,.碌节Chll=Fpz

13、.Chl2=Fp2.Chl3=F4.盼Chl4=C4.碌节ChlB=CPz. . .ChlO=P4.碌节Chl7=POz. . . 19.VChl8=PO8. . .Chl9=O2.碌Ch20=FCz, r . 借助sloreta对差异区进行定位第一部分：准备工作两个文件夹和一个txt文档。每个文件夹包含了一种条件下被试的数据，数据的格式是txto 一个txt文档是编写好的电极顺序，电极顺序必须与导出的ERP数据电极顺序一致为了确保 一致，可以对照vhdr文件中的电极信息顺序，将其写在txt文档中，用于分析，如图所示。1生成坐标电极。打开sloreta,点击主莱单Main Utilities然

14、后选择Electrode names to coordinates将含有电极顺序的txt文本拖入第一个空格，转换为以sxyz为后缀的电极 矩阵，这个肘候有三个坐标值（X-Y-Z）,是相应电极所对应的脑区位置坐标。2.转换矩阵。点 Electrode coordinates to transformat io n matrix。将含有 sxyz 后缀 文件拖入，生成spinv为后缀的文件。该文件含有电极信息，这一步就将电极信息映射到一 个脑模型里。3将从ERP里导出的txt数据转换成sloreta数据，选择EEG/ERPs to sLORETA ,注意在右侧方框将数据类型改为txt,还可将其上面

15、的复选框勾选上，然后，依次将两个文件夹直接 拖入大方框即可。Number of lines里面显示你拖入的文件数目。最后将spinv为后缀的文 件拖入到第二个空格中。（根据你的目的来选择统计分析的方法） ROUI Fosloreta的统计分析有两种分析思路：以独立样本t检验为例，第一种：就是用ERP操 作方法，选择Convert all time frameso前提是：我们不知道哪个成分有显著差异，只是比较被试在两种条件下激活的差异脑区，例如：向被试分别呈现花朵和灰屏，我们想看一下 被试在两种不同刺激条件下所澈活的脑区是否有显著差异，就可以用这种方法。在某个特定 时间段找到找到两组差异最大的时

16、间点；基于这个最大的时间点，然后在其前后3-5个时间 点围取一时间段（比如：我们找到最大的时间点是75,因此我们的时间窗口可设置72-78以， 之后的sloreta的统计检验就是在此时间段进行）。笫二种：就是直接选取某个时间段进行 sloreta的统计检验，选择Specify time frames （比如：我们对在花朵剌激条件下和灰屏 刺激条件下的平均波幅进行统计分析，发现两种条件在N1成分有差异，其N1时间窗为 80-120ms,我们现在相对某一种脑电成分溯源，就可以采用第二种方法。在这种条件下，我 们的工作就是将ERP数据（只转换某一脑电成分的ERP数据，不需转换整个时间轴的所有数 据。

17、例如：NL,只需将80-120ms之间的ERP数据转换成slor格式即可）直接转换成slor 格式，之后直接进行sloreta的统计分析。第二部分：eLORETA and sLORETA or sLORETAeLORETA and sLORETA （比较两种不用条件下激活的差异脑区，分别对应的是 “When”的工作和“Where”的工作eLORETA : “When” 的工作1. 在这一过程中，我们使用的是从analyser里导出的ERP数据，也就是txt格式的数据。 点击sLORETA主菜单Statistics,进入下图这个工作界面，选择ERP,然后点击NEXT。2. 进入下个界面，如图所示

18、，选择average reference。Number of electrodes里填写上 你选择的电极数目。Number of time frames中填写上你的数据一共包含有多少个时间点。 例如：你截取的时间是600ms,你的基线时间是50ms,原始数据的采样率是1000HZ, 一共 有650个时间点，所以在这个方框里填上650,然后点击NEXT。（特别提示：所有被试原始 数据的采样率必须保持一致，否则导出每个被试的数据大小不一致，在Statistics中将无 法进行统计检验，提示数据提取出错。）3进入下一个界面，根据研究目的选择进行配对样本t检验，还是独立样本t检验，然后点 击 NEXT

19、。4进入下一个界面。选择None,点击Next。| Test details I 5/13|p None厂 Use baseine correction5进入下一个界面。将数据导入，此阶段使用的数据是从analyser里面导出的数据，后缀 为txt。在All files下拉菜单中选择txt,所有txt文件就被选出，然后将其中一个文件 夹中的txt数据拖入到ListA中，将另一个文件夹中的数据拖入到ListB中（或勾选左上 Indude subfolders前面的复选框，直接拖动相关文件夹），这样更方便。点击Nexto6. 进入到下一个界面，选择第三个All tests for each Tim

20、eFrames/Frequencies进行 点对点的分析。点击NextoonesLORETA/eLORETA SnPMQpcn 5ove 八 tlcarAll&Restart Exit HelpHELPTestdetoisBOnX X tests for ol T mcFrorcsrcoucnocs1test for Bvcrge r a second TF rterve Al tests for ccch TimcFrcrc/Frequency in aI :6ficd htervd关于每个复选框的解释如下：1、All tests for al 1 TimeFrame/Frequencies

21、表示对所有文件中的所有时间点（例如： 时间轴是1秒的，time frame是1000,那在此步骤就是对两种条件下的每一个电极点 下的1000点进行平均，之后进行统计检验以获得差异最大的时间点）进行的eloreta 分析，得出的是1000个点的数据。2、1 single test for averge in a specified TF interval 表示对某个时间窗口中（例 如：我们要对N1成分进行定位，那么N1成分的时间窗是80-120ms00003. 320690E-00003. 260374E-00003. 194393E-00003. 12S354E-00003. 067124L0

22、0003.014611E+00002. 973122E+00002. 948212E十00002. 923686E*00002.911722E*00002. 903594E-00002. 89490血00002. 8810S1E*00002. 857934E*00002. 822121E-00002. 771633E*00001234S1011111Z34111156781922O122223425ZO27283壬2+。 3. 4-lH：OOi 3金OE；CC 3 38ET+OO 33 丁圧十OO 3 SSE-OQ 3 32E-OO 3 NTH：十OO 3 ZOE-OQ 3 19E-*-OO

23、3. i 3E：十aa 3 07E700 3 01E-*-OO N 97EZ+OO 2 9 4-E：4-OOi 2 9NR：十OO N. 91E：4-OO22.Z.2-2.2.2.N2-Z.90E4-OO:3 41740? .cic)5 332023 32009.2e)TS5S 2003714393.1.N83S 斗3.3.2.2.2.OOTL24 014011 973122 4-3212911722 03594-9 0 12 3z 3 3 3 3N2.2勺圧十aa SSE-OO 86E-*-OO azHz-oo 77E-OO 了 IE：亠 OO 63EOO DD+OO 5 4-E-*-OO

24、5 4-E4-OO 4-E：4-OO 53H：OO 88LO8L 85TO34 .822121 77LO33 700151 D4TOO1 54342 N 54319 D3ODOH 623C1722.2.2.53E4-OO N 53116952E4-OO 2 52027SlHzOO N 51OB):5. 3494. 290One-Tailed(A00. 3417470. 3592000. 6834930.4820001. 0252400. 5990001. 3669870. 5652001. 7OS7340. 5052002.0504800.3404002. 3922270. 2576002.

25、7339740. 2432003. 0757210. 2050003.4174670. 186000T (0. Oo)丈柯go)歐叭 IgttxH)I .187 10.3.3.4.ExtremeP0. 186200. 853200. 27720Thrsh(lTailel0)Prob(lTa订edThrsh(2Tailed)-0. 1935610. 2374000. 341747-0. 3S71220.3116000.683493-0. 5806830. 3500001. 025240-0. 7742430. 3916001.366987-0. 967S040. 398SOO1.708734-1

26、. 16136&0. 3588002. O5O4SO-1.3549260 6032002. 392227-1.54S4S70. 7076002. 733974 1. 7420480. 8040003. 075721-1.9356090. 8532003. 417467PrsLORETA : “Where” 的工作12. 回到Statistics,回到上次离开时的界面，点击Clearall & Restart,会出现Are you sure you want to clear all dataM 然后选择yes”。13回到下图，这一做统计分析用到的是slor数掲。选择Sloreta(slor)

27、t点击Next。Cp” JFirit DoS hpe 1/1C WsLCftETA-WZ C.ster3)Cross soecrumMeasures deoendence or time-varyrg log-4pecta14.选择 None,点击 Next。15根据实验目的，选择统计方法，例如：配对样本t检验，独立样本t检验，点击NextoNextBackTest details I 5/13tJjoneJC Use baselne correction17.将数据导入，此阶段使用的数据是经过转换的带有slor后缀的数据。在All f订es下 拉菜单中选择slor,所有slor文件就被选出，

28、然后将其中一个文件夹中的slor数据拖入 到ListA中，将另一个文件夹中的数据拖入到ListB中。点击Next。jCOfiETA tier)LfetAOur Sort PMiU-UWljojnjr Jor一 C31_goanwvr jlorJ二 OBqm.wjlor034jon_wrJor3 JtoSjoan.w Jot R 0350 an_*r star 3038oanjr Jtor 038_gom_r JorJ039jO8jr Jor3WOjoan.ir JorJ CHOjooan.wJorJ(H58JT4tof. (HlqmjrjfarJ 35*Vr Jfer1S4MBXMMB*414

29、5 第呦1S.4MI4MB15佃1S.4M8I*I5.4W4M0154MBSt文件 SLOW文件 迪文件 SUX1文件 SLOR 5ft 51(X3 Xft SLOR文井 SLOW文件 SLOR 刪 SLOW文件 心文件 SIOW文件 5LOR文件 98】文件 SLCRJtt2016 WQS19 2016,/2710.2016 附2519 2016S . 2016,/2519. 2016.W 10. 20I6W. XX6.W10.20跡応归 20!6咖10二20 血 W2519 2010CfUsersftoie _Orr CewtDp trietes 护护jtCrthersffcDielOie

30、-Orittopt*iete5C：lbNcote.CherCXDdeMltHet 90如mVW5 jo a C：ltes ffcole.Oe-OesctDpcretflflNetej Qo.an.Mi&Xjo a UUwNcoie.OOMWnodaeto 旨 oo_arW*j C:UsesfSe.ggooaaetmiess 3.m.i5 j)o d C：Uf f kde.Oe* VefntooX*etoNee yv*019o. UUWb 他C*V021OGtberificole.OvOwtopUetoNdo W23jo C：XJK3qvrgtopaaa；ii3c 护尹.30_90 *20咖19

31、 2016怡/2710 2016,6/2519. UgrsMcW-ChenmDpHoretaVxdMy go.an.M0： C：VJs,t*cote-Chene*iDpcLoret3prflbiary 90.an.wrT：C:gyMcote.6engktopMoresmry go.m j*Q. CzVJsefkote.CheHectopetoretordrary po.an.wV C：Ustkole_ChenC*Wet(yetArdrarF 84J*O C：U$Ncote.6engWtymgwy go.an.w C：yNggnZgSSgwyC:grs 仙8.an.*O叙9尹 MOC:V5cn,

32、Nco*e_ChentwaDp,fiioreaccdw go_anXO C：V5cr5f*frdbwy 卩.*。 C：V5er5t*co*e_OwOAttp,filoretaVx,dry go.an_r,|X C:VserjNcoieChenedopcU/etaprflbian 90.6_wg18选择 I single test for average in a specified IF interval v 点击 Next。sLORETA/eLORCTA SrtPM19在“when这一步骤中，我们寻找出差异最大的点，比如说这个最大的点是161,在这 个值前后3-5个值，形成一个区间，例如：1

33、61-167。愷sLORETA/eLORETA SnPMQpen barAll3 Restart xh Help20按照下列标准进行选择，然后最下方灰色的方框，设置输出文件的文件名点击Next。21选择Perform tests,此阶段需要一定时间，请耐心等待。22此过程结束后，如果你在步骤9中设定的文件名称为t-test,目标文件夹中会生成5个 文件，如图所示。雇xStatistics中会生成一个最大值，将这个最大值与Threshold &ExtremePs中的值对比，看是否具有显著性差异。(0 t-test.slort-test.TestSetupt-test-MaxStatistics.

34、txt2016/10/15 22:092016/10/15 22:082016/10/15 22:09SLOR文件TESTSETUP 文件文本曲25 KB5 KB1 KBj 1 t-test-TestSetup.txt_j t-test-Thresholds&Extre?mePstxt7 KB2 KB2016/10/1 5 22:08文本如2016/10/1 5 22:09文本文栏查看定位结果23打开 sloieta,选择 Viewer/Exploreo 定义电极数目（选择 Number of Electrodes,例 如：用于分析的有20个电极，就输入20）,定义采样率（Sampling R

35、ateHz。例如:1000HZ/S. 就输1000）o24. Data Type中选择Sloreta time domain,然后红框就变成下面这样。UORtT A/etORETA wcpta w (iLORE WelORHA time do-nan)24选择 FlieExplore,打开 Folders and files,选择 t-test. slor,双击。然后再 Slice Viewer中会显示出激活的脑区，如图所示。主要激活区域是Brodmann area 10,上额回合和中额回最前测的部分。 IX 1吕trHat*Y.ZX5.W.l5)hH (97C 1) LORKTAa p(21

36、IRIZI00 w5-5iri .5D-5inert-56(H|a rr IFind Max/Min直接进行sLORETA：对某种已知差异显著的脑电成分进行定位。1.选择Specify time frames,输入某一成分的时间窗加基线时间的数据。例如：对N1进 行定位，71的时间是80-120mst基线时长是200ms,那么First time frame应该输入280,2直接进行sLORETA。步骤同13-24 致。第三部分：单样本t检验做组间比较的时候并没有定位到我所需要的脑区，那么可以试图通过单样本 t检验来检测组是否可以出现我所需要看到的定位脑区，并且查看与0比较，是 否有显著激活。

37、1 进行单样本t检验所需要的数据格式为Slor3o拖入一种条件下所需要的数据，数据格式 是txt。选择Specify time frames输入某一成分的时间窗加基线时间的数据。例如：对 N1进行定位，N1的时间是80-120ms,基线时长是200ms,那么First time frame应该输入280, Last time frame 应该输入 320。选择 Compute sLORETA-xyz。Electrodeto coordinates I Regte” e-sl etevod I Gletrodw se、e M/oorQUon mMrs II eeg/ERPb to BLCrrAE

38、EGo tobogcaoectrvm | Croce aoectro to dORTA | Svaby| PM terrm electrodesUOS| 口O】r S LCxrty Tm-vwynQ cfo|、*w tm.-ywyfo cro“ pxtrtJii Tme-vorvrg crotw-Bpedro to pLORCTA I Coriectfvitv2 I Mcrostate segmentatton【ndependent component gnalvBio | CRL bLCRCTA norms6EG/ERP5 to SLORETAIC mjyuantt2 LO_AQQdo-tt

39、2 C *mjyun vrhanQcydcrxj ktt2 3P3rzhengdoy tt2 L8_A.tt2 L9_Aersg txt C *njyuan Vrbangcdor-iQ (tt2 2O_Azr”. txti 93/3n550Qdoy tt23_Ave8e txt I !euruan Uhengador- MtZ4_Averwje txt f fnjyuan hanQcydor-cj tt26_Avara .txt 3/3n5ecg1ytt28_Ag8亡 xxt 3.txt3.txtr at IIFonnt convert I ROI mo匕&rl vnZP Include bu

40、rttranslormattsn matrocOccult2.点击sLORETA主菜单Statistics,进入下图这个工作界面，选择sLORETA-xyz(slor3), 然后点击NEXTosLORETA/eLORETA SnPMOpe n Save - Cl ea r Al !8Re&ta rt - Exit - HelpHELPData type 1/13C ERPC sLORETA (.dor) sLORETA-XYZ (x.sbr3)C Cross spectrumQ Measures of dependesnee or time-varyi3选择None。点击NextsLORETA

41、/eLORETA SnPMOpen Save ClearAll&RestartExit HelpHELPDoto ncrmolizotion【3/13Kone(TimeFrame/FreqL*ney wise Subject-wse(Electrode/oxel) wtse (relative power4.选择 Single group, zero mean testA二0, 点击 NextsLORETA/eLORETA SnPMOpen Save 一 ClearAII&Restart Exit - HelpHELP1List AslORETAfiics (x.sfor3)夫小项目类型1览改

42、餉一 tt23_Averaoe.sior36.49SLOR3文件2D16/1U17 2|ct24_Aver 邓Bor36.49 TOSLOR3文件20L5/1L/17Ztt26_AvTege.dof36.49 rBSLOR3文件2016/11/17 2_ tt28_Averege.dor36.49 MBSLOR3文件2016/11/17；.tt42_Averege3.d.6.49 TOSL0R3文件2016/11/17 2_ tt_AvcrDac3.d.6.49 WSLOR3文件3DL6/1L/17 2tt48_Averaqe3.d.6.49 TOSLOR3文件3DLO/11/17 2.tt_

43、Avreqe3.d.6.49 rBSLOR3袁件3DL6/1L/17Z1 O.AH 姜.2.6.49 r/BSLOR3文件2016/11/17；_ tt212-ApoKV 6.49 MBSL0R3文件2016/11/17 2_ tt21S_A7erag46.49 roSLOR3文件3DL6/1U17 2_ tt216A/eraoEt 6.49SLOR3文件2DL6/1U17 2tt218_A7erageaD.49T/BSLOR3夬件3D 15/11/17219*/39=-.$1.6.49 r*BSLOR3文件2016/11/17；_ “220_A-/erag?.sl.6.49 r*SLOR3文

44、件2016/11/17；_ tt4KLA7eeg3.6 你佗SL0R3文件20 11/17 2_ tt413_A7eraos36.49 k8SLOR3文件2DL6/1U17 2匚gl4_A7erat35.49SLOR3夬件3DL5/1L/17Ciear Sort MJTLres-iaF ：爛鼻疋tt210_Averege.or3F ：渤淖近电.tt212_Avercge.知 3F:明直 Etejyuonbxgq亡op 僭215-A、b?oc 03F:瑚童疋 0Ysrf少crcqUcrc :t216_AvcrMC sior3 F：t5fi：S3J?uan2ercqconatt218_versoe.dor3F：BS3J?uan2erK;qcorcn2Z)_Aver=qe.dor3 尸：耳8滇亡U：心an2hen;qCo为沁3_几已a;e hr3 F： % 憑曰pjanVher创Vongt24.Av臼時皿3 F:湖:皋EZyxn计E%qdorgW26AuEege stor 3 F:斓瘪EbysnVrqdorg皿28二旧心也 stor 3 F ：谢: