正在阅读《第四节 马尔可夫链蒙特卡洛方法》
第(2/3)页


如何产生收敛于某个平稳分布π(x)的随机变量序列(马尔可夫链),是实现马尔可夫链蒙特卡洛方法的关键。根据马尔可夫链的建立原则,链中的各个观测值的产生过程应该是依次进行的,后一个观测值的产生只依赖于前一个观测值,而与更早产生的观测值无关。也就是说,在建立马尔可夫链过程中的任一时刻t(t≥0),序列中下一时刻t+1处的取值X(t+1)由条件分布P(x|X(t))决定,它只依赖于时刻t处的状态,而与时刻t以前的历史状态无关。因此,马尔可夫链的建立过程就是从一个时刻状态向下一个时刻状态转移的过程,其一步(次)转移概率函数为P(X(t+1)=x′|X(t)=x),该转移概率函数可以记为p(x,x′),称为该马尔可夫链的转移核,通过该转移核产生的且收敛的随机变量序列的分布即为平稳分布π(x)。因此,转移核的构造就成了实现马尔可夫链蒙特卡洛方法的关键,不同的马尔可夫链蒙特卡洛方法,其实质往往也就是构造转移核的方法不同。而在现在的实际应用研究中,大多数马尔可夫链蒙特卡洛方法在构造转移核时,关键就是要写出模型参数(变量)的满条件分布(fullconditionals)。所谓满条件分布,就是在条件分布中,所有的模型参数(变量)应该全部出现,或出现在条件中,或出现在变元中。在关于项目反应理论模型的相关研究文献中一般会提供各相关参数的满条件分布。

由于模型参数的特性是通过建立的马尔可夫链的样本进行估计的,而该马尔可夫链应该收敛到某个平稳分布π(x)上,因此,理论上,该平稳分布应该要服从模型参数的后验分布(posteriordistribution)。

另外,任何一个马尔可夫链的建立都需要有一个起始点,即初始值X(0)。虽然经过研究发现非常不同的初始值最终也能够收敛到平稳分布,但选择一个更加合理的初始值应该能够提高收敛的效果。接着的另外一个问题是,什么时候才能认为抽样过程已经收敛了呢?这虽然没有统一的标准检验判断方法,但是,从经验出发,如果产生的观测值足够多,如6000或更多,那么,它应该是收敛的。另外,也可以通过以下经简化的方法来进行验证,在建立马尔可夫链的过程中,每隔一段计算该段所有观测值的平均数,当所有平均数趋于稳定时,可以认为抽样过程收敛。

马尔可夫链的具体建立方法主要有以下几种:Gibbs抽样方法、Metropolis-Hastings抽样方法、Gibbs结合Metropolis-Hastings方法。

在项目反应理论研究和应用中,Gibbs方法就是从模型参数的满条件分布中直接进行反复迭代抽样,这个过程与项目反应理论中的联合极大似然估计过程类似,即一个参数估计(产生)值作为估计(产生)另一个参数的已知值,如此循环进行迭代更新,因此,Gibbs抽样是单元素抽样,因为每一步只从单个参数的满条件分布中进行抽样。然而,该方法的实现通常涉及计算正态化常数,使样本难以抽取,这个步骤使得该方法显得不够简洁便利。

三、马尔可夫链的使用

建立了服从某个平稳分布的马尔可夫链之后,通常以这个马尔可夫链的样本平均数作为参数的估计值。当然也可以通过这个马尔可夫链的样本来估计参数估计的标准差等。在实现马尔可夫链的过程中,开头的部分观测值会显得非常不稳定,因此,在实践研究当中,通常是以删除了马尔可夫链开头的部分观测值(burn-in)之后剩余的观测值作为研究样本。

下面还是以DINA模型为例,简单介绍利用马尔可夫链蒙特卡洛方法估计项目参数和被试属性掌握模式的过程。这里需要估计的参数为项目的猜测参数s、失误参数g,以及被试的属性掌握模式α。这些参数的联合后验分布可以表示如下:

相关小说

绝代仙医 都市 / 完结
绝代仙医
神话
张铁柱英雄救美,在水井里获得了龙血传承,从此医术逆天,成为了村里的致富带头人,过上了幸福的生活。
434.86万字 2 小时前
楼下的房客 推理 / 完结
楼下的房客
九把刀
一部悬疑而又诡异的出租屋的故事,详细描述了每个人的外型、特征、癖好,运用一个思想怪异的房东把他们串联成一个个恐怖的故事。讲述了继承出租公寓的无业游民,在东海大学附近贴起招租布告,他不要家世清白的乖乖牌学生、不要收支稳定的上班族、不要勤俭质朴的小家庭,...
9.32万字 4 小时前
逍遥小地主 历史 / 完结
逍遥小地主
堵上西楼
有幸穿越了,还是生在地主家,此生不缺吃穿却也不想混吃等死,所以傅小官随意的做了些事情,没料到产生的影响如此巨大。皇帝要让他官居一品,公主要招他为驸马,尚书府的千金非他不嫁,荒人要他的头,夷国要他的命,樊国要他的钱……可是,傅小官就想当个大地主啊!
303.7万字 3 小时前
山野村情 玄幻 / 完结
山野村情
水漾胭脂
桃酒村的酒很香,这没错。可是李阳实在不明白,为什么他的酒,就只醉女人!隔壁家的嫂子,城里来的美人村长,放假回来的俏丽村花……但凡跟李阳一起喝酒,就逢酒必醉,这些美女的酒量都这么差吗?而且醉就醉了,为啥还会往人身上扑呢?直到很久之后,李阳才懂得一个道理:“女人三分醉,演到你流泪!”
215.76万字 1 小时更新
荣华锦医 都市 / 完结
荣华锦医
苏羽宣
她救了他,他却亲手杀了她父亲。他说杀父之仇不共戴天,可现在,是她与他不共戴天!他独独留了她的性命,他说,他不杀手无缚鸡之力之人。他嫌她弱,那她就变强!她身怀医中瑰宝,一路来到皇城,入到七王爷府。皇权争斗,风起云涌,她卷身其中,更是查出,当年杀父一事,...
169.17万字 1 小时更新
暖婚蜜爱:景少宠妻至上 都市 / 完结
暖婚蜜爱:景少宠妻至上
壹元硬币
景珩说温尔是个地痞流氓,从见她第一天开始,他就在生气和悸动之间来回。 温尔说,景珩这张脸,想让人包养,藏起来。 外界不解,说景少是一朵鲜花插在了温尔这坨牛粪上,殊不知,某男正在百度怎么取悦女人…… 景先生把所有的温柔都给了温小姐,而温小姐把所有的撩汉技术都用在了景先生身上。 ...
73.3万字 1 小时更新