Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Evaluating the Maximal Violation of the Original Bell Inequality by Two-Qudit States Exhibiting Perfect Correlations/Anticorrelations
Previous Article in Journal
Nonlocal Means Two Dimensional Histogram-Based Image Segmentation via Minimizing Relative Entropy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geometric Average Asian Option Pricing with Paying Dividend Yield under Non-Extensive Statistical Mechanics for Time-Varying Model †

College of Mathematics and Information Science, Henan Normal University, Xinxiang 453007, Henan Province, China
*
Author to whom correspondence should be addressed.
This work is supported by the National Natural Science Foundation(NNSF) of China (No. U1504701), Doctoral research startup project of Henan Normal University (No. qd15184).
Entropy 2018, 20(11), 828; https://doi.org/10.3390/e20110828
Submission received: 23 September 2018 / Revised: 20 October 2018 / Accepted: 24 October 2018 / Published: 28 October 2018
(This article belongs to the Section Statistical Physics)

Abstract

:
This paper is dedicated to the study of the geometric average Asian call option pricing under non-extensive statistical mechanics for a time-varying coefficient diffusion model. We employed the non-extensive Tsallis entropy distribution, which can describe the leptokurtosis and fat-tail characteristics of returns, to model the motion of the underlying asset price. Considering that economic variables change over time, we allowed the drift and diffusion terms in our model to be time-varying functions. We used the I t o ^ formula, Feynman–Kac formula, and P a d e ´ ansatz to obtain a closed-form solution of geometric average Asian option pricing with a paying dividend yield for a time-varying model. Moreover, the simulation study shows that the results obtained by our method fit the simulation data better than that of Zhao et al. From the analysis of real data, we identify the best value for q which can fit the real stock data, and the result shows that investors underestimate the risk using the Black–Scholes model compared to our model.

1. Introduction

In financial markets, the movement of asset price is the foundation of the pricing of financial assets and derivatives. Many researchers are interested in option valuation, where Black and Scholes [1] carried out milestone work. In 1973, Black and Scholes [1] constructed the Black–Scholes model to price the option, and they assumed that the empirical distributions of returns are log-normal distributions. Thereafter, many researchers continued to study option pricing based on the Black–Scholes model. For example, Wang [2] studied the Black–Scholes stock option pricing model based on dynamic investment strategy, deriving new option pricing models based on the Black–Scholes option pricing theory. Glazyrina [3] showed how the normal approximation of the binomial distribution leads to an alternative derivation of the Black–Scholes formula from a binomial option pricing model. Ulyah et al. [4] considered short-dated foreign equity options and proposed a new model based on the Black–Scholes theory for their pricing.
In the above studies, the financial markets were Brownian motion-driven models. The assumption that the price of a risky asset follows geometric Brownian motion implies that the price changes are independent. However, many papers have shown that the distribution of empirical returns does not follow log-normal distribution and has the characteristic of fat tails [5,6,7]. In subsequent papers, a number of researchers discovered that the implied volatility calculated from the Black–Scholes model shows a volatility smile. Therefore, some scholars modified the Black–Scholes model to correct for the volatility smile. For example, Hull and White [8] introduced a stochastic volatility model to correct for the smile effect. Merton [9] presented jumps to characterize the intermittent fluctuations in price changes, namely, the jump-diffusion model. Hubalek et al. [10] presented some results on geometric Asian option valuation for affine stochastic volatility models with jumps. Necula [11], Xiao [12], and Gu [13] used fractional Brownian motion with self-similarity and long-term correlation to study option pricing. However, these approaches are very complicated and cannot achieve manageable closed-form solutions.
In 1988, the Brazilian physicist Tsallis [14] proposed the non-extensive Tsallis entropy theory. Tsallis theory regards the financial market as a complex system, defines the price process of assets as an abnormal diffusion process, and obtains a simple form of a distribution function that can describe complex systems with nonlinear, long-range interactions and long-term memory effects. Recently, Tsallis theory has been widely used in the financial field. Farmer and Geanakoplos [15] pointed out that there are complex systems with asymptotic power-law behaviors in finance and economics, indicating that non-extensive statistics can be used in finance. Borland [16] used the Tsallis entropy distribution to study the European pricing issue. Ferrari [17] applied the Tsallis entropy theory to the research of value at risk, return, volatility, and asset portfolios. Based on the data from the Chinese financial market, Li [18] analyzed Tsallis entropy of the financial market. Wang et al. [19] obtained the pricing formulas of power European options based on Tsallis entropy distribution. Devi [20] obtained a good fit to a Tsallis q-Gaussian distribution for the distributions of all the returns using the method of maximum likelihood estimation. Borland [21] compared the Tsallis distribution model with the L e ´ v y model. Sosa-Correa et al. [22] found that the Tsallis distribution model performs better than the Black–Scholes model in about one-third of the option chains analyzed in the Brazilian option market.
The Asian option is one of the most active exotic options in the financial derivatives market today. In 1987, it was first introduced by the Bankers Trust in Tokyo, Japan. The difference from the standard option is that, when determining the option income on the maturity date, instead of using the current market price of the underlying asset, the Asian option uses the average of the asset price over a certain period of time during the option contract period. As Asian options are widely traded, they have received considerable attention in the financial literature. For example, Kemna and Vorst [23] proposed an expression for the geometric Asian option. Rogers and Shi [24] approached the problem of computing the price of an Asian option with the finite difference method. Benhamou and Duguet [25] presented an efficient method for pricing discrete Asian options in the presence of smile and non-proportional dividends. For geometric Asian options, Fusai and Meucci [26] studied the pricing of Asian options under the L e ´ v y process, and they provided closed-form solutions in terms of the Fourier transform. Kirkby [27] developed a method for efficiently inverting analytic characteristic functions using frame projection for geometric Asian pricing under general Levy models. Cai and Kou [28] obtained a closed-form solution for the double-Laplace transform of Asian options under the hyper-exponential jump-diffusion model. Cui et al. [29] proposed a general framework for the valuation of options in stochastic local volatility (SLV) models with a general correlation structure, and provided single Laplace transform formulae for arithmetic Asian options, as well as occupation time derivatives.
Considering that the return distribution of the underlying stock has a peak and fat tails in actual financial markets, in this study, we used the non-extensive Tsallis entropy distribution with long-term interaction and historical memory characteristics to replace the normal distribution without historical memory and modeled the motion of the underlying asset price. This model can depict the leptokurtosis and fat-tail characteristics of the distribution of returns. Furthermore, we used the I t o ^ formula, Feynman–Kac formula, and P a d e ´ ansatz to obtain a closed-form solution of geometric average Asian option pricing. Compared to Zhao’s paper [30], our method can better fit the simulation data.
The rest of this paper is organized as follows. In Section 2, we introduce the Tsallis distribution to model the price of a risky asset. In Section 3, we investigate the geometric average Asian call option and derive the pricing formula. In Section 4, we provide several simulation studies. In Section 5, we select the daily returns of a stock and perform an analysis of real data. In Section 6, we summarize our paper.
Our work differs from the existing approach of Zhao et al. [30] in the following important respects:
  • In our model, we allow the coefficients of the drift term and diffusion term to be time-varying functions. In some sense, we extend the model of Zhao et al. [30].
  • In this paper, we propose the geometric average Asian option pricing with paying dividend yield. So, our results are more widely applicable than that of Zhao et al. [30].
  • When we derive the price formula of geometric average Asian option, the most important part is calculating the integral which is defined by Equation (13). Zhao et al. [30] got their results by using Lemma 1 in their paper. We use the Feynman–Kac formula and P a d e ´ ansatz to deal with the integral. A simulation study shows that our method can better fit the simulation data than that of Zhao et al. [30].

2. Model Setting

In this section, we suppose that there are stocks and bonds in the continuous financial market. W ( t ) is the bond price which satisfies the equation below:
d W ( t ) = r ( t ) W ( t ) d t , W ( 0 ) = 1 ,
where r ( t ) is the risk-free interest rate.
The stock price X ( t ) follows the time-varying coefficient stochastic process,
d X ( t ) = μ ( t ) X ( t ) d t + σ ( t ) X ( t ) d Ω ( t ) ,
where X ( 0 ) = X 0 , and μ ( t ) and σ ( t ) denote the drift term and the diffusion term, which both depend on time t, respectively. Also, μ ( t ) = r ( t ) d ( t ) , where d ( t ) is dividend yield.
For the probability space ( Ω , F , P ) , we assume that Ω ( · ) follows the statistical feedback process
d Ω ( t ) = P ( Ω ( t ) ) 1 q 2 d B ( t ) ,
where B ( t ) is the standard Brownian motion. The probability distribution of Ω ( t ) derives from the nonlinear Fokker–Planck equation,
t P ( Ω ( t ) , t ) = 1 2 2 Ω ( t ) 2 P 2 q ( Ω ( t ) , t ) .
The probability density function P ( Ω ( t ) , t ) is called the Tsallis distribution (see [30]), which is given by
P ( Ω ( t ) , t ) = 1 Z ( t ) 1 β ( t ) ( 1 q ) Ω ( t ) 2 1 1 q ,
with
β ( t ) = k 1 q 3 q ( 2 q ) ( 3 q ) t 2 3 q , Z ( t ) = + 1 ( 1 q ) β ( t ) Ω ( t ) 2 1 1 q d Ω ( t ) = ( 2 q ) ( 3 q ) k t 1 3 q , k = π q 1 Γ 2 1 q 1 1 2 Γ 2 1 q 1 ,
where Γ ( · ) is the gamma function.
It is easy to calculate that its mean is zero, and its variance is
V a r ( Ω ( t ) ) = E [ Ω ( t ) 2 ] = 1 ( 5 3 q ) β ( t ) .
According to Equation (4), the variance diverges for q > 5 3 . Thus, we only consider 5 3 > q > 1 , for which the variance is limited.
Then, we define an equivalent martingale measure Q . Let
θ ( t ) = μ ( t ) r ( t ) σ ( t ) P 1 q 2 , B ˜ ( t ) = B ( t ) + 0 t θ ( s ) d s ,
where P 1 q 2 = P ( Ω ( t ) ) 1 q 2 and E [ exp ( 1 2 0 t θ 2 ( s ) d s ) ] < (see [31]). We construct the Radon–Nikodým derivative between the equivalent martingale measures Q and P as follows
d Q d P = exp 0 t θ ( s ) d B ( s ) 1 2 0 t θ 2 ( s ) d s .
Under the probability measure Q , B ˜ ( t ) is the standard Brownian motion, which can be derived from Girsanov theorem.
In order to calculate the stock price, we need to prove that the discounted stock price process
X ( t ) = e 0 t r ( s ) d s X ( t )
is a martingale under the probability measure Q .
In fact, applying Equation (5) to Equation (2), we obtain
d X ( t ) = μ ( t ) X ( t ) d t + σ ( t ) X ( t ) P 1 q 2 d B ( t ) = μ ( t ) X ( t ) d t + σ ( t ) X ( t ) P 1 q 2 d B ˜ ( t ) θ ( t ) d t = σ ( t ) X ( t ) P 1 q 2 d B ˜ ( t ) + r ( t ) X ( t ) d t .
Using the I t o ^ formula on Equation (6), we get
d X ( t ) = e 0 t r ( s ) d s d X ( t ) + r ( t ) e 0 t r ( s ) d s X ( t ) d t .
Substituting Equation (7),
d X ( t ) = e 0 t r ( s ) d s σ ( t ) X ( t ) P 1 q 2 d B ˜ ( t ) + r ( t ) X ( t ) d t       + r ( t ) e 0 t r ( s ) d s X ( t ) d t = e 0 t r ( s ) d s σ ( t ) X ( t ) P 1 q 2 d B ˜ ( t ) = X ( t ) σ ( t ) X ( t ) P 1 q 2 d B ˜ ( t ) .
Based on the above discussion, X ( t ) is a martingale under the probability measure Q .
Theorem 1 below gives the solution of Equation (2).
Theorem 1.
Under the measure Q , we obtain the approximate value of the solution of Equation (2) as below
X ( t ) = X ( 0 ) exp [ 0 t r ( s ) d s + σ ( t ) Ω ( t ) 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s + 1 2 ( 1 q ) g 0 ( t ) + g 2 ( t ) Ω ( t ) 2 ] .
where
g 0 ( t ) = 0 t Z ( s ) q 1 g 2 ( s ) d s g 2 ( t ) = t 5 q 9 ( 2 q ) ( 3 q ) ( 2 q ) ( 3 q ) 0 t σ ( s ) 2 s 1 5 q 9 ( 2 q ) ( 3 q ) d s .
Proof. 
According to Equation (2), we can get
d X ( t ) = μ ( t ) d t + σ ( t ) P 1 q 2 d B ( t ) X ( t ) .
Applying the It o ^ formula to ln X ( t ) , we obtain
d ln X ( t ) = 1 X ( t ) d X ( t ) 1 2 σ 2 ( t ) P 1 q d t .
Substituting Equation (7), we have
d ln X ( t ) = σ ( t ) P 1 q 2 d B ˜ ( t ) + r ( t ) 1 2 σ 2 ( t ) P 1 q d t .
Integrating both sides of Equation (12) simultaneously in [ 0 , t ] ,
ln X ( t ) ln X ( 0 ) = 0 t σ ( s ) P 1 q 2 d B ˜ ( s ) + 0 t r ( s ) 1 2 σ 2 ( s ) P 1 q d s .
Consequently, we infer that
X ( t ) = X ( 0 ) exp 0 t σ ( s ) P 1 q 2 d B ˜ ( s ) + 0 t r ( s ) 1 2 σ 2 ( s ) P 1 q d s = X ( 0 ) exp [ σ ( t ) Ω ( t ) 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 1 β ( s ) ( 1 q ) Ω ( s ) 2 d s + 0 t r ( s ) d s ] = X ( 0 ) exp [ 0 t r ( s ) d s + σ ( t ) Ω ( t ) 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s + 1 2 ( 1 q ) 0 t σ ( s ) 2 Z ( s ) q 1 β ( s ) Ω ( s ) 2 d s ] = X ( 0 ) exp [ 0 t r ( s ) d s + σ ( t ) Ω ( t ) 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s + 1 2 ( 1 q ) I ( t ) ] .
where
I ( t ) = 0 t σ ( s ) 2 Z ( s ) q 1 β ( s ) Ω ( s ) 2 d s .
It is hard to calculate the integral I ( t ) , which contains the term of the type 0 t F ( Ω ( s ) ) d s . Thus, we use the Feynman–Kac formula and the P a d e ´ ansatz, which are valid according to [31], to evaluate I ( t ) .
First, we define U ( Ω ( T ) , T ) = i Ω ( T ) ω i exp { ε 0 T F ( Ω ( t ) ) d t } , where ω i is the weight of a path i ending at Ω ( T ) at time T, and ε is a small parameter. Then, we expand the polynomial to the second order in Ω ,
U ( Ω ( T ) , T ) = i Ω ( T ) ω i exp ε 0 T j = 1 2 f j ( t ) Ω ( t ) j d t .
We establish a generalized Feynman–Kac equation as follows:
U T = 1 2 2 U 0 1 q U Ω 2 + ε j f j ( T ) Ω ( T ) j .
When we insert the ansatz U = U 0 [ 1 + ε ( g 0 ( T ) + g 1 ( T ) Ω ( T ) + g 2 ( T ) Ω ( T ) 2 ) ] into Equation (15), the coefficients g i must satisfy
d g 0 d t = Z ( t ) q 1 g 2 ( t ) d g 1 d t = 2 ( q 2 ) Z ( t ) q 1 β ( t ) g 1 ( t ) + h 1 ( t ) d g 2 d t = ( 5 q 9 ) Z ( t ) q 1 β ( t ) g 2 ( t ) + h 2 ( t ) .
For I ( t ) , we have h 1 = 0 , h 2 = σ ( t ) 2 Z ( t ) q 1 β ( t ) . The coefficients g i are calculated from Equation (16):
g 0 ( t ) = 0 t Z ( s ) q 1 g 2 ( s ) d s g 2 ( t ) = t 5 q 9 ( 2 q ) ( 3 q ) ( 2 q ) ( 3 q ) 0 t σ ( s ) 2 s 1 5 q 9 ( 2 q ) ( 3 q ) d s .
Next, we can get I ( t ) = g 0 ( t ) + g 2 ( t ) Ω ( t ) 2 . Finally, we substitute the terms of X ( t ) for I ( t ) ,
X ( t ) = X ( 0 ) exp [ 0 t r ( s ) d s + σ ( t ) Ω ( t ) 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s + 1 2 ( 1 q ) g 0 ( t ) + g 2 ( t ) Ω ( t ) 2 ] .
This completes the proof of Theorem 1. □

3. Geometric Average Asian Option Pricing Formula

In this section, we discuss the geometric average Asian call option pricing. Suppose that the maturity time is T, the strike price is K, and J T = exp 1 T 0 T ln X ( t ) d t , which is the average price. It is well known that a geometric average Asian call option price can be written as below (see [32]):
C = e 0 T r ( t ) d t E [ ( J T K ) + ] .
Theorem 2.
Given a geometric average Asian call option, it has the following payoff at the maturity time T:
( J T K ) +
where K is the strike price, and J T = exp 1 T 0 T ln X ( t ) d t . Then, the approximation of the option price in the risk-neutral world is given as below:
C = e 0 T ( r ( t ) d ( t ) ) d t X ( 0 ) ω 1 ω 2 e c T + g 3 ( t ) T Ω ( t ) + 1 q 2 T g 5 ( t ) Ω ( t ) 2 P ( Ω , T ) d Ω e 0 T r ( t ) d t K ω 1 ω 2 P ( Ω , T ) d Ω ,
where
c = 0 T 0 t r ( s ) d s d t 0 T 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s d t + 1 2 ( 1 q ) g 4 ( t ) + 0 T 1 2 ( 1 q ) g 0 ( t ) d t g 0 = 0 t Z ( s ) q 1 g 2 ( s ) d s g 2 = t 5 q 9 ( 2 q ) ( 3 q ) ( 2 q ) ( 3 q ) 0 t σ ( s ) 2 s 1 5 q 9 ( 2 q ) ( 3 q ) d s g 3 = t 2 3 q 0 t σ ( s ) s 2 3 q d s g 4 = 1 ( 2 q ) ( 3 q ) 0 t Z ( s ) q 1 g 2 ( s ) d s g 5 = t 5 q 9 ( 2 q ) ( 3 q ) ( 2 q ) ( 3 q ) 0 t 0 u σ ( s ) 2 s 1 5 q 9 ( 2 q ) ( 3 q ) d s d u ω 1 = g 3 ( t ) Δ ( 1 q ) g 5 ( t ) ω 2 = g 3 ( t ) + Δ ( 1 q ) g 5 ( t ) Δ = g 3 ( t ) 2 4 1 2 ( 1 q ) g 5 ( t ) c + T ln X ( 0 ) K .
Proof. 
According to Equation (19), we can obtain
C = e 0 T r ( t ) d t E [ J T ] { J T > K } e 0 T r ( t ) d t E [ K ] { J T > K } = M N ,
where
M : = e 0 T r ( t ) d t E [ J T ] { J T > K } N : = e 0 T r ( t ) d t E [ K ] { J T > K } .
Then, we should calculate the domain of inequality { J T > K } , which is equal to exp 1 T 0 T ln X ( t ) d t > K . It follows that
0 T ln X ( t ) d t > T ln K .
Using Theorem 1, we plug X ( t ) back into the above inequality,
T ln X ( 0 ) + 0 T 0 t r ( s ) d s d t + 0 T σ ( t ) Ω ( t ) d t + 0 T 1 2 ( 1 q ) g 0 ( t ) d t 0 T 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s d t + 0 T 1 2 ( 1 q ) g 2 ( t ) Ω ( t ) 2 d t > T ln K .
After the computation, we get
T ln X ( 0 ) + 0 T 0 t r ( s ) d s d t 0 T 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s d t T ln K + 0 T 1 2 ( 1 q ) g 0 ( t ) d t + 0 T σ ( t ) Ω ( t ) d t + 1 2 ( 1 q ) 0 T g 2 ( t ) Ω ( t ) 2 d t > 0 .
For the convenience of calculation, we denote the following:
A = 0 T σ ( t ) Ω ( t ) d t B = 0 T g 2 ( t ) Ω ( t ) 2 d t .
We use the Feynman–Kac approach and the P a d e ´ ansatz to evaluate A and B, for A and B both have the terms of 0 t F ( Ω ( s ) ) d s .
For A, the coefficients g i are calculated from Equation (16) with h 1 = σ ( s ) :
g 1 ( t ) = t 2 3 q 0 t σ ( s ) s 2 3 q d s .
To avoid confusion, we let g 3 = g 1 and get
A = g 3 ( t ) Ω ( t ) .
For B, we have h 0 = 0 , h 2 = t 5 q 9 ( 2 q ) ( 3 q ) ( 2 q ) ( 3 q ) 0 t σ ( s ) 2 s 1 5 q 9 ( 2 q ) ( 3 q ) d s . Then, we can calculate the coefficients g i from Equation (16), as below:
g 0 ( t ) = 1 ( 2 q ) ( 3 q ) 0 t Z ( s ) q 1 g 2 ( s ) d s g 2 ( t ) = t 5 q 9 ( 2 q ) ( 3 q ) ( 2 q ) ( 3 q ) 0 t 0 u σ ( s ) 2 s 1 5 q 9 ( 2 q ) ( 3 q ) d s d u .
In order to make the expression clearer, we allow g 4 = g 0 , g 5 = g 2 and get
B = g 4 ( t ) + g 5 ( t ) Ω ( t ) 2 .
According to the above approximations, we replace the terms of Equation (25) with A and B,
1 2 ( 1 q ) g 5 ( t ) Ω ( t ) 2 + g 3 ( t ) Ω ( t ) + c + T ln X ( 0 ) K > 0 ,
where
c = 0 T 0 t r ( s ) d s d t 0 T 0 t 1 2 σ ( s ) 2 Z ( s ) q 1 d s d t + 1 2 ( 1 q ) g 4 ( t ) + 0 T 1 2 ( 1 q ) g 0 ( t ) d t .
Because 5 3 > q > 1 , it is clear that 1 2 ( 1 q ) g 5 ( t ) < 0 .
Let K < exp ( ln X ( 0 ) 1 T ( g 3 ( t ) 2 2 ( 1 q ) g 5 ( t ) c ) ) ; then, we can obtain
Δ = g 3 ( t ) 2 4 1 2 ( 1 q ) g 5 ( t ) c + T ln X ( 0 ) K > 0 .
Therefore, the quadratic inequality’s corresponding quadratic equation has two roots,
ω 1 = g 3 ( t ) Δ ( 1 q ) g 5 ( t ) ,
and
ω 2 = g 3 ( t ) + Δ ( 1 q ) g 5 ( t ) .
Hence, we get the solution set of { J T > K } as ω ( ω 1 , ω 2 ) . Then, we can calculate M and N in the risk-neutral world as follows
M = e 0 T r ( t ) d t E [ J T ] { J T > K } ω 1 ω 2 d Ω = X ( 0 ) e 0 T ( r ( t ) d ( t ) ) d t ω 1 ω 2 e c T + g 3 ( t ) T Ω ( t ) + 1 q 2 T g 5 ( t ) Ω ( t ) 2 P ( Ω , T ) d Ω .
Similar to M, we obtain
N = e 0 T r ( t ) d t E [ K ] { J T > K } ω 1 ω 2 d Ω = e 0 T r ( t ) d t K ω 1 ω 2 P ( Ω , T ) d Ω .
Hence, we get the option price as follows:
C = M N = e 0 T ( r ( t ) d ( t ) ) d t X ( 0 ) ω 1 ω 2 e c T + g 3 ( t ) T Ω ( t ) + 1 q 2 T g 5 ( t ) Ω ( t ) 2 P ( Ω , T ) d Ω e 0 T r ( t ) d t K ω 1 ω 2 P ( Ω , T ) d Ω .
This completes the proof of Theorem 2. □

4. Simulation Study

In this section, the results of a simulation study are presented to show the difference between our method and the method of Zhao [30]. In order to compare with Zhao, we assume that r = 0.5 , σ = 0.25 . Without loss of generality, we can take μ = 0 .
When we calculate the geometric average Asian option price, it is hard to solve the terms of 0 T β ( t ) Z ( t ) q 1 Ω ( t ) 2 d t since the Ω ( t ) in the integral function is a stochastic process. We used the Feynman–Kac formula and P a d e ´ ansatz to solve them, while Zhao et al. [30] introduced the transformation Ω ( t ) = β ( T ) β ( t ) Ω ( T ) . In order to compare the two methods, we paint a picture of 0 T β ( t ) Z ( t ) ( q 1 ) Ω ( t ) 2 d t versus Ω ( T ) in the following three ways:
(1)
Monte-Carlo simulation;
(2)
Ω ( t ) = β ( T ) β ( t ) Ω ( T ) in Zhao et al. [30];
(3)
Feynman–Kac formula and P a d e ´ ansatz method in our paper.
Assuming T = 0.5 , q = 1.5 , then we get the following figure. It is also discussed in the appendix of [31]. In Figure 1, it is obvious that the black line fits the discrete points better than the red one. In other words, the data obtained by our method are closer to the simulation data of the real market. The black line is extremely close to those points obtained by Monte-Carlo simulation. From Figure 1, the data calculated by the Feynman–Kac formula and P a d e ´ ansatz are more responsive to the real market data than the method obtained by Zhao et al. [30].
Next, we describe some simulation studies to show the difference between our model and the Black–Scholes model. According to Michael and Johnson’s work [33], the data generated by the Tsallis distribution can describe the real market data accurately, so we used the Tsallis distribution to generate the simulation data. Some steps of the numerical simulation are given as follows:
Step 1:
We assume the value of q , K , r , σ , T , etc.
Step 2:
We generate 1000 random numbers of the Tsallis distribution.
Step 3:
According to step 2, we can calculate the price of a risky asset by Equation (8).
Step 4:
We derive the geometric average Asian call option price by Equation (21) and generate the resulting figures.
Before examining the numerical simulation, we first give the pricing formula of the Black–Scholes model, which is
P = X ( 0 ) e 1 2 0 T ( r ( t ) + 1 6 σ ( t ) 2 ) d t N ( d 1 ) K e 0 T r ( t ) d t N ( d 2 )
where
d 1 = ln X ( 0 ) K + 1 2 0 T ( r ( t ) + 1 6 σ ( t ) 2 ) d t 3 T 0 T σ ( t ) 2 d t , d 2 = d 1 3 T 0 T σ ( t ) 2 d t
Let X ( 0 ) = 100 . Without loss of generality, we used d ( t ) = 0 . Figure 2 depicts the difference between our model and the Black–Scholes model.
In Figure 2, the call option prices are becoming lower as the strike price K becomes larger. Moreover, the call option price calculated by our model is lower than that from the Black–Scholes model. This suggests that investors underestimate the risk using the Black–Scholes model. According to Figure 2, we can see that our model is a better fit than the Black–Scholes model.

5. Analysis of Real Data

In this section, we use actual cases to test the model. We selected the daily closing prices of a stock with a code of 601318 (hereinafter called 601318 ) in the Chinese stock market as real data. The stock is issued by the Ping An Insurance (Group) Company of China. The time period is from 1 March 2007 to 1 March 2018 and the sample size is 2612.
In Table 1, “J-B” is the value of the Jarque–Bera test, which is a test for normality comparing sample skewness and kurtosis against the theoretical values for a normal distribution. “P” is the p-value associated with this test. From the basic statistical values in Table 1, it can be seen that the daily returns of “601318” have obvious characteristics of leptokurtosis and of fat tail. The kurtosis of the empirical data is 200.9483 , and it is well known that the kurtosis of a normal distribution is 3. Comparing these values of kurtosis, 200.9483 is much larger than 3. The estimated kurtosis of 200.9483 demonstrates that the empirical tails are much heavier than those of a normal distribution, which has a kurtosis of 3. Therefore, it is not appropriate to describe real data with a normal distribution. The p-value corresponding the Jarque–Berra test of normality leads us to reject the null hypothesis of normality for this sample.
In Figure 3, we compare the fitting of the empirical distribution of the daily returns for the empirical distribution, normal distribution, and q-Gaussian distribution with q = 1.55 . We can see that the Gaussian distribution cannot describe the peak of the empirical data. According to the calculation results and several experiments, we discover that a Tsallis distribution with the parameter q = 1.55 can fit the empirical density distribution of daily returns more accurately than the normal distribution.
Figure 4 shows the difference in the geometric average Asian call option price between our model and the Black–Scholes model. The data are the stock prices of “601318”. In the graph, we can observe that the price calculated by the Black–Scholes model is higher than that of our model. Theoretically speaking, it is because the tail of the distribution of the stock price is heavy. This suggests that the investors overestimate the option price using the Black–Scholes model.

6. Summary

This paper mainly examines the geometric average Asian call option pricing under the time-varying coefficient diffusion model. The underlying asset price is modeled by using the non-extensive Tsallis entropy distribution. Considering that economic variables change from time to time, we allow both drift and diffusion terms in our model to be time-varying functions. We obtain the closed-form solution of geometric average Asian option pricing with paying dividend yield for the time-varying model by using the I t o ^ formula, Feynman–Kac formula, and P a d e ´ ansatz. Moreover, the simulation studies show that the results obtained by our method fit the real data better than the method of Zhao et al. [30]. Also, investors estimate the risk more reasonably using our model than with the Black–Scholes model. Based on the results obtained through the analysis of real data and reported in our paper, we determine that q = 1.55 is the optimal q value to fit the real stock data.

Author Contributions

All authors contributed equally to the manuscript, read and approved the final manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. U1504701), Doctoral research startup project of Henan Normal University(Grant No. qd15184).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Black, F.; Scholes, M. The pricing of options and corporate liabilities. J. Polit. Econ. 1973, 81, 637–654. [Google Scholar] [CrossRef]
  2. Wang, X.; Wang, L. Study on Black–Scholes stock option pricing model based on dynamic investment strategy. Int. J. Innov. Comput. Inf. Control 2007, 3, 1755–1780. [Google Scholar]
  3. Glazyrina, A.; Melnikov, A. Bernstein’s inequalities and their extensions for getting the Black–Scholes option pricing formula. Stat. Probab. Lett. 2016, 111, 86–92. [Google Scholar] [CrossRef]
  4. Ulyah, S.M.; Lin, X.C.; Miao, D.W. Pricing short-dated foreign equity options with a bivariate jump-diffusion model with correlated fat-tailed jumps. Financ. Res. Lett. 2017, 24, 1–16. [Google Scholar] [CrossRef]
  5. Thomas, L. The stable paretian hypothesis and the frequency of large returns: An examination of major German stocks. Appl. Financ. Econ. 1996, 6, 463–475. [Google Scholar]
  6. Kittiakarasakun, J.; Tse, Y. Modeling the fat tails in asian stock markets. Int. Rev. Econ. Financ. 2011, 20, 430–440. [Google Scholar] [CrossRef]
  7. Bonato, M. Modeling fat tails in stock returns: A multivariate stable-garch approach. Comput. Stat. 2012, 27, 499–521. [Google Scholar] [CrossRef]
  8. Hull, J.C. Options, Futures, and Other Derivatives; Tsinghua University Press: Beijing, China, 2006. [Google Scholar]
  9. Merton, R.C. Option pricing when underlying returns are discontinuous. J. Financ. Econ. 1976, 3, 125–144. [Google Scholar] [CrossRef]
  10. Hubalek, F.; Keller-Ressel, M.; Sgarra, C. Geometric asian option pricing in general affine stochastic volatility models with jumps. Quantit. Financ. 2017, 17, 873–888. [Google Scholar] [CrossRef]
  11. Necula, C. Option Pricing in a Fractional Brownian Motion Environment. 2002. Available online: http://dx.doi.org/10.2139/ssrn.1286833 (accessed on 26 October 2018).[Green Version]
  12. Xiao, W.; Zhang, W.; Zhang, X.; Wang, Y.L. Pricing currency options in a fractional Brownian motion with jumps. Econ. Model. 2010, 27, 935–942. [Google Scholar] [CrossRef]
  13. Gu, H.; Liang, J.; Zhang, Y. Time-changed geometric fractional Brownian motion and option pricing with transaction costs. Phys. A Stat. Mech. Appl. 2012, 391, 3971–3977. [Google Scholar] [CrossRef]
  14. Tsallis, C. Possible Generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  15. Farmer, J.D.; Geanakoplosy, J. Power Laws in Economics and Elsewhere; Santa Fe Institute: Santa Fe, NM, USA, 2008. [Google Scholar]
  16. Borland, L. A theory of non-Gaussian option pricing. Quant. Financ. 2002, 2, 415–431. [Google Scholar] [Green Version]
  17. Ferrari, D.; Paterlini, S. Efficient and Robust Estimation for Financial Returns: An Approach Based on Q-Entropy. 2010. Available online: http://dx.doi.org/10.2139/ssrn.1906819 (accessed on 26 October 2018).
  18. Li, S.; He, J.; Song, K. Network entropies of the chinese financial market. Entropy 2016, 18, 331. [Google Scholar] [CrossRef]
  19. Wang, Y.; Li, D.; Wei, J. Pricing of Power European Options Based on Tsallis Entropy and O-U Process under Stochastic Interest Rate. J. Zhengzhou Univ. 2017, 49, 1–4. [Google Scholar] [CrossRef]
  20. Devi, S. Financial market dynamics: Superdiffusive or not? J. Stat. Mech. Theory Exp. 2017, 2017, 083207. [Google Scholar] [CrossRef]
  21. Borland, L. Exploring the dynamics of financial markets: From stock prices to strategy returns. Chaos Solitons Fractals 2016, 88, 59–74. [Google Scholar] [CrossRef]
  22. Sosa-Correa, W.O.; Ramos, A.M.T.; Vasconcelos, G.L. Investigation of non-Gaussian effects in the Brazilian option market. Phys. A Stat. Mech. Appl. 2018, 496, 525–539. [Google Scholar] [CrossRef]
  23. Kemna, A.G.Z.; Vorst, A.C.F. A pricing method for options based on average asset values. J. Bank. Financ. 2014, 14, 113–129. [Google Scholar] [CrossRef]
  24. Rogers, L.C.G.; Shi, Z. The value of an Asian option. J. Appl. Probab. 1995, 32, 1077–1088. [Google Scholar] [CrossRef]
  25. Benhamou, E.; Duguet, A. Small dimension pde for discrete asian options. J. Econ. Dyn. Control 2010, 27, 2095–2114. [Google Scholar] [CrossRef]
  26. Fusai, G.; Meucci, A. Pricing discretely monitored Asian options under lvy processes. J. Bank. Financ. 2008, 32, 2076–2088. [Google Scholar] [CrossRef]
  27. Kirkby, J.L. Efficient option pricing by frame duality with the fast Fourier transform. Soc. Sci. Electron. Publ. 2014, 6, 713–747. [Google Scholar] [CrossRef]
  28. Cai, N.; Kou, S. Pricing Asian Options Under a Hyper-Exponential Jump Diffusion Model. Oper. Res. 2012, 60, 64–77. [Google Scholar] [CrossRef] [Green Version]
  29. Cui, Z.; Kirkby, J.L.; Nguyen, D. A general valuation framework for SABR and stochastic local volatility models. SIAM J. Financ. Math. 2018, 9, 520–563. [Google Scholar] [CrossRef]
  30. Zhao, P.; Zhou, B.; Wang, J. Non-Gaussian closed form solutions for geometric average Asian options in the framework of non-extensive statistical mechanics. Entropy 2018, 20, 71. [Google Scholar] [CrossRef]
  31. Borland, L.; Bouchaud, J. A non-Gaussian option pricing model with skew. Quant. Financ. 2004, 4, 499–514. [Google Scholar] [CrossRef] [Green Version]
  32. Devreesea, J.P.A.; Lemmensa, D.; Temperea, J. Path integral approach to Asian options in the Black–Scholes model. Phys. A Stat. Mech. Appl. 2011, 389, 780–788. [Google Scholar] [CrossRef]
  33. Michael, F.; Johnson, M.D. Financial market dynamics. Phys. A Stat. Mech. Appl. 2003, 320, 525–534. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The x-axis is Ω ( T ) , and the y-axis is 0 T β ( t ) Z ( t ) ( q 1 ) Ω ( t ) 2 d t . The points were obtained by Monte-Carlo simulation, the red line is calculated in Zhao’s paper, and the black line is obtained by the Feynman–Kac formula and P a d e ´ ansatz in our paper.
Figure 1. The x-axis is Ω ( T ) , and the y-axis is 0 T β ( t ) Z ( t ) ( q 1 ) Ω ( t ) 2 d t . The points were obtained by Monte-Carlo simulation, the red line is calculated in Zhao’s paper, and the black line is obtained by the Feynman–Kac formula and P a d e ´ ansatz in our paper.
Entropy 20 00828 g001
Figure 2. The x-axis is the strike price K, and the y-axis is the geometric average Asian call option price C.
Figure 2. The x-axis is the strike price K, and the y-axis is the geometric average Asian call option price C.
Entropy 20 00828 g002
Figure 3. Comparison of the fitting of the empirical distribution of the daily returns for the empirical distribution, normal distribution, and q-Gaussian distribution.
Figure 3. Comparison of the fitting of the empirical distribution of the daily returns for the empirical distribution, normal distribution, and q-Gaussian distribution.
Entropy 20 00828 g003
Figure 4. Comparison of the geometric average Asian call option price for our model and the Black–Scholes model.
Figure 4. Comparison of the geometric average Asian call option price for our model and the Black–Scholes model.
Entropy 20 00828 g004
Table 1. The basic statistics of daily returns of 601318 .
Table 1. The basic statistics of daily returns of 601318 .
Sample SizeMeanStdKurtosisJ–BP
26121.4838 × 10 4 8.7055 × 10 4 200.94834.2867 × 10 6 0

Share and Cite

MDPI and ACS Style

Wang, J.; Zhang, Y. Geometric Average Asian Option Pricing with Paying Dividend Yield under Non-Extensive Statistical Mechanics for Time-Varying Model. Entropy 2018, 20, 828. https://doi.org/10.3390/e20110828

AMA Style

Wang J, Zhang Y. Geometric Average Asian Option Pricing with Paying Dividend Yield under Non-Extensive Statistical Mechanics for Time-Varying Model. Entropy. 2018; 20(11):828. https://doi.org/10.3390/e20110828

Chicago/Turabian Style

Wang, Jixia, and Yameng Zhang. 2018. "Geometric Average Asian Option Pricing with Paying Dividend Yield under Non-Extensive Statistical Mechanics for Time-Varying Model" Entropy 20, no. 11: 828. https://doi.org/10.3390/e20110828

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop