脑成像的元分析研究：基于脑成像的元分析研究：基于坐标点坐标点的元分析方法的元分析方法 报告人：封春亮报告人：封春亮 单位：北京师范大学，脑与认知科学研究院单位：北京师范大学，脑与认知科学研究院 脑成像元分析的背景简介 脑成像数据元分析的发展 元分析 (Glass, 1976) “将现有研究整合并对其结果进行综合分析的统计方法” 脑成像数据元分析的发展 (胡传鹏等, 2015) 脑成像研究的数量 脑成像数据元分析的发展 元分析对脑成像研究的意义 (Wager et al., 2007; Eickhoff et al., 2009) v 单个脑成像研究的被试数量普遍较少，结果不够稳定； v 单个脑成像的结果很可能受到特定的实验操作（如，扫描参数）的影响； v 单个脑成像研究对某个脑区功能的解释往往局限于使用的单个或几个任务。 (e.g., vmPFC=self-reflection?) -单个脑成像研究的局限性 (Wager et al., 2007) 脑成像数据元分析的发展 脑成像技术：where? (位置信息) 坐标信息 (x, y, z) 最大激活点的坐标 (x, y, z) 脑成像数据元分析的发展 脑成像的元分析：对不同研究中位置信息的整合 早期脑成像数据的 “元分析”：对不同研究中报告的 位置信息的主观描绘 (Peyron, 2000)(Joseph, 2000) Activation likelihood estimation (激活似然性估计，ALE, Turkeltaub et al., 2002) 目标： 最大化地量化研究间位置信息的一致性；尽可能减小分析方法的主观 性。 脑成像数据元分析的发展 主观的描述 客观的统计 ALE的计算 (Turkeltaub et al., 2002) 概念：把每个激活点拟合成一个概率分布 也即这个激活点落在大脑每个voxel上的概率大小：与激活点越近的voxel，概率越大。 脑成像数据元分析的发展 文献中报告的激活点 一个激活点落在大脑每个voxel上的概率。 （d为每个voxel与激活点之间的距离；为分布的方差） ALE的计算 (Turkeltaub et al., 2002) 脑成像数据元分析的发展 ALE：至少有一个激活点落在某个voxel内的概率。 ALE= 第i个激活点落在某个voxel内的概率。 一个激活点落在大脑每个voxel上的概 率。 相关研究中报告的所有激活点 ALE的统计检验 (Turkeltaub et al., 2002) Permutation test (置换检验) v 零假设：激活点没有重合，所有重合都是随机因素导致的 v 假设元分析采集到了100个坐标（激活点）；则产生100个随机坐标，然后根据这些随机 坐标计算出大脑中每个voxel在零假设下的ALE值（如，重复1000次，然后求均值） 产生零假设分布； 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 具备显著性检验，得到的结论有统计上的支持 脑成像数据元分析的发展 早期ALE算法的缺点： 以激活点为单位，在每个voxel上，求出所有激活点对应概率的联合概率（即ALE值）。 v 忽略了研究间的变异，属于固定效应的统计推断 (fixed-effects inference): 元分析的结论不能 推广到该元分析未包含的研究 (胡传鹏等, 2015; Eickhoff et al., 2009)。 v 报告了较多激活点的研究对元分析的影响更大，甚至主导元分析的结果 (Wager et al., 2007) 。 脑成像数据元分析的发展 ALE的改进算法 (Eickhoff et al., 2009) 脑成像数据元分析的发展 ALE/MA= Modeled activation (MA) map 某个voxel在特定一个实验中激活的可能性 第一步：在每个实验的内部先做一次整合 类比：把每个实验看做一个“被试”；每个voxel看做一个实验条件；每个实 验报告出来的激活点看做一个采集到的数据点。 第一步的整合 “” 对每个被试各个条件下的所有数据点求和。 第一步整合的结果：每个实验在每个voxel上，都有只一个MA/ALE值。 脑成像数据元分析的发展 Modeled activation (MA) map 跨实验的整合 第二步：以实验 (“被试”)为单位，计算每个voxel上的总ALE值，并做统计检 验。 结果：每个voxel的总ALE值 统计检验 ：零假设：元分析中各实验的MA map之间没有重合，所有重合均是 随机因素导致的。 置换检验： 零分布：从每个实验的MA图中各自随机抽一个坐标点和对应的 MA值，根据这些MA值算出一个总ALE值。 重复1011次，每次能得到一个总ALE值。构建出零分布。根据不同 ALE值在零分布中出现的概率，确定显著性水平。 ALE的改进算法 (Eickhoff et al., 2009) ALE算法进一步的优化： 脑成像数据元分析的发展 Eickhoff et al., 2012 (NeuroImage): 以MA值为单位代替以voxel为单位，减少了迭代次数，加快了运算速度； Turkeltaub et al., 2012 (Human Brain Mapping): 取最大值代替求和的方法，进一步减少单个实验对结果的影响。 脑成像元分析的步骤 脑成像元分析的步骤 脑成像元分析的步骤 1. GingerALE: https:/www.brainmap.org/ale/ (邓沁丽, 2015) 脑成像元分析的步骤 使用到的软件： 脑成像元分析的步骤 使用到的软件： 2. Mango: rii.uthscsa.edu/mango 用于结果呈现 第一步：选好主题，查找和筛选文献 1.1. 选好元分析的主题 v 工作记忆 (Nee et al., 2013); v Go/No-go 任务 (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); “共性共性” 脑成像元分析的步骤 单个任务对应的神经机制：多个任务共同的神经机制 ： 第一步：选好主题，查找和筛选文献 1.2. 查找和筛选文献 v 查找途径： google学术 ( http:/scholar.glgoo.org/) PubMed (http:/www.ncbi.nlm.nih.gov/pubmed/) ISI Web of Science ( 已查找文献的引用列表 相关领域综述的引用列表 v 文献筛选： 被试群体：如，正常成人被试； 是否使用fMRI技术； 使用的是全脑(whole-brain)分析，而不是感兴趣区(ROI)分析； 报告了感兴趣的统计结果； 激活坐标报告在标准空间里 (如MNI，Talairach). 脑成像元分析的步骤 v 查找内容： 与主题相关的几个关键词 fMRI/functional fairness/fair/ultimatum game 第二步：采集坐标点 v 需要提供的信息（规定的格式）： 标准空间的类型（第一行）：MNI或者Talairach; 文献信息： / 作者名，年份等； / Subjects=N (被试量的信息) 坐标点。每个点一行 不同文献之间需要有一个空行； 注意MNI与Talairach空间之间的转换。 脑成像元分析的步骤 第三步：使用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. 矫正结果 关于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, magnetic resonance imaging, neuroimaging 文献编码 ： 基于坐标的元分析 v 任务一：UG_unfair_vs_fiar_main_effect.txt v 任务二：TG_trustee_decision_making.txt P 0.05, FDR-corrected at the voxel level; Mini. Cluster size = 300 mm3 基于坐标的元分析 结果： 写文章： 基于坐标的元分析