论文部分内容阅读
序贯蒙特卡罗(简称SMC)方法是一种估计非线性非高斯状态空间模型的常用方法,它通过逐维地产生高维蒙特卡罗样本来近似状态变量的后验分布并进行统计推断。随着计算机性能的提升,该方法已经被广泛地应用于计算机视觉,信号处理,生物信息学,金融等领域。与另一常用的采样方法—马尔科夫链蒙特卡罗方法(简称MCMC)相比,SMC具有不易陷入局部模式(mode)、可进行实时(online)估计和易于并行计算等优势。重抽样步骤是SMC方法中的一个重要环节,它用概率的手段复制权重大的样本,去掉权重小的样本,从而解决了序贯重要性抽样(简称SIS)方法中的权重退化问题。但是在模型信噪比很高的情况下,比如在动态随机一般均衡(简称DSGE)模型中,状态方程扰动项的方差通常远大于观测方程扰动项的方差,这种情况下使用一般的重抽样方法会导致样本的快速退化,从而降低估计精度。Fearnhead和Clifford(2003)指出在状态空间是离散的情形下,重复的样本会导致重复计算以及样本多样性退化的问题。因此他们提出一种新的重抽样方法,保证相同的粒子最多被抽取一次。但在状态空间是连续的情况下,由于新产生的样本不会完全相同,上述关于重复样本的结论并不成立。针对状态空间是连续的情形,我们提出在重抽样中限制每个粒子最多被抽取R次(R>1)以兼顾被抽取样本的多样性需求(即拥有大量不同的样本)和聚焦需求(即保证权重大的样本多出现)。数值实验表明在高信噪比的情况下,新的重抽样方法可以减小状态变量的估计误差以及提高似然函数的估计精度。本文还将新的重抽样方法应用到粒子马尔科夫链蒙特卡罗(简称PMCMC)算法中,对一个DSGE模型的参数进行贝叶斯估计,数值结果表明使用新的重抽样方法可以提高参数的估计效率。