该检验可表述为标准似然比检验，其中对数似然比渐近分布为χ（1）LRUC=-2 lnhLUC（p；I（δ）（si；p），I（δ）nobs（si；p））/LUC（π；I（δ）（si；p），I（δ）nobs（si；p））I~ χ(1).我们在EUR3M曲线上执行相同的测试。为了产生样本外测试程序所需的置信区间，我们使用两种方法：第2.3.1节概述的高斯扩散模型和第2.3.2节解释的bootstrapmethodology。然后我们比较这两种设置的测试结果。对于高斯扩散模型，表4、6和8中报告了EONIA瞬时前向曲线的无条件覆盖测试结果，而表5、7和9包含通过自举程序获得的结果。符号(*)及(**)分别对应于95%和99%的统计显著性。在图15和图16中，我们展示了三个月预测置信区间（蓝线）与实际利率（黑线）的时间演变，以定位负异常和正异常（分别为红点和交叉点）。我们选择短期到期利率f（t，1m）和长期到期利率f（t，5y）。在第一行中，我们使用95%的覆盖概率来比较高斯扩散模型（左面板）和自举法（右面板），而第二行显示了99%覆盖概率的结果，再次比较高斯扩散（左面板）和自举法（右面板）。对于EUR3M FRA期限结构，高斯扩散模型的无条件覆盖测试结果见表10、12和14，而表11、13和15包含通过自举程序获得的结果。

最后，在图17和图18中，我们报告了短期和长期EUR3M曲线（即FRA3M（t，3m）FRA3M（t，5y））的三个月预测置信区间的时间演变。就EONIA曲线而言，在短期预测期（一周）内，测试结果非常令人满意，两种方法（多变量和bootstrap）的覆盖率均为95%。当考虑更高的覆盖概率值时，即p=99%，引导法进一步提高了测试质量。这一事实也适用于更长的预测期——三个月零一年——即使结果总体恶化。通过观察已实现的时间序列与预测分布平均值之间的关系，可以猜测这种行为背后的原因。例如，请参阅图15的左上面板。从2008年10月到2009年4月，欧洲央行将基准短期利率下调了五倍。EONIA曲线的短期部分始终跟随利率水平的快速下降。然而，尽管历史时间序列从显著上升的趋势突然转变为显著下降的趋势，但预测率不会立即遵循相同的行为。预测需要一段时间才能恢复漂移，即风险溢价上升到较大的正值，并提供强烈的负修正。实际上，从图6中我们可以看出，短期风险溢价在2009年最后一个季度之前不会增加。由于利率动态的突然变化，风险溢价的延迟调整是大多数预测失败的原因。可以很容易理解的是，预测期越长，预测与实际利率之间的不匹配就越严重。

因此，我们得出结论，如果一个模型没有明确考虑到由于欧洲央行货币政策而导致利率突然变动的可能性，就无法成功地生成可靠的长期预测。EUR3M曲线的预测结果与观察到的EONIA曲线一致，在短期预测期（一周）内令人满意，预测率为95%。然而，对于三个月的期限曲线，基于高斯模型的结果比使用bootstrapforecasting方法得到的结果稍差。同样在这种情况下，当考虑覆盖概率的更高值时，bootstrap方法提高了测试的质量。尽管如此，针对EONIA案例得出的关于更长预测期的相同结论也适用于三个月期限结构。4结论在本文中，我们提出了一种新的利率期限结构预测方法，该方法适合于多收益率曲线的后危机世界。与基于因子分析、过滤历史模拟或流行的RiskMetricsTM方法的生成技术的单曲线法不同，目前利率风险管理的质量和效果取决于正确描述EONIA期限结构以及三个月、六个月和一年曲线的能力。据我们所知，这是首次尝试捕捉历史时间序列的波动相关性结构，该时间序列是为多收益率曲线而设计的。我们提出了aHJM模型框架，其中我们描述了即时远期利率的贴现曲线，而远期利率的动态决定了更长期限利率结构的演变。

我们展示了如何通过有限维的向量自回归过程来近似HJM设置中有限维SDE的连续时间。简化为离散时间模型显著简化了模型参数的估计。通过迭代最大似然法，我们估计了波动率、相关系数和风险溢价。特别是，我们展示了调查期间（2005-2013年）风险溢价的演变，这说明了在利率期限结构中观察到的变化。显著正的风险溢价表明，市场期限结构隐含的远期利率需要进行强烈的负修正，以描述即期利率的实际演变。然后，我们进行样本外数值测试，最终证明我们的方法在一周的预测期内的可靠性。如果在预测过程中采用简单的自举方法，我们的高斯模型可以进一步改进以捕获尾部事件。随着预测时间的延长，模型的性能会恶化。实证分析有力地支持了这样一个假设，即绝大多数故障都是随着速率动态的突然变化而累积的。这些事件——似乎是由欧洲央行货币政策的变化驱动的——在当前的建模方法中很难适应。作为当前工作的未来展望，我们计划在一定程度上适应体制转变，以增强我们的估计程序，包括基于调查的短期收益率预测。

就波动性时间序列中随机性的证据而言，我们目前正在研究离散时间设置的修改版本，能够捕捉瞬时远期和远期利率波动的异方差性质和可能的不对称性。感谢弗拉维亚·巴尔索蒂、安德烈亚·贝尔塔尼亚、托马索·科洛扎、富尔维奥·科尔西、尼科洛·科蒂尼、洛伦佐·利希、斯特凡诺·马米、阿尔多·纳西格、安德烈亚·帕拉维奇尼和罗伯托·雷诺进行了许多令人振奋的讨论。我们也感谢Andrea Sillari提供的历史数据。CS的研究活动得到了联合信贷银行（UniCredit S.p.A.）第1300240/2013号拨款的支持，该拨款由SCOLA NORMAL Superiore管理。英国广播公司感谢特斯拉师范高等教育基金SNS14 B BORMETTI的研究支持。贝塞尔三次样条每次我们需要插值在一组有限点上定义的曲线时，我们都会使用贝塞尔三次样条方法。让我们考虑通过以下K点{（s，g（s）），…，（sK，g（sK））}定义的一般曲线g（u）。三次样条插值函数由一组K- 1三阶多项式样条线（x）=ah+bh（x）-sh）+ch（x-sh）+dh（x-sh），对于sh≤ 十、≤ sh+1和h=1，K-1.每个多项式有4个系数需要确定，因此需要施加的约束总数为4K- 4.在[53]中，作者讨论了约束条件，并详细说明了后续结果的推导。虽然ai仅等于si处插值函数的值，即ai=g（si），但系数的解析表达式与tob=s对应- ss+s- 2ss- s（g（s）- g（s））-s- 党卫军- s（g（s）- g（s））,bi=si+1- 硅-1.si+1- 西西- 硅-1（g（si）- g（si）-1） ）+si- 硅-1si+1- si（g（si+1）- g（si））对于i=2。

K- 1，bK=-sK- sK-2.sK- sK-1sK-1.- sK-2（g）sK-1) - g（sK）-2)) -2sK- sK-1.- sK-2sK- sK-1（f（tk，sK）- f（tk，sK）-1)),上述方程可以改写为TKW处三个相邻桶的瞬时远期利率的线性叠加，其系数仅取决于这三个桶的到期时间SB=A（s，s，s）g（s）+B（s，s）g（s）+C（s，s，s）g（s），bi=Ai（si）-1，si，si+1）g（si-1） +Bi（si）-1，si，si+1）g（si）+Ci（si）-1，si，si+1）g（si+1），对于i=2，K- 1，bK=AK（sK-2，sK-1，sK）g（sK）-2） +BK（sK）-2，sK-1，sK）g（sK）-1） +CK（sK）-2，sK-1，sK）g（sK），其中A（s，s，s）=2s-s-s（s）-s） （s）-s） ，B（s，s，s）=s-s（s）-s） （s）-s） ，C（s，s，s）=s-s（s）-s） （s）-s） ,，艾未未-1，si，si+1）=si-si+1（si-硅-1） （si+1）-硅-1） ，Bi（si）-1，si，si+1）=si-1.-2si+si+1（si-硅-1） （si+1）-硅-1） ，Ci（si）-1，si，si+1）=si-硅-1（si+1）-si）（si+1-硅-1） ，i=2，K- 1.AK（sK）-2，sK-1，sK）=sK-sK-1（sK）-1.-sK-2） （sK）-sK-2） ，BK（sK）-2，sK-1，sK）=sK-2.-sK（sK-1.-sK-2） （sK）-sK-1） ，CK（sK-2，sK-1，sK）=2sK-sK-1.-sK-2（sK）-sK-1） （sK）-sK-2） ，同样的推理也适用于系数chand dhc=s- sg（s）- g（s）s- s-g（s）- g（s）s- s,d=0ci=si- si+12（g（si）- g（si）-1） ）（是的- 硅-1)-2（g（si+1）- g（si）-1） ）si+1- 硅-1.-g（si+2）- g（si）si+2- si+g（si+2）- g（si+1）si+2- si+1,di=（si）- si+1）g（si）- g（si）-1） （是的- 硅-1)-g（si+1）- g（si）-1） si+1- 硅-1.-g（si+2）- g（si）si+2- si+g（si+2）- g（si+1）si+2- si+1,对于i=2，K-1.因此，两个chand Dh都可以写成四个相邻瞬时远期利率值的线性叠加，系数完全取决于到期时间向量C=A（s，s，s）g（s）+B（s，s）g（s）+C（s，s）g（s），ci=Ai（si）-1，si，si+1，si+2）g（si-1） +Bi（si）-1，si，si+1，si+2）g（si）+Ci（si）-1，si，si+1，si+2）g（si+1），+Di（si）-1，si，si+1，si+2）g（si+2），对于i=2。

K- 2，cK-1=AK-1（sK）-2，sK-1，sK）g（sK）-2） +BK-1（sK）-2，sK-1，sK）g（sK）-1） +CK-1（sK）-2，sK-1，sK）g（sK），其中A（s，s，s）=（s-s） （s）-s） ，B（s，s，s）=-（s）-s） （s）-s） ，C（s，s，s）=（s-s） （s）-s） ,，艾未未-1，si，si+1，si+2）=（si-硅-1） （si+1）-硅-1） ，Bi（si）-1，si，si+1，si+2）=si-1+si-2si+2（si）-硅-1） （si+1）-si）（si+2-si），Ci（si）-1，si，si+1，si+2）=-硅-1+si+1-2si+2（si+1-硅-1） （si+1）-si）（si+2-si+1），Di（si）-1，si，si+1）=-（si+2）-si）（si+2-si+1），i=2，K- 2.AK-1（sK）-2，sK-1，sK）=（sK-1.-sK-2） （sK）-sK-2） ，BK-1（sK）-2，sK-1，sK）=-（sK）-1.-sK-2） （sK）-sK-1） ，CK-1（sK）-2，sK-1，sK）=（sK-sK-2） （sK）-sK-1） ，d=0，di=Ai（si-1，si，si+1，si+2）g（si-1） +Bi（si）-1，si，si+1，si+2）g（si）+Ci（si）-1，si，si+1，si+2）g（si+1），+Di（si）-1，si，si+1，si+2）g（si+2），对于i=2，K- 2dk，K-1=0，其中艾未未-1，si，si+1，si+2）=（si-硅-1） （si+1）-硅-1） （si+1）-si），Bi（si）-1，si，si+1，si+2）=si+2-硅-1（si）-硅-1） （si+2）-si）（si+1-si），Ci（si）-1，si，si+1，si+2）=si+2-硅-1（si+1）-硅-1） （si+2）-si+1）（si+1）-si），Di（si）-1，si，si+1）=（si+1-si）（si+2-si）（si+2-si+1），i=2，K- 总之，所有系数都可以表示为矩阵向量乘积ah=[g]h，h=1，K- 1，bh=[M（s）g]h，h=1，K、 ch=[M（s）g]h，h=1，K- 1，dh=[M（s）g]h，h=1，K- 1，（32）式中，g=[g（s），…，g（sK）]|，M（s）=ABC00。0ABC0 0 0。00 ABC00。0......0 0 0 . . . 0 AK-1BK-1CK-10 0 0 . . . 0 AKBKCK,男(s)=ABC00。0ABCD0。00 ABCD0。0......0 0 . . . 0 AK-2BK-2CK-2DK-20 0 . . . 0 AK-1BK-1CK-1.,男(s)=0 0 0 0 0 0 . . . 0ABCD0。00 ABCD0。0......0 0 . . . 0 AK-2BK-2CK-2DK-20 0 . . .

0 0 0 0 0.函数g（u）的样条插值形式在插值点处的导数将减少到系数bixgspline（x）x=xi=bi=M（s）g]i，i=1，K.至于积分，我们需要计算以下数量zSidu gspline（u），我们已经放弃了对Ai，Bi。为便于记谱，请在s上单击。其中Si是插值点之一。我们假设s>0，就像我们分析中的情况一样，并选择将函数fl从sdown外推到0，以便zsidu gspline（u）=Zsdu a+i-1Xh=1Zsh+1shduah+bh（u）- sh）+ch（u- sh）+dh（u- （上海）= as+i-1Xh=1啊（sh+1）- sh）+bh（sh+1- sh+ch（sh+1）- sh）+dh（sh+1- （上海）.现在我们可以使用系数ah，bh，chand dhas报告的不等式（32）Zsidu gspline（u）=g（s）s+i的显式表达式-1Xh=1生长激素（sh+1）- sh）+[Mg]h（sh+1- [sh（sh）+1毫克- sh）+[Mg]h（sh+1- （上海）= g（s）s+KXj=1i-1Xh=1δjh（sh+1）- sh）+Mhj（sh+1- sh）+Mhj（sh+1- sh）+Mhj（sh+1- （上海）g（sj）≡KXj=1P（s）ijg（sj）=[P（s）g]i，其中P（s）ij=sδj1+i-1Xh=1δjh（sh+1）- sh）+Mhj（sh+1- sh）+Mhj（sh+1- sh）+Mhj（sh+1- （上海）,所以P（s）=沙P1j=0，j>1，P（s）2j=sδj1+M1j（s- s） +M1j（s）- s） +M1j（s- （s）=> P2j（s）=0j>3，为了便于记谱，我们放弃了M和Mon的依赖关系。对于较大的i和j，依此类推。最后，矩阵P的形式由P（s）给出=P（s）0 0。0P（s）P（s）P（s）0 0。0P（s）P（s）P（s）P（s）0。0......P（s）K 1P（s）K 1P（s）K 2。P（s）K.我们得到了方程（9）和（16）中给出的结果，用EONIA瞬时正向曲线fint（tk，x）和FRA项结构f RA的时间观测代替gspline（x）,分别是int（tk，x）。用σf，int（tk，x）分量的时间tk观测代替gspline（x），可以很容易地导出方程（11）。参考文献[1]R.B.Litterman和J。

Scheinkman，“影响债券收益的共同因素”，《固定收益杂志》，第一卷，第一期，第54-61页，1991年。[2] P.J.Knez，R.Litterman和J.Scheinkman，“探索货币市场回报的解释因素”，《金融杂志》，第49卷，第5期，第1861-1882页，1994年。[3] F.Jamshidian和Y.Zhu，“情景模拟：理论和方法”，金融与随机，第一卷，第1期，第43-67页，1996年。[4] R.利率期权模型。威利金融工程系列，威利，第二版，1998年5月。[5] K.P.Scherer和M.Avellaneda，“人人为一……人人为一？1994年至2000年拉丁美洲布雷迪债券债务的主要成分分析”，《国际理论和应用金融杂志》，第05卷，第01期，第79-106页，2002年。[6] J.Driessen，P.Klassen和B.Melenberg，“定价和对冲上限与互换期权的多因素期限结构模型的表现”，《金融与定量分析杂志》，第38卷，第635-672页，2003年9月。[7] C.R.Nelson和A.F.Siegel，“收益率曲线的简约建模”，《商业杂志》，第60卷，第4期，第473-489页，1987年。[8] F.X.Diebold和C.Li，“预测政府债券收益率的期限结构”，《经济计量学杂志》，第130卷，第2期，第337-364页，2006年。[9] R.R.Bliss，“利率期限结构的变动”，《经济评论》，第四季度，第16-33页，1997年。[10] Q.Dai和K.J.Singleton，“有效期限结构模型的具体分析”，《金融杂志》，第55卷，1943-1978页，2000年10月。[11] F.de Jong和P.Santa Clara，“远期利率曲线的动态：带有状态变量的公式化”，《金融与定量分析杂志》，第34卷，第131-157页，1999年3月。[12] F.de Jong，“有效期限结构模型中的时间序列和横截面信息”，CEPR讨论论文2065，C.E.P.R.讨论论文，1999年2月。[13] G。

R.Duffee，“有效模型中的期限溢价和利率预测”，《金融杂志》，第57卷，第405-443页，2002年2月。[14] C.Bernadell，J.Coche和K.Nyholm，“战略投资者的收益率曲线预测”，欧洲央行工作文件系列第472号，2005年4月。[15] L.Coroneo，K.Nyholm和R.Vidova Koleva，“Nelson-Siegelmodel的无套利性如何？”《实证金融杂志》，第18卷，第3期，第393-407页，2011年。[16] Monfort，Alain和Pegoraro，Fulvio，“切换VARMA期限结构模型-扩展版”，法国银行第191号工作文件。，2007年[17]C.Gourieroux、A.Monfort、F.Pegoraro和J.-P.Renne，“体制转换和债券定价”，法国银行第456号工作文件。，2013年[18]G.Barone Adesi，K.Giannopoulos和L.Vosper，“衍生证券组合的无相关性VaR”，期货市场杂志，第19卷，第5期，第583-602页，1999年。[19] G.Barone Adesi、K.Giannopoulos和L.Vosper，“使用过滤历史模拟（FHS）对衍生品组合进行反向测试”，欧洲金融管理，第8卷，第1期，第31-58页，2002年。[20] F.Audrino和F.特洛伊，“使用函数梯度下降生成精确收益率曲线场景”，法国银行第456号工作文件。，2004年[21]J.Teichman和M.V.W–uthrich，“一致产量曲线预测”，预印本，苏黎世，2013年。[22]O.Vasicek，“期限结构的均衡表征”，《金融经济学杂志》，第5卷，第2期，第177-1881977页。[23]D.Heath，R.Jarrow和A.Morton，“债券定价和利率期限结构：未定权益估值的新方法”，经济计量学，第60卷，第1期，第77-105页，1992年。[24]Y.Ait-Sahalia和R.L.Kimmel，“利用闭式似然展开估计一个有效的多因素期限结构模型”，《金融经济学杂志》，第98卷，第1期，第113-144页，2010年。[25]米。

亨拉德，《衍生品贴现中的讽刺》，威尔莫特杂志，2007年7月。[26]M.Morini，“解决利率市场的难题，”http://ssrn.com/abstract=1506046, 2009.[27]F.Mercurio，“利率和信贷紧缩：新公式和市场模型”，彭博投资组合研究论文第2010-01-FRONTIERS号。，2009年[28]F.Mercurio，“随机基础的LIBOR市场模型”，风险杂志，2010年12月，第84-89页，2010年。[29]F.Mercurio，“现代LIBOR市场模型：使用不同曲线预测利率和贴现”，国际理论与应用金融杂志，第13卷，第01期，第113-137页，2010年。[30]A.Pallavicini和M.Tarenghi，“多收益率曲线的利率建模”，arXiv:1006.4767v1[q-fin.PR]，2010年。[31]M.Fujii、Y.Shimada和A.Takahashi，“存在抵押品和多种货币的动态基础利差利率市场模型”，Wilmott，2011年第54卷，第61-73页，2011年。[32]S.Cr\'epey，Z.Grbac和H.-N.Nguyen，“银行间风险的多曲线HJM模型”，数学与金融经济学，第6卷，第3期，第155-190页，2012年。[33]S.Crpey，Z.Grbac，N.Ngor和D.Skovmand，“应用于cva计算的l`evy hjm多曲线模型，”http://ssrn.com/abstract=2334865, 2013.[34]N.Moreni和A.Pallavicini，“多重收益率曲线动力学的简约HJM建模”，量化金融，第14卷，第2期，第199-210页，2014年。[35]C.Cuchiero、C.Fontana和A.Gnoatto，“多产量曲线建模的通用hjm框架，”http://ssrn.com/abstract=2454003, 2014.[36]M.Henrard，“衍生品贴现中的讽刺第二部分：危机”，《威尔莫特期刊》，第2卷，第6期，第301-316页，2010年。[37]M.P.A.Henrard，“随机价差的多曲线框架：刺激期货及其期权的连贯方法”，OpenGamma定量研究，第号。

2013年3月11日。，2013年[38]M.Bianchetti，“两条曲线，一个价格”，风险杂志，第23卷，第8期，第66-72页，2010年。[39]M.Bianchetti和M.Carlicchi，“信贷危机后的利率：多重CurveVilla衍生品和SABR”，arXiv预印本arXiv:1103.25672011。[40]M.Kijima，K.Tanaka和T.Wong，“利率的多质量模型”，量化金融，第9卷，第2期，第133-145页，2009年。[41]C.Kenyon，“冲击后短期利率定价”，风险杂志，2010年11月，第83-872010页。[42]M.Grasselli和G.Miglietta，“柔性点多曲线模型，”http://ssrn.com/abstract=2424242, 2014.[43]L.Morino和W.J.Ruggaldier，“关于期限结构的多曲线模型，”http://arxiv.org/abs/1401.5431, 2014.[44]M.Henrard，多曲线框架下的利率建模。应用定量金融，帕尔格雷夫·麦克米伦，2014年版，2014年5月。[45]P.H.Kupiec，“验证风险度量模型准确性的技术”，《衍生工具杂志》，第3卷，第2期，第73-84页，1995年。[46]P.F.Christo Offersen，“评估区间预测”，《国际经济评论》，第39卷，第4期，第841-862页，1998年。[47]J.-P.Renne，“具有离散政策利率的欧元区收益率曲线模型”，可查阅SSRN 20152362014。[48]D.H.Kim和A.Orphanides，“利率预测调查数据的期限结构估计”，《金融与定量分析杂志》，第47卷，第01期，第241-272页，2012年。[49]M.Musiela，“随机偏微分方程和期限结构模型”，预印本，1993年。[50]A.Brace和M.Musiela，“Heath、Jarrow和Morton的多因素Gauss-Markov实现”，数学金融，第4卷，第3期，第259-283页，1994年。理论上，连续时间套利。牛津大学出版社，2004年。[52]F.M.Ametrano和M。

Bianchetti，“引导非流动性：市场一致性远期利率估计的多收益率曲线构造”，利率建模，法比奥·默库里奥主编，风险书籍，锐思传媒，2009年。[53]P.S.Hagan和G.West，“曲线构造的插值方法”，应用数学金融，第13卷，第2期，第89-129页，2006年。[54]T.W.Anderson，“主成分分析的渐近理论”，《数学统计年鉴》，第34卷，第122-148页，1963年3月。[55]B.Flury，多元统计的第一门课程。《统计学中的斯普林格文本》，斯普林格出版社，1997年版，1997年8月。[56]B.Efron，“引导方法：另一种刀切法”，安。统计学家。，第7卷，第1-26页，1979年1月1日。图1：根据期限曲线计算的两年连续复合收益率之间的基点差异 = 3M、6M和1Y以及EONIA两年利率。（蓝色→ 300万欧元，橙色→ 600万欧元紫色→ 1欧元。）5.05.25.45.65.86.0窗口结束日期F`EΦHFLr 95%。考虑到EONIA曲线的K=12个桶，得出了结果。每个点指的是三年的数据，每周采样一次。x轴上报告的日期对应于三年期的结束日期，然后将其提前五天以获得连续点。ccccccccccccóóóóóóóóóóóóíííííííííííí~o~o~o~o~o~o~o~o~o~o~o~o-0.8%-0.6%-0.4%-0.2%0.%0.2%0.4%s HyearsLWeekly modified volatility Figure 3：首次在单曲线框架中修改了EONIA期限结构的波动率函数。在x轴上，我们每年报告到期时间。这些点对应于K=12个桶。蓝色→ w（f），红色→ w（f），绿色→ w（f），橙色→ w（f），紫色→ w（f）。

数据指的是2010年1月5日至2013年1月5日这三年期间的数据，采样频率为每周。窗口结束日期F`EΦHFLr 95%图4：要保留的协方差矩阵^C主成分的最小数量，以保持原始方差的一部分大于95%。考虑到EONIA曲线的K=12个铲斗和EUR3M TERM结构的K3M=10个铲斗，得出的结果总共为D=22个部件。每个点指的是三年的数据，每周采样频率。x轴上报告的日期对应于三年期的结束日期，然后将其提前五天获得连续点。ccccccccccccóóóóóóóóóóóóíííííííííííí~o~o~o~o~o~o~o~o~o~o~o~o-0.8%-0.6%-0.4%-0.2%0.%0.2%0.4%s Hyears每周修正波动率ccccccccccóóóóóóóóóóíííííííííí~o~o~o~o~o~o~o~o~o~o-0.8%-0.6%-0.4%-0.2%0.%0.2%0.4%s HyearsLWeekly modified volatility Figure 5：在多曲线框架中，首次对EONIA（K=12）和EUR3M（K3M=10）曲线的波动率函数进行了五次修改。在x轴上，我们每年报告一次成熟度。左面板：w（f）m=1，5.右侧面板：w（3M）mwithm=1，5.颜色如图3所示。数据指的是三年期间2010年1月5日至2013年1月5日，每周抽样一次。0.00.51.01.5窗口结束日期∧`s图6：黑线：在周间隔重叠收益率的三年样本上估计模型得到的λs值。红线：bootstrap估计器的中心值。黑色虚线：黑色虚线±一个标准偏差，用自举法计算-0.3-0.2-0.10.00.10.2窗口结束日期∧`L图7：黑线：基于周间隔重叠收益率的三年样本估计模型的λL值。

红色和黑色虚线如图6.0.00.10.20.30.40.5窗口结束日期∧\'3 M图8：黑线：根据三年间隔的重叠收益率样本估计模型得出的λ3值。图6.0.0030.0040.0050.006窗口结束日期标题中的红色和黑色虚线Ohm`H1ML图9：黑线：估计模型在三年的周间隔重叠收益率样本上获得的EONIA曲线ω（1m）值。图6.0.0070.0080.0090.0100.0110.012窗口结束日期中的红色和黑色虚线Ohm`H10YL图10：黑线：根据每周间隔重叠收益的三年样本估计模型获得的EONIA曲线ω（10y）值。图6.0.0000.0020.0040.0060.008窗口结束日期中的红色和黑色虚线Ohm`H3ML图11：黑线：根据每周间隔重叠收益的三年样本估计模型获得的EUR3M曲线的ω（3m）值。图6.0.0080.0090.0100.0110.012窗口结束日期中的红色和黑色虚线Ohm`H10YL图12：黑线：根据每周间隔重叠收益的三年样本估计模型获得的EUR3M曲线的ω（10y）值。红色和黑色虚线如图6所示-0.15-0.10-0.050.000.050.100.150.20窗口结束日期g`H1m-10YL图13：黑线：Γ（1m）的值- 10y）对于获得的EONIA曲线，根据每周间隔重叠收益的三年样本估计模型。图6标题中的红色和黑色虚线-0.20-0.15-0.10-0.050.000.050.10窗口结束日期G`H3m-10YL图14：黑线：Γ（3m）的值- 10 y）对于获得的EUR3M曲线，根据周间隔重叠收益的三年样本估算模型。

图6标题中的红色和黑色虚线。本周三（23日）的霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍窗口结束日期FHTK，1mL窗口结束日期FHTK，1mLaeaeaeaeaeaeaeaeaeaeaeae''''0.000.010.020.030.04窗口结束日期FHTK，1mL图15：蓝色实线：f（t，1m）利率三个月前预测分布的平均值。蓝色虚线：预测未来三个月f（t，1m）率分布的平均值±两个标准差。黑色实线：实现的f（t，1m）速率。红点（十字）：第一行和第二行的覆盖率分别为95%和99%的负（正）异常。左栏中的曲线图指的是高斯扩散模型，右栏中的曲线图指的是自举预测方法。3.0.0.020.020.020.020.04.05窗口收尾日期FHTK、5 yl-3.0 0.0.0.0.0.0.0.0.0.0.0.0.0.0 0.0.0.0.0.0.0.0.0.0.0.0.0.0.0 0.0.0.0.0.0.0.0.0.0.0.0.0.05窗口收尾日期FTTK、5窗口收收尾日期FHTTTK、5窗口收尾日期、5基日期、5基基、5基基、5日、5基方方方、基方、基方、基方方、基方、基方方、基方、基方方、基方、基方、3.3.3.3.3.3.3.0 0 0无无无托托托托托托托托托托托托托托基、基、无无无无无无0.05窗口结束日期FHTK，5YL图16：蓝色实线：未来三个月f（t，5y）利率预测分布的平均值。蓝色虚线：预测未来三个月f（t，5y）率分布的平均值±两个标准差。黑色实线：实现的f（t，5y）速率。红点（十字）：第一行和第二行的覆盖率分别为95%和99%的负（正）异常。

左栏和右栏如图15所示。本周三（23日）的霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍aeaeae'''''''''''''''''0.000.010.020.030.040.05窗口结束日期FRA3 MHtk，本周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的方方方（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三（周三）的周三）的周三（周三）的日日日（周三）的）的日日（周三）的日日（周三）的日日）的日（周三））的日（周三）的日（日日））的日日（日）的日（日日））的日日（日）））的日（日日日（日）））的日（日（日日）的）））的日（日（日日）））的日（日（日日））））的日（日日日日（日日）））的日（日））的日（日（日）的）直线：FRA3M（t，3m）利率三个月前预测分布的平均值。蓝色虚线：FRA3M（t，3m）利率三个月前预测分布的平均值±两个标准偏差。黑色实线：3米（t，3米）速率。红点（十字）：负（正）异常，第一行和第二行的覆盖概率分别为95%和99%。左栏和右栏如图15所示。（2300）托托托霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍霍奈奈尔，，目前目前的他们他们他们的两两0 0 0 0 0 0 0 0 0 0 0 0 0.0 0 0 0 0.10.10.10.0.0 0.0.0.0 0.0.0 0.0.0.0.0 0.0.0.0 0.0.0.0.0.0.0 0.0.0.0.0.0 0.0.0 0.0.0.0 0.0.0.0.0 0.0.0.0.0.0.0 0 0.0 0.0.0.0 0.0 0 0.010.020.030.040.05窗口结束日期FRA3 MHtk，5YL图18：蓝色实线：FRA3M（t，5y）利率三个月前预测分布的平均值。蓝色虚线：FRA3M（t，5y）利率三个月预测分布的平均值±两个标准偏差。黑色实线：3米（t，5年）速率。红点（十字）：负（正）异常，第一行和第二行的覆盖概率分别为95%和99%。

左栏和右栏如图15所示。期限开始日期结束日期2013年12月27日2005年8月2日2013年12月27日2013年2月8日EUR3M表1：历史时间序列的开始和结束日期。到期期限EONIA 1d 7d 1m 2m 3m 6m 9m 1y 1y 6m 2y 3y 4Y 5Y 6y 7y 8y 10y 12y 15y 20y 25y 30Y 3m 1d 2m 3m 6m 1y 1y 6m 2y 3y 4Y 5Y 6y 7y 8y 10y 12y 15y 25y 30Y表2：到期时间网格。期限到期时间EONIA 1m 2m 3m 6m 9m 1y 5y 10y 15y 20y 25y 30yEUR3M 6m 9m 1y 5y 10y 15y 20y 25y 30Y表3：EONIA和EUR3M曲线实证分析中使用的到期时间网格。4.5.76 12 16.24（*）0.012米20 2.01 15.58 10 10 10.78（*）0.15 15.58 10 10 10.78（*）0.13米15.15.10 10 10.78（*）0.13米15 0.15 15 15 15.10 10 10.10 10 10.78（10）10.10 10 10 10 10.10 10 10 10 10 10.78（10）10.15 15 15 15 15 15 15 0.15 0.15 15 0.10.10 10 10 10.78（10.78（10.78（10.78）10.78（*）10.78（*）10（10（10（10）10）7）10.78（7）10（10）10（10）10）10.78）10）10.78（7）10.78（7）10）10.78（10（10）10）10.78（7）10（10）10）10）10.78（7）7）10（10.78（7）10）10））3.9825y 19 1.38 24.04 7 4.22（*）3.9830y 17 0.45 50.26 12 16.24（**）0.01表4：使用高斯扩散模型获得的EONIAcurve一周前预测的无条件覆盖率测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value 1M 19 1.38 24.04 7 4.22（*）3.982m 17 0.45 50.26 4 0.38 53.513m 11 0.94 33.19 3 0。94.856m 120.46 49.63 30。94.859m 9 2.49 11.49 2 0.31 57.751y 11 0.94 33.19 3 0。94.855y 18 0.85 35.53 7 4.22（*）3.9810y 18 0.85 35.53 2 0.31 57.7515y 16 0.17 68.07 2 0.31 57.7520y 14 0.01 90.29 3 0。94.8525y14 0.01 90.29 5 1.28 25.8430y 18 0.85 35.53 7 4.22（*）3.98表5：使用自举方法获得的EONIAcurve一周前预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value1m 51 63.87（**）0。2567.23（**）0.2m 5369.77（**）0。2567.23（**）0.3M4855.38（**）0。2567.23（**）0.6m 4547.34（**）0。

2358.28（**）0.9m 40 35.03（**）0。22 53.95（**）0.1y 37 28.33（**）0。20 45.59（**）0.5 Y 25 7.62（**）0.58 13 20.05（**）0.10 Y 17 0.68 40.9 4 0.48 49.0215 Y 22 4.25（*）3.92 7 4.55（*）3.2920 Y 35 24.17（**）0。12 16.97（**）0.25y 34 22.19（**）0。14 23.29（**）0.30y 21 3.32 6.83 15 26.68（**）0。表6：使用高斯扩散模型获得的EONIAcurve三个月前预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value1m 32 18.44（**）0。21.49.72（*）0.21 49.72（*）0.21 49.72（*）0.21.72（7）49.72（*）0.21 21.72（7）21.72（7）21.72（7）21.72（7）21.72（7）21.72（*）0.21.72（7）21.21.72（7）21.72（7）7（7）0.30（3）19（7）19（7）10（7）10（7）10（7）10）10.19（7（7）7）7（7）7（7（7）7）7）7（7（7（7）7）7）7（7）7（7）7（7（7）7）7）7（7（7（7）7）7（7（7）7）7）7）7（7（7（7）7）7）7）7）7（7（7（7）7）7）7（7（7）7）7*）0.30y 15 0.09 76.49 12 16.97（**）0。表7：使用bootstrap方法获得的EONIAcurve三个月前预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value1m 93 253.6（**）0。70359.74（**）0.2M88229.04（**）0。69352.32（**）0.3M88229.04（**）0。63308.68（**）0.6m 78 182.65（**）0。51226.17（**）0.9m 64124.19（**）0。2051.26（**）0.1比3634.17（**）0。7.5.95（*）1.475y 30 20.77（**）0。9 10.89（**）0.110y 27 15.07（**）0.01 30.15 69.7915y 46 61.62（**）0。2051.26（**）0.20y 71 152.44（**）0。27 84.58（**）0.25y 67 136.05（**）0。31105.53（**）0.30y 55 90.98（**）0。21 55.72 (**) 0.表8：使用高斯扩散模型获得的EONIA曲线的1年前预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value1m 35 31.76（**）0。2051.26（**）0.2m 3224.95（**）0。2051.26（**）0.3m 3020.77（**）0。2051.26（**）0.6m 2511.69（**）0.061219.99（**）0.9m9 0.8136.830 4.78（*）2.871y 0 24.42（**）0。

0 4.78（*）2.875y 17 2.04 15.29 0 4.78（*）2.8710y 0 24.42（**）0。0.4.78（*）2.8715y 2 13.09（**）0.03 0 4.78（*）2.8720y 3 9.88（**）0.17 0 4.78（*）2.8725y 12 0。97.63 1.03 30.9330y 15 0.79 37.47 7 5.95（*）1.47表9：使用自举方法获得的EONIA曲线1年前预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value 3M 10 1.61 20.46 7 4.22（*）3.986m 26 7.94（**）0.48 13 19.24（**）0.9m 35 22.4（**）0。19 40.27（**）0.1 Y 26 7.94（**）0.48 11 13.42（**）0.025y 25 6.72（**）0.95 7 4.22（*）3.9810y 24 5.59（*）1.81 6 2.58 10.8215y 33 18.69（**）0。17 32.73（**）0.20y 19 1.38 24.04 8 6.16（*）1.325y 18 0.85 35.53 9 8.36（**）0.3830y 16 0.17 68.07 8 6.16（*）1.3表10：使用高斯扩散模型获得的EUR3M曲线提前一周预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value 3M 6 6.61（*）1.01 4 0.38 53.516m 22 3.6 5.76 8 6.16（*）1.39m 34 20.51（**）0。14 22.4（**）0.1 Y 19 1.38 24.04 8 6.16（*）1.35 Y 21 2.76 9.68 3 0。94.8510y 17 0.45 50.26 3 0。94.8515y 25 6.72（**）0.95 15 25.7（**）0.20y 18 0.85 35.53 6 2.58 10.8225y 18 0.85 35.53 7 4.22（*）3.9830y 15 0.02 88.27 6 2.58 10.82表11：使用自举方法获得的EUR3M曲线的提前一周预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value 3M 57 82.12（**）0。2986.14（**）0.6m 5885.32（**）0。33 106.29（**）0.9m 64 105.42（**）0。36122.14（**）0.1比4547.34（**）0。29 86.14（**）0.5y 37 28.33（**）0。19 41.57（**）0.10y 21 3.32 6.83 10 11.35（**）0.0815y 35 24.17（**）0。21 49.72（**）0.20y 32 18.44（**）0。12 16.97（**）0.25y 39 32.74（**）0。16 30.21（**）0.30y 38 30.5（**）0。

16 30.21 (**) 0.表12：使用高斯扩散模型获得的EUR3M曲线三个月前预测的无条件覆盖测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value 3M 35 24.17（**）0。21 49.72（**）0.6m 38 30.5（**）0。20 45.59（**）0.9m 40 35.03（**）0。26 71.83（**）0.1 Y 26 8.93（**）0.28 20 45.59（**）0.5 Y 28 11.78（**）0.06 9 8.85（**）0.2910y 6 5.95（*）1.47 0 5.59（*）1.8115y 23 5.28（*）2.15 5 1.45 22.8920y 15 0.09 76.49 13 20.05（**）0.25 Y 19 1.78 18.26 12 16.97（**）0.30 Y 16 0.32 57.21 16.97（**）0。表13：使用bootstrap方法获得的EUR3M曲线三个月前预测的无条件覆盖率测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value 3M 97 273.88（**）0。78420.52（**）0.6m 85214.73（**）0。69352.32（**）0.9m 80191.63（**）0。56 259.73（**）0.1 Y 75 169.47（**）0。38 145. （**）0.5y 78 182.65（**）0。43175.1（**）0.10y 112354.86（**）0。43175.1（**）0.15y 4764.67（**）0。17 38.53（**）0.20y 45 58.61（**）0。8.29（**）0.425y 65 128.1（**）0。2574.59（**）0.30y 75169.47（**）0。39 150.9 (**) 0.表14：使用高斯扩散模型获得的EUR3M曲线1年前预测的无条件覆盖率测试结果。sip=0.95 p=0.99nLRUCp-value nLRUCp-value 3M 57 98.04（**）0。35127.68（**）0.6m 3531.76（**）0。19 46.9（**）0.9m 21 6.03（*）1.41 17 38.53（**）0.1 Y 17 2.04 15.29 13 23.39（**）0.5 Y 2 13.09（**）0.03 0 4.78（*）2.8710y 1 17.36（**）0。0.4.78（*）2.8715y 9 0.81 36.83 1 1.03 30.9320y 3 9.88（**）0.17 0 4.78（*）2.8725y 13 0.1 74.7 2 0.06 79.9130y 16 1.35 24.56 7 5.95（*）1.47表15：使用自举方法获得的EUR3M曲线的1年前预测的无条件覆盖率测试结果。