Special Issue

## Stochastic Process Theory and Its Applications

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 1603509 | https://doi.org/10.1155/2020/1603509

Yan Chen, Ya Cai, Chengli Zheng, "Estimation of Tail Risk and Moments Using Option Prices with a Novel Pricing Model under a Distorted Lognormal Distribution", Mathematical Problems in Engineering, vol. 2020, Article ID 1603509, 25 pages, 2020. https://doi.org/10.1155/2020/1603509

# Estimation of Tail Risk and Moments Using Option Prices with a Novel Pricing Model under a Distorted Lognormal Distribution

Guest Editor: Wenguang Yu
Accepted15 Jun 2020
Published20 Jul 2020

#### Abstract

Risk measures based on the trading option prices in the market are forward-looking, such as VIX. We propose a new method combining distorted lognormal distribution with interpolation to price options accurately and then estimate tail risk. Our method can price the option of any strikes between the maximum and the minimum value of strikes in the real market, which reduces the instability and inaccuracy of using the limited option to measure the risk. In addition, our novel method treats the underlying asset price as a stochastic indicator rather than a fixed indicator as described in previous research studies for risk measurement. Moreover, even if the available sample size is very small, we can measure the risk stably and precisely after interpolation. Finally, the empirical test results of SP500 market show that this method has good performance, especially for the option markets with sparse strikes.

#### 1. Introduction

Risk measure is always being regarded as one of the most important parts for risk management. To quantify the risk well, the future distribution of the risk must be characterized, which involves forecasting and is difficult. The traditional method is based on the assumption that the history has a habit of repeating itself. So, the historical time series are used to estimate the risk measures. This method is backward-looking, which will never contain information about future.

To solve this problem, people turn to options market. Derivative markets are viewed as forward-looking markets because the prices of the derivative instruments traded in those markets incorporate the market participants’ perception of future market changes, which is extra information compared to those from stock markets. As a publicly traded market, the options market including stock options and index options can bring all the opinions about the future information of investors together. So, the options market can be exploited to compute the risk measures of their underlying assets. The standard approach is to determine the volatility of a stock or a stock index that is implied by today’s market prices of traded options.

The VIX volatility index is one of the application methods of option price. It is disseminated by the Chicago Board Options Exchange (CBOE) following the results of Whaley  and has attracted much attention in recent years. The index is computed from current price on a wide range of out-of-the-money European style call and put options written on the SP500 equity index, and it is built to serve as a model-free option-implied return volatility measure for the SP500 index over the coming 30 days, expressed in annualized percentage terms. VIX captures attention from the real market, and it is even labeled “fear gauge” and routinely cited by the media when describing current market conditions or “investor sentiment.” Besides the popularity in practice, VIX is increasingly used in financial economics, especially it provides a model-free measurement of the expected value of the SP500 return variation under the risk-neutral (Q probability) measure. Hence, the VIX embodies a model-free market volatility forecast. The literature about VIX can be classified in three main directions. One direction is about the forecast power of variation in the VIX for future realized return volatility. For example, Jiang and Tian  thought that the information content of the VIX volatility forecast is superior to alternative implied volatility measures as well as forecasts based on historical volatility. The second direction emphasizes the fact that the VIX measure is constructed directly from observed option prices so it is bound to incorporate any pricing of variance risk that may be embedded in the market prices. Thus, VIX does not represent a pure estimator of future return volatility for the underlying asset because it includes compensation for variance risk as well. Instead, the wedge between regular time series forecasts for volatility, developed under the actual or objective (P probability) measure, and the VIX, representing the risk-neutral return volatility forecast, may be interpreted as a market volatility risk premium (see, e.g., Bondarenko ; Bollerslev et al. ; and Carr and Wu  and their extensions Todorov  and Todorov and Bollerslev ). The third one is about seeking to draw inference regarding the stochastic properties of the underlying market return volatility process directly from the high-frequency behavior of the VIX (see Todorov and Tauchen  and Jiang and Tian ). Andersen et al.  considered a novel Corridor Implied Volatility index (CX) based on a range of strikes with high-frequency data including jumps and asymmetries. For more information about model-free volatility indices, see Gonzalez-Perez .

Besides research on VIX, there are some others using the ideas behind this kind of model-free estimation from options. Linders et al.,  proposed a novel group of herd behavior indices (HIX), which are model-free and risk-neutral, derived from available option data from the market. Their numerical illustration showed that the HIX based on the Dow Jones is identical to the CIX from Dhaene et al., , which further demonstrated that model-free estimation is an appropriate approach to calculate the HIX. Related empirical research studies are limited and need to be further explored.

This idea of model-free estimation based on options data is very important and contributes to reduce model error dramatically. However, their applications are based on the options market, where the number of options traded are limited. Usually, the strikes of options are located mainly close to the strike of at-the-money option, and the options with strikes far from that are less, even sparse. However, when this model-free method is applied, the options used mainly are those out-of-the-money options with strikes far from that of at-the-money option. Moreover, there are not enough options traded in some markets, such as China; the size of samples of data is very small. Both of the above cases will lead to non-negligible computation errors of risk measurement.

The second problem of the current method to compute risk measures is attributable to range estimation of risk measures. Andersen et al.,  thought that some extreme options value with strikes on the tails might be noises and will make the estimation of risk to be not correct and not stable. To make it robust, the corridor method is developed, discarding the data of tails and only considering a range of strikes. Dhaene et al.,  have the similar idea of range when building their downside heard behavior indices (DHIX). In fact, because the trading options cannot cover all ranges of strikes from zero to infinity, current risk measures based on market data are all within the strike range from the minimum to the maximum in the real market. The ranges of those measures are consistent with the fixed range of strikes, including Andersen et al.  and Linders et al. . But the real range of risk measures must be built on the range of underlying assets, which is a stochastic process, so the range must be random, too. Apparently, it is not suitable to describe risk features with a fixed range and there are differences between the two ranges of fixed and random. Therefore, the current method (proposed by VIX, Andersen et al.  and Linders et al. ) cannot deal with this problem well because they only used the option prices traded in the market and measure risk on a fixed range of strikes.

To make some improvements on above problems, we propose a new method based on distorted lognormal distribution and interpolation in this paper. Firstly, according to the characteristics of market price of options, we introduce a distorted variable to traditional lognormal distribution of underlying asset. This distorted variable is the function of the strikes, and it makes the options price of market to fit very well. Through market data, the distorted variables of options traded in the market can be estimated, and then using interpolation, all the distorted variables of arbitrary strikes between the minimum strike to maximum strike in the real market can be obtained, and then the prices of all the options (with or without options traded in the real market) can be computed well. This method can enlarge the size of samples of data from limited traded options, which will make the computation of risk measures more precise and robust. To test validation of our method, we decrease the size of samples gradually. And we present the whole equations for risk measures under the stochastic range using this method; it shows that our new method can solve the second problem very well. The empirical results show that the differences between the two ranges of fixed and random can be very big.

Besides the traditional risk measure such as volatility which is computed on the whole range from zero to infinity and their corridor versions, we also consider the computations of tail risk, including left tail and right tail, and their range versions, discarding the extreme points of the tails. The empirical results from SP500 market show that our novel method of distorted lognormal distribution combined with interpolation is very good, especially for the option markets with sparse strikes.

The rest of the paper is arranged as follows. In Section 2, the idea of distorted lognormal distribution will be introduced, and the empirical estimations of distorted variables and their interpolations will be presented. In Section 3, we will use this method to compute moments and tail risk estimations, compared with the current traditional method; especially, the stochastic ranges are treated for three different kinds of range moments and tail risk. Section 4 concludes the paper.

#### 2. Distorted Lognormal Distributions

Precise option price is the basis of computing risk measurement with the model-free estimation. But many empirical results (see, e.g., MacBeth and Merville ; Gultekin et al. ; and Long and Officer ) show that the traditional B-S-M model mispriced the options, sometimes undervalued options and sometimes overvalued options out of the money or in the money. So, we think this is one kind of investing psychology or behavior. The investors prefer to adjust subjective probability according to the moneyness. This makes the final equilibrium market prices vary with the moneyness, which may be the source of the puzzle of volatility smiles or smirks. So, we describe this adjustment of probability by a parameter , which varies with the moneyness . We can view that the market distorts the probability with the moneyness. The idea of distorting the normal distribution has been accepted and developed gradually to fit the “spike and fat tail” distribution better (see, e.g., Hamada and Sherris ; Godin et al. ; Labuschagne and Offwood ; Gerber and Shiu ; and Xiao-nan et al. ). Specifically, we assume that the distortion function is as follows:where is the price of underlying asset at the time and means expectation. Then, we can use this distortion function to options pricing.

##### 2.1. Theoretical Approach

In this section, we will explain how to define distorted variables of for a call option and for a put option. Given the price of a call and of a put as a function of , , , , , and . Now we consider the option pricing model based on the distortion function as equation (1). For simplicity, we only consider to distort the lognormal distribution here.

Under the traditional assumptions, an underlying stock price follows the stochastic process as follows:where is a Brownian motion under risk-neutral probability measure , is the risk-free rate, and is the constant volatility of . So, we havewhere with a probability distribution as

Here, we just distort the probability of directly with weight function as

is omitted as , briefly. So, the probability is transformed to , and is distorted to as

Obviously, it still follows a normal distribution. Based on this normal distribution, we can get the closed solution for European options easily.

For European call:where

Apparently, when , equation (7) is transformed into the traditional B-S-M model. We can estimate from the real market data of call options, and then we can price for call options by (7).

Here, it is further illustrated that this distorted method can solve the problem of volatility smiles and smirks, which is described as follows:where is the price based on the traditional B-S-M model, is the price in the real market, and is the implied volatility. For the same underlying assets with the same , but different strike price , the implied volatility varies with strike price , which means the traditional B-S-M model does not work in the real market. Here, we show that our new model (equation (7)) can value the options well and solve the problem of volatility smiles and smirks. We suppose that our price is the same as the real market price:and guess the shape of through the shape of volatility smiles and smirks, and then we compute the implied volatility as follows:

According to the shapes of volatility smiles and smirks reported, there are some kinds of distortion function that can be chosen, such as logistic, exponential, quadric, and their composite shapes. Here, we just show an example with quadric function as

Equation (13) means is zero at and is asymmetrical around . To show the volatility smiles and smirks, we use an example. We suppose that , and we select proper value for parameters and in equation (13). Figure 1 shows the example of ; Figures 2 and 3 show the corresponding volatility smile and smirk from equation (12).

For European put, we denote distortion as because we think it is different from that of call options. At this situation, we can get the closed solution for European put options easily.where the definitions of are the same as equations (8) and (9)fd9.

Apparently, when , equation (14) is transformed into the traditional B-S-M model.

Similarly, we can estimate from the real market data of put options, and then we can price the put options by equation (14). And we show volatility smiles and smirks with put option similarly. We still use quadric function as

Equation (16) means is zero at and is asymmetrical around . To show the volatility smiles and smirks, we still use the example with , and we select proper value for parameter and in equation (16). Figure 4 shows the example of ; Figures 5 and 6 show the corresponding volatility smile and smirk from equation (12). To show a little difference, here we show left smirk in Figure 6, while in Figure 3, we show a right smirk.

So, we show that our distorted method can explain the puzzle of volatility smiles and smirks well with simulation examples theoretically. To show that this method can be good in practice, we will use the real market data to fit in next section.

##### 2.2. Practice

We now consider a date with options available on the SP500 index and all the maturity shorter than 3 years.

We collect all data about options whose underlying asset is SP500 index. The whole data are collected from the Reuters database. The data period starts from the date of Feb 7, 2017, so the Current date is Feb 7, 2017. We have 2666 types of prices of European options, half of which are calls and the others are puts. The maturity of options listed in Table 1 has 11 types and each maturity corresponds to a different risk-free rate (%). On that day, the close price of SP500 index is 2293, which means , and the volatility .

 0.03 0.1 0.2 0.28 0.35 0.6 0.85 0.95 1.35 1.87 2.87 0.50 0.51 0.53 0.53 0.57 0.66 0.75 0.78 0.926 1.11 1.42 0.125 2293

The scatter graph of options market price with different moneyness is shown in Figure 7. The basic descriptive statistics are listed in Table 2. From Table 2, we can see that the moneyness in the real market changes around in the interval of . Its mean and median values are lower than 1 and the kurtosis can be nearly approximate to that of normal distribution. Both the market prices of calls and puts are right-skewness distributions, which means the market is more likely to be optimistic, especially for the puts market due to its significantly high kurtosis. Surely, the basic descriptive statistics show a little difference among different maturities; here we omit them. In Table 2, means a mixture of all the put and call options.

 Max Min Mean Median Std Skewness Kurtosis 1.53 0.04 0.79 0.81 0.26 −0.27 3.05 2136.45 0.05 311.49 92.63 431.86 1.65 5.27 2136.45 0.05 536.67 440.35 489.98 0.90 3.19 1225.40 0.05 86.32 7.60 177.71 3.30 15.71

Now, we represent the theoretical (for calls) and (for puts); for each strike, we solve for the exact value of so that the market price of the option matches the model price. We assume that our option pricing model based on distorted lognormal distributions is reasonable and equals the real market price:and then by using equation (17), we can estimate for calls as

Similarly, we can estimate for puts as

We call estimated from real market implied and implied . Then, we get all the implied and at all the 11 types of maturities . Figures 8 and 9 show the scatter graphs for all the implied and with different moneyness at all 11 types of maturities ; the basic descriptive statistics are listed in Table 3. From Figures 8 and 9, it can be seen that the implied and at different maturities varying with the moneyness have the similar trend, which means that the relationship between these two variables can be described by the same kind of function. We can find that the shapes of q are almost the same as we guessed in Figure 4; we can find a function similar to equation (13) to fit them well. But the shapes of m are not the same as we guessed in Figure 1, only similar; they are more complicated, with a big convex around (see Figures 9 and 10). Figures 10 and 11 show the relationship between implied (or ) and the moneyness at maturities . If we can find suitable functions to fit them very well, we can use the method of distorted lognormal distribution we proposed to all the prices of options with all the strikes from zero to (or big enough), whether they are traded in the real market or not. Then, we can use these options to compute the moments and tail risk measures more precisely. However, we cannot find a good function to fit them well, which may be our work in future. So, we use another method—interpolation—to extend the number of strikes. Though the interpolation cannot extend the range of strikes outside of the range of real market, it can enlarge the number of strikes inside that range, where there are no options traded in the real market. As we know, the number of options traded in the market is limited.

 T Call (m) Put () Max Min Mean Std Skewness Kurtosis Max Min Mean Std Skewness Kurtosis 0.03 5.15 −0.51 −0.17 0.58 5.46 39.50 −0.27 −71.41 −10.4 11.5 −2.02 8.69 0.1 4.82 −0.65 −0.17 0.54 6.99 56.60 −0.09 −48.17 −4.79 6.64 −3.09 15.90 0.2 0.30 −0.67 −0.21 0.17 −1.48 4.78 −0.05 −24.93 −3.34 4.34 −2.37 9.67 0.28 1.55 −0.66 −0.16 0.33 2.79 15.12 −0.06 −9.05 −2.62 2.65 −0.80 2.41 0.35 3.11 −0.66 −0.01 0.63 3.28 14.05 −0.08 −24.90 −4.01 4.98 −1.84 6.79 0.6 1.17 −0.62 −0.14 0.25 2.30 14.18 −0.11 −18.97 −3.02 3.70 −1.95 7.40 0.85 0.65 −0.54 −0.16 0.15 0.87 11.16 −0.13 −24.92 −3.30 4.22 −2.47 10.97 0.95 1.32 −0.58 −0.11 0.29 3.21 15.62 −0.15 −23.90 −3.01 4.04 −2.58 11.36 1.35 0.02 −0.48 −0.14 0.12 −0.71 3.47 −0.17 −19.52 −2.55 3.29 −2.82 12.63 1.87 0.03 −0.37 −0.13 0.10 0.01 1.74 −0.19 −16.26 −2.41 2.82 −2.43 10.20 2.87 0.02 −0.31 −0.13 0.11 −0.23 1.60 −0.27 −12.72 −1.84 2.04 −2.87 13.29

The method of interpolation will be useful for the markets where there are not enough options. To verify this, we narrow the number of samples and then make interpolation treatment. Figures 10 and 11 show the original data and their interpolations. The interpolation method here chosen is spline. The blue dot is the original data, there are 63 points for maturity , and the red curve denotes interpolation. We can see that even for SP500 options market, the points are not enough, especially for the right side where is high. For some option markets, such as some single stock-based option market and Chinese options market, there are not enough options traded, which will be similar to Figure 12. In Figure 12, we decrease the number of samples by deleting part of points. From Figures 12(a)12(d), the number of samples decreases and becomes more and more sparse. Here, we do not delete the first and the last original points to hold the same range of strikes. In the real market, that will not be the truth. For the traditional method, to replicate risk measures by options, enough samples are very important, especially for tail risk. In case of Figure 12(d), there are only 12 points; if we compute the tail risk, there are only one or two points that can be used, and the result may be very bad. Our interpolation method can overcome this problem.

In the next section, we will use this method to price the options and apply these prices to compute moments and tail risks.

#### 3. Applications to Moments and Tail Risk Estimation

In this section, we will use the methods suggested by Andersen et al.  and LDS (2014, 2015) to construct moments and tail risk measures via options with a range of strike prices, and for comparison, we will use the pricing model we proposed to compute them.

##### 3.1. Theory

Firstly, we will explain how to compute with twice differential as a function of the prices of the options and approximation for the finite number of strikes available.

We will use the same approach as the paper on corridor VIX  or on HIX [23, 24]. This approach is from the idea of Carr and Madan .

From Carr and Madan , for any twice continuously differentiable function , it can be expressed as follows:

Linders et al.  used this to express the swap rate with , as follows:where and . Equation (21) means that for any twice continuously differentiable function , its expectation can be expressed as the integral of European puts and calls at different strikes. Because the options market is looking forward, they can react with the information more actively for the future than the underlying assets markets, so we can use this method to construct some kind of measures in a way of model-free estimation to forecast the market more precisely.

Apparently, equation (21) requires that there must exist options with all the strikes from 0 to continuously. However, only some of them exist in the real options market with limited discrete strikes. To deal with this problem, Linders et al.  used the composite trapezoidal rule; they approximated formula (21) as follows:where

Throughout the whole paper, we will assume that there are only a finite number of European options with maturity . In particular, for the underlying asset (SP500 index), the strikes of the traded puts and calls are denoted by , withwhere is assumed to be finite. In reality, the underlying asset and call option based in it have unknown upward potential. For discrete approximation, here a finite upper bound can be chosen arbitrarily large. Note that as long as there is at least one strike for which the prices and are traded, the forward rate can be computed in model-free way using the put-call parity. In practical situations, we follow the methodology proposed in Chicago Board Options Exchange  to determine the forward rate of the SP500 index:with

Now we use this method to estimate some kinds of measures.

Example 1. Estimation of .
The first is to estimate moments of .
When , can be estimated.
When , we haveWhen , we have

Example 2. Estimation of range of variance .
In the real market, the extreme situations are not stable; this means that the options price with extreme strikes may not be good to estimate the variance, so these extreme strike situations must be discarded. Andersen et al.  used this method to replicate implied volatility; they called it Corridor Implied Volatility Index, or CX, and through the empirical test, they proved that this CX is more robust than the ordinary Implied Volatility Index (VIX). In fact, the basic theory they considered is as follows:where and ; apparently, and are not options in the real market, so they cannot use equation (22) to approximate with the real market prices of options. The real reason of this problem is that is stochastic. To solve this problem, they substitute it with , which is not stochastic. This is an approximation method. Then, equation (29) transforms toUsing equation (21), we havewhere and . If , then equation (31) will beIt only considers the range without the left and right tails and can be replicated both by call and put options which are out of the money. The Corridor Implied Volatility Index of Andersen et al.  is one of the special cases of this situation.
In most of situations, people only consider the downside risk, especially for the left tail; many risk measures are created to quantify the risk of this tail, such as Value at Risk () and Average Value at Risk (). Similarly, their calculations need to estimate the risk-neutral density (RND), which will be not forward-looking. Because the second derivatives for and are zero, we cannot replicate them by options prices as the method we proposed here. We can consider the moments higher than first order, such as second and third moments of tail, to quantify the risk of this tail. Using equation (32), moments of tails can be quantified by replicating with options.
If , then equation (32) will beIt only considers the left tail and can be replicated only by put options which are out of the money.
If , then equation (31) will beIt only considers the right tail and can be replicated only by call options which are out of the money.
For these three kinds of risk measures, we can compute arbitrary moments of them, except the first-order moment.
For the same reason, to calculate the moments through real options market, approximation by discretization must be done. Then using equation(22), equation (31) can be approximated by discretization as follows:where and . , and ; we only need to choose proper and to make the results stable or robust. That is a problem, not only because of selection of and , but also because of the calculation of the quantile of , which means we need to estimate the risk-neutral density (RND) of , which is not easy. How to solve this problem? Firstly, we try to determine by a fixed proportion of under the standard lognormal assumption simply, such as . And can be chosen based onwhere refers to geometrical Brownian motion.
But, we find it is not good. Figure 13 shows with the confidence . From this figure, we can find that even when confidence is approaching zero, is still higher than 0.8. From Table 2, we can see that the median value of , which means this method will lose a large part of information of the options market.
To avoid this difficulty, Andersen et al.  proposed a good method. They adopted the approach of Andersen and Bonarenko . It builds on the fact that option prices reflect tail moments. The left and right tail moments of a positive random variable, with strictly positive density for all , are given byand then they define the ratio statistic, , as an indicator of how far in the tail a given point, , is located within the support of : is akin to a cumulative density function or CDF as it is increasing on with and . So, for a given percentile, , in the range of the function, we can get the quotient, , asLet denote the risk-neutral density for SP500 forward price at maturity , or one month, and be the strike price of European style put and call options; we haveThe ratio statistic is computed only from put and call prices. Hence, if the range is located where option prices may be extracted reliably, can be computed without estimation of risk-neutral density. Here, we can borrow this method to determine the . After considering the market liquidity, we can settle the truncation level for the tails. That is,Utilizing this method, we can compute the moments via replicating options prices in the real market. However, the options traded are limited, and these computations may not be so good. And all of the range moments are just approximations as shown in equation (30); it will have big error. To solve these two problems, we propose the alternative approach in the next section.

##### 3.2. Alternative Approach

As mentioned above, the options traded in the real market are limited; especially, for some extreme strikes, it will be sparse. Using the limited options to replicate the risk measures will lead to errors. And when calculating the range moments, the method of approximations as shown in equation (30) is not convincing. In this section, we will propose some methods to solve these two problems.

In the first section, we propose a novel method to evaluate the options well. It shows that the weight function can be fitted very well. So, for the first problem, we can use the method of distorted lognormal distribution to get all and , from the limited options traded in the real market. Then, we use interpolation to get enough points of and between the minimum and maximum of . Based on this, we can get all the prices of call and put options with these strikes (). Then, we can use equation (21) to compute the risk measures precisely.

In fact, in equation (21), the options are supposed to be existing in the real market with price , and by using the implied volatility equation, we havewhere is for call options and for put options. Then, equation (20) transforms into

This means that equation (43) is the function of implied volatility of different options with all the strikes. Because we do not have the real traded prices for options with all the strikes, we can approximate equation (43) well. Here, we remember that if and can be obtained, we can use them to price all the options with all the strikes. Namely, we assume that the price of our new price model can fit the market prices well, so we can offer the option prices for strikes that are not traded in the real market:

So, equation (43) can be rewritten as

and are from our new pricing model equations (7) and (14). Using this method, we can compute equation (45) smoothly with integral. In fact, if the points by interpolation are enough, equation (22) is still a good approximation to (45):

Moreover, denote ; then, we havewhereand the definitions of are

We can see that equations (48) and (49) are the functions of , and .

For the second problem, we can still solve it using our novel method of distorted lognormal distribution. In equation (29),where and .

Because and are not options in the real market, we cannot use equation (22) to approximate with the real market prices of options. The papers before [14, 24] use a nonstochastic indicator to take the place of the stochastic indicator . Using our method, we can calculate them easily:wherewhere and are from equations (7) and (14), and

As we show that our novel distorted method can price the options very well, all , , , and at arbitrary strikes between of the real market (which corresponds to ) can be calculated precisely. We do not need to use the approximation equation (30) any more. Similarly, we can get different situations of range moments, including tail moments as equations (32) and (34):

If , then equation (52) transforms to

This situation is to measure the risk of the left tail; it can be replicated by a series of put options out of the money.

If , then equation (53) transforms to

This situation is to measure the risk of the middle without left and right tails; it can be replicated by a series of put and call options out of the money.

If , then equation (53) transforms to

This situation is to measure the risk of the right tail; it can be replicated by a series of call options out of the money.

Compared with equations (32) and (34), the results with stochastic indicator are different. The first difference is that the integral terms are different. The integral terms with stochastic indicator have some adjustment items. The second difference is the integral interval. For the tail moments, the intervals with stochastic indicator are bigger than those with nonstochastic indicator. So, we can expect that the final results of numerical calculation of moments may be different.

We can compute the integral of these equations by MATLAB. Surely, if there are enough interpolation points (here we have 3000 points), we can use approximation method of equation (22). Similarly, equation (52) can be approximated by discretization as follows:whereand and other variables have the same definitions as that of Section 3.1.

In next section, we will use the real data of SP options market to calculate the moments and tail risk and compare these different methods mentioned here.

##### 3.3. Comparison

In this section, we will utilize the real data of SP options market to calculate the moments and tail risk and compare these different methods mentioned here. We compare three kinds of approaches to compute two moments of risks with 3 pairs of ranges.

The three kinds of approaches are as follows: the first approach is to replicate risk measures by limited options traded in the market, used by current papers (see [14, 24]); the second approach is to replicate risk measures by enough interpolation points through our novel distorted lognormal distribution; these two methods are based on the nonstochastic indicator . Here, we call the first method NSO (nonstochastic indicator with original data) and call the second method NSDI (nonstochastic indicator with distorted and interpolation data). And the third approach is to replicate risk measures by enough interpolation points through our novel distorted lognormal distribution based on the stochastic indicator . We call the third method SDI (stochastic indicator with distorted and interpolation data).

The two moments are for with . We choose three pairs of range through changing. Here, means calculations will be done by all the options traded in the market without discarding the extreme strikes of the tails. This situation can be viewed as a special case of range risk measures. The range is the maximum range of strikes for the options traded in the real market. means calculations will be done by some of the options with discarding of the extreme strikes of the both tails. According to Andersen et al. , this range will make the risk measures more stable.

The second pair of range is for the left tail risk without and with discarding of the extreme strikes of the left tail.

The third pair of range is for the right tail risk without and with discarding of the extreme strikes of the right tail.

We want to check if discarding of the extreme strikes can make the risk measure to be more stable. will be replicated by some of the put options out of the money with discarding of the extreme strikes of the left tail to make it stable. The right tail risk will be replicated by some of the call options out of the money with discarding of the extreme strikes of the right tail to make it stable.

Besides these situations, to show that our interpolation methods are good for markets with sparse options, we decrease the number of the samples and then compute those risk measures again. The method to decrease is as follows: . For example, there are samples; means only 6 samples are left, and they are ; means only 2 samples are left, and they are . Here, we will take . We can see that the number of samples decreases to of original data. With increasing, when or , it will be very sparse, and there are not enough samples to finish the calculations, especially for the tail risk.

The data are the same as those in Section 2.2, and we use the method of Section 2 to estimate and , and then we use the interpolation method to get all the points of and between and . We will get 3000 points, which is enough.

Firstly, we compute with . To be convenient, we denote

We compute all with different situations for all the 11 maturities. Here, we only show results of as an example because the results for the other maturities are similar.

Table 4 shows the results of and its truncated range with three approaches and 6 kinds of number of samples (). It shows that . is expected because the former is coarser than the latter. can be explained too. They have the same integral interval, but there is a negative adjustment of integral term for because of the truncated tails. For , there are still truncated tails with the range of . For , the range is smaller, so difference between and is bigger (theoretically, with the narrowing of range, this difference will increase first and then decrease to zero from equation (52)). This means that the risks of this kind (with ) can be overvalued by fixed indicator methods ( and ), relative to the stochastic indicator method (). From Table 4, it shows that the distorted and interpolation methods ( and ) are more stable than the traditional method () through different samples from to . And the range moments are more stable than .

 k NSO NSDI SDI NSO NSDI SDI 1 0.0762 0.0761 0.0760 0.0713 0.0707 0.0574 2 0.0765 0.0759 0.0757 0.0720 0.0706 0.0574 3 0.0773 0.0761 0.0759 0.0725 0.0707 0.0575 4 0.0780 0.0754 0.0748 0.0743 0.0704 0.0575 5 0.0907 0.0877 0.0877 0.0766 0.0709 0.0579 6 0.0806 0.0755 0.0744 0.0777 0.0707 0.0575

Table 5 shows the results of