本方案描述了使用部分小波变换相干性(pWTC)来计算人际神经同步(INS)的时间滞后模式,以推断社交互动过程中信息传递的方向和时间模式。两项实验验证了pWTC在消除信号自相关对INS干扰方面的有效性。
方法文章
* These authors contributed equally
本方案描述了使用部分小波变换相干性(pWTC)来计算人际神经同步(INS)的时间滞后模式,以推断社交互动过程中信息传递的方向和时间模式。两项实验验证了pWTC在消除信号自相关对INS干扰方面的有效性。
社会互动对人类至关重要。尽管超扫描方法已被广泛用于研究社会互动过程中的人际神经同步(INS),但由于功能近红外光谱(fNIRS)具有相对较高的空间分辨率、良好的解剖定位能力以及对运动伪影极高的耐受性,它已成为超扫描自然情境下社会互动的最受欢迎技术之一。以往基于fNIRS的超扫描研究通常采用小波变换相干性(WTC)计算时间滞后的INS,以描述个体间信息流动的方向和时间模式。然而,该方法的结果可能受到每个个体fNIRS信号自相关效应的干扰。为解决这一问题,研究引入了一种称为部分小波变换相干性(pWTC)的方法,旨在消除自相关效应,同时保持fNIRS信号的高时间-频谱分辨率。本研究首先通过模拟实验展示了pWTC在消除自相关对INS影响方面的有效性;随后,基于一项社会互动实验的fNIRS数据集,提供了pWTC操作的逐步指导。此外,本文还比较了pWTC方法与传统WTC方法、以及pWTC方法与格兰杰因果关系(GC)方法之间的差异。结果表明,pWTC可用于确定不同实验条件下INS的差异,以及自然社会互动中个体间INS的方向性和时间模式。同时,与传统WTC相比,pWTC具有更优的时间和频率分辨率;与GC方法相比,则具有更高的灵活性。因此,pWTC是推断自然社会互动中个体间信息流方向与时间模式的有力候选方法。
社会互动对人类至关重要1,2。为了理解社会互动的双脑神经认知机制,近年来研究者广泛采用超扫描(hyperscanning)方法,发现人际神经同步(interpersonal neural synchronization, INS)的模式能够很好地表征社会互动过程3,4,5,6,7,8,9,10,11,12,13,14。在近期研究中,一个有趣的发现是,成对个体在互动中的角色差异可能导致INS出现时间滞后模式,即某一个体的脑活动比另一个体滞后数秒,例如从听者到说话者5,9、从领导者到追随者4、从教师到学生8、从母亲到儿童13,15,以及在情侣关系中从女性到男性6。尤为重要的是,时间滞后的INS间隔与社会互动行为的时间间隔具有良好对应关系,例如教师提问与学生回答之间8,或母亲育儿行为与儿童顺从行为之间15。因此,时间滞后的INS可能反映了一种从一个个体向另一个个体传递的定向信息流,正如近期提出的人际言语交流层级模型所假设的那样16。
以往,由于功能近红外光谱(fNIRS)信号具有相对较高的空间分辨率、良好的解剖定位能力,以及在研究自然情境下的社会互动时对运动伪影具有极高的耐受性17,时间滞后的脑间同步性(INS)主要基于fNIRS信号进行计算。此外,为了精确刻画社会互动过程中神经时间滞后与行为时间滞后之间的对应关系,有必要获取每个时间滞后点(例如,从无时间滞后到10秒的时间滞后)的INS强度。为此,以往研究广泛采用小波变换相干性(WTC)方法,在将一个个体的脑信号相对于另一个体的脑信号进行前移或后移之后进行分析5,6,18。然而,在将这种传统的WTC方法应用于fNIRS信号时,可能存在一个潜在挑战,即观测到的时间滞后INS可能受到个体fNIRS信号自相关效应的干扰19,20,21。例如,在二元社会互动过程中,参与者A在时间点t的信号可能与参与者B在同一时间点的信号同步;同时,由于自相关效应,参与者A在时间点t的信号也可能与其自身在稍后时间点t+1的信号同步。因此,参与者A在时间点t的信号与参与者B在时间点t+1的信号之间可能出现虚假的时间滞后脑间同步性。
Mihanović 及其同事22首次提出了一种称为部分小波变换相干性(pWTC)的方法,并随后将其应用于海洋科学领域23,24。该方法最初的目的是在估计两个信号的相干性时控制外源性混杂噪声。本文中,为解决fNIRS超扫描数据中的自相关问题,将pWTC方法扩展用于计算fNIRS信号上的时间滞后型脑间同步(INS)。具体而言,可利用以下公式(公式1)23计算从参与者A到参与者B的时间滞后型INS(以及方向性信息流)。

此处假设有两个信号, A 和 B, 分别来自参与者 A 和 B。信号的发生情况 B 总是先于信号 A 存在时间延迟 n,其中 WTC (At, Bt+n是传统的时滞WTC。 WTC (At, At+n) 是参与者 A 的自相关小波相干性(autocorrelated WTC)。 WTC (At, Bt) 是时间点对齐的小波相干性(WTC) t 在参与者 A 和 B 之间. * 是复共轭算符(图1A).

图 1:pWTC 概述。(A)pWTC 的逻辑。在一个二元组中有两个信号 A 和 B,信号 A 的出现总是滞后于信号 B,滞后时间为 n。灰色方框表示在某一时间点 t 或 t+n 的小波窗口。根据图中所示的 pWTC 公式,需要计算三个 WTC:信号 At+n 与 Bt 的时间滞后 WTC;参与者 A 内部 At 与 At+n 的自相关 WTC;以及在时间点 t 处对齐的 At 与 Bt 的时间对齐 WTC。(B)光极探头组的布局。根据国际 10-20 系统27,28,CH11 置于 T3 位置,CH25 置于 T4 位置。请点击此处查看该图的放大版本。
本方案首先引入了一个模拟实验,以展示pWTC在解决自相关挑战方面的效果。随后,基于自然情境下社会互动的实证实验,逐步说明了如何实施pWTC。本文采用交流情境来介绍该方法,原因是此前时间滞后的INS通常在自然交流情境中进行计算3,4,6,8,13,15,18。此外,还进行了pWTC与传统WTC的比较,并通过格兰杰因果关系(Granger causality, GC)检验进行了验证。
访问受限。请登录或开始试用以查看此内容。
本研究的人体实验方案已获得北京师范大学认知神经科学与学习国家重点实验室机构审查委员会和伦理委员会的批准。所有参与者在实验开始前均签署了书面知情同意书。
1. 模拟实验


2. 实证实验



访问受限。请登录或开始试用以查看此内容。
模拟结果
结果表明,具有自相关的时间滞后 INSWTC 显著高于无自相关的时间滞后 INSWTC(t(1998) = 4.696,p < 0.001),也显著高于时间滞后 INSpWTC(t(1998) = 5.098,p < 0.001)。此外,无自相关的时间滞后 INSWTC 与 INSpWTC 之间无显著差异(t(1998) = 1.573,p = 0.114,图 2A)。这些结果表明,pWTC 能有效消除自相关效应对 INS 的影响。此外,当 WTC 值接近 0 或 1 时,只要 WTC 值偏离 0 或 1,时间滞后 INSpWTC 仍能呈现可靠的结果(补充图 2)。
实证实验结果
访问受限。请登录或开始试用以查看此内容。
在超扫描研究中,通常需要描述个体之间信息传递的方向性和时间模式。以往大多数基于功能性近红外光谱(fNIRS)的超扫描研究采用传统的WTC25方法,通过计算时间滞后的脑间同步性(INS)来推断这些特征。然而,fNIRS信号的一个固有特性是存在自相关效应20,21,这可能会混淆时间滞后的INS结果。为解决这一问题,本方案引入了一种称为pWTC的方法22。该方法在部分剔除自相关效应后估计时间滞后的INS,同时保留了WTC方法的优势。本方案提供了pWTC操作的逐步指导,并通过将pWTC结果与传统WTC和格兰杰因果关系(GC)检验结果进行比较,验证其有效性。
本方案展示了在基于功能性近红外光谱(fNIRS)的超扫描数据中应用部分小波相干性(pWTC)的关键步骤。具体而言,首先需要基于时间滞后的fNIRS时间序列计算时间滞后的小波相干性(WTC)、自相关WTC以及时间对齐的WTC。接...
访问受限。请登录或开始试用以查看此内容。
作者声明不存在竞争性财务利益。
本工作得到国家自然科学基金(61977008)和万人计划青年拔尖人才项目的资助。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| fNIRS拓扑系统 | Shimadzu Corporation | Shimadzu LABNIRS systen | LABNIRS系统包含40个发射器和40个探测器,用于fNIRS信号的测量。在本方案中,我们利用这些发射器和探测器构建了两个定制的26通道探头组,并根据10-20系统将它们分别固定在两个帽具上。此外,LABNIRS系统还内置了图形用户界面(GUI)软件,可用于数据质量检查、数据转换和数据导出。 |
| MATLAB | The MathWorks, Inc. | MATLAB 2019a | 在本方案中,使用了MATLAB中的多个内置工具箱和函数: 使用SPM12工具箱通过其图形用户界面(GUI)对有效的MRI数据进行标准化处理。 使用NIRS_SPM工具箱通过其图形用户界面(GUI)将探头的MNI坐标投影到AAL模板上。 使用Homer3工具箱中的hmrMotionCorrectWavelet函数(默认参数)去除运动伪影。 使用小波工具箱中的wcoherence函数计算WTC和pWTC。 |
| MRI扫描仪 | Siemens Healthineers | TRIO 3-Tesla扫描仪 | 在本方案中,使用MRI扫描仪获取每个通道和光极的MNI坐标。扫描参数在正文中有详细描述。 |
| 定制帽具 | 在本方案中,我们首先在两个尼龙帽上标记出10-20系统的位置。然后,制备了两个26通道的定制光极探头组。最后,将探头组按照解剖标志对齐并固定在帽具上。 |
访问受限。请登录或开始试用以查看此内容。
申请许可以重复使用本 JoVE 文章的文本或图表
申请许可