22构建静态利率期限结构模型.pptVIP

标准5：有异常值时要剔除 按照一定的标准选择2003年7月21日的国债样本，到期收益率散点图可以看到，在靠近3年和靠近10年的位置都出现了异常点，像这样的债券就应该从债券样本中予以剔除。 标准5：剔除异常值后得到合理的期限结构 拟合利率期限结构既是一门科学，也是一门艺术 利率期限结构模型不再像那些债券基本指标计算一样，只有一个正确答案，而是需要包含一定的主观判断； 模型的建立者和使用者，都需要用自己的经验、对数据的直观判断，以及金融的理论与实务，来评价模型的质量。 22.4 NSS模型搜索算法 由于Nelsen-Siegel模型涉及到非线性最优化，实际应用起来比较麻烦。这里给出一种比较简单的近似算法—数据搜索算法。首先确实参数的范围，通过循环计算，最小化目标函数，确定出参数的估计值。由于在一定时期内，曲线形状不会发生很大变化，因此，可以确定相近日期的参数范围，从而减少计算量。 第22章 构建静态利率期限结构模型 清华大学经管学院 朱世武 Zhushw@ Resdat样本数据： SAS论坛： 22.1 银行间债券利率期限结构拟合 计算环境 2005年1月31日作为计算时点指标。 从2005年1月31日之前发行、2005年1月31日之后到期的固定利率政策性金融债券中选择样本，拟合政策性金融债券利率期限结构。 计算数据集：2005年1月31日固定利率金融债样本债券ResDat.SampFbd050131。 数据集的变量说明如下页。 ResDat.SampFbd050131变量说明 变量名 中文全称 Date 日期 Bondcls 债券种类 Bdid 债券标识 Bdcd 债券代码 Resbdid 锐思债券标识 Price 债券价格（收盘全价） Par 债券面值 Intmd 计息方式 Couprt 票面利率 Freq 年付息频率 Intdt 起息日 Matdt 到期日 Maturity 到期期限 Modifdur 修正久期 Yrstmat 剩余期限 Yield 到期收益率 /* 产生样本债券SampFbd050131未来现金流与对应的时刻 */ data SampFbd050131; set ResDat.SampFbd050131; if freq=0 then _t=Yrstmat+1e-12; else _t=mod(yrdif(date,matdt,act/act),1/Freq); /*_t为当前日到下一个付息的时间 */ if mod(_t,1)=0 then do; _t=_t+1e-12; /* 正好付息日时，因为_t为整数时，包含_t的函数没有结果 */ end; run; data tbond_info; set SampFbd050131; if freq=0 then do; t=_t; output; end; else do t=_t to Yrstmat by 1/freq; /* 产生现金流对应的时刻 */ output; end; proc sort; by resbdid; data tbond_info; set tbond_info; if Intmd=2 and last.resbdid=0 then CF=couprt*par/freq; /*到期日前附息债产生的现金流*/ if Intmd=2 and last.resbdid=1 then CF=couprt*par/freq+par; /*到期日附息债产生的现金流*/ if Intmd=0 and last.resbdid=0 then CF=0; /*到期日前贴现债现金流*/ if Intmd=0 and last.resbdid=1 then CF=par; /*到期日贴现债现金流*/ if Intmd=1 and last.resbdid=0 then CF=0; /*到期日前零息票债现金流*/ if Intmd=1 and last.resbdid=1 then CF=par+couprt*par*maturity; /*到期日贴现债现金流*/ by resbdid; run; /* 数据集tbond_info增加了3个变量： _t为当前日到下一个付息的时间 t：现金流发生时间 CF：现金流量 */ 多项式样条法 采用上一章得到的简化模型： /* 多项式样条法 */ data tbond_info1; set tbond_info; p1=CF; pb1=CF*t; pc1=CF*t**2; pd1=CF*((0=t=5)*t**3+(5t=30)*(t**3-(t-5)**3)); pd2=CF*((5t=10)*(t-5)**3+(10t=30)*((t-5)**3-(t-10