MAXIMUM LIKELIHOOD ESTIMATION OF THE
COX-INGERSOLL-ROSS PROCESS: THE MATLAB
IMPLEMENTATION
Kamil Klad´ıvko
1
Department of Statistics and Probability Calculus, University of Economics, Prague
and Debt Management Department, Ministry of Finance of the Czech Republic
kladivk@vse.cz or kamil.kladivko@mfcr.cz
Abstract
The square root diffusion process is widely used for modeling interest rates
behaviour. It is an underlying process of the well-known Cox-Ingersoll-Ross
term structure model (1985). We investigate maximum likelihood estimation
of the square root process (CIR process) for interest rate time series. The
MATLAB implementation of the estimation routine is provided and tested on
the PRIBOR 3M time series.
1 CIR Process for Interest Rate Modeling
A continuous-time model in finance typically rest on one or more stationary diffusion processes
{X
t
, t ≥ 0}, with dynamics represented by stochastic differential equations:
dX
t
= µ(X
t
)dt + σ(X
t
)dW
t
, (1)
where {W
t
, t ≥ 0} is a standard Brownian motion. The functions µ(·) and σ
2
(·) are, respectively,
the drift and the diffusion functions of the process.
The fundamental process in interest rate modeling is the square root process given by the
following stochastic differential equation:
dr
t
= α(µ − r
t
)dt +
√
r
t
σdW
t
, (2)
where r
t
is the interest rate and θ ≡ (α, µ, σ) are model parameters. The drift function µ(r
t
, θ) =
α(µ − r
t
) is linear and possess a mean reverting property, i.e. interest rate r
t
moves in the
direction of its mean µ at sp ee d α. The diffusion function σ
2
(r
t
, θ) = r
t
σ
2
is proportional to the
interest rate r
t
and ensures that the process stays on a positive domain.
The square ro ot process (2) is the basis for the Cox, Ingersoll, and Ross short-term interest
rate model [1] and therefore often denoted as the CIR process in the financial literature.
1.1 CIR process densities
If α, µ, σ are all positive and 2αµ ≥ σ
2
holds, the CIR process is well-defined and has a steady
state (marginal) distribution. The marginal density is gamma distributed.
For maximum likelihood estimation of the parameter vector θ ≡ (α, µ, σ) transition den-
sities are required. The CIR process is one of few cases, among the diffusion processes, where
the transition density has a closed form expression. We follow the notation given in [1] on page
391. Given r
t
at time t the density of r
t+∆t
at time t + ∆t is
p(r
t+∆t
|r
t
; θ, ∆t) = ce
−u−v
(
v
u
)
q
2
I
q
(2
√
uv), (3)
1
This work was supported by grant of IGA VSE nr. IG410046
咚咚跳
- 粉丝: 1
- 资源: 1
我的内容管理
收起
- 我的资源
快来上传第一个资源
我的收益 登录查看自己的收益
我的积分
登录查看自己的积分
我的C币
登录后查看C币余额
我的收藏 
我的下载 
下载帮助
会员权益专享
最新资源
- 保险服务门店新年工作计划PPT.pptx
- 车辆安全工作计划PPT.pptx
- ipqc工作总结PPT.pptx
- 车间员工上半年工作总结PPT.pptx
- 保险公司员工的工作总结PPT.pptx
- 报价工作总结PPT.pptx
- 冲压车间实习工作总结PPT.pptx
- ktv周工作总结PPT.pptx
- 保育院总务工作计划PPT.pptx
- xx年度现代教育技术工作总结PPT.pptx
- 出纳的年终总结PPT.pptx
- 贝贝班班级工作计划PPT.pptx
- 变电值班员技术个人工作总结PPT.pptx
- 大学生读书活动策划书PPT.pptx
- 财务出纳月工作总结PPT.pptx
- 大学生“三支一扶”服务期满工作总结(2)PPT.pptx
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
