Advances in Mathematical Physics
Volume 2017 (2017), Article ID 7150203, 7 pages
https://doi.org/10.1155/2017/7150203

Research Article

Stochastic Volatility Effects on Correlated Log-Normal Random Variables

Yong-Ki Ma

Department of Applied Mathematics, Kongju National University, Chungcheongnam-do 32588, Republic of Korea

Correspondence should be addressed to Yong-Ki Ma; ykma@kongju.ac.kr

Received 12 July 2017; Revised 25 September 2017; Accepted 9 December 2017; Published 28 December 2017

Academic Editor: Antonio Scarfone

Copyright © 2017 Yong-Ki Ma. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The transition density function plays an important role in understanding and explaining the dynamics of the stochastic process. In this paper, we incorporate an ergodic process displaying fast moving fluctuation into constant volatility models to express volatility clustering over time. We obtain an analytic approximation of the transition density function under our stochastic process model. Using perturbation theory based on Lie–Trotter operator splitting method, we compute the leading-order term and the first-order correction term and then present the left and right skew scenarios through numerical study.

1. Introduction

The multivariate log-normal distribution is a widely used stochastic model in social sciences. What is the probability of the sum or difference of log-normal random variables? The solution to this question has wide applications in many fields such as finance [1, 2], actuarial science [3, 4], and physics [5]. Especially, in physics, examples include wave transmittance in random media, dissipation rate of turbulence energy, and temporal fluctuations of some nonlinear systems. However, to the best of the author of this paper’s knowledge, almost nothing is known about the question yet. So, finding analytical approximations to it is as important as an alternative approach (see [6] for details). The main advantage of the analytical approximation approaches compared to other numerical methods is that in general the first ones are much faster and precise at least under certain model parameter regime. In addition, analytic approximation formulas retain qualitative model information and preserve an explicit dependence of the results on the underlying parameters. Many approximate methods can be categorized by their scope into three classes: generally correlated [7–10], independent [11–14], and independent and identically distributed [15, 16]. The probability density function of the correlated log-normal especially is carried out to propose various approximations of the sum distribution. These approximations can be divided into two categories. The first method is approximated by another log-normal random variable. For example, the log-normal parameters are obtained by moment matching [7], the mean and the variance are computed recursively [17], and the shifted log-normal distribution [18]. A second method of these methods is introduced to construct efficient numerical integrations techniques such as a steepest-descent integration [19] and a Smolyak’s algorithm [20].

In option pricing theory, it is well-known that constant volatility for the stock price in the Black-Scholes model [21] can hardly capture the accumulated empirical evidence in financial markets since the parameters contained in the model actually change over time. A main drawback in the assumption of this model lies in the flat implied volatility surface, which is contradictory to empirical results that the implied volatilities of the equity options exhibit the smile or skew curve, in particular, from the 1987 crash onward. Among those several methods of overcoming the above drawback and relaxing the assumption of the Black-Scholes constant volatility model, a variety of stochastic volatility models have been suggested to incorporate the volatility skew. Fast mean reversion especially is one of the notable features of volatility. For example, Fouque et al. [22] estimated that its volatility reverts to the mean with characteristic time about days in high-frequency S&P 500 index data. Also, Engle [23] and Bollerslev [24] introduced the family of autoregressive conditional heteroskedasticity and generalized autoregressive conditional heteroskedasticity, respectively, to describe the evolution of the volatility of the stock price in discrete time setting and showed that econometric tests of these models reject the assumption of constant volatility and find evidence of volatility clustering over time.

Based on these observations together with the renowned contribution of stochastic volatility formula to option pricing, we propose the following fast mean-reverting (FMR) model in which asset prices are conditionally log-normal and the volatility process is a positive increasing function of a FMR process:

where and are log-normal random variables and , , and are standard Brownian motions correlated as follows:

with . Here, the parameter measures the rate at which the ergodic Markov process reaches its long-term mean value , is a function, is a constant, and the correlations and control the slope of the skew of each underlying. We assume that is large so that volatility is sufficiently fast mean-reversion when looked over the time scale of options Fouque et al. [22]. We do not specify the concrete forms of , , and since they will not play an essential role in the perturbation theory performed in this paper but the functions must satisfy a sufficient growth condition to avoid some kind of bad cases such as nonexistence of the moments of . This generalization is a clear advantage of this model in that the FMR model encompasses the models which reflect stylized facts such as a feedback effect between volatility and volatility of volatility.

2. Asymptotic Analysis

2.1. Problem Formulation

Now we will denote and as and , respectively, for simplicity. The probability distribution of the sum or difference of the two correlated log-normal distributions can be obtained by calculating the integral

where is the joint probability distribution of the two log-normal random variables and is the Dirac delta function. Unfortunately, a closed-form representation for this probability distribution still does not exist.

FMR volatility enables us to make good use of perturbation theory. To end this, we use the small positive parameter , which denotes the inverse of the rate of mean-reversion (which is assumed to be large). We suppose , where the variance of the invariant distribution of is a constant with respect to . If rewriting (2) in terms of , we have

This gives a singular perturbation problem with respect to the small parameter , which furnishes us with analytic tractability. We define the scaled FMR (SFMR) model as the FMR model used with an first-order asymptotic method in this paper.

If we compare the Heston model [25], one of representative stochastic volatility models, with the SFMR model, we can get some strengths. The Heston model has five parameters required to be estimated from market data while the SFMR model has two. In addition, the Heston model has certain restrictions on the volatility functions , , and whereas the SFMR model does not have any.

2.2. Perturbation Theory Based on Lie–Trotter Operator Splitting Method

Perturbation theory as developed by Fouque et al. [22] is a methodology utilized to find an approximated solution when the original problem is difficult to solve by separating it into more easily solvable, simple parts. If we apply the Feynman-Kac formula, we find that satisfies the following Kolmogorov backward equation:

where

Here, is the infinitesimal generator of the SFMR process . contains the mixed partial derivatives due to the correlations of the two Brownian motions and and and , respectively. is the operator of a generalized version of the two-dimensional standard Brownian motion at the volatility levels and instead of constant volatilities, respectively.

Before we solve problem (6), we write a useful lemma about the centering (or solvability) condition on the Poisson equation related to the operator as follows.

Lemma 1. If which is a solution to the Poisson equation

exists, then the centering condition must be satisfied, where the notation is the average (or expectation) with respect to the invariant distribution (namely, ) of . Then, solutions of (10) are given by the form

Proof. See Fouque et al. [22].

Using the resultant partial differential equation (PDE) (6) and expanding in powers of , one can approximate to the sum of the leading term and the first correction term as follows:

Since we focus on the first-order correction term for , we reset (12) with respect to and denote it by . Using perturbation theory based on Lie–Trotter operator splitting method, and must satisfy the following PDEs with boundary conditions, respectively:

and then one obtains the solutions and of the PDEs, respectively,

where

and the constant parameters , and are defined as follows:

See the Appendix for the relevant technique.

As a result, the distributions of the sum and difference of correlated log-normal random variables under SFMR model are shown to follow combining a shifted log-normal distribution with mixed partial derivatives of it. Also, all the original parameters are absorbed in the group parameters , , , and and the present level of the hidden process driving the fast time-scale volatility needs not be specified in the present approximation.

2.3. Numerical Experiment

In this subsection, we illustrate the effectiveness of our result (12) by showing numerical results.

Several parameter sets are obtained from Lo except for group parameters. Figure 1(a) depicts log-normal distribution (leading-order term or Lo’s result) and log-normal distributions under SFMR model with positive or negative group parameters (leading-order term plus first-order correction term), respectively. The solid line denotes the shifted log-normal distribution and the dash lines show combining the shifted log-normal distribution and the first correction term, respectively. Figure 1(b) depicts the errors calculated by subtracting the log-normal distributions with positive or negative group parameters from the log-normal distribution, respectively. We also apply the same result to the difference of the probability density function in Figure 2. Our numerical results show that the left and right skew scenarios are presented through first-order correction term and the major discrepancies appear around the peak of the probability density function. These pictures are sensitive to the choice of the involved parameters and give a lot of flexibility to the shape of the transition densities.

(a)

(b)

Figure 1: Sum of probability density function: ,  ,  ,  ,  , and .

Figure 2: Difference of probability density function: ,  ,  ,  ,  , and  .

3. Final Remarks

Stochastic processes are popular in modeling various economics and financial variables. The transition density function especially plays a key role in the analysis of continuous-time diffusion models. In this paper, we obtained an analytic approximation of correlated log-normal random variables under SFMR model.

This paper offers various possible directions for further development. Our result can be applied to pricing and hedging spread options and is to incorporate a slowly varying volatility-driving process into the SFMR model. Also, this result can provide a very useful guide for credit risk management (see [26, 27]). We leave these issues as future research topics.

Appendix

Perturbative Analysis

Now, we delineate the derivations of the PDEs (13) and (15).

We expand in powers of in order to apply the perturbation theory to the PDE problem (6):

Substituting (A.1) into (6), we reorganize the formula as follows:

For (A.2) to hold for any , each term of the equation must be zero. From the -order term, we obtain

Since the operator is the infinitesimal generator of the SFMR process , the solution of the PDE must be a constant with respect to the variable ; . From -order term,

Since each term of the operator contains -derivative, the -independence of yields . The PDE (A.4) then reduces to , so that is also independent of ; . Namely, the two terms and do not depend on the current level of the process driving the fast scale volatility. One can continue to eliminate the terms of order , , , and so on. From the constant order terms, we have

This PDE becomes

due to the -independence of as seen above. The PDE (A.6) can be considered as a Poisson equation (see Lemma 1). From the centering condition with , the leading-order has to satisfy

where

Then solves the PDE (A.7).

Next, we derive the first correction term . From the -order term, we have

Applying the centering condition with respect to to (A.9), we have

From (9) and (A.7), we get

Substituting (A.12) into (A.6), we reorganize the form as follows:

for arbitrary finite-valued functions , , and where the functions , , and are solutions of the Poisson equations, respectively:

From (A.11) and (A.13), we get

Plugging (A.16) into (A.10), a PDE for is given by

Since we focus on the first-order correction term for , we reset (A.17) with respect to and denote it by so that the fast scale correction satisfies the following PDE:

Here, the operator is expressed as

where the constant parameters , , , and are defined as follows:

It can be checked directly that is given by

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (MSIP) (no. 2015R1C1A1A01054663).

References

M. A. Milevsky and S. E. Posner, “Asian options, the sum of lognormals, and the reciprocal gamma distribution,” Journal of Financial and Quantitative Analysis, vol. 33, no. 3, pp. 409–422, 1998.

R. Carmona and V. Durrleman, “Pricing and hedging spread options,” SIAM Review, vol. 45, no. 4, pp. 627–685, 2003.

J. Dhaene, M. Denuit, M. J. Goovaerts, R. Kaas, and D. Vyncke, “The concept of comonotonicity in actuarial science and finance: applications,” Insurance: Mathematics and Economics, vol. 31, no. 2, pp. 133–161, 2002.

S. Vanduffel, T. Hoedemakers, and J. Dhaene, “Comparing approximations for risk measures of sums of nonindependent lognormal random variables,” North American Actuarial Journal, vol. 9, no. 4, pp. 71–82, 2005.

M. Romeo, V. Da Costa, and F. Bardou, “Broad distribution effects in sums of lognormal random variables,” The European Physical Journal B, vol. 32, no. 4, pp. 513–525, 2003.

S. Nadarajah, “A review of results on sums of random variables,” Acta Applicandae Mathematicae, vol. 103, no. 2, pp. 131–140, 2008.

L. F. Fenton, “The sum of log-normal probability distributions in scatter transmission systems,” IEEE Transactions on Communications Systems, vol. 8, no. 1, pp. 57–67, 1960.

F. Berggren and S. B. Slimane, “A simple bound on the outage probability with lognormally distributed interferers,” IEEE Communications Letters, vol. 8, no. 5, pp. 271–273, 2004.

C. Lam and . Tho Le-Ngoc, “Outage probability with correlated lognormal interferers using log shifted gamma approximation,” in Proceedings of the 2005 5th International Conference on Information Communications & Signal Processing, pp. 618–622, Bangkok, Thailand.

N. B. Mehta, A. F. Molisch, J. Wu, and J. Zhang, “Approximating the sum of correlated lognormal or lognormal-rice random variables,” in Proceedings of the 2006 IEEE International Conference on Communications, ICC 2006, pp. 1605–1610, Turkey, July 2006.

S. B. Slimane, “Bounds on the distribution of a sum of independent lognormal random variables,” IEEE Transactions on Communications, vol. 49, no. 6, pp. 975–978, 2001.

N. C. Beaulieu and Q. Xie, “An optimal lognormal approximation to lognormal sum distributions,” IEEE Transactions on Vehicular Technology, vol. 53, no. 2, pp. 479–489, 2004.

J. C. S. Santos Filho, M. D. Yacoub, and P. Cardieri, “Highly accurate range-adaptive lognormal approximation to lognormal sum distributions,” IEEE Electronics Letters, vol. 42, no. 6, pp. 361–363, 2006.

Z. Liu, J. Almhana, F. Wang, and R. McGorman, “Mixture lognormal approximations to lognormal sum distributions,” IEEE Communications Letters, vol. 11, no. 9, pp. 711–713, 2007.

N. C. Beaulieu and F. Rajwani, “Highly accurate simple closed-form approximations to lognormal sum distributions and densities,” IEEE Communications Letters, vol. 8, no. 12, pp. 709–711, 2004.

L. Zhao and J. Ding, “A strict approach to approximating lognormal sum distributions,” in Proceedings of the 2006 Canadian Conference on Electrical and Computer Engineering, pp. 916–919, Ottawa, ON, Canada, May 2006.

S. C. Schwartz and Y. S. Yeh, “On the Distribution Function and Moments of Power Sums With Log‐Normal Components,” Bell System Technical Journal, vol. 61, no. 7, pp. 1441–1462, 1982.

C. F. Lo, “The sum and difference of two lognormal random variables,” Journal of Applied Mathematics, Art. ID 838397, 13 pages, 2012.

D. Senaratne and C. Tellambura, “Numerical computation of the lognormal sum distribution,” in Proceedings of the 2009 IEEE Global Telecommunications Conference, GLOBECOM 2009, December 2009.

M. Di Renzo, L. Imbriglio, F. Graziosi, and F. Santucci, “Smolyak's algorithm: A simple and accurate framework for the analysis of correlated log-normal power-sums,” IEEE Communications Letters, vol. 13, no. 9, pp. 673–675, 2009.

F. Black and M. Scholes, “The pricing of options corporate liabilities,” Journal of Political Economy, vol. 81, pp. 637–659, 1973.

J.-P. Fouque, G. Papanicolaou, and K. Ronnie Sircar, Derivatives in Financial Markets with Stochastic Volatility, Cambridge University Press, 2000.

R. F. Engle, “Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation,” Econometrica, vol. 50, no. 4, pp. 987–1007, 1982.

T. Bollerslev, “Generalized autoregressive conditional heteroskedasticity,” Journal of Econometrics, vol. 31, no. 3, pp. 307–327, 1986.

S. L. Heston, “A closed-form solution for options with stochastic volatility with applications to bond and currency options,” The Review of Financial Studies, vol. 6, no. 2, pp. 327–343, 1993.

J.-H. Kim, Y.-K. Ma, and C. Y. Park, “Joint survival probability via truncated invariant copula,” Chaos, Solitons and Fractals, vol. 85, pp. 68–76, 2016.

Y.-K. Ma, “Modeling the dependency structure of integrated intensity processes,” PLoS ONE, vol. 10, no. 8, 10 pages, 2015.