我们详细介绍了如何使用基于种子点的被试图像置换元分析方法(Seed-based d Mapping with Permutation of Subject Images, SDM-PSI)对基于体素的神经影像学研究进行荟萃分析。
我们详细介绍了如何使用基于种子点的被试图像置换元分析方法(Seed-based d Mapping with Permutation of Subject Images, SDM-PSI)对基于体素的神经影像学研究进行荟萃分析。
目前大多数用于进行基于体素的神经影像研究荟萃分析的方法,并非评估效应是否不为零,而是评估统计显著性峰值是否存在空间聚集性,并将证据评估简化为仅基于 p 值的二分类(即体素只能被判定为“具有统计学显著性”或“不具有统计学显著性”)。本文详细介绍了如何使用基于种子点的 d 映射结合被试图像置换法(Seed-based d Mapping with Permutation of Subject Images, SDM-PSI)开展荟萃分析,该新方法采用标准的置换检验来评估效应是否不为零。我们还展示了如何根据一系列标准对证据强度进行分级,这些标准涵盖不同水平的统计显著性(从较宽松到较严格)、数据量大小以及潜在偏倚的检测(如小样本效应和显著性过高现象)。为示范该流程,本文详细描述了一项针对强迫症的基于体素的形态学研究的荟萃分析,并提供了从已发表文献中提取的全部数据,以便读者能够轻松复现该分析。SDM-PSI 还可用于功能磁共振成像、弥散张量成像、正电子发射断层扫描以及基于表面的形态学研究的荟萃分析。
自磁共振成像技术问世以来,神经影像学领域已发表了数千项关于心理功能和神经精神障碍神经基础的研究。为了总结这些发现,研究人员已开发出多种方法1,2,3,4,5,6。原始的基于体素的神经影像学研究通常报告统计显著性峰值的坐标(例如,患者组与对照组之间灰质体积比较的结果),而元分析方法则通常评估这些峰值在特定脑区内是否存在空间聚集。
然而,我们先前已表明,这些用于检验峰收敛性的方法依赖于一些微妙的假设,这些假设可能影响元分析结果的模式及其统计学显著性7。具体而言,这些检验假设体素之间相互独立,且每个体素出现“假”峰的概率相同;而在真实的灰质中,体素与其邻近体素存在相关性,且某个体素出现“假”峰的概率取决于其组织成分。此外,这些方法还包含一些悖论,例如当真实效应较少时统计效能反而增加,而当存在多个真实效应时统计效能却下降。
为解决这些问题,我们开发了一种方法,该方法对每项研究的统计效应脑图进行填补,然后进行标准的随机效应元分析,以正式检验这些效应是否与零存在显著差异。该方法称为“基于种子点的d映射结合被试图像置换法”(Seed-based d Mapping with Permutation of Subject Images, SDM-PSI)8,其主要特点包括:
我们已在其他地方详细描述并完全验证了SDM方法4,8,10,13,14。
此外,我们建议不要依赖基于统计显著性水平(显著与不显著)对体素进行二分类,而是应使用一组标准来评估证据强度22。二分类的统计显著性简化方法会导致对假阳性率和假阴性率的控制不佳15,而采用的标准则利用统计显著性水平的范围,并考虑数据量或潜在偏差。SDM-PSI 软件可提供进行此类分类所需的各项要素8,因此可用于实现对证据强度更精细的分类。
本文展示了如何使用 SDM-PSI 对基于体素的神经影像学研究进行元分析。为示范该操作流程,我们采用了一项已发表的关于强迫症(OCD)患者灰质异常的基于体素的形态学研究元分析中的数据4。然而,我们并未采用该早期元分析中所使用的方法,而是采用上述最先进的分析流程。读者可从我们的网站(http://www.sdmproject.com/)下载软件及这些数据,以重复本分析过程。
所有旨在对基于体素的神经影像学研究进行荟萃分析的研究人员均可遵循本方案。该方法适用于功能磁共振成像(fMRI,例如对刺激的BOLD反应)16、基于体素的形态学分析(VBM,例如灰质体积)17、扩散张量成像(DTI,例如分数各向异性)18、正电子发射断层扫描(PET,例如受体占有率)19以及基于表面的形态学分析(SBM,例如皮层厚度)等研究或数据集。
1. SDM-PSI 的安装
2. Meta分析计划
3. 穷举搜索
4. 数据采集
5. 将数据导入 SDM-PSI
6. 预处理
7. 主要分析
8. 异质性、发表偏倚与分级
如在 MRICron 中打开的主分析阈值图所示(图2,步骤 7.6),强迫症患者的背侧前扣带/内侧额叶皮层灰质体积在统计学上显著减小。配套网页显示,该脑区簇体积中等偏小(143 个体素),主要位于 Brodmann 32 区,其峰值位于 MNI 坐标 [2, 32, 32],z 值为 -4.97,经 FWER 校正后的 p 值为 0.01。
在步骤 8.1 和 8.2 获取的网页中,较低的 I2 统计量(1.5%)表明异质性极小,漏斗图未显示不对称性(图 3),小研究效应检验和显著性过剩检验结果均为阴性。然而,背侧前扣带回皮层灰质体积较小的证据强度中等偏弱,如在 MRICron 中打开地图进行证据分级时所示,这尤其与数据量有限有关。
当主要分析采用较宽松的统计显著性水平进行阈值设定时(步骤 7.3),患者在纹状体和上顶叶皮层也显示出统计学上显著更大的灰质体积(未经校正的 p 值分别为 0.00006 和 0.0002),但这些异常的证据应被视为较弱。

图1:SDP-PSI 图形用户界面在平均执行过程中的主窗口。 请点击此处查看此图的放大版本。

图 2:与匹配的健康对照相比,强迫症患者灰质体积显著较小的脑区。
该显著性簇包含 143 个体素,峰值位于 MNI [2,32,32],主要涉及背侧前扣带回/内侧额叶皮层,Brodmann 32 区。 请点击此处查看该图的放大版本。

图3:背侧前扣带回皮层灰质体积显著减小区簇峰值的漏斗图。 请点击此处查看该图的高清版本。
| 软件包 | SDM 中的编码 |
| 统计参数映射(SPM) | spm |
| FMRIB 软件库(FSL) | fsl |
| 其他软件包 | other |
| 立体定向空间 | SDM 中的编码 |
| 蒙特利尔神经学研究所(MNI) | mni |
| 原始 Talairach 空间 | tal |
| 使用 Brett 变换将 MNI 转换为 Talairach 空间 | brett |
表1:SDM-PSI 支持的软件包与立体定向空间列表。
| 当以下情况时,t 值必须为正值: | 当以下情况时,t 值必须为负值: | |
| 单样本 fMRI 研究 | 任务 > 基线(激活) | 任务 < 基线(去激活) |
| 双样本 fMRI 研究 | 患者 > 对照组,任务 > 基线(过度激活) | 患者 < 对照组,任务 > 基线(低度激活) |
| 患者 < 对照组,任务 < 基线(去激活失败) | 患者 > 对照组,任务 < 基线(过度去激活) | |
| 双样本 VBM / FA 研究 | 患者 > 对照组(体积 / FA 更大) | 患者 < 对照组(体积 / FA 更小) |
表2:峰的t值符号。
如前所述,大多数基于体素的元分析方法采用峰值汇聚检验,该检验存在一些局限性,然后仅基于 p 值对证据进行二分类。
在本方案中,我们详细介绍了如何使用SDM-PSI进行基于体素的元分析,该方法具有多项优势,包括采用标准置换检验来评估效应的统计学显著性。此外,我们还展示了如何利用一组标准对证据强度进行分级,这些标准超越了仅依赖单一统计学显著性水平的二元分类方法。
为了便于复制本示例中的荟萃分析,我们提供了从先前荟萃分析中已从文献中提取出的数据。有趣的是,在该荟萃分析的论文中,证据“似乎”比我们使用更新方法所得到的证据更强。因此,我们建议对以往基于体素的荟萃分析中非系统性评估的证据应谨慎对待。
我们希望在遵循本方案后,神经影像学荟萃分析能够对神经影像学研究结果的证据提供更丰富、更精细的描述。
作者无任何利益冲突需要披露。
本工作得到了 Miguel Servet 研究合同 MS14/00041 和 2013–2016 年国家 I+D+i 计划的研究项目 PI14/00292、Carlos III 健康研究所科研评估与促进总局、欧洲区域发展基金(FEDER)以及 PFIS 博士前合同 FI16/00311 的资助。资助方在研究的设计与实施、数据的收集、管理、分析与解释、稿件的撰写、审阅或批准,以及决定提交稿件发表等过程中均未发挥任何作用。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 计算机 | 32 位或 64 位 Windows、64 位 Mac OSX 或 64 位 Linux |