脑成像数据的元分析-ALE原理与操作.pptx

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