# Gaussian Estimation of One-Factor Mean Reversion Processes

Authors: Freddy H. Marín Sánchez; J. Sebastian Palacio
Journal: Journal of Probability and Statistics
 (2013)

---

## Abstract

We propose a new alternative method to estimate the parameters in one-factor mean reversion processes based on the maximum likelihood technique. This approach makes use of Euler-Maruyama scheme to approximate the continuous-time model and build a new process discretized. The closed formulas for the estimators are obtained. Using simulated data series, we compare the results obtained with the results published by other authors.


---

## Body

## 1. Introduction

Different applications of stochastic differential equations have attracted the attention of many researchers in relation to the theoretical advances needed for modeling and simulation of the dynamic behavior of the main variables affecting the evolution of some phenomena. Applications in fields of economics, finance, physics, insurance, and others have been drivers of these developments. Particularly, through diffusion processes have been carried out modeling interest rates and commodity prices, where they used in their modeling processes mean reversion. Some trends of these studies are based on some a priori knowledge to determine the structure of the model to be used, while other alternatives are part of a general structure, and statistical procedures through the particular structure are determined to represent the phenomenon. Whichever approach is used, in general, in the specification of the model, it is necessary to estimate parameters to achieve a proper fit of the model.Although some parameter estimation methods have been extensively studied in diffusion processes, this has a serious problem, which is that the specification of the models is done in continuous time but the data that is available for adjusting the parameters are discretely sampled in time. A common solution to this problem is to discretize the continuous-time model based on Euler-Maruyama scheme and perform procedures on the theoretical formulation resulting in discrete time. From this scheme alternatives can be proposed for parameter estimation approach among which stands out the maximum likelihood estimation.Following the classification proposed by Yu and Phillips [1], different parameter estimation procedures can be classified globally into three categories according to how to proceed from the model formulation: discretization of continuous process and estimation of the parameters with the resulting discretized process, estimation of the transition probability function for the continuous-time model, and independent estimations of the trend and volatility components by kernel methods.Thus, among the methods that are based on the discretization of the continuous-time model is Nowman [2], who considers a stochastic differential equation with two unknown variables in the drift term and two in the volatility, elasticity allowing the variance. Similarly, Yu and Phillips [1] present an improvement to the method proposed by Nowman [2], which is based on a nonuniform sampling of the observations, and the estimation process makes use of the equivalence of the methods of maximum likelihood and least squares under the conditions laid down in the model. Alternatively, Yu [3] proposes a new scheme that is based on the idea of transforming a stochastic differential equation (EDE) in its discrete-time equivalent, an AR(1) to determine the formulas that approximate the bias error present in the estimation of one of the parameters.Some of the works that could be classified as methods in which it seeks to find the transition probability function and use Martingale properties to approximate the likelihood function are Durham and Gallant [4], Valdivieso et al. [5], and Elerian et al. [6]. Durham and Gallant [4] propose a scheme based on simulation. Establishing different model modifications, Pedersen [7] improves computational performance, and characteristics of the model are studied using simulated synthetic series. Koulis and Thavaneswaran [8] study combined martingale estimation functions for discretely observed interest rates models. They derive closed-form expressions for the first four conditional moments of the CIR model via Ito’s formula and for extended versions of the CIR model using the Milstein’s approximation and construct the optimal estimating functions. They also show that combined estimating functions have more information when the conditional mean and variance of the observed process depend on the parameter of interest. Valdivieso et al. [5] proposed a method for estimating the transfer function through the inversion of the characteristic function of the process and the use of Fourier transform. This methodology is developed for some particular cases of D-OU processes (Ornstein-Uhlenbeck processes with broken car distribution D). Elerian et al. [6] used simulation-based techniques for estimating the transition probability depending on the sample to incorporate auxiliary dummy observations. The parameters are obtained by simulation and optimization using numerical algorithms.Methods based on kernel techniques approximate density functions of the components and spreading tendency of the process. The works of Jiang and Knight [9] and Florens-Zmirou [10] are within this framework.In this paper, we propose an alternative scheme for parameter estimation in one-factor mean reversion models, using the Euler-Maruyama discretization. Using the properties of discretized error normality, the maximum likelihood technique to obtain closed formulas for the estimators applies. We present some performance measures of the maximum likelihood estimators obtained with the proposed method in simulated processes, and we compare them with the results published by Nowman [2], Yu and Phillips [1], and Durham and Gallant [4].This paper is organized as follows. In Section2, a brief description of mean reversion processes of one factor with constant parameters. Some Gaussian estimation techniques are presented in Section 3. Section 4 develops the proposed method, and Section 5 presents some experimental numerical results with different sampling rates.


## 2. Mean Reversion Processes

Mean reversion processes have been considered as a natural choice to describe the dynamic behavior of interest rates in Vasicek [11], Brennan and Schwartz [12], Cox et al. [13], and Chan et al. [14]; daily prices of some commodities such as oil and natural gas in Lari-Lavassani et al. [15] have a fundamental role in modeling spot prices when they exhibited long-term trends, such as electrical energy (Pilipovic [16]).Mean reversion processes of one factor with constant parameters can be written as(1)dXt=α(μ-Xt)dt+σXtγdBt;t∈[0,T]
with initial condition X0=x, where α>0, σ>0, μ∈ℝ, and γ∈[0,3/2] are constants and {Bt}t≥0 is a unidimensional standard Brownian motion defined on a complete probability space (Ω,ℱ,ℙ).Theα parameter is called reversion rate, μ is the mean reversion level or trend of long-run equilibrium, σ is the parameter associated with the volatility, and γ determines the sensitivity of the variance to the level of Xt.The model (1) is known as CKLS proposed by Chan et al. [14] and is a generalization of the models of interest rate in the short term with constant parameters in which γ = 0, 1/2, and 1 or γ=3/2 with α=0.In general, mean reversion processes are stationary. For example, note that (1) can be written in integral form as
(2)Xt-X0=α∫0t(μ-Xs)ds+σ∫0tXsγdBs.Taking the expected value,E[Xt]-E[X0]=α∫0t(μ-E[Xs])ds and thus we obtain the differential equation
(3)m˙(t)=α(μ-m(t)),m(t)=E[Xt].The solution of this equation ism(t)=μ+(X0-μ)e-αt, so limt→∞E[Xt]=limt→∞m(t)=μ.It can be concluded then thatm(t)=μ is the equilibrium level for the linear differential equation. In other words, the expected value of the process m(t) returns in the long term to μ. In a finite-time horizon, μ plays the role of attractor in the sense that when Xt>μ the trend term α(μ-Xt)<0, and therefore, Xt decreases and when Xt<μ a similar argument establishes that Xt grows.The stationarity condition is the key to perform the Gaussian estimation of parametersα,σ, and μ, assuming that the elasticity (γ=0,1/2,1,3/2) is given a priori.Although the general model (1) derives various interest rate models studied in the literature, our interest is focused primarily on those models that are strictly stationary.For example, whenγ=0, it is said that the process has additive noise so that we get the linear stochastic differential equation
(4)dXt=α(μ-Xt)dt+σdBt;X0=x,t∈[0,T].This model is also known as Ornstein-Uhlenbeck mean reversion process type used by Vasicek [11] in modeling the equilibrium price of the coupon. This Gaussian process has been widely used in the valuation of options on bonds, futures, options on futures, and other financial derivative instruments, while in Lari-Lavassani et al. [15] it has been used to study the dynamic behavior of the prices of oil and natural gas.Although the model of Vasicek [11] is analytically tractable, it has some drawbacks. Since the interest rate in the short term is normally distributed, for each t∈[0,T] exists a probability that Xt is negative, and this is not reasonable from an economic standpoint. Another drawback of this model is that conditional volatility assumes changes in the interest rate that are constants, independent of the level of Xt.Ifγ=1 it is said that the process has proportional noise and gives the nonhomogeneous linear stochastic differential equation given by
(5)dXt=α(μ-Xt)dt+σXtdBt;X0=x,t∈[0,T].This model is used by Brennan and Schwartz [12] to obtain a numerical model of convertible bond prices. This process is also used by Courtadon [17] in developing a model for option pricing on discount bonds.Whenγ=1/2, we obtain the nonlinear stochastic differential equation
(6)0dXt=α(μ-Xt)dt+σXt1/2dBt,2αμ≥σ2,X0=x,t∈[0,T].This process is naturally linked to theχ centered-square distribution and is the (square root) version of the model of Brennan and Schwartz [12], proposed by Cox et al. [18], which produces types of instantaneous short-term interest strictly positive.This model has also been used extensively in the development of pricing models for derivative instruments sensitive to interest rates. Examples include pricing models for mortgage-backed in Dunn and McConnell [19], models of options for discount bonds in Cox et al. [18], future and option pricing models for futures in Ramaswamy and Sundaresan [20], the swap valuation models in Sundaresan [21], and option pricing models for the convenience yield in Longstaff [22].Finally, for the special case in whichγ=3/2 and α=0, equation (1) takes the form
(7)dXt=σXt3/2dBt,X0=x,t∈[0,T].
This model is introduced by Cox et al. [13] in their study of variable-rate (VR) securities. A similar model is also used by Constantinides and Ingersoll [23] for value bonds in the presence of taxes.Table1 shows the summary of the models associated with (1) that are of interest in this work (Vasicek, Brennan-Schwartz, CIR-SR,CIR-VR, and CKLS).Table 1
Summary of the models.
Model
α
μ
σ
γ
Vasicek
0
Brennan-Schwartz
1
CIR SR
1
/
2
CIR VR
0
3
/
2
CKLS
[0,3/2]


## 3. Gaussian Estimation of Diffusion Models

Many theoretical models in economics and finance are formulated in continuous time as either a diffusion or diffusion systems. For this reason, identification, estimation, and, therefore, the investigation of the asymptotic properties of samples of the estimators of these processes turn out to be a difficult task. These problems arise mainly due to the lack of availability of continuous sampling observations. In fact, the estimation of diffusion models has been considered in the literature for many years, and in most scientific journals it only refers to the estimate considering continuous sampling observations, although the data that these models described can only be sampled at discrete points in time (See, e.g., Jiang and Knight [9] and references therein.)There are various techniques and methods for estimation of models of mean reversion of a single factor; with different structural forms in the trend and distribution functions, most simulation techniques-based on statistical and quite rigorous arguments. The work of Jiang and Knight [9] proposes a parametric identification and estimation procedure for diffusion processes based on discrete sampling observations. The nonparametric kernel estimator for diffusion functions developed in this work deals with general diffusion processes and avoids any functional form specification for either the trend function or distribution function. They show that under certain regularity conditions, the nonparametric estimator function is consistent, and timely dissemination also asymptotically follows a mixed normal distribution. In Jiang and Knight [24], there are some numerical experimental results obtained by a Monte Carlo study using the results of Jiang and Knight [9], which include parameter estimation for Ornstein-Uhlenbeck process. In this sense, the lack of analytical solutions has hampered empirical evidence and validation of alternative hypotheses about these continuous-time models.Although the maximum likelihood method is a widely used parametric method, the analytical expressions for the likelihood functions exist only for processes of mean reversion with additive noise. In fact, in most of the work on estimation of diffusion models, processes become more complex processes of Ornstein-Uhlenbeck type. For example, Giet and Lubrano [25] consider a diffusion process in which the elasticity constant of the nonlinear term of diffusion can be freely chosen, achieving after reducing it a transformation to an Ornstein-Uhlenbeck process. The major drawback is that there is no degree of freedom in the form of the trend function which is completely determined by the shape of the distribution function and the choice of the underlying linear model. Valdivieso et al. [5] proposed a methodology to estimate maximum likelihood parameters of a one-dimensional stationary process of Ornstein-Uhlenbeck type, which is driven by general Levy process and which is constructed through a distribution D self-decomposed. This approach is based on the inversion of the characteristic function and using the discrete fractional fast Fourier transform. The maximization of the likelihood function is performed using numerical methods.The framework proposed in our work has some features in common with Nowman [2], Yu and Phillips [1], and Elerian et al. [6]. The similarity of our approach lies in the use of Euler numerical scheme, using discretely sampled data and maximum likelihood estimation, but the main difference is that in our proposal does not need any numerical optimization because it is possible to obtain a closed formula for the estimators.For example, in Elerian et al. [6] developed a technique based on Bayesian estimation of nonlinear EDEs approximating the Markov transition densities of broadcasts, including Ornstein-Uhlenbeck processes and elasticity models for interest rates, when observations are discretely sampled. The frame of the estimate is based on the introduction of auxiliary data to complete the diffusion latent missing between each pair of measurements. Combine methods of Monte Carlo Markov Chain algorithm based on the Metropolis-Hastings, in conjunction with the numerical scheme of Euler-Maruyama, to sample the posterior distribution of latent data and model parameters.Another alternative for Gaussian estimation was proposed by Nowman [2]. In this paper, an estimate using the exact discretization of the model is proposed,
(8)dXt=α(μ-Xt)dt+σX(sΔ)γdBt,sΔ≤t<(s+1)Δ,
assuming that volatility stays constant on the intervals [sΔ,(s+1)Δ); s=0,1,….In this way, we obtain an equation of the form(9)Xt=e-αΔX(t-Δ)+μ(1-e-αΔ)+ηt,
where the conditional distribution ηt∣Ft-1B≈N(0,(σ2/2α)(1-e-2αΔ){Xt-1}2γ).However, despite the obtainable approximate conditional distribution, this approach does not correspond to the Euler discretization applied to (1). This estimation technique requires the discrete model corresponding to (2) which is given by
(10)Xt=e-αX(t-1)+μ(1-e-α)+ηt,(t=0,1,…,T),
where ηt, (t=0,1,…,T) satisfies the conditions E[ηsηt]=0; (s≠t) and
(11)E[ηt2]=∫t-1te-2(t-τ)ασ2{X(t-1)}2γdτ=σ22α(1-e-2α){X(t-1)}2γ=mtt2.
This states that the logarithm of the Gaussian likelihood function may be written as
(12)L(θ)=∑t=1T‍[2logmtt+{Xt-e-αX(t-1)+μ(1-e-α)mtt}2]=∑t=1T‍[2logmtt+εt2],
where ε=[ε1,…,εT] is a vector whose elements can be calculated.The optimization ofL(θ) is carried out using the methods developed by Bergstrom [26–30] and Nowman [2].Yu and Phillips [1] propose a Gaussian estimation using random time changes which is based on the idea that any continuous time martingale can be written as a Brownian motion after a suitable time change. It considers the explicit solution (1) to obtain the expression
(13)X(t+h)=μ(1-e-αh)+e-αhX(t)+∫0hσe-α(h-τ){X(t+τ)}γdBτX(t+h)=μ(1-e-αh)+e-αhX(t)+∫0hσe.e-α(h-τ)forh>0.LetM(h)=σ∫0he-α(h-τ){X(t+τ)}γdBτ.  M(h) is a continuous martingale (In Yu and Phillips [1],   M(h) was assumed as a continuous martingale, but this assumption is generally not warranted, and [M]h does not induce a Brownian motion as Yu and Phillips showed in [31]. In addition, they establish that a simple decomposition splits the error process into a trend component and a continuous martingale process.) with quadratic variation
(14)[M]h=σ2∫0he-2(h-τ)α{X(t+τ)}2γdτ.They introduce a sequence of positive numbers{hj} which deliver the required time changes. For any fixed constant a>0, let
(15)hj+1=inf{s∣[Mj]s≥a}=inf{s∣σ2∫0se-2(s-τ)α{X(tj+τ)}2γdτ≥a}.Then a succession of points{tj} is constructed using the iterations tj+1=tj+hj+1 with t1=0 to obtain
(16)X(tj+1)=μ(1-e-αhj+1)+e-αhj+1X(tj)+Mhj+1,
where Mhj+1=B(a)≈N(0,a).Note that according to (16), (12) takes the form
(17)L(θ)=1N∑j‍[{{X(tj+1)-μ(1-e-αhj+1)-e-αhj+1X(tj)a}2}2loga+{X(tj+1)-μ(1-e-αhj+1)-e-αhj+1X(tj)a}2].Since the data are observed in discrete time, the time change formula cannot be applied directly and therefore must be approximated by(18)hj+1=Δmin{s∣∑i=1s‍σ2e-2α(s-i)Δ{X(tj+iΔ)}2γ≥a}.This should be selected so thata is the volatility of the error term unconditional in the equation
(19)X(t+Δ)=μ(1-e-αΔ)+e-αΔX(t)+ε.The implementation of this method is based on the estimate given by Nowman [2] to find the parameters estimates σ and γ. Furthermore, the estimation of the remaining parameters μ and α is carried out using the numerical optimization algorithm by Powell [32].


## 4. Gaussian Estimation of One-Factor Mean Reversion Processes with Constant Parameters

The development of this section will focus on the description of the proposed technique for the Gaussian estimation of mean reversion models of a single factor, with linear and nonlinear functions of diffusion constant parameters and linear trend functions with constant parameters.As the density function for most diffusion processes is not known analytically, the proposed idea is not to find or approximate the transition density function or the marginal density function of the process. Otherwise, we use the numerical approximation of Euler on the EDE to build a new discrete process, Huzurbazar [33], and use the fact that in this approximation the normalized residuals resulting from Brownian process are independent standard normal random variables. As pointed out by Pedersen [7], Kohatsu-Higa and Ogawa [34], and Bally and Talay [35], the density function of the discrete process obtained with Euler scheme converges to the density function of the continuous process, and as discussed by Zellner [36], these residues can be interpreted as a vector of unknown parameters with the advantage of knowing their distribution. So we build a new variable for the models defined in (1) so that it can adapt to the observed data and later can be evaluated by examining the normalized residuals that will be useful in the estimation process. By knowing the density function of this new process we can find discretized likelihood function and use the conventional method of maximum likelihood estimation to find a closed formula for the parameters.Consider the stochastic differential equation (1), with the initial condition X0=x. The terms α(μ-Xt) and σXtγ being the trend and diffusion functions, respectively, that depend on Xt and a vector of unknown parameters θ=[α,μ,σ,γ], where {Bt}t≥0 is a unidimensional standard Brownian motion defined on a complete probability space (Ω,ℱ,ℙ). The trend function describes the movement of the process due to changes in time, while the diffusion function measures the magnitude of the random fluctuations around the trend function. Suppose that the conditions under which we guarantee the existence and uniqueness of the solution of Mao [37] are satisfied. Assume further the process observations in discrete time Xt=X(τt) at times {τ0,τ1,…,τT} where Δt=τt-τt-1≥0, t=1,2,…,T. The objective is to estimate the unknown parameter vector θ given the historical information of the sample X=(X0,X1,…,XT)′ with equally spaced data. Note that the explicit solution of (1) is given by
(20)Xt=μ+(X0-μ)e-αt+∫0tσe-α(t-s)XsγdBs.In the context of the likelihood, the estimation ofθ is based on the likelihood function
(21)logL(X1,X2,…,XT∣X0,θ)=∑t=1T‍[logg(Xt+1∣Xt,θ)],
where g(Xt+1∣Xt,θ) are the Markov transition densities, or L(θ∣X1,X2,…,XT)=f(X1,X2,…,XT,θ) and f(X1,X2,…,XT,θ) is the joint density function.Although the differential equation (1) can be solved analytically on [0,T], the functions g(Xt+1∣Xt,θ) and f(X1,X2,…,XT,θ) are not available in closed form, and, therefore, the estimate can not be accomplished, except for the case where γ=0.Note that the exact discretized solution of (1) takes the form
(22)X(τt)=μ+(X(τt-1)-μ)e-α(τt-τt-1)+∫τt-1τσe-α(τt-s)X(s)γdBs
which can be approximated by
(23)X(τt)=μ+(X(τt-1)-μ)e-αΔt+σe-αΔtXγ(τt-1)ΔBt,
where ΔBt=Bτt-Bτt-1.Thus,(24)X(τt)=X(τt-1)+α(μ-X(τt-1))Δt+Xγ(τt-1)ετt,
where ετt~N(0,σ2Δt), corresponding to the numerical approximation of Euler.Now define the new variable(25)Yi=Xi-Xi-1-α(μ-Xi-1)ΔXi-1γ=εi,Δconstant.Note thatYi, i=1,2,3,…,T are independent random variables distributed normal, N(0,σ2Δ).In terms of the density function can be written:(26)f(Y1:θ)=(12πσ2Δ)1/2×exp[-12σ2Δ(X1-X0-α(μ-X0)ΔX0γ)2]f(Y2:θ)=(12πσ2Δ)1/2×exp[-12σ2Δ(X2-X1-α(μ-X1)ΔX1γ)2]f(Y2:θ).⋮×exp[-12σ2Δ(X2-X1-α(μ-X1)ΔX1γ)2]f(YT:θ)=(12πσ2Δ)1/2×exp[-12σ2Δ(XT-XT-1-α(μ-XT-1)ΔXT-1γ)2].Since eachYi is independent can be expressed the joint density function as their marginal densities product, that is,
(27)f(Y1,Y2,…,YT)=f(Y1:θ)f(Y2:θ)⋯f(YT:θ).Consequently, the likelihood function is given by(28)L(θ∣YT)=(12πσ2Δ)T/2×exp[-12σ2Δ∑i=1T‍(Xi-Xi-1-α(μ-Xi-1)ΔXi-1γ)2].Taking the natural logarithm on both sides, it follows that(29)logL=-T2log(2πσ2Δ)-12σ2Δ[∑i=1T‍(Xi-Xi-1-α(μ-Xi-1)ΔXi-1γ)2].Now consider the problem of maximizing the likelihood function given by(30)∂log(L)∂θ=0→,i.e.,θ^=Argθmax(log(L)).The estimation process leads to the estimated parametersα and μ, are the solutions of which the next equation system, where it is assumed that the value of γ is known. Consider
(31)A-E(1-αΔ)-αμΔB=0C-B(1-αΔ)-αμΔD=0,
where
(32)A=∑i=1T‍XiXi-1Xi-12γ,B=∑i=1T‍Xi-1Xi-12γ,C=∑i=1T‍XiXi-12γ,D=∑i=1T‍1Xi-12γ,E=∑i=1T‍(Xi-1Xi-1γ)2.
Thus, the estimated parameters are given by
(33)α^=ED-B2-AD+BC(ED-B2)Δ,μ^=A-E(1-α^Δ)α^BΔ.Additionally, the parameterσ is determined by an equation independent of the system of (31) and is completely determined by the estimated parameters:
(34)σ^=1TΔ∑i=1T‍(Xi-Xi-1-α^(μ^-Xi-1)ΔXi-1γ)2.This estimation scheme has several important qualities. First, you come to a closed formula for Gaussian estimation without numerical optimization methods. Second, the computational implementation is simple and does not require sophisticated routines to perform calculations very extensive.


## 5. Simulations

This section will present comparisons of the results with those published by other authors using the techniques described above. Comparisons are made with the estimate based on simulated data sets to analyze statistical properties of the proposed estimators.Yu and Phillips [1] present the results of a simulation study with the Monte Carlo technique in comparing their method with Nowman’s [2]. To make this comparison are suggested three scenarios: to resolve data monthly, weekly, and daily. In all cases the basic pattern of square root is defined as
(35)000dXt=(λ+βXt)dt+σXtγdBt;withλ>0,β<0,σ>0,γ=0.5.Note that with appropriate substitution of this model parameters are equivalent to those described in (1) (The equivalence of models (1) and (35) is satisfied by having λ=αμ and β=-α. Originally in Yu and Phillips [1], the parameter λ in (35) is defined as α but was changed to avoid confusion with the model (1).)The specific parameters used in each stage are shown in Table2.Table 2
Monte Carlo study parameters.
Monthly
Weekly
Daily
Size
500
1000
2000
λ
0.72
3
6
β
−0.12
−0.5
−1
σ
0.6
0.35
0.25
γ
0.5
0.5
0.5To perform the first comparison parameters are set as specified in Table2 and obtained a series simulated. Then you set the value to γ=0.5 and proceed to the estimation by the method of Section 4 to obtain λ^, β^, and σ^. This procedure was repeated 1000 times (as in Yu and Phillips [1]). Some statistical estimators are presented in Tables 3, 4, and 5.Table 3
Comparison  of monthly data series with Nowman, Yu-Phillips, and proposed methods.
Nowman
Yu, Phillips
Nowman-Yu, Phillips
Proposed
λ
β
λ
β
σ
γ
λ
β
σ
Mean
1,3440
−0,2332
1,2330
−0,2275
0,6173
0,4919
1,2797
−0,2200
0,5989
Bias
0,6240
−0,1132
0,5130
−0,1075
0,0173
−0,0081
0,5597
−0,1000
−0,0011
Var
0,5897
0,1713
0,5732
0,1589
0,0099
0,0071
0,5406
0,0149
0,0004
MSE
0,9791
0,1841
0,8363
0,1705
0,0102
0,0071
0,8533
0,0249
0,0004
The results of the methods of Nowman and Yu, Phillips are taken into account, Yu and Phillips [1]. Our results were obtained with 1000 simulations with 500 observations each. The real values of the parameters are λ=0.72, β=-0.12, σ=0.6, and γ=0.5.Table 4
Comparison of weekly data series with Nowman, Yu-Phillips, and proposed methods.
Nowman
Yu, Phillips
Nowman-Yu, Phillips
Proposed
λ
β
λ
β
σ
γ
λ
β
σ
Mean
4,4090
−0,7320
4,1650
−0,7011
0,3762
0,4925
4,2038
−0,7022
0,3499
Bias
1,4090
−0,2320
1,1650
−0,2011
0,0262
−0,0075
1,2038
−0,2022
−0,0001
Var
4,0663
0,1109
3,8608
0,1030
0,0172
0,0357
3,5427
0,0992
0,0001
MSE
6,0855
0,1647
5,2180
0,1435
0,0179
0,0358
4,9883
0,1400
0,0001
The results of the methods of Nowman and Yu, Phillips are taken into account, Yu and Phillips [1]. Our results were obtained with 1000 simulations with 1000 observations each. The real values of the parameters are λ=3, β=-0.5, σ=0.35, and γ=0.5.Table 5
Comparison  of daily data series with Nowman, Yu-Phillips, and proposed methods.
Nowman
Yu, Phillips
Nowman-Yu, Phillips
Proposed
λ
β
λ
β
σ
γ
λ
β
σ
Mean
9,8250
−1,5360
8,8440
−1,4750
0,2521
0,4970
8,9542
−1,4933
0,2500
Bias
3,8250
−0,5360
2,8440
−0,4750
0,0021
−0,0030
2,9542
−0,4933
0,0000
Var
19,1959
0,5219
17,8220
0,4798
0,0244
0,0746
15,5818
0,4334
0,0000
MSE
33,8266
0,8092
25,9100
0,7054
0,0244
0,0746
24,2938
0,6763
0,0000
The results of the methods of Nowman and Yu, Phillips are taken into account, Yu and Phillips [1]. Our results were obtained with 1000 simulations with 1000 observations each. The real values of the parameters are λ=6.0, β=-1.0, σ=0.25, and γ=0.5.Tables3, 4, and 5 show that the method proposed by Nowman [2] gives good estimates for the parameters σ and γ, as manifested by Yu and Phillips [1] it becomes clear that a starting point for the method is defined by the latter. Furthermore, through the resulting estimate of the parameter σ, our method becomes much more accurate than Nowman [2], introducing a bias less than 0.2% for the monthly data series and a small value for cases of weekly and daily series.In the case ofλ and β are parameters observed inefficient performance (as discussed Yu and Phillips [1, page 220]) in the method of Nowman, λ yielding bias of 87%, 47% and 64% for data monthly, weekly and daily, respectively. Also, for β the bias present in the estimate is 94%, 46%, and 54% for the respective monthly, weekly, and daily series. The method proposed by Yu and Phillips has a better performance than Nowman, measured from the bias in the parameters, reducing it by 15%, 8%, and 16% for the case of λ and 5%, 6%, and 6% in the case of β for the respective monthly, weekly, and daily series. Our method has a performance similar to that of Yu and Phillips for daily data series, reducing the bias regarding Nowman method in 15% for λ and 4% for β. For weekly data series for reducing λ bias, regarding the method of Nowman, is 7% whereas for β is 6%, as reduction Yu and Phillips method.An important point to note in the results obtained by our method is the condition that the parameterγ is known a priori. This is necessary in order to find closed formulas for the estimators.In the approach to the model proposed by Nowman and Yu-Phillips,dXt=(λ+βXt)dt+σXtγdBt is implied to have particular difficulty in estimating the parameter of speed of reversion. To make this condition clear, Table 6 shows the results obtained by the method of Durham and Gallant [4] and our method to estimate the parameters of the model dXt=α(μ+Xt)dt+σXtγdBt with α>0, μ>0, σ>0, γ=0.5. The results are the average of 512 estimations of monthly series with 1000 observations.Table 6
Comparison of Durham and Gallant, and proposed methods.
Parameter(*)
Real
Bias(a)
Bias(b)
Real
Bias(a)
Bias(b)
α
0,5
0,048900
0,037700
0,5
0,059730
0,021400
μ
0,06
0,000610
0,000127
0,06
−0,000310
−0,000021
σ
0,15
0,000200
0,000032
0,22
0,000120
−0,000046
α
0,5
0,043780
0,032400
0,4
0,045890
0,032200
μ
0,06
0,000050
0,000016
0,06
0,000770
0,000046
σ
0,03
0,000010
−0,000003
0,15
0,000080
−0,000025
The data are obtained with 512 simulations monthly series (dt = 1/12) with 1000 observations each, and the parameters are listed as real.
(a)Results of the method proposed by Durham and Gallant.
(b)Results of the method proposed in this paper.
(
*
)In Durham and Gallant [4]: μ=θ1, α=θ2, and σ=θ3.The results shown in Table6 show that the errors in estimation of μ and σ are small (less than 1%) for both methods. On the other hand, the estimated speed of reversion has an error of about 10% resulting in Durham and Gallant, while our method has an error close to 7%. From this analysis, it is clear that the parameter α has a higher error in estimating than the other parameter estimates.In this regard, the estimation error ofα in Yu [3] provides an analysis of the bias in the parameter estimation inherent reversion total observation time which is sampled in the continuous process, for the special case, where the long-term average μ is zero. Reference is made to the fact that the bias in parameter depends asymptotically on total observation time which of the process and not the number of observations made in this period.Figures1, 2, and 3 show the results of the estimation of the parameters for different observation horizons with different sampling frequencies. (The results are obtained as the average of the individual estimates on 1000 simulated series in 4 sampling scenarios: (a) bimonthly series with 500 observations, (b) monthly series with 1000 observations; (c) weekly series with 4333 observations and (d) daily series with 21,000 observations. The reference model for the simulations was dXt=α(μ-Xt)dt+σXtγdBt, with α=1.5, μ=0.5, γ=0.5, and σ=0.4.)Figure 1
α estimation results. dt = 1/6, 1/12, 1/52, and 1/252, n = 500, 1000, 4333, and 21000, respectively.Figure 2
μ estimation results. dt = 1/6, 1/12, 1/52, and 1/252, n = 500, 1000, 4333, and 21000, respectively.Figure 3
σ estimation results. dt = 1/6, 1/12, 1/52, and 1/252, n = 500, 1000, 4333, and 21000, respectively.Figure1 shows the evolution of the estimate of α for observations from 10 to 83 years. Behavior of the estimates in the four scenarios can be established empirically by the dependence of the bias on the estimate and the horizon of observation as they have more years of observation bias decreases, with no evidence of a significant difference between the estimates different frequencies.From Figure2 it is observed that after few observations, virtually independent of the year of observation, estimations for μ are obtained with a bias of less than 1%, showing the efficiency of the estimation of the long-term average reversion processes average. In the case of estimation of σ, Figure 3 shows that the bias in the parameter estimates is more dependent on the number of observations which is the total time of observation of the process. This is reflected in the fact that the estimates are in series with a higher sampling rate reaches a lower bias for the same observation horizon than lower sample rate and therefore have fewer observations. On the other hand, as for the case of the estimate of μ, the bias present in the estimates of σ is practically negligible (less than 1%) in most scenarios.


## 6. Conclusions and Comments

This paper presents a new alternative for Gaussian estimation models of mean reversion of one factor under the Euler numerical approximation that includes, among others, some models of interest rates in the short term.Monte Carlo simulations show a good approximation of the diffusion term parameter but less approximation in a trend term parameter. The estimation results have almost the same order as that in Nowman’s [2] and in some cases better than that in Durham and Gallant’s [4]. It should be noted that the proposed methodology outperforms the estimate of σ compared to the other alternatives presented. These results become more relevant since in practice the estimation of volatility is extremely useful in the evaluation of financial derivatives, in which the trend term loses significance when it is absorbed by the processes of change measurement and valuation risk neutral.The estimate based on this new method, though does not eliminate the bias in estimation of the reversion rate parameter for a small sample, is in some cases comparable with the results in Yu and Phillips [1], if samples have longer time horizons.This estimation scheme has several important qualities. First, we come to a closed formula for Gaussian estimation without numerical optimization methods. Second, the computational implementation is simpler and does not require sophisticated routines to perform extensive calculations.



---

*Source: 239384-2013-10-20.xml*