关于“sde”程序包计算随机值时结果有问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

使用程序包中的rcCIR和rcOU函数估计CIR和Vasicek随机利率模型的利率。
程序包中CIR和Vasicek的形式是：

和 ；
而论文中使用的一般是
和 。
程序包把括号外的回归速度参数代入到括号里面了，于是我在进行运行时也这么做了。

1. k<-rcOU(n=1000, Dt=0.003, x0=0.03, theta=c(0.15,5,0.005))

复制代码

假设计算三年内的利率变化，划分为1000份，Dt就是单位变化，近似连续变化，此时3/1000记为0.003；x0为初始利率，theta为程序包中对应的值。此时回归速度=5，远期利率均值为0.03，波动率为0.005.

1. > max(k)
   
2. [1] 0.03077587
   
3. > min(k)
   
4. [1] 0.02912205

复制代码

这时候有点不对劲，于是把远期利率均值改为0.05.

1. k<-rcOU(n=1000, Dt=0.003, x0=0.03, theta=c(0.25,5,0.005))

复制代码

1. > max(k)
   
2. [1] 0.03113812
   
3. > min(k)
   
4. [1] 0.02949922

复制代码

远期利率为0.05，但是最大的随机数还是这么点，按理说应该趋向0.05才对，结果最大值这么小。（实际上使用rcCIR也是这样）

但是论文中实际的结果是这样。
究竟是哪里出现了问题，是我的思路错了？
如果不能用这个程序，求推荐程序能模拟随机利率的；另外这样的利率变化图具体应该怎么画。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：SDE 程序包 Vasicek ASIC 远期利率 程序