When AR(1) Calibration Does Not Fit an Ornstein Uhlenbeck Process
Summary
The note presents the standard mapping from an exact discrete time Ornstein Uhlenbeck process to an AR(1) regression. In that mapping, the AR coefficient determines the continuous time mean reversion rate through its logarithm, while the intercept and residual variance determine the long run mean and diffusion volatility. The question arises because an ordinary least squares fit returns a negative AR coefficient, making the stated logarithm based formulas unusable; it also asks whether taking the absolute value is appropriate for estimating half life.
The replies caution that a time series should not be assumed to follow AR(1) merely because that model is convenient. The data may require a different ARIMA specification, and fitting a desired OU process calls for methods suited to diffusion parameter estimation. The answers cite research on estimation bias and bias correction but do not provide a resolution for negative coefficients or a worked calibration. The formulas therefore depend on model fit and their assumptions, rather than supporting an automatic absolute value adjustment.
Key ideas
- The OU to AR(1) parameter mapping assumes the discrete time coefficient is compatible with the OU model.
- A negative fitted AR coefficient makes the logarithm based mean reversion formulas invalid.
- The replies advise checking the time series model fit instead of presuming it is AR(1).
- Alternative ARIMA specifications may be more appropriate for the data.
- Diffusion parameter estimation can have bias, and the cited papers address bias correction.
Tags
Full text
# Calibrating OU parameters using AR(1)
# Calibrating OU parameters using AR(1)
I have a mean reverting time series and want to find the Ornstein-Uhlenbeck (OU) parameters of it. I researched the internet and found that we can calibrate the model as a simple AR(1) process, $$\text dS_{t} = \lambda(\mu-S_t)\text dt+\sigma \text dW_t,$$ where $\lambda$ is the mean reversion rate, $\mu$ the mean and $\sigma$ the volatility.
The exact solution of the above SDE is \begin{align*} S_{i+1} = S_i e^{-\lambda\delta} + \mu(1-e^{-\lambda\delta}) + \sigma \sqrt{\frac{(1-e^{-2\lambda\delta})}{2\lambda}}N_{0,1}, \tag{1} \end{align*} where $\delta$ is a small time increment.
An AR(1) process is \begin{align*} S_{i+1} = aS_i+b+\varepsilon. \tag{2} \end{align*} Comparing the AR(1) process with the exact solution of the SDE, we can get the following relations \begin{align*} \lambda &= -\frac{\ln a}{\delta} \tag{3} \\ \mu&=\frac{b}{1-a} \tag{4} \\ \sigma &= \text{stdev}(\epsilon) \sqrt{\frac{-2\ln a}{\delta(1-a^2)}} \tag{5} \end{align*}
I fitted a simple OLS model for (2) and $a$ turns out to be negative (e.g., $a=-0.03$). We can neither obtain $\lambda$ nor $\sigma$ as they have $\ln(a)$ and we cannot take log of a negative value.
My question is very similar to link. I looked into these supporting stack exchange links (1, 2, 3) but none of the links could give a solution to my issue
Also I understand when estimating half life using AR(1) we should use $-\frac{\ln(2)}{\ln(|a|)}$ and half life of OU is $\frac{\ln(2)}{\lambda}$. Associated link. Should I take absolute value of $a$ in Equations (3) and (5)?
## Answer by eruiz (score 1)
https://quant.stackexchange.com/a/61641
You cannot assume you can fit your time series to AR(1). You should fit your data and see what ARIMA coefficients it gives you. If the data you're using is not an AR(1) model then obviously it wouldn't work. ARIMA(1,1,0) is common for example but just last month I ran an ARIMA on 3M and got like ARIMA(3,5,0) so play around with it.
## Answer by Andrea Di Iura (score 0)
https://quant.stackexchange.com/a/66383
As stated by eruiz in his answer you cannot assume that an AR(1) model can fit your data. Nevertheless, if you are interested in a OU process it might be useful to look at this paper:
Cheng Yong Tang and Song Xi Chen. Parameter estimation and bias correction for diffusion processes. Journal of Econometrics, 149(1):65–81, 2009.
This is in the same vein of mark leeds comment about bias and variance reduction.
Also, this paper might be useful:
Zi-Yi Guo. Out-of-sample performance of bias-corrected estimators for diffusion processes. Journal of Forecasting, 40(2):243–268, 2021.Shown in full with attribution under the source's licence. Licence: CC BY-SA 4.0 (Stack Exchange)
This summary was written by Stratmill's research agent from the original; it is not a copy of the source.