原文链接:http://tecdat.cn/?p=11680
这篇文章的目的是指导读者逐步使用R编程语言实现Nelson-Siegel模型的步骤。您可能已经知道,估计利率期限结构是任何资产定价的关键,因此对投资者和政策制定者起着重要的作用 ( 点击文末“阅读原文”获取完整代码数据 )。
想法是使一条连续曲线适合现有数据。就是说,给定可获取的利率和相应的到期日(通过彭博社或任何其他数据提供商),可以使用Nelson-Siegel方法得出利率的期限结构。
- R或RStudio
- LIBOR / OIS利率和相应的到期日(通过彭博社或其他数据提供商)
一点理论…
在开始执行模型之前,让我们回顾一下基础知识。
**
**
Nelson-Siegel零息债券的到期收益率:
零息债券价格:
优化问题:
这个想法是两个价格应该相等。因此,我们想找到使两个价格之间的平方差平方和最小的Nelson-Siegel因素。
** **
实施模型的步骤
第1步:数据导入和变量定义
- 我们导入LIBOR / OIS利率和相应的到期日,对其进行过滤,然后将其存储在数据框中。
- 然后,我们使用导入的LIBOR / OIS汇率计算每个到期日的零息票价格。
- 我们最终计算出最长(50年)和最短(1个月)到期的到期收益率(YTM)。
-
- 可以用50年的YTM来近似。
- 可以通过50年期和1个月YTM之间的差异来近似估算。
- 为了简化步骤2中描述的网格搜索和步骤3中描述的优化问题,我们将近似 并使用50年和1个月的YTM。
步骤2:对目标函数进行编程
- 我们对函数进行编程,该函数计算LIBOR / OIS利率给出的零息债券价格与Nelson-Siegel模型给出的零息债券价格之间的平方偏差的平方和。
点击标题查阅往期内容
R语言用神经网络改进Nelson-Siegel模型拟合收益率曲线分析
左右滑动查看更多
01
02
03
04
步骤3:网格搜索
- 我们定义为我们的参数范围 :
-
我们创建一个包含所有可能的组合矩阵 。
-
我们选择产生最小平方偏差总和的组合作为我们的初始参数,以插入优化问题。
步骤4:优化问题
- ,我们可以解决非线性优化问题:
-
- 初始参数(x0)是在网格搜索中找到的参数。
- 目标函数(eval_f)是在步骤2中编程的目标函数。
- 上限和下限(ub和lb)定义如下:
步骤5:调整模型
此时获得完美契合的机会非常渺茫。必须重复步骤3和4调整模型。进一步来说:
-
我们必须围绕步骤5中获得的结果执行第二次网格搜索,搜索范围较窄,然后重新运行优化问题。
-
- 您可能还想尝试使用不同的参数组合,得出平方偏差的第二,第三或第四最小和。
-
我们也将要执行的最后一个网格搜索 在第二轮优化得到的数值。
-
技巧 –在模型中尝试不同的初始参数时,针对LIBOR / OIS Bloomberg数据点绘制通过求解参数获得的最终收益曲线,以了解其拟合程度。没有完美的方法可以完成–这是一个反复试验的过程。
步骤6:绘制估计的收益曲线
现在,我们有了Nelson-Siegel因子,可以估算出收益率曲线:
结果与结论
下图以红色显示了所获得的LIBOR收益率曲线。彭博LIBOR / OIS汇率以绿色显示。
尽管曲线可以很好地拟合数据点,但是我们可以看到拟合并不完美。Nelson-Siegel模型的替代方案是Svensson模型,该模型增加了两个参数以实现更好的拟合。
**
**
Svensson模型–零票息收益率:
实施Svensson模型的步骤与实施Nelson-Siegel模型的步骤相同。
点击文末 “阅读原文”
获取全文完整代码数据资料。
本文选自《R语言中的Nelson-Siegel模型在汇率预测的应用》。
点击标题查阅往期内容
R语言用神经网络改进Nelson-Siegel模型拟合收益率曲线分析
分解商业周期时间序列:线性滤波器、HP滤波器、Baxter滤波器、Beveridge Nelson分解等去趋势法
用R语言用Nelson Siegel和线性插值模型对债券价格和收益率建模
R语言用神经网络改进Nelson-Siegel模型拟合收益率曲线分析
R语言和QuantLib中Nelson-Siegel模型收益曲线建模分析
R语言使用随机技术差分进化算法优化的Nelson-Siegel-Svensson模型
用R语言用Nelson Siegel和线性插值模型对债券价格和收益率建模R语言用神经网络改进Nelson-Siegel模型拟合收益率曲线分析
R语言中的Nelson-Siegel模型在汇率预测的应用
python使用LASSO回归预测股票收益
R语言数据的收益率和波动性交易
R语言用线性模型进行预测:加权泊松回归,普通最小二乘,加权负二项式模型,多重插补缺失值
使用SAS,Stata,HLM,R,SPSS和Mplus的分层线性模型HLM
R语言用线性回归模型预测空气质量臭氧数据