小波变换相干性(WTC)是功能近红外光谱(fNIRS)超扫描研究中用于评估信号间耦合的常用方法。本文介绍了一种用于评估信号相互作用方向性的工具箱。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
小波变换相干性(WTC)是功能近红外光谱(fNIRS)超扫描研究中用于评估信号间耦合的常用方法。本文介绍了一种用于评估信号相互作用方向性的工具箱。
尽管功能性近红外光谱(fNIRS)超扫描研究日益增多,但利用小波变换相干性(WTC)评估两个神经信号之间耦合性的方法似乎忽略了相互作用的方向性。目前该领域缺乏一个能够帮助研究人员判断通过WTC函数获得的高相干值究竟反映的是同相位同步(即双人组中两名成员同时出现神经激活)、滞后同步(即双人组中一名成员的神经激活先于另一名成员)还是反相位同步(即双人组中一名成员神经激活增强,而另一名成员神经激活减弱)的分析框架。为满足这一需求,本研究提出了一种互补且更敏感的双神经信号相位相干性分析方法。该工具箱允许研究者通过将传统WTC获得的相位角值分类为同相位同步、滞后同步和反相位同步,来估计耦合的方向性。此外,该工具箱还使研究人员能够评估交互动态在任务过程中如何发展和变化。通过在fNIRS超扫描研究中应用这一新型WTC方法及配套工具箱,将有助于深化我们对复杂社会互动机制的理解。
近年来,为理解社会行为的神经基础而开展的研究类型发生了转变1,2。传统上,社会神经科学研究主要关注个体在执行与社会相关任务时,单个孤立大脑中的神经活动。然而,神经影像技术的进步使得研究者能够在“真实生活”情境下,同步考察一个或多个个体在社会互动过程中的脑内神经激活情况3。在“真实生活”情境中,个体能够自由活动,随着信息的交换以及社会互动双方相互接收反馈,脑激活模式也可能随之动态变化4。
超扫描是一种通过同时测量两个或更多个体的脑活动来评估这种双向信息交换的方法5近年来,越来越多的研究采用功能性近红外光谱成像(fNIRS)这一非侵入性神经影像技术,相较于其他神经影像技术,该方法对运动伪影的敏感性较低。6超扫描 通过 fNIRS 可在真实生活场景中评估互动伙伴间的脑间同步(IBS),同时允许参与者自由、自然地活动。这对于婴幼儿研究尤为重要,因为婴幼儿通常活动较为频繁。已有研究报道,IBS 能够反映互动双方之间的相互理解,这种相互理解是有效社会互动与交流的基础,并介导了共享意图性1,7,8.
评估两个大脑间 IBS 的方法有多种,包括时间序列相关性分析,如互相关和皮尔逊相关系数9,10(参见 Scholkmann 等人10 的综述)。其他方法则涉及在频域中评估耦合强度,例如相位锁定值(PLV)和相位相干性(参见 Czeszumski 等人11 的综述)。在 fNIRS 研究中,最常用的方法之一是小波变换相干性(WTC),即衡量两个时间序列在频率和时间上的互相关性10。
WTC 使用相关分析来计算两个时间序列在时频域中的相干性和相位滞后。功能性近红外光谱(fNIRS)超扫描研究已采用 WTC 在多种功能领域中估计个体间的脑同步性(IBS),包括动作监测。12,合作与竞争行为5,13,14,15,模仿16母婴问题解决17以及教学行为18,19,20,21通常,在双人同步神经成像研究中,会通过小波时频相干性(WTC)测量实验任务期间的跨脑相干性,并将其与控制任务期间的跨脑相干性进行比较。这些结果通常以WTC“热图”形式呈现,该图显示两个大脑在各个时间点和频率上的相干性水平(参见 图1).
正如 Czesumaski 等人11所建议的,小波相干性分析(WTC)已成为分析 fNIRS 超扫描的标准化分析方法。WTC 分析是一种灵活且“与工具无关”的数据可视化和解释方法22。其主要优势在于相干系数热图,该图提供了一种叙述性分析形式,便于识别同步或异步行为的时间段以及在完成任务期间大脑活动的强度,这使得 WTC 成为应用研究中的有力工具22。WTC 相较于相关性分析技术具有优势。相关性分析对血流动力学响应函数(HRF)的形状敏感,而 HRF 被认为在个体之间(尤其是年龄方面)以及不同脑区之间存在差异。相比之下,WTC 不受脑区间 HRF 变化的影响23。研究人员已使用小波方法研究 fMRI 时间序列。Zhang 等人24比较了常用的几种功能连接度量方法,包括皮尔逊相关、偏相关、互信息以及小波相干性变换(WTC)。他们利用来自静息态 fMRI 数据和观看视频的自然刺激 fMRI 数据中提取的大规模功能连接模式进行了分类实验。研究结果表明,WTC 在分类性能(特异性、敏感性和准确性)方面表现最佳,提示 WTC 至少在分类应用中是研究功能脑网络更优的功能连接度量方法24。

图 1:小波变换相干性(WTC)。 WTC 显示了两个时间序列在时间(x 轴)和频率(y 轴)上的相干性和相位角。图中红色表示相干性增强,图中的小箭头表示两个时间序列的相位角。箭头向右表示同相同步;箭头向下和向上分别表示滞后同步;箭头向左表示反相同步30。该图改编自 Pan 等人19。请点击此处查看此图的放大版本。
最近,Hamilton25 阐述了在 fNIRS 超扫描研究中解释跨脑相干性数据的若干局限性。Hamilton 的主要关注点之一在于,相干性度量(例如,WTC)仅能报告对称性效应(即两个大脑相关,表现出相同的变化模式)。然而,许多社会互动具有不对称性(例如,说话者与听者之间的信息流动),其中两名参与者可能扮演不同角色,而 WTC 是否能够捕捉此类信息尚不明确。本文通过一种新框架解决了这一问题,该框架利用交叉小波相位来检测方向性,从而对交叉小波功率进行直接解释。该框架还将有助于考察互动动态过程如何在任务过程中发展和变化。
尽管WTC和相关性方法用于评估功能连接性,但其他方法则用于评估有效连接性,旨在提取一个神经单元对另一个神经单元的因果影响。转移熵是信息论领域中一种用于描述联合依赖过程之间信息传递的度量方法26。另一种相关方法是格兰杰因果分析(GCA),该方法已被认为在本质上等同于转移熵26。
在现有的fNIRS超扫描研究文献中,格兰杰因果分析(Granger causality analysis, GCA)已被广泛用于估计在多种不同任务过程中获取的fNIRS时间序列数据之间的耦合方向性,例如合作5、教学19和模仿16。GCA采用向量自回归模型来评估脑数据中时间序列之间耦合的方向性。格兰杰因果性基于预测和时序优先性:“如果变量X的过去值包含有助于预测变量Y未来值的信息,且该信息超出Y自身过去值所包含的信息,则称变量X‘格兰杰导致’(G-cause)变量Y”27。因此,格兰杰因果性在两个方向上进行分析:1)从被试A到被试B;2)从被试B到被试A。
尽管GCA分析可作为补充分析手段,用于判断通过WTC函数获得的高相干值反映的是IBS还是滞后同步(即一个信号领先于另一个信号),但它无法确定是否发生了反相位同步。在传统的神经影像学研究中,仅对一名参与者进行扫描(即“单脑”方法),反相位模式意味着一个脑区的活动增强,而另一个脑区的活动减弱28。在超扫描研究文献中,反相位同步的存在可能表明,一名受试者的神经激活增强,而与此同时,另一名受试者的神经激活减弱。因此,有必要建立一个能够检测方向性的综合模型。更具体地说,该模型除了能够检测同相位同步和滞后同步外,还能够检测反相位同步(即个体的活动方向与其同伴相反)。
为了解决WTC仅显示对称效应(即两个大脑表现出相同变化模式)的担忧25,本文提出了一种新方法,通过考察同步的相位(即同相、滞后或反相)来识别相互作用的类型(见图2)。为此,开发了一个基于WTC方法对不同类型相互作用进行分类的工具箱。该分类通过使用交叉小波变换分析得到的相对相位数据实现。

图2简单正弦波不同相位关系的示意图。 (A)当两个信号,信号1(蓝线s)和信号2(橙色线s),若它们在同一时间点达到各自的最大值、最小值和零值,则称其表现出同相位同步32. (B当一个信号达到其最大值而另一个信号同时达到零值时,称为表现出滞后同步(其中一个信号超前) 90°)32,33,34. (C当两个时间序列向相反方向移动,即一个信号达到最大值的同时另一个信号达到最小值,这种情况被称为反相同步28. (D-P) 在两个时间序列的所有其他相位关系中,一个信号领先于另一个信号。在所有正相位中,信号2领先于信号1(例如,图中各面板 E, F, M,以及 N),而在所有负相位中,信号1超前于信号2(例如,图示面板 D, G, H, O,以及 P)。值得注意的是,当相位的绝对值较大时,哪一个时间序列领先于另一个序列的特征变得更加明显(例如,在图示面板中,领先的特征更为显著 J 比图中面板 I,以及在面板中 K,该领导作用比在面板中更为显著 L). 请点击此处查看此图的放大版本。
访问受限。请登录或开始试用以查看此内容。
该研究在佛罗里达大西洋大学(FAU)进行,并获得了FAU机构审查委员会(IRB)的批准。
1. 使用Homer3软件(材料表)对fNIRS超扫描数据进行预处理
注意:Homer3 是一个 MATLAB 应用程序,用于分析 fNIRS 数据以获得脑激活的估计值和激活图谱29。Homer3 可通过以下链接下载并安装(https://openfnirs.org/software/homer/)。
2. 开始使用 LeaderFollowerByPhase 工具箱

图3:工作流程概览。(A)母子配对在自由游戏过程中,同步采集超扫描fNIRS数据。(B)母子时间序列示意图。(C)使用Homer3对时间序列进行预处理。(D、E)利用工具箱分析不同类型的互动,如同相位同步、反相位同步和滞后同步。请点击此处查看该图的高清版本。

图 4:基于相位的四种不同类型相互作用的分类。 该图表示两个神经时间序列在360°模量下的相位差。相位差可被视为两个值之间的时间滞后,以度、弧度或波长的分数来度量。此处,360°模量被划分为四个不同的区间,代表四种不同的相互作用相位:(A)信号1超前(约90°区间,介于45°至135°之间),(B)信号1与信号2反相同步(约180°或−180°区间,介于135°至−135°之间),(C)信号2超前(介于−135°至−45°之间),(D)同相同步(约0°区间,介于−45°至45°之间)。此划分方式为默认方法(每个点周围45°);然而,该工具箱允许用户自定义不同的划分方式。尽管其他配置可能不覆盖完整的360°,但可能对每种相互作用类型提供更精确的定义。请点击此处查看该图的放大版本。
访问受限。请登录或开始试用以查看此内容。
本节演示了使用该工具箱可执行的分析类型(工具箱可在 https://www.ariel.ac.il/wp/sns/download/ 或 https://github.com/Minisharmaa/Leader-Follower-By-Phase. 下载)。这些分析使用了一组少量婴儿-父母配对采集的 fNIRS 数据。共招募了六对母婴配对,采用一种经过验证的行为任务——自由游戏任务31进行测试,该任务尽可能接近真实的婴儿-母亲互动情境。实验前,为婴儿和父母佩戴了定制的光学探头组,用于采集 fNIRS 数据。本研究中使用的光学探头组包含 8 个光源(红点)和 8 个探测器(蓝点),配置形成 18 个通道,双侧覆盖前额叶和颞顶区域(见图 5)。使用 NIRScout 系统通过两个波长采集光学成像数据:760 nm(对脱氧血红蛋白 HbR 更敏感)和 850 nm(对氧合血红蛋白 HbO 更敏感)。所有父母均为女性(年龄范围 = 26–36 岁),婴儿均为健康足月儿(两名女性,四名男性,年龄范围 = 1–2 岁),无已知发育迟缓。配对...
访问受限。请登录或开始试用以查看此内容。
在功能性近红外光谱(fNIRS)研究中,最常用的方法之一是小波变换相干性(wavelet transform coherence, WTC),它用于衡量两个时间序列在频率和时间维度上的互相关程度10。WTC通过相关性分析计算两个时间序列之间的相干性和相位滞后(补充文件1)。在fNIRS超扫描研究中,WTC已被用于估计多种功能领域的脑间同步性(IBS),包括动作监控12、合作与竞争行为5,13,14,15、模仿16、母婴问题解决17以及教学-学习行为18,19...
访问受限。请登录或开始试用以查看此内容。
作者声明,本研究是在不存在任何可能被解释为潜在利益冲突的商业或财务关系的情况下进行的。
感谢国家自然科学基金(编号:62207025)、教育部人文社会科学研究项目(编号:22YJC190017)以及中央高校基本科研业务费对潘亚峰的资助。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| NIRScout | NIRx Medical Technologies, LLC | n.a. | 8 个光源,8 个探测器 |
| MATLAB | The Mathworks, Inc. | Matlab 2022a | 本实验方案中使用了多个 MATLAB 工具箱及内置函数:采用 HOMER3 工具箱将光强信号转换为光密度(OD),并通过其函数 hmrMotionCorrectWavelet(使用默认参数)去除运动伪影,再将 OD 转换为浓度(Conc);使用小波工具箱(Wavelet Toolbox)计算小波相干性(WTC)。 |
访问受限。请登录或开始试用以查看此内容。