Cox–Ingersoll–Ross model

Three trajectories of CIR Processes

In mathematical finance, the CoxIngersollRoss model (or CIR model) describes the evolution of interest rates. It is a type of "one factor model" (short rate model) as it describes interest rate movements as driven by only one source of market risk. The model can be used in the valuation of interest rate derivatives. It was introduced in 1985 by John C. Cox, Jonathan E. Ingersoll and Stephen A. Ross as an extension of the Vasicek model.

The model

CIR process

The CIR model specifies that the instantaneous interest rate follows the stochastic differential equation, also named the CIR Process:

dr_t = a(b-r_t)\, dt + \sigma\sqrt{r_t}\, dW_t

where W_t is a Wiener process (modelling the random market risk factor) and  a ,  b , and  \sigma\, are the parameters. The parameter  a corresponds to the speed of adjustment,  b to the mean and  \sigma\, to volatility. The drift factor,  a(b-r_t), is exactly the same as in the Vasicek model. It ensures mean reversion of the interest rate towards the long run value b, with speed of adjustment governed by the strictly positive parameter a.

The standard deviation factor, \sigma \sqrt{r_t}, avoids the possibility of negative interest rates for all positive values of a and b. An interest rate of zero is also precluded if the condition

2 a b \geq \sigma^2 \,

is met. More generally, when the rate (r_t) is close to zero, the standard deviation (\sigma \sqrt{r_t}) also becomes very small, which dampens the effect of the random shock on the rate. Consequently, when the rate gets close to zero, its evolution becomes dominated by the drift factor, which pushes the rate upwards (towards equilibrium).

This process can be defined as a sum of squared Ornstein–Uhlenbeck process. The CIR is an ergodic process, and possesses a stationary distribution. The same process is used in the Heston model to model stochastic volatility.

Distribution

The distribution of future values of a CIR process can be computed in closed form:

r_{t+T} = \frac{Y}{2c},

where c=\frac{2a}{(1 - e^{-aT})\sigma^2}, and Y is a non-central Chi-Squared distribution with \frac{4ab}{\sigma^2} degrees of freedom and non-centrality parameter 2 c r_te^{-aT}. Formally the probability density function is:

f(r_{t+T};r_t,a,b,\sigma)=c\,e^{-u-v} \left (\frac{v}{u}\right)^{q/2} I_{q}(2\sqrt{uv}),

where q = \frac{2ab}{\sigma^2}-1, u = c r_t e^{-aT}, v = c r_{t+T}, and I_{q}(2\sqrt{uv}) is modified Bessel function of the first kind of order q.

Due to mean reversion, as time becomes large, the distribution of r_{\infty} will approach a gamma distribution with the probability density of:

f(r_{\infty};a,b,\sigma)=\frac{\omega^\nu}{\Gamma(\nu)}r_\infty^{\nu-1}e^{-\omega r_\infty},

where \omega = 2a/\sigma^2 and  \nu = 2ab/\sigma^2 .

Properties

Calibration

The continuous SDE can be discretized as follows

 r_{t+\Delta t}-r_t = a (b-r_t)\,\Delta t + \sigma\, \sqrt{r_t \Delta t} \epsilon_t ,

which is equivalent to

 \frac{r_{t+\Delta t}-r_t}{\sqrt r_t} =\frac{ab\Delta t}{\sqrt r_t}-a \sqrt r_t\Delta t + \sigma\, \sqrt{\Delta t} \epsilon_t , provided \epsilon_t is n.i.i.d. (0,1). This equation can be used for a linear regression.

Simulation

Stochastic simulation of the CIR process can be achieved using two variants:

Bond pricing

Under the no-arbitrage assumption, a bond may be priced using this interest rate process. The bond price is exponential affine in the interest rate:

P(t,T) = A(t,T) \exp(-B(t,T) r_t)\!

where

A(t,T) = \left(\frac{2h \exp((a+h)(T-t)/2)}{2h + (a+h)(\exp((T-t)h) -1)}\right)^{2ab/\sigma^2}
B(t,T) = \frac{2(\exp((T-t)h)-1)}{2h+(a+h)(\exp((T-t)h)-1)}
h = \sqrt{a^2+2\sigma^2}

Extensions

Time varying functions replacing coefficients can be introduced in the model in order to make it consistent with a pre-assigned term structure of interest rates and possibly volatilities. The most general approach is in Maghsoodi (1996). A more tractable approach is in Brigo and Mercurio (2001b) where an external time-dependent shift is added to the model for consistency with an input term structure of rates. A significant extension of the CIR model to the case of stochastic mean and stochastic volatility is given by Lin Chen (1996) and is known as Chen model. A CIR process is a special case of a basic affine jump diffusion, which still permits a closed-form expression for bond prices.

See also

References

This article is issued from Wikipedia - version of the Tuesday, March 22, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.