评估非侵入性脑刺激中的大脑-行为耦合：来自可靠变化指数与蒙特卡洛模拟的见解

《Cortex》：Assessing Brain-Behaviour Coupling in Non-invasive Brain Stimulation: Insights from Reliable Change Indices and Monte Carlo Simulations

　　

非侵入性脑刺激（NiBS）研究经常报告应用NiBS后的探索性大脑-行为相关性。然而，这些分析通常功效不足，并依赖于具有已知统计学局限的基于比值的变化分数。研究人员通过汇集三项独立研究的个体数据（N = 69），并将实证分析与蒙特卡洛模拟框架相结合，解决了这些局限。所有研究均在停止信号任务中施加20 Hz经颅交流电刺激（tACS），靶向前辅助运动区（preSMA）和右侧额下回（rIFG）。研究人员使用可靠变化指数（RCI）量化preSMA–rIFG脑电图（EEG）虚部相干性（ImCoh）和停止信号反应时（SSRT）的变化，并在研究内进行z标准化，以进行合并混合效应回归。研究人员发现tACS诱导的ImCoh变化与SSRT变化之间无有意义的关联（r = .013，边际R² = .004），项目层面的相关性较小、不显著且方向不一致。虽然使用基于比值的变化分数的敏感性分析收敛于同样的零结果（r = .014），但模拟框架揭示，基于比值的分数在所有合理的效应量下几乎完全无法恢复真实的大脑-行为关联。相反，基于RCI的分数保留了方向准确性，并提供了校准良好的I型错误控制，尽管在现实噪声条件下存在中度衰减。这些结果凸显了在20 Hz tACS后检测个体水平大脑-行为耦合的严重统计局限。研究人员的结果并非排除真实的生物学关联，而是表明标准样本量和典型测量噪声造成了对现实效应量的实质性“盲区”，强调了在NiBS研究中需要更可靠的指标和充分功效的设计。

非侵入性脑刺激（NiBS）包括经颅磁刺激（TMS）和经颅电刺激（tES）等，被认为是调节神经生理活动与行为的潜力工具。然而，NiBS诱导的神经生理调节（如运动诱发电位）能否可靠预测相应行为结果（如反应时改善）仍是该领域的关键问题。确立可靠的大脑-行为关系，对于识别刺激功效的机制标志物以及指导NiBS干预的优化与转化至关重要。既往研究多报告施加NiBS后的探索性大脑-行为相关分析，但这些分析通常功效不足，并依赖具有已知统计学缺陷的基于比值的变化分数。比值类指标会增加抽样变异、违反正态性假设，并可能引入虚假相关。相反，可靠变化指数（RCI）将前后差异相对于测量的预期变异和重测信度进行标准化，为量化刺激诱导的个体变化提供更合理的替代方案。基于此，研究人员开展本研究，旨在评估20 Hz经颅交流电刺激（tACS）诱导的神经生理与行为变化之间的个体水平关联，并系统比较不同变化分数指标的统计性能。

研究人员汇集了三项独立研究的个体数据（总N = 69），所有研究均采用20 Hz tACS靶向前辅助运动区（preSMA）和右侧额下回（rIFG），并在停止信号任务中采集脑电图（EEG）虚部相干性（ImCoh）和停止信号反应时（SSRT）。研究人员使用RCI量化个体变化，并将RCI分数在研究内z标准化后进行混合效应回归（mixed-effects regression），以检验ImCoh变化与SSRT变化的关联。同时，通过敏感性分析考察比值变化分数，并运用蒙特卡洛模拟（Monte Carlo simulation）在已知真实相关的情况下对比RCI与比值方法的恢复能力。实证结果显示，ImCoh RCI与SSRT RCI之间无显著关联；模拟结果揭示RCI方法能保留方向性并控制I型错误，而比值方法几乎无法恢复任何真实关联。这些结果强调，当前样本量和测量噪声条件下，标准设计对现实效应量存在实质性“盲区”。

主要技术方法：数据来自三项独立研究（Fujiyama等人，2023；Tan等人，2025；Fujiyama等人，2025）的69名健康年轻参与者；使用RCI量化变化，标准公式假设重测信度为.80；将RCI分数在各研究内z标准化，以项目为固定效应、参与者为随机截距进行混合效应回归；使用post/pre比值进行敏感性分析；蒙特卡洛模拟生成与样本匹配的合成数据，在6个真实相关水平（0–0.5）、3个均值偏移和3个噪声水平下全交叉，每种条件迭代1000次，比较RCI与比值方法的平均估计、标准差、偏倚和显著性比例。分析使用R 4.5.0的lme4、lmerTest、MASS等包。

研究结果部分如下。

变化分数总结：描述性统计显示，总体ImCoh RCI均值为0.20（SD = 2.03），SSRT RCI均值为0.38（SD = 1.28），表明平均而言参与者表现出轻微改善；但研究间均值差异大（ImCoh RCI范围从-0.13到0.63），个体间变异性显著。

大脑-行为关联：混合效应回归分析显示，ImCoh RCI不是SSRT RCI的显著预测因子（β = -0.05，SE = 0.10，t = -0.51，p = 0.99），固定效应仅解释0.4%的方差（边际R² = .004），总体相关r = 0.013（p = 0.92），表明不存在可检测的个体水平关联。

项目层面关联：三个项目的描述性相关均较小且不显著（Fujiyama等人，2023：r = 0.04，p = 0.82；Tan等人，2025：r = -0.20，p = 0.43；Fujiyama等人，2025：r = 0.05，p = 0.86），方向不一致，不支持系统性大脑-行为耦合。

敏感性分析：使用比值分数重新计算，总体相关仍接近零（r = 0.014），混合效应模型不显著（β ≈ 0，p = 0.91），与RCI结果一致；但比值分数分布严重偏离正态（SSRT比值的Shapiro-Wilk W = 0.16，p < .001），而RCI接近正态（W = 0.975，p = .183）。

模拟结果：模拟框架显示，RCI估计相关随真实相关增加呈线性上升，但存在系统性衰减（例如真实r = 0.5时，恢复为r = 0.373，偏倚-0.127）；比值估计在所有真实相关水平下均接近零，偏倚最高达-0.500。真实效应的恢复：RCI方法能保留方向准确性，而比值方法几乎无敏感性。偏倚与变异性：两种方法变异性相近，但RCI的负偏倚较小且随效应量增大，比值方法的偏倚严重且随效应量线性增大。H₀拒绝率（α = .05；I型错误与统计功效）：零相关下两种方法I型错误率接近名义水平（RCI：0.050；比值：0.051）；在非零相关下，RCI功效随效应量增加而提高（真实r = 0.3时功效为0.467，r = 0.5时为0.831），而比值方法的功效始终保持在0.051–0.058之间，无法检测任何真实效应。

讨论部分总结：研究人员指出，实证零结果不应被解释为生物学上不存在大脑-行为耦合，而是反映了当前统计功效和测量噪声的限制。即便真实存在中等关联（r ≈ 0.3），在N = 69和典型噪声条件下检测概率仅约46.7%，因此零结果更像是对效应量上限的约束。RCI方法比比值方法更可靠，且表明组水平效应与个体水平关联可以分离。研究局限性包括RCI计算假定重测信度为0.80、三项研究在刺激参数和任务细节上存在差异，以及静态线性连接指标可能无法捕捉动态非线性耦合。未来需要更敏感的神经标志物、整合多模态成像和闭环刺激方案。

研究结论翻译：总之，尽管既往工作已报道20 Hz tACS后preSMA–rIFG连接性的组水平增强，但这些神经生理变化转化为个体水平反应抑制改善（以SSRT测量）仍然微弱且不一致。即使在合并和协调的数据集中，这种微弱的大脑-行为耦合也凸显了tACS反应中巨大的个体间和个体内变异性。除对解释设计用于捕获前后变化的研究中的事后相关性提供警示外，本研究还提供了方法论基准：结合RCI变化分数和混合效应建模提供了一种结构化途径，以更好地评估个体水平变化并提高研究间可比性。同时，结果凸显了需要识别更敏感的神经生物标志物、整合多模态成像，并探索状态依赖性（即闭环）刺激方案，以改善行为结局并推进对NiBS神经生理反应机制的理解。