基于bold信号的事件相关设计.docxVIP

基于bold信号的事件相关设计 近年来，功能磁共振成像在方法学中取得了重大进展。换句话说，与活动相关的项目的不断发展和广泛使用。同传统的区块设计相比,事件相关设计在灵活安排实验刺激,准确描述血氧水平依赖(blood oxygen level dependent, BOLD)信号特征和同其他手段结合等方面都具有很明显的方法学优势,从而在功能磁共振成像的研究中得到了日益广泛的运用,并产生了很多具有重要意义的发现。在本文中,我们拟就事件相关功能磁共振成像的基本原理及其设计优化做一个初步的探讨。 1 事件相关设计的基本类型 迄今为止,虽然事件相关设计已经发展出来多个子类型,但总的说来,从其涉及的内在原理的差异上讲,事件相关设计可以大致分为慢速呈现事件相关设计和快速呈现事件相关设计两大类型。下面我们将分别阐述这两类设计的基本原理。 1.1 bold的成像特点 典型的慢速呈现事件相关设计可以用Fig 1a表示,这种设计是早期事件相关设计所采用的主要形式,在构思上与ERPs有很多类似的地方。慢速呈现设计的基本特征是刺激呈现的事件间隔(SOA)相对较长,通常在12 s或者以上。由于当前功能磁共振成像研究主要采用BOLD的成像方式,而单个事件所诱发的典型的BOLD信号持续的时间大概在12~16 s。因此,在慢速事件相关设计里,可以粗略地忽略前后刺激所引发的信号的重叠,也不考虑BOLD信号的饱和与叠加的线性等问题,而可以直接采用时间锁定(time-lock)的方法,分离出单个刺激所导致的信号变化。这样的设计方式相对比较直观,也比较容易理解。 1.2 间隔相对较短 典型的快速呈现事件相关设计的方式可以用Fig 1b表示。其特点就是刺激呈现的事件间隔相对较短(不需要间隔12 s或者以上),可以在2~3 s,甚至1 s以下。因此,在快速呈现刺激的条件下,必然要考虑前后两个刺激引发的BOLD信号的叠加和相互干扰问题,这就构成了快速呈现事件相关设计和慢速呈现事件相关设计在原理上的本质差异。 1.3 加速事件相关设计的bold信号检验方法 在统计思路上,事件相关实验设计也与区块设计有着本质的区别。区块设计主要是通过检验实验任务区块和控制区块的信号的差异,以确定大脑的激活区域。在这里面,统计的元素是一个区块所引起的信号之和。而事件相关设计分析的是单个刺激的信号的特征。对于慢速呈现事件相关设计来说,首先要根据刺激的类型选择性地进行平均,以此代表一类刺激所引起BOLD信号。然后通过对分离出来的BOLD信号进行统计检验,通过其与基线,或者参考方程(如泊松分布方程)的差异来确定显著激活的脑区和激活水平。同时,也可以直接比较两类刺激诱发的BOLD信号的差异来确定它们在某个区域的激活是否有显著差异。 在对快速呈现事件相关设计所得的数据进行统计分析时,第一步就是要从混合的BOLD信号中分离出单个刺激诱发的信号。当前一个常用的方法是在一般线性模型(GLM)的基础上采用去卷积(deconvolve)算法来实现。显然,这种方法的一个前提是BOLD信号具有线性叠加性。分离出单个事件所对应的BOLD信号以后,其余的统计分析思路就与慢速呈现事件相关设计非常类似。 以上所介绍的方法主要是为了更准确地描述刺激所诱发的BOLD信号的特征,这也是事件相关设计的一个重要目的和优势。另一个方面,如果我们只是关心哪些地方有显著的激活,则可以采用不同的统计方法。目前使用较广的一个方法是假设一个理想的BOLD信号,然后在一般线性模型的基础上与刺激序列做卷积运算,形成一个“理想波”。最后计算实际信号和理想波的拟合程度,从而得到某个区域的激活强度。这种方法可以同时运用于区块设计,慢速事件相关设计和快速事件相关设计。其中的差异在于对于区块设计来说,处理的基本单元是一个区块,而事件相关设计的基本单元是单个事件。同时,快速事件相关设计同样也是以线性叠加为前提的。 2 bold信号基本特征 从事件相关的基本原理可以看到,无论是快速呈现设计还是慢速呈现设计,都是建立在对单个事件所诱发的BOLD信号的分析的基础之上的。它们也自然包括了对BOLD信号的一些基本特征的了解,包括BOLD信号的时间进程、稳定性、一致性和线性叠加性等等。而这些特征都将对实验设计的效果产生非常重要的影响。因此,在讨论事件相关设计的优化之前,有必要对血液动力学反应和BOLD信号的基本特征做一个简要的分析。 2.1 bold信号的时间进程特征 当前研究发现,同神经元电冲动信号相比,BOLD信号是一个相对很慢的反应过程。单个刺激所诱发的BOLD信号大约在1~2 s后逐渐上升,在5~6 s达到顶点。然后在10 s左右回到基线水平。同时,还有研究者发现在信号上升之前会有较小幅度的低于基线水平的下冲(pre-undershoot)。在回到基线水平后,也会存