Next Article in Journal
Multiscale Approximate Entropy for Gait Analysis in Patients with Neurodegenerative Diseases
Next Article in Special Issue
Correction: Contreras-Reyes, J.E.; Cortés, D.D. Bounds on Rényi and Shannon Entropies for Finite Mixtures of Multivariate Skew-Normal Distributions: Application to Swordfish (Xiphias gladius Linnaeus). Entropy 2016, 18, 382
Previous Article in Journal
Segmented Thermoelectric Generator under Variable Pulsed Heat Input Power
Previous Article in Special Issue
Asymptotic Behavior of Memristive Circuits
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

European Option Based on Least-Squares Method under Non-Extensive Statistical Mechanics

College of Mathematics and Information Science, Henan Normal University, Xinxiang 453007, China
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(10), 933; https://doi.org/10.3390/e21100933
Submission received: 18 August 2019 / Revised: 15 September 2019 / Accepted: 23 September 2019 / Published: 25 September 2019
(This article belongs to the Collection Advances in Applied Statistical Mechanics)

Abstract

:
This paper is devoted to the study of the pricing of European options under a non-Gaussian model. This model follows a non-extensive statistical mechanics which can better describe the fractal characteristics of price movement in the financial market. Moreover, we present a simple but precise least-square method for approximation and obtain a closed-form solution of the price of European options. The advantages of this technique are illustrated by numerical simulation, which shows that the least-squares method is better compared with Borland’s two methods in 2002 and 2004.

1. Introduction

Since the famous Black–Scholes [1] model was put forward in 1973, many scholars have begun to deal with financial investment problems in terms of the Black–Scholes model. Wang [2] studied the Black–Scholes option pricing model based on the dynamic investment strategy. Glazyrina [3] showed another derivation of how the normal approximation of the binomial distribution leads to the Black–Scholes formula from the binomial option pricing model. Ulyah et al., [4] proposed a new pricing model for short-term foreign stock options by Black–Scholes theory.
It is well-known that the Black–Scholes model assumed that stock prices follow a logarithmic normal distribution and have independent increments. However, Lo [5] and Lux [6] have confirmed that the returns of financial assets have some fractal properties, such as the sharp peak, fat tail, and long-term memory. Therefore, the hypothesis in Black–Scholes model is inconsistent with the empirical results. In addition, the underestimation of the options based on Black–Scholes model results in the volatility smile curve. Hence, many scholars have begun to correct the Black–Scholes model. Merton [7] proposed a jump model to describe the intermittent fluctuation of price. Hull [8] introduced the stochastic volatility model. Necula [9] established a new price model by using fractional Brownian motion. Since these models can well depict the characteristics of long memory or sharp peak and fat tail of the return distribution, Hubalek et al., [10], Xiao [11] and Gu [12] popularized the application of these models in the financial field.
In 1988, Tsallis [13] put forward the theory of non-extensive Tsallis entropy. Then Borland [14] studied the fractal characteristics of Tsallis distribution from the perspective of microdynamics. Liu [15] systematically discussed the self-similarity, non-Markovian and long dependence of this distribution. Since the non-extensive Tsallis has the fractal characteristics, it can well model the distributions of asset returns of many financial assets. In 2002, Borland [16] began to apply Tsallis theory to the pricing problem and achieved a quiet better result. Later, some scholars used this theory to asset portfolio, option pricing and other aspects and achieved better results [17,18,19,20,21].
In the process of option pricing, the key issue is to solve the integral 0 T f ( Ω ( t ) , t ) d t . To address this problem, Borland transformed this integral into a function in terms of Ω ( T ) and T by substituting Ω ( t ) with β ( T ) β ( t ) Ω ( T ) for the reason that β ( T ) β ( t ) Ω ( T ) and Ω ( t ) have identical distribution. Zhao [22] extended this method to the average geometric Asian option pricing problem. In 2004, Borland [23] used the Feynman–Kac equation and P a d e ´ ansatz operator to reconsider the integral approximation and got the price formula of European call option. Then Wang [24] extended this method to the average geometric Asian option pricing problem. However, the above scholars did not discuss the effectiveness of the two approximations of Borland. In this paper, we used a new least-square method to reconsider the approximation problem and discussed the effectiveness of our method, by which to promote the wide application of Tsallis theory in the financial field.
The rest of the paper is organized as follows. In Section 2, we establish the asset price model by non-extensive statistics theory. In Section 3, we use the least-square method to solve the integral approximation problem. In Section 4, we deduce the closed solution of European option which is similar to the Black–Scholes formula. In Section 5, the effectiveness of our method is discussed. In Section 6, we summarize our paper.

2. Price Model Based on the Non-Extensive Statistics

This section considers a price model which is same as Borland’s model in [16,23]. The characteristic of this model is that the noise of stock price fluctuation is a non-Gaussian process.
Suppose that there are two kinds of assets in the market. One is the risk-free bond with a price satisfying the equation below:
d B ( t ) = r B ( t ) d t , B ( 0 ) = 1 ,
where r is risk-free rate, B ( t ) is the bond price. The other asset in the market is the stock, its price S ( t ) follows
d S ( t ) = μ S ( t ) d t + σ S ( t ) d Ω ( t ) ,
with
d Ω ( t ) = P ( Ω ( t ) , t ) 1 q 2 d W ( t ) ,
where Ω ( t ) is a non-Gaussian noise source, W ( t ) is a standard Brownian motion defined on the probability space ( Ω , { F t } t 0 , P ) . P ( Ω ( t ) , t ) comes from the non-extensive entropy theory in statistical physics, which refers to the probability distribution when Tsallis entropy is maximized under the constraints of regularization and normalization.
P ( Ω ( t ) , t ) satisfy the following nonlinear Fokker-Planck equation:
P ( Ω ( t ) , t | Ω ( t ) , t ) t = 1 2 2 Ω ( t ) 2 P 2 q ( Ω ( t ) , t | Ω ( t ) , t ) ,
where q is the non-extensive parameter. Plastino and Plastino [25] have proved that the conditional probability in Equation (4) has a form as
P ( Ω ( t ) , t | Ω ( t ) , t ) = 1 Z ( t ) ( 1 β ( t ) ( 1 q ) ( Ω ( t ) Ω ( t ) ) 2 ) 1 1 q ,
with
β ( t ) = c 1 q 3 q [ ( 2 q ) ( 3 q ) ( t t ) ] 2 3 q ,
Z ( t ) = [ ( 2 q ) ( 3 q ) c ( t t ) ] 1 3 q ,
c = π q 1 Γ 2 ( 1 q 1 1 2 ) Γ 2 ( 1 q 1 ) ,
where Γ ( · ) is the Gamma function. In addition, as q 1 , Ω ( t ) becomes a Gaussian process. Let t = 0 , Ω ( 0 ) = 0 , we can get the following Tsallis distribution [26]:
P ( Ω ( t ) , t ) = 1 Z ( t ) [ 1 β ( t ) ( 1 q ) Ω ( t ) 2 ] 1 1 q .
From Equation (6), it is easily to verify that the mean function of Ω ( t ) is zero and the variance function is
V a r [ Ω ( t ) ] = 1 ( 5 3 q ) β ( t ) .
Since we are only interested in the process with finite variance, it is reasonable to assume 1 q < 5 3 .
When there is no arbitrage opportunity in the market, for a risk-neutral investor, the expected return of the asset in the future must be the risk-free rate. In other words, the fair market value should discount accordingly in the risk-neutral framework at the risk-free rate. The following Theorem deduces the formula of stock price based on Tsallis distribution in the fair market, which is the reorganization and restatement of Borland [16].
Theorem 1.
In the risk-neutral market, the stock price S ( t ) satisfying the stochastic differential Equation (2) is
S ( t ) = S ( 0 ) exp { σ Ω ( t ) + r t σ 2 2 δ ( Ω ( t ) , t ) } ,
where S ( 0 ) is the price at time 0, and
δ ( Ω ( t ) , t ) = 0 t 1 Z ( s ) 1 q [ 1 β ( s ) ( 1 q ) Ω ( s ) 2 ] d s .
Proof. 
Define the discounted stock price J ( t ) is
J ( t ) = e r t S ( t ) .
Since there is generally no arbitrage opportunity in a risk-neutral financial market, Equation (10) is required to be a martingale process. According to Girasanov’s theorem, we can find an equivalent measure Q corresponding to another noise term W ˜ ( t ) that transforms Equation (10) into a martingale. Let
d W ˜ ( t ) = κ d t + d W ( t ) ,
and
κ = μ r σ P 1 q 2 ,
where
P 1 q 2 = P ( Ω ( t ) , t ) 1 q 2 .
Then the new measure Q is related to original P by the following Radon–Nikodym derivative
d Q d P | F t = exp 0 t κ d W ( t ) 1 2 0 t κ 2 d s .
Using I t o ^ formula on Equation (10), we obtain
d J ( t ) = ( μ r ) J ( t ) d t + σ P 1 q 2 J ( t ) d W ( t ) .
Inserting (11) into (12) yields
d J ( t ) = ( μ r ) J ( t ) d t + σ P 1 q 2 J ( t ) d W ( t ) = σ J ( t ) P 1 q 2 ( μ r σ P 1 q 2 d t + d W ( t ) ) = σ J ( t ) P 1 q 2 d W ˜ ( t ) .
Thus, J ( t ) becomes a martingale process.
Using I t o ^ formula we get
d ln S ( t ) = r d t 1 2 σ 2 P 1 q d t + σ P 1 q 2 d W ˜ ( t ) ,
which implies that
S ( t ) = S ( 0 ) exp { σ Ω ( t ) + r t σ 2 2 δ ( Ω ( t ) , t ) } ,
where δ ( Ω ( t ) , t ) is defined by (9).  ☐

3. Approximate Solution Based on Least-Square Method

Since integral δ ( Ω ( t ) , t ) has the form of 0 t f ( Ω ( s ) , s ) d s , it cannot be solved directly. In this section, we use the least-square method (LSM) to find a simple regression equation to approximate this integral.
As mentioned in Equation (7), the stochastic process Ω ( t ) has the mean of 0 and the variance of 1 ( 5 3 q ) β ( t ) . Suppose the standardized variable of Ω ( t ) is Ω * expressed as
Ω * = Ω ( t ) 1 ( 5 3 q ) β ( t ) ,
then Ω * follows a standardized Tsallis distribution with the mean of 0 and the variance of 1. Similarly, for any time s, it is easily to obtain
Ω ˜ * = Ω ( s ) 1 ( 5 3 q ) β ( s ) .
This means that Ω * and Ω ˜ * are identically distributed, i.e.
Ω ( s ) = d β ( t ) β ( s ) Ω ( t ) .
It can be clearly seen that if we use the Ω ( t ) to substitute the Ω ( s ) in Equation (9), then the form of δ ( Ω ( t ) , t ) can be transformed more simple so that we can solve. Following this intuition, inserting Equation (15) into Equation (9) yields
δ ( Ω ( t ) , t ) = 1 2 ( 3 q ) ( ( 2 q ) ( 3 q ) ε ) q 1 3 q t 2 3 q ( 1 ( 1 q ) β ( t ) Ω ( t ) 2 ) ,
By observing Equation (16), we find δ is a function about t 2 3 q and Ω ( t ) 2 . Hence, the basis functions are t 2 3 q , Ω ( t ) 2 and t 2 3 q · Ω ( t ) 2 . Suppose the approximate expression of δ ( Ω ( t ) , t ) is
δ ^ ( Ω ( t ) , t ) = d 1 ( q ) · t 2 3 q + d 2 ( q ) · Ω ( t ) 2 + d 3 ( q ) · t 2 3 q Ω ( t ) 2 + d 4 ( q ) .
Definition 1.
The estimated error E between δ ^ and δ is
E = E P ( δ ( Ω ( t ) , t ) δ ^ ( Ω ( t ) , t ) ) 2 ,
where E P [ · ] represents the mean under measure P .
Borland [16] in 2002 mapped Ω ( s ) onto the Ω ( t ) by Equation (15) and derived an approximate expression for δ ( Ω ( t ) , t ) .
Remark 1.
The approximation of δ ( Ω ( t ) , t ) derived by Borland in 2002 is
δ 1 ^ ( Ω ( t ) , t ) = 1 2 ( 3 q ) ( ( 2 q ) ( 3 q ) c ) q 1 3 q t 2 3 q ( 3 q ) ( 1 q ) 2 β ( t ) ( ( 2 q ) ( 3 q ) c ) q 1 3 q t 2 3 q Ω ( t ) 2 .
By comparing with Equation (17), it can be seen that d j in Equation (17) take the following form respectively,
d 1 = 1 2 ( 3 q ) ( ( 2 q ) ( 3 q ) c ) q 1 3 q , d 2 = 0 , d 3 = ( 3 q ) ( 1 q ) 2 β ( t ) ( ( 2 q ) ( 3 q ) c ) q 1 3 q , d 4 = 0 .
However, the disadvantage of the method (15) is that Ω ( t ) is a random process not a simple distribution, the substitution of (15) will cause a large deviation of the approximation of δ ( Ω ( t ) , t ) . Then Borland proposed a new method in 2004 taken in [23], in which he used the Feynman–Kac equation and P a d e ´ ansatz to reconsider the approximate of δ ( Ω ( t ) , t ) and got another approximation result.
Remark 2.
The approximation of δ ( Ω ( t ) , t ) derived by Borland in 2004 is
δ 2 ^ ( Ω ( t ) , t ) = ( ( 2 q ) ( 3 q ) c ) q 1 3 q 3 q 2 ( 1 q ) ( 3 q ) 2 ( 9 5 q ) t 2 3 q 1 q 9 5 q Ω ( t ) 2
Similar comparison with Equation (17), it can be seen that the d j in (17) are
d 1 = ( ( 2 q ) ( 3 q ) c ) q 1 3 q 3 q 2 ( 1 q ) ( 3 q ) 2 ( 9 5 q ) , d 2 = 1 q 9 5 q , d 3 = 0 , d 4 = 0 .

4. European Option Price

In this section, the approximate function δ ^ ( Ω ( t ) , t ) is used to solve the pricing problem of European options.
We first consider the pricing of European call options. A European call option has a payoff of
C T = max [ S ( T ) K , 0 ] ,
where S ( T ) represents the price of a risky asset at maturity T and K is the strike price. Under the measure Q , the price of this option is
C 0 = E Q [ e r T C ] = E Q [ e r T S ( T ) I { S ( T ) > K } ] E Q [ e r T K I { S ( T ) > K } ] = A 1 A 2 ,
with
A 1 = E Q [ e r T S ( T ) I { S ( T ) > K } ] , A 2 = E Q [ e r T K I { S ( T ) > K } ] ,
where E Q [ · ] represents the mean value under measure Q and I { S ( T ) > K } is an indicator function.
Theorem 2.
The price of the European call option is given by
C 0 = S ( 0 ) M q ( γ 1 , γ 2 ) K e r T N q ( γ 1 , γ 2 ) ,
where
M q ( γ 1 , γ 2 ) = e r T γ 1 γ 2 exp { σ Ω ( T ) σ 2 2 δ ^ ( Ω ( T ) , T ) } P ( Ω ( T ) , T ) d Ω ( T ) , N q ( γ 1 , γ 2 ) = γ 1 γ 2 P ( Ω ( T ) , T ) d Ω ( T ) , P ( Ω ( T ) , T ) = 1 Z ( T ) [ 1 β ( T ) ( 1 q ) Ω ( T ) 2 ] 1 1 q , δ ^ ( Ω ( T ) , T ) = d 1 ( q ) · T 2 3 q + d 2 ( q ) · Ω ( T ) 2 + d 3 ( q ) · T 2 3 q Ω ( T ) + d 4 ( q ) , γ 1 , 2 = a 2 ± a 2 2 4 a 1 · a 3 2 a 1 a 2 , a 1 = σ 2 2 ( d 2 ( q ) + d 3 ( q ) T 2 3 q ) , a 2 = σ , a 3 = r T σ 2 2 d 1 ( q ) T 2 3 q σ 2 2 d 4 ( q ) ln K S ( 0 ) .
Proof. 
By Theorem 1, the approximate expression of the stock price is
S ( T ) = S ( 0 ) exp { σ Ω ( T ) + r T σ 2 2 δ ^ ( Ω ( T ) , T ) } = S ( 0 ) exp { σ Ω ( T ) + r T σ 2 2 ( d 1 ( q ) · T 2 3 q + d 2 ( q ) · Ω ( T ) 2 + d 3 ( q ) · T 2 3 q Ω ( T ) + d 4 ( q ) ) } .
To calculate A 1 and A 2 in Equation (22), we should first solve the inequality { S ( T ) > K } . Using Equation (24) yields
S ( 0 ) exp { σ Ω ( T ) + r T σ 2 2 ( d 1 ( q ) · T 2 3 q + d 2 ( q ) · Ω ( T ) 2 + d 3 ( q ) · T 2 3 q Ω ( T ) + d 4 ( q ) ) } > K .
Taking logarithms of both sides of Equation (25) obtain
a 1 Ω ( T ) 2 + a 2 Ω ( T ) + a 3 > 0 ,
where
a 1 = σ 2 2 ( d 2 ( q ) + d 3 ( q ) T 2 3 q ) , a 2 = σ , a 3 = r T σ 2 2 d 1 ( q ) T 2 3 q σ 2 2 d 4 ( q ) ln K S ( 0 ) .
The quadratic equation
a 1 Ω ( T ) 2 + a 2 Ω ( T ) + a 3 = 0
has the discriminant as
= a 2 2 4 a 1 · a 3 = σ 2 + 2 × ( σ 2 ( d 2 ( q ) + d 3 ( q ) T 2 3 q ) ) × ( r T σ 2 2 d 1 ( q ) T 2 3 q σ 2 2 d 4 ( q ) ln K S ( 0 ) ) .
Assuming that > 0 , there are two roots of quadratic Equation (27)
γ 1 , 2 = a 2 ± a 2 2 4 a 1 · a 3 2 a 1 a 2 .
Therefore, we get the solution set of { S ( T ) > K } as Ω ( T ) ( γ 1 , γ 2 ) . Substituting the roots into A 1 and A 2 , we finally get
A 1 = E Q [ e r T S ( T ) · I S ( T ) > K ] = γ 1 γ 2 e r T · S ( T ) · P ( Ω ( T ) , T ) d Ω ( T ) = e r T S ( 0 ) ( γ 1 γ 2 exp { σ Ω ( T ) + r T σ 2 2 δ ^ ( Ω ( T ) , T ) } P ( Ω ( T ) , T ) d Ω ( T ) ) = S ( 0 ) M q ( γ 1 , γ 2 ) . A 2 = E Q [ e r T K · I S ( T ) > K ] = γ 1 γ 2 e r T · K · P ( Ω ( T ) , T ) d Ω ( T ) = K e r T N q ( γ 1 , γ 2 ) .
where
M q ( γ 1 , γ 2 ) = e r T γ 1 γ 2 exp { σ Ω ( T ) σ 2 2 δ ^ ( Ω ( T ) , T ) } P ( Ω ( T ) , T ) d Ω ( T ) , N q ( γ 1 , γ 2 ) = γ 1 γ 2 P ( Ω ( T ) , T ) d Ω ( T ) .
Hence the value for a European call option is
C 0 = A 1 A 2 = S ( 0 ) M q ( γ 1 , γ 2 ) K e r T N q ( γ 1 , γ 2 ) .
 ☐
Corollary 1.
According to the parity formula of European call-put option, the price of European put option is
P 0 = S ( 0 ) ( M q ( γ 1 , γ 2 ) 1 ) K e r T ( N q ( γ 1 , γ 2 ) 1 ) .

5. Numerical Results

In this section, we use numerical simulation to compare LSM method and Borland’s methods from two aspects, one is the values of integral approximation, the other is the option prices.
To compare the three methods, we first generate the paths of this process by using the following Euler’s iterative formula:
Ω ( t i ) Ω ( t i 1 ) = 1 Z ( t i 1 ) 1 q 2 ( 1 β ( t i 1 ) ( 1 q ) Ω ( t i 1 ) 2 ) 1 2 ( W ( t i ) W ( t i 1 ) ) .
Moreover, the values of δ ( Ω ( t ) , t ) is calculate by rectangle integral principle that are expressed as
δ ( Ω ( t ) , t ) = i = 0 n 1 { 1 2 ( Z ( t i ) ( q 1 ) [ 1 β ( t i ) ( 1 q ) Ω ( t i ) 2 ] + Z ( t i + 1 ) ( q 1 ) [ 1 β ( t i + 1 ) ( 1 q ) Ω ( t i + 1 ) 2 ] ) × ( t i + 1 t i ) } ,
where n in Equations (28) and (29) represents the number of nodes in the path. t i refers to the time point corresponding to each node.
The steps of the numerical simulation mainly include three aspects represented as follows.
Step 1. Use the Euler’s iterative formula (28) to generate the path of Ω ( t ) .
Step 2. Use the rectangle integral principle (29) to generate δ ( Ω ( t ) , t ) .
Step 3. Perform the multiple nonlinear regression based on LSM via (17) to get the regression parameters.
Use 1000 paths, and each path has 100 nodes. The following Table 1 reports the corresponding simulation results of regression parameter values. As can be seen from Table 1, the value of each parameter d j varies with the difference of q.
Assuming q = 1.3 , T = 1 , we get the following Figure 1 which depicts the quadratic relationship between δ ^ and Ω ( T ) . Scatter points represent Monte Carlo simulation results. It can be seen that the approximate functions obtained by the three methods are close to each other in the place with dense scatter points, but differ greatly in the place with few scatter points. The approximation based on LSM method can better fit the sparse points distributed at both ends implying that our method is much closer to the results of Monte Carlo.
The following Table 2 is the results of errors of δ calculated by LSM and the other two methods, respectively. From Table 2, we can see that the values of error vary with the values of q. It is obviously to find that the δ ^ calculated with the LSM has the smallest errors in the three methods.
Before examining the numerical simulation of the call option price by three method, we first give the price formula of European call option deduced by Borland in 2002 and 2004, respectively.
Remark 3.
The price of European call option based on Borland’s method in 2002 is
C 0 = S ( 0 ) ξ 1 ξ 2 exp { σ Ω ( T ) σ 2 2 δ 1 ^ ( Ω ( T ) , T ) P ( Ω ( T ) , T ) d Ω ( T ) } K e r T ξ 1 ξ 2 P ( Ω ( T ) , T ) d Ω ( T ) ,
where δ 1 ^ is same as Equation (19) and ξ 1 , 2 = s 1 , 2 σ 1 ( 5 3 q ) β ( T ) , s 1 , 2 take the form as
s 1 , 2 = 1 ρ ( T ) ( 1 q ) σ β ( T ) ± 1 ρ ( T ) T 2 3 q ( 1 q ) 2 σ 2 β ( T ) 2 2 ( 1 q ) ρ ( T ) σ 2 β ( T ) × ( r T + ln S ( 0 ) K σ 2 2 ρ ( T ) ) ,
with
ρ ( t ) = 1 2 ( 3 q ) ( ( 2 q ) ( 3 q ) ε ) q 1 3 q · t 2 3 q .
Remark 4.
The price based on Borland’s method in 2004 is
C 0 = S ( 0 ) π 1 exp { σ Ω ( T ) σ 2 3 q ( 1 ( q 1 ) ( b 1 ( T ) + b 2 ( T Ω ( T ) 2 ) ) ) P ( Ω ( T ) , T ) d Ω ( T ) } K e r T π 1 P ( Ω ( T ) , T ) d Ω ( T ) ,
where π 1 = σ T 2 and b 1 , b 2 are defined by
b 1 ( t ) = ( ( 2 q ) ( 3 q ) ε ) q 1 3 q · 3 q 2 ( 9 5 q ) · t 2 3 q ,
and
b 2 ( t ) = 1 9 5 q .
Letting S ( 0 ) = 50 , r = 0.04 , σ = 0.2 , q = 1.3 and T = 0.8 , we use the price Formulas (23), (30) and (32) to calculate the option prices and absolute errors based on three method respectively. The absolute errors represent the absolute values of the difference of prices between each method and the Monte Carlo simulation. In Table 3, as can be seen, the price based on LSM is closest to the result of Monte Carlo simulation. Therefore, the obtained option price based on LSM is better than the two methods of Borland’s.
Figure 2 below is an implied volatility curve based on LSM model, which is plotted as a function of strike price K. We substitute the prices calculated by Equation (23) ( q = 1.3 , σ = 0.2 , S ( 0 ) = 50 , T = 0.8 and r = 0.04 ) into Blake-Scholes model to back out the implicit volatility. Obviously, these implied fluctuations form a smiling shape, very similar to the shape implied by real market data. Moreover, the downward sloping smile curve reproduces well-known systematic features of the volatility smile that appears when using the standard Black–Scholes to price real options.

6. Conclusions

This paper mainly solves the pricing problem of European options with a non-Gaussian model. Considering the characteristics of abnormal diffusion of financial asset prices, we apply the non-extended Tsallis entropy theory with this characteristic to the price movement model of assets. Since the integral 0 t f ( Ω ( s ) , s ) d s in the price formula cannot be calculated, we use LSM to evaluate. Finally, the numerical simulation results show that the least-square method is better than that of Borland’s in 2002 and 2004.

Author Contributions

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

Funding

This research was funded by the Henan Science and Technology Tackling Key Project (Grant No. 182102210364), Project of Humanities and Social Sciences Research Program of Ministry of Education (Grant No. 16YJA790063), National Social Science Fund Project (Grant No. 18BJY247).

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. Lo, A. Long term memory in stock market prices. Econometrica 1991, 59, 1279–1313. [Google Scholar] [CrossRef]
  6. Lux, T. The stable Paretian hypothesis and the frequency of large returns: An examination of major German stocks. Appl. Financ. Econ. 1996, 6, 463–485. [Google Scholar] [CrossRef]
  7. Merton, R.C. Option pricing when underlying returns are discontinuous. J. Financ. Econ. 1976, 3, 125–144. [Google Scholar] [CrossRef]
  8. Hull, J.C. Options, Futures, and Other Derivatives; Tsinghua University Press: Beijing, China, 2006. [Google Scholar]
  9. Necula, C. Option Pricing in a Fractional Brownian Motion Environment. 2002. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=1286833 (accessed on 24 September 2019).
  10. Hubalek, F.; Keller-Ressel, M.; Sgarra, C. Geometric asian option pricing in general affine stochastic volatility models with jumps. Quant. Financ. 2017, 17, 873–888. [Google Scholar] [CrossRef]
  11. 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]
  12. 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]
  13. Tsallis, C. Possible Generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  14. Borland, L. Microscopic dynamics of the nonlinear fokker-planck equation: A phenomenological model. Phys. Rev. E 1998, 57, 6634–6642. [Google Scholar] [CrossRef]
  15. Liu, L.; Cui, Y.; Xu, J.; Li, C.; Gao, Q. The Non-Markovian Property of q-Gaussian Process. Unpublished work. 2019. [Google Scholar]
  16. Borland, L. A theory of non-Gaussian option pricing. Quant. Financ. 2002, 2, 415–431. [Google Scholar]
  17. Ferrari, D.; Paterlini, S. Efficient and Robust Estimation for Financial Returns: An Approach Based on Q-Entropy. 2010. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=1906819 (accessed on 24 September 2019).
  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]
  20. Devi, S. Financial market dynamics: Superdiffusive or not? J. Stat. Mech. Theory Exp. 2017, 2017, 083207. [Google Scholar] [CrossRef]
  21. Liu, L.; Zhang, L.; Fan, S. Optimal investment problem under non-extensive statistical mechanics. Comput. Math. Appl. 2018, 75, 3549–3557. [Google Scholar] [CrossRef]
  22. 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]
  23. Borland, L.; Bouchaud, J. A non-Gaussian option pricing model with skew. Quant. Financ. 2004, 4, 499–514. [Google Scholar] [CrossRef]
  24. 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. [Google Scholar] [CrossRef]
  25. Plastino, A.R.; Plastino, A. Non-extensive statistical mechanics and generalized Fokker-planck equation. Phys. Stat. Mech. Its Appl. 1995, 222, 347–354. [Google Scholar] [CrossRef]
  26. Tsallis, C.; Bukman, D.J. Anomalous diffusion in the presence of external forces: Exact time-dependent solutions and their thermostatistical basis. Phys. Rev. Stat. Phys. Plasmas Fluids Relat. Interdiscip. Top. 1996, 54, R2197. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The horizontal axis is Ω ( T ) , and the vertical axis is the corresponding values of δ . Scatter points are obtained by Monte Carlo simulation.
Figure 1. The horizontal axis is Ω ( T ) , and the vertical axis is the corresponding values of δ . Scatter points are obtained by Monte Carlo simulation.
Entropy 21 00933 g001
Figure 2. The horizontal axis is K, and the vertical axis is the corresponding values of implied volatilities.
Figure 2. The horizontal axis is K, and the vertical axis is the corresponding values of implied volatilities.
Entropy 21 00933 g002
Table 1. Parameter values at different values of q.
Table 1. Parameter values at different values of q.
q Value d 1 ( q ) d 2 ( q ) d 3 ( q ) d 4 ( q )
q = 1.1 1.07230.02500.00670.0015
q = 1.2 1.55560.06140.02340.0059
q = 1.3 1.24010.10490.09040.0234
q = 1.4 1.19630.05760.44630.1468
q = 1.5 0.4591−0.26231.51901.0991
Table 2. Error analysis of the δ for three methods.
Table 2. Error analysis of the δ for three methods.
q Value Borland ( 2002 ) Borland ( 2004 ) LSM
q = 1.1 0.02740.01560.0154
q = 1.2 0.20580.15140.1433
q = 1.3 1.43251.42361.2196
q = 1.4 25.196829.820520.6674
q = 1.5 2179.12334.41316.6
Table 3. Error of option price for three methods.
Table 3. Error of option price for three methods.
Strike Monte Prices of Prices of Prices of Errors of Errors of Errors of
Price Carlo Borland ( 2002 ) Borland ( 2004 ) L S M Borland ( 2002 ) Borland ( 2004 ) L S M
457.45467.49007.49167.45220.03540.03710.0023
466.67496.70786.71086.67220.03290.03590.0027
475.93255.96285.96775.92980.03040.03520.0026
485.23315.26065.26775.23070.02750.03450.0025
494.58204.60614.61574.57960.02420.03370.0024
503.98604.00664.01913.98400.02060.03310.0020
513.45063.46733.48293.44900.01670.03230.0017
522.97992.99263.01132.97860.01270.03140.0013
532.57292.58192.60362.57210.00900.03080.0008
542.22502.23022.25482.22440.00520.02980.0006
551.92861.93061.95761.92820.00200.02900.0004

Share and Cite

MDPI and ACS Style

Liu, L.; Cui, Y. European Option Based on Least-Squares Method under Non-Extensive Statistical Mechanics. Entropy 2019, 21, 933. https://doi.org/10.3390/e21100933

AMA Style

Liu L, Cui Y. European Option Based on Least-Squares Method under Non-Extensive Statistical Mechanics. Entropy. 2019; 21(10):933. https://doi.org/10.3390/e21100933

Chicago/Turabian Style

Liu, Limin, and Yingying Cui. 2019. "European Option Based on Least-Squares Method under Non-Extensive Statistical Mechanics" Entropy 21, no. 10: 933. https://doi.org/10.3390/e21100933

APA Style

Liu, L., & Cui, Y. (2019). European Option Based on Least-Squares Method under Non-Extensive Statistical Mechanics. Entropy, 21(10), 933. https://doi.org/10.3390/e21100933

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