序贯重要采样和重采样在动态投资组合信用风险中.doc

《序贯重要采样和重采样在动态投资组合信用风险中.doc》由会员分享，可在线阅读，更多相关《序贯重要采样和重采样在动态投资组合信用风险中.doc（8页珍藏版）》请在三一办公上搜索。

1、序贯重要采样和重采样在动态投资组合信用风险中的应用摘要：我们提出一种序贯蒙特卡洛方法，用于估计在动态、基于强度点过程模型中的投资组合信用风险的极少事件概率。这种方法是基于测度的改变，并涉及一种重采样机制。我们确定重采样权重，使得在技术条件下，能够得到对巨额组合损失的概率的对数有效模拟估计。通过一种数值分析说明了这种序贯蒙特卡洛方法的特征，并与其他最新用于研究组合信用风险的极少事件方法进行了比较，这些极少事件方法包括交互粒子方法和重要采样方法。1、Introduction组合信用风险是一种由于像贷款和公司债券这样信用敏感资产组合的违约而导致的财务损失的分布。蒙特卡洛模拟方法就被广泛用于估计组合损

2、失的分布。该方法几乎应用于任何与违约期限相关的模型中并且相对容易实现；但另一方面，精确估计组合巨额损失概率的计算工作量可能是巨大的。这个概率在风险管理应用中处于中心地位；比如，在风险价值中风险度量的估计。本文描述和分析了一种序贯蒙特卡洛方法，对组合巨额风险和其他极少事件的概率做出有效和无偏估计。该方法应用于违约期限相关的动态点过程模型中，在这些广泛应用的模型中，组合中的一种资产的违约受随机强度过程控制。强度过程与组合中的各种资产相关联，以反映组合的违约依附结构。类似于传统重要采样方法，序贯蒙特卡洛方法涉及概率测度的转变。违约事件是按照不同于标准测度的概率测度序贯采样的。另外，连续产生的样本路径

3、需要使用一组状态依赖权重重新采样。重采样机制区分了序贯重要采样和重采样方法和现存的序贯重要采样方法。我们按渐进最优方式选择重采样权重，提供条件保证按SISR方法产生的组合巨额损失的概率的估计为对数有效的。For personal use only in study and research; not for commercial useSISR方法涉及到由Moral和Garnier（2005）研究的交互粒子系统方法，该方法被Carmona和Crepey(2010)、Carmona (2009)等人用于估计组合巨额损失的概率。IPS方法在标准测度下按序贯进行并包含了一个重采样机制。所谓重采样机制

4、就是在路径空间上测度转变。重采样机制的有效性及由IPS方法产生的估计量的性能的高低，很大程度上取决于指定的重采样权重的参数的明智选择上，这个参数的渐进最优选择方法在文献中尚未解决。因此，点对点方法被用于决定这个参数。我们所采用的SISR方法排除了选择一个参数的需要，而是产生一个可证明有效的极少事件估计量。而且，SISR方法允许人们使用一种不同于标准测度的采样测度，这种方法能使采样更方便并可能导致额外的方差减少。数值经验证明了应用于组合信用风险的自激模型的SISR算法的性能良好，组合中资产的违约强度追随相关跳跃扩散过程。在给定计算预算条件下，SISR算法比IPS方法能更精确的估计巨额损失的概率。

5、而且，SISR方法对非常小的概率产生有意义的估计。如果资产组合是相对多样化的，SISR方法也能超越对数效率IS方案。尽管SISR方法很早就用于发生在非线性滤波中的复杂多维的分布的样本中，但是该方法在极少事件仿真中的应用也只在Chan and Lai(2011)的文章中出现过。他们给出了一个针对一般SISR估计量和渐进方差的一致估计的CLT（计算机语言翻译程序）。针对经典的大偏差设定，Chan and Lai(2011)展示如何选择重采样权重，以便获得伴随有限短随机游走过程中的某种极少事件概率的对数有效估计量。在本文中，我们在多变量、基于强度点过程的设定下提出重采样权重用于构建相关事件到来的模型

6、。建立在Chan and Lai(2011)的文章中有关重采样权重讨论的基础上，我们发展了条件保证这些权重能产生巨额损失概率的对数有效估计量。For personal use only in study and research; not for commercial use文章其他部分的安排如下：第二部分阐述投资组合信用风险问题。第三部分描述了一个基本的SISR算法，第四部分分析了重采样权重的渐进最优选择。第五部分描述了伴随偶尔重采样的一个延伸的SISR算法。第六部分提供了数值结果。第七部分总结。2.动态投资组合信用风险考虑一个拥有n家公司股票的投资组合，这些公司易遭受违约风险。这些公司的随

7、机违约时间以几乎完全不同的到期时间被引入到模型中，这些时间变量被定义在右连续和完全信息过滤的完全测度空间上（f,P）。在风险管理应用中，P是统计概率，而在衍生品定价应用中，P是风险中性定价测度。与相关的的是过程指标，其中=()，在这里是一个事件的指示函数（数学中，指示函数是定义在某集合X上的函数，表示其中有哪些元素属于某一子集A）。对每一个，都会有一个严格正的、可积的、逐步衡量的过程使得随机变量 （1）形成一个鞅。过程代表条件违约率或者公司I的强度。被认为是给定的相关的随机过程。之间的相关性反映了投资组合个资产之间的违约依赖结构。大量的关于的详述在文献中都能看到，比如，见谁谁。与投资组合相联系

8、的的信用风险被描述成组合损失，这里为违约指标向量，是损失头寸向量。我们首要目标是损失分布的尾部分布，这就代表了非典型巨额违约损失的概率。这些极少时间概率在投资组合风险管理和其他应用中处在中心的位置。比如，这些概率被用来估计风险组合测度，像高置信水平下的VaR。在固定水平T0上计算和的分布问题可被认为是马尔科夫链问题（该过程中，在给定当前知识或信息的情况下，过去（即当期以前的历史状态）对于预测将来（即当期以后的未来状态）是无关的。P(Xn+1=J |X0=I0, X1=I1, Xn=In) = P(Xn+1=J| Xn=In)。Giesecke et al(2010)中的命题3陈述：存在一个连续

9、时间马尔科夫链使得对固定的和所有的都有成立。模拟链自身各要素之间没有联合过渡，同时在0处开始并且转换密度为，这里, （2）期望方程2可以被的很多标准模型计算出来，参见Giesecke et al(2010)。马尔科夫链的存在将计算分布问题归纳到计算上来。同样的，在每一个从固定分布中得出的独立于的变量已知的假设条件下，的存在将计算分布的的问题归纳到计算的分布上，这里=。这种归纳转换时非常重要的，它允许我们分析一个马尔科夫链模型，而不是分析一个潜在复杂的点过程模型（通过把点过程转换到马尔科夫过程，大大降低运算的难度）。跳跃过程本身不是一个马尔科夫链，因为跳跃时间有一个形式的强度。然而，的分布可以由

10、的分布获得，的分布能通过解前向柯尔莫哥洛夫方程获得。然而这种方式是不切实际的，这是因为的状态空间的维度在实际应用中将趋向于非常大。在定价问题中，投资组合在100到125之间，在风险管理设置中，可能是更大的。我们提出使用一种蒙特卡洛模拟的方法估计的分布，的转换可以使用稀释方法来抽样，这能导致分布的无偏模拟估计量。然而，为了获得精确的尾部分布的统计，（这是我们关注得中心），我们需要做大量的复制工作（即多次重复抽样）。我们发展了一种序贯模拟方法，该方法能极大的减少因获得精确的尾部分布的概率而需要的大量模拟试验。我们集中考虑的情况，后面我们会怎么对待一般的情况。3.序贯重要采样和重采样这一部分描述了用

11、来有效估计分布尾部概率的SISR方法。3.1回顾让，这里是跳跃过程的第个到达的时间，=。并且，让,。是在上的离散时间马尔科夫链。我们用来表示的P转换概率。假设对一些合适的设置，感兴趣的极少事件采用的形式，这里及接下来的叙述中，我们让 。我们估计 的算法结合了序贯重要采样和重采样。估计量采用 （3）的乘积形式，在这里是一个重要测量，是一个采样测量 ，两者都定义在流域。假设在条件下,是一个伴有跃迁密度的马尔科夫链。那么（3）第一个似然比为 （4） 按照 顺序取样而改变。已经产生的样本路径使用状态依赖权重被重新采样；重采样机制努力使（3）中的似然比 有效。尽管不能排除的情况，但是在一个不同于的条件下

12、采样可能是方便的，同时可能导致额外方差减少。我们会在6用数字情况的背景来说明。3.2基本算法我们用一般重采样权重函数来描述基本SISR算法，不一定试图模仿方程（3）中的 。在算法1中总结了操作步骤。对于 ，让，这里是一个整数。 对每一个阶段，我们从概率密度中抽取个独立变量，形成 个样本路径或粒子。从这些过程中，我们给出粒子并使用定义在标准化的权重的概率，该权重为 （5） 对我们在阶段停止重采样，的基本的SISR估计量 如下 （6）在这里， 并且 （7） = （8） 算法1是一个整数。开始 到 for doFor alldo Generate from .set .compute end for

13、ifthen Resample using weight to draw end forreturn estimator 在估计量（6）中， 这一项是标准测度和一个粒子的采样测度之间的似然比。它解释了我们按测度而不是来抽样的事实。 这一项解释了在阶段 的重采样。重采样权重被选择的目的是使 接近于目标重要测度的似然比。在这种情况下，的分母接近于目标重要测度 对的似然比，而 的分子接近于这些权重所构建的样本均值这一标准化常数。具体来说，定义权重 ，由（9）式递归 （9） 当在和下是一个马尔科夫链，（9）的右端归纳为 ，其中 是在下的的转换密度， 是在 的转换密度。重采样机制消除了直接从采样的需要。

14、使用接近这个似然比，我们有 （10） （11）于是，SISR被用来接近重要测度 ，这种测度是经由的序贯IS方法得出，这对应于（10）中的 ，而重采样对应于 。注意在阶段 处没有重采样。3.3中心极限定理及方差估计量因为在每一个阶段进行重采样，SISR粒子不在是独立的，意识到这一点很重要。因此，针对独立试验的CLT不在有效，Chan and Lai(2011)提出一个CLT用于由SISR方法产生的估计量和渐进方差的一致估计量。这些结果应用于任何权重函数的选择，同时能方便SISR估计量（6）中的分析。一个粒子的“原始祖先”扮演了一个重要角色。在重采样之前，粒子的第一代是。如果第一个同等级的是，比如

15、如果是的原型，就设=。用表示从中产生的的复制数以形成第代的粒子。还有，用表示标准化增量的权重（式（5）中的第代的第个粒子。Chan and Lai(2011)表明了了的SISR估计量满足 当， （12）同时渐进方差可由（13）式得出一致估计式（12）和（13）的导出是基于的鞅表示，这在附录中对的情况进行概述。这个鞅表示与4中的渐进有效性的分析相关。同时也暗含了是无偏的。最后，这个鞅表示表明每一个重采样对的方差贡献一个时期。因此，应避免不必要的重采样，比如，只有当重采样是有效的时候才进行。直觉是党（5）中相关权重的直方图是相对平坦的，比如当这些权重是非常小的时候，重采样是无效的。在那种情况下，重

16、采样产生额外的方差而没有重大改变估计量。我们将在5中回到这个话题。3.4SISR 与IPSSISR方法是与由Del Moral 和Garnier(2005)发展的IPS方法相关联的。尽管从不同的视角设计，但是在抽样测量和相等以及重采样权重采取特殊的方式的情况下，IPS方法可以被看做SISR方法。为了估计，人们经常用作为某个参数。尽管有Del Moral 和Garnier(2005)发展的IPS方法在原则上能用在其他方面，但是这个构想被关于投资组合信用风险的文献广泛采用。IPS倾斜参数在整个阶段中保持固定，为了确保算法的良好性能，必须采取的形式。的渐进最优选择在文献中并没有给出。因此需要用特设的

17、方法来决定。一种方式是选择使在给定的条件下，估计量的相对错位最小化的。相应的最优化就是基于大量的辅助估计的经验，即测试一系列的建议值。而且，所获得的的值有不稳定的趋势，因为它们受到估计量方差的影响。因为这个原因，获得极端尾部分布概率的良好性能的估计将会非常困难。在下一部分，我们提出针对SISR算法1的重采样权重，SISR算法能导致被证明是有效的的估计量，表达式如上式（6）。这些权重函数不要求选择调整参数。仅供个人用于学习、研究；不得用于商业用途。For personal use only in study and research; not for commercial use.Nur fr den persnlichen fr Studien, Forschung, zu kommerziellen Zwecken verwendet werden.Pour l tude et la recherche uniquement des fins personnelles; pas des fins commerciales. , , . 以下无正文