脑成像数据的元ALE原理与操作

1、脑成像元分析的背景简介脑成像数据元分析的发展元分析元分析 (Glass, 1976)“将现有研究整合并对其结果进行综合分析的统计方法统计方法”脑成像数据元分析的发展(胡传鹏等, 2015)脑成像研究的数量脑成像研究的数量脑成像数据元分析的发展元分析对脑成像研究的意义元分析对脑成像研究的意义 (Wager et al., 2007; Eickhoff et al., 2009)v 单个脑成像研究的被试数量普遍较少，结果不够稳定；v 单个脑成像的结果很可能受到特定的实验操作（如，扫描参数）的影响；v 单个脑成像研究对某个脑区功能的解释往往局限于使用的单个或几个任务。 (e.g., vmPFC=se

2、lf-reflection?)-单个脑成像研究的局限性单个脑成像研究的局限性(Wager et al., 2007)脑成像数据元分析的发展脑成像技术：脑成像技术：where? (位置信息位置信息) 坐标信息 (x, y, z)最大激活点的坐标(x, y, z)脑成像数据元分析的发展脑成像的元分析：对不同研究中位置信息的整合脑成像的元分析：对不同研究中位置信息的整合早期脑成像数据的早期脑成像数据的 “元分析元分析”：对不同研究中报告的：对不同研究中报告的位置信息的主观描绘位置信息的主观描绘(Peyron, 2000)(Joseph, 2000)Activation likelihood esti

3、mation (激活似然性估计，激活似然性估计，ALE, Turkeltaub et al., 2002)目标：目标：最大化地量化量化研究间位置信息的一致性；尽可能减小分析方法的主观性。脑成像数据元分析的发展主观的描述主观的描述客观的统计客观的统计ALE的计算的计算 (Turkeltaub et al., 2002)概念：概念：把每个激活点拟合成一个概率分布 也即这个激活点落在大脑每个voxel上的概率大小：与激活点越近的voxel，概率越大。脑成像数据元分析的发展文献中报告的激活点一个激活点落在大脑每个voxel上的概率。（d为每个voxel与激活点之间的距离；为分布的方差）ALE的计算的计

4、算 (Turkeltaub et al., 2002)脑成像数据元分析的发展ALE：至少有一个激活点落在某个voxel内的概率。ALE=第i个激活点落在某个voxel内的概率。一个激活点落在大脑每个voxel上的概率。相关研究中报告的所有激活点ALE的统计检验的统计检验 (Turkeltaub et al., 2002)Permutation test (置换检验)v 零假设：激活点没有重合，所有重合都是随机因素导致的v 假设元分析采集到了100个坐标（激活点）；则产生100个随机坐标，然后根据这些随机坐标计算出大脑中每个voxel在零假设下的ALE值（如，重复1000次，然后求均值）产生零假

5、设分布；v 通过比较实际ALE值与随机坐标产生的ALE值，来确定显著性。脑成像数据元分析的发展ALE的统计检验的统计检验 (Turkeltaub et al., 2002)Permutation test (置换检验)在H0分布下，ALE为5.7210-3的概率为0.0001。因此，如果设置=0.0001为显著水平，那ALE值大于5.7210-3的voxel都视为显著（也即至少有一个激活点落在这个voxel内）。脑成像数据元分析的发展ALE的优点的优点 (Turkeltaub et al., 2002)v 自动化的分析v 除了位置信息外，也有对一致性水平的量化 (ALE值)v 具备显著性检验，

6、得到的结论有统计上的支持脑成像数据元分析的发展早期早期ALE算法的缺点：算法的缺点：以激活点为单位，在每个以激活点为单位，在每个voxel上，求出所有激活点对应概率上，求出所有激活点对应概率的的联合概率（即联合概率（即ALE值）。值）。v 忽略了研究间的变异，属于固定效应的统计推断 (fixed-effects inference): 元分析的结论不能推广到该元分析未包含的研究 (胡传鹏等, 2015; Eickhoff et al., 2009)。v 报告了较多激活点的研究对元分析的影响更大，甚至主导元分析的结果 (Wager et al., 2007)。脑成像数据元分析的发展ALE的改进算

7、法的改进算法 (Eickhoff et al., 2009)脑成像数据元分析的发展ALE/MA=Modeled activation (MA) map某个某个voxel在特定一个实验中激活的可能在特定一个实验中激活的可能性性第一第一步：在每个实验的内部先做一次整合步：在每个实验的内部先做一次整合类比类比：把每个实验看做一个“被试”；每个voxel看做一个实验条件；每个实验报告出来的激活点看做一个采集到的数据点。第一步的整合 “” 对每个被试各个条件下的所有数据点求和。第一步整合的结果：每个实验在每个第一步整合的结果：每个实验在每个voxel上，都有只一个上，都有只一个MA/ALE值。值。脑成像

8、数据元分析的发展Modeled activation (MA) map跨实验的整合跨实验的整合第二步：以实验第二步：以实验 (“被试被试”)为单位，计算每个为单位，计算每个voxel上的总上的总ALE值，并做统计检验。值，并做统计检验。结果：每个结果：每个voxel的总的总ALE值值统计检验：统计检验：零假设：元分析中各实验的零假设：元分析中各实验的MA map之间没有重合，所有重合均之间没有重合，所有重合均是随机因素导致的。是随机因素导致的。置换检验：置换检验：零零分布分布：从每个实验的：从每个实验的MA图中各自随机抽一个坐标点和对应的图中各自随机抽一个坐标点和对应的MA值，根据这些值，根据

9、这些MA值算出一个总值算出一个总ALE值。值。重复重复1011次，每次能得到一个总次，每次能得到一个总ALE值。构建出零分布。根据不同值。构建出零分布。根据不同ALE值在零分布中出现的概率，确定显著性水平。值在零分布中出现的概率，确定显著性水平。ALE的改进算法的改进算法 (Eickhoff et al., 2009)ALE算法进一步的优化：算法进一步的优化：脑成像数据元分析的发展Eickhoff et al., 2012 (NeuroImage): 以MA值为单位代替以voxel为单位，减少了迭代次数，加快了运算速度；Turkeltaub et al., 2012 (Human Brain

10、Mapping): 取最大值代替求和的方法，进一步减少单个实验对结果的影响。 脑成像元分析的步骤脑成像元分析的步骤脑成像元分析的步骤1. GingerALE: /ale/ (邓沁丽, 2015)脑成像元分析的步骤使用到的软件：使用到的软件：脑成像元分析的步骤使用到的软件：使用到的软件：2. Mango: /mango用于结果呈现用于结果呈现第一步：选好主题，查找和筛选文献第一步：选好主题，查找和筛选文献1.1. 选好元分析的主题选好元分析的主题v 工作记忆 (Nee et al., 2013);v Go/No-go 任

11、务 (2008);v 奖赏加工 (Bartra et al., 2013);v 社会决策 (Feng et al., 2015).v 共情与自身疼痛感受是否 基于共同的神经机制 (Lamm et al., 2011);v 自传体记忆、展望未来、心理理论以及默认网络是否基于共同的神经机制 (Spreng et al., 2009);v 初级和次级奖赏的加工是否基于共同的神经机制 (Sescousse et al., 2013)v 自己接受奖赏与观看别人接受奖赏是否基于共同的神经机制 (Morelli et al., 2015);脑成像元分析的步骤单个任务对应的神经机制：单个任务对应的神经机制：多

12、个任务共同的神经机制：多个任务共同的神经机制：第一步：选好主题，查找和筛选文献第一步：选好主题，查找和筛选文献1.2. 查找和筛选文献查找和筛选文献v 查找途径：查找途径： google学术 (https:/ /) PubMed (/pubmed/) ISI Web of Science (http:/ 已查找文献的引用列表 相关领域综述的引用列表v 文献筛选：文献筛选： 被试群体：如，正常成人被试； 是否使用fMRI技术； 使用的是全脑(whole-brain)分析，而不是感兴趣区(ROI)分析

13、； 报告了感兴趣的统计结果； 激活坐标报告在标准空间里 (如MNI，Talairach).脑成像元分析的步骤v 查找查找内容内容： 与主题相关的几个关键词 fMRI/functional fairness/fair/ultimatum game第二步：采集坐标点第二步：采集坐标点v 需要提供的信息（规定的格式）：需要提供的信息（规定的格式）： 标准空间的类型（第一行）：MNI或者Talairach; 文献信息： / 作者名，年份等； / Subjects=N (被试量的信息) 坐标点。每个点一行 不同文献之间需要有一个空行； 注意MNI与Talairach空间之间的转换。脑成像元分析的步骤第三

14、步：使用第三步：使用GingerALE进行元分析进行元分析3.1. 计算计算ALE (activation-likelihood-estimation)脑成像元分析的步骤v Eickhoff et al., 2012 (NeuroImage): v Turkeltaub et al., 2012 (Human Brain Mapping): 内置在GingerALE软件里的算法第三步：使用第三步：使用GingerALE进行元分析进行元分析3.2. 矫正结果：矫正结果： Settings脑成像元分析的步骤第三步：使用第三步：使用GingerALE进行元分析进行元分析3.2. 矫正结果矫正结果关于

15、FDR的bug (在2.3.3以后的版本中更正)FDR矫正的原理：矫正的原理： 假设有N个比较，对应有N个p值。 i). 所有p值从小到大排序：p(1), p(2), p(3).p(N)；ii). 若想控制FDR不能超过q，则只需找到最大的正整数i，使得 p(i) 公平对待公平对待 信任游戏信任游戏 (trust game)的被投资者的被投资者: 返还决定时的脑活动返还决定时的脑活动文献查找：文献查找： Google学术学术, PubMed等等基于坐标的元分析关键词：关键词：normative decision making, fair, altruistic punishment, ultimatum game, trust, trust game, investment game, trustor, investor, trustee, trustworthiness, reciprocity, fMRI