Recent Trends in Statistical Methodologies within Communications and Systems Engineering
View this Special IssueResearch Article  Open Access
Modeling and Risk Analysis Using Parametric Distributions with an Application in EquityLinked Securities
Abstract
In this study, we model the returns of a stock index using various parametric distribution models. There are four indices used in this study: HSCEI, KOSPI 200, S&P 500, and EURO STOXX 50. We applied 12 distributions to the data of these stock indices—Cauchy, Laplace, normal, Student’s t, skew normal, skew Cauchy, skew Laplace, skew Student’s t, hyperbolic, normal inverse Gaussian, variance gamma, and general hyperbolic—for the parametric distribution model. In order to choose the bestfit distribution for describing the stock index, we used the information criteria, goodnessoffit test, and graphical tail test for each stock index. We estimated the valueatrisk (VaR), one of the most popular management concepts in the area of risk management, for the return of stock indices. Furthermore, we applied the parametric distributions to the risk analysis of equitylinked securities (ELS) as they are a very popular financial product on the Korean financial market. Relevant risk measures, such as VaR and conditional tail expectation, are calculated using various distributions. For calculating the risk measures, we used Monte Carlo simulations under the bestfit distribution. According to the empirical results, investing in ELS is more risky than investing in securities, and the risk measure of the ELS heavily depends on the type of security.
1. Introduction
The normal or Gaussian distribution is a widespread distribution for modeling in finance. However, in general, financial asset return distributions are not normal, which is one of the stylized facts of stock returns (see [1]). Much empirical research has shown that real data for stock price returns are generally characterized by skewness, kurtosis, and fat tails (see [2–4]). Therefore, multiple distributions have been used as an alternative to the normal one.
First, by adjusting the number of degrees of freedom, the Student’s t distribution has a fat tail compared to the normal distribution. The Student’s t distribution is used in financial engineering, such as option pricing and risk management (see [5–7]). Second, skew distributions allow us to take advantage of the skewness value. The skew normal (see [8–10]) and skew Student’s t (see [11–13]) are typical skew distributions. Third, a normal variancemean mixture with generalized inverse Gaussian distribution can generate, for example, the generalized hyperbolic distribution introduced by BarndorffNielsen [14]. These types of distributions can be both symmetric and skewed, and their tails are heavier than those of the Gaussian distribution. The generalized hyperbolic distribution was used in many studies to fit a series of stock index returns (see [15–18]). Eberlein and Keller [19] and Eberlein [20], especially, demonstrated that a mediumtailed generalized hyperbolic family of distributions produces a more suitable fit to stock returns observed in the stock market. In addition, the variance gamma distribution is a subclass of the generalized hyperbolic distribution. Its skewness and kurtosis together describe the shape of the distribution. In the financial literature, Madan and Seneta [21] initially introduced the variance gamma distribution, and it has been used in various fields, for example, option pricing and synthetic CDO (collateralized debt obligations) pricing (see [22–24]).
In this study, we utilized twelve parametric distributions—Cauchy, Laplace, normal, Student’s t, skew normal, skew Cauchy, skew Laplace, skew Student’s t, hyperbolic, normal inverse gamma, variance gamma, and generalized hyperbolic distributions—to fit the distribution of the stock index return. We also used four different stock indices: HSCEI, KOSPI 200, S&P 500, and EURO STOXX 50. Furthermore, we applied the fitting results to the risk management because the return distribution plays an important role in risk measurement. These research approaches can be found in other studies (see [2, 12, 13, 25, 26]). Furthermore, Vernic [27], Bolance et al. [28], and Eling [29, 30] analyzed the skew normal and skewstudent as favorable models for describing actuarial loss of data and the investment returns of insurance companies by fitting various parametric distributions. The difference from the previous studies is in the applications of various parametric distributions to the model risk of equitylinked securities (ELS). Although many financial products have been sold on the Korean market, we focus on ELS because it is one of the bestselling financial derivatives. ELS are a hybrid debt security whose investment return is connected to an underlying equity, such as stock indices—a group of stocks. ELS worth approximately 45.9 trillion won were issued in 2013 in Korea. Furthermore, risk management for ELS has become critical in the last several years. For example, only recently, the financial supervisor warned that sales of ELS tracking the Hang Seng China Enterprises Index (HSCEI) will be restricted as there is an increasing possibility that if such sales continue, it could trigger a knockin option in most ELS products.
Since the Basel Committee on Banking Supervision introduced the valueatrisk (VaR) in 1996, it has been widely used as a risk measure in the risk management industry. VaR specifies the maximum amount an investment may lose, within a given probability, in a specified period of time. It is well known that VaR provides several benefits in terms of risk management. Above all, it captures an important aspect of risk in a single number, and it is useful to compare different assets and portfolios (see [31, 32]).
Measuring risk is crucial to control the risk and expected losses for both banks and companies. In addition, Hu and Kercheval [33] said that the choice of risk measure is less important for portfolio management than the choice of distribution family. In other words, the best choice of return distribution is a critical issue for portfolio management. However, there are two typical measures for risk, VaR (see [34]) and conditional tail expectation (CTE) (see [35]). They both have some pros and cons. For example, in the case of valueatrisk (VaR), it is easy to understand the VaR figures, but those figures are volatile. If we refer to CTE, it has subadditivity (the portfolio effect) and it is less volatile than the VaR. The use of risk measures—VaR and CTE—is required by the Basel Committee on Banking Supervision in determining a bank’s risk profile.
There are two different approaches to the estimation of VaR: the nonparametric method and the fully parametric method. A short review of these methods is as follows. Historical simulation is one of the most commonly used nonparametric approaches, and it is the easiest way to estimate VaR. One of the most popular parametric models is the deltanormal model. This model assumes that the return of a portfolio follows a normal distribution. However, this assumption is not supported by results from the real market; it has been shown that many financial assets have return distributions with “fatter tails” than present in the normal distribution (i.e., values that are further from the mean have higher probabilities than a normal distribution would suggest). In consequence, a great deal of effort has been put into developing parametric models to describe the fattailed distribution.
In this study, we use VaR and CTE in order to measure the risk of a stock index and financial derivative. Thus, we estimate VaR for the return of stock indices and evaluate the VaR backtest. In order to measure the risk of ELS, we implement several steps as follows: first, we investigate whether these models are appropriate for describing the stock indices—that is, we fit the return to the parametric distributions and find the bestfit distribution for each stock index. Second, we calculate VaR and CTE for a single stock index, portfolios consisting of stock indices and ELS, using the historical data. Finally, we use the Monte Carlo method to simulate VaR and CTE for the ELS.
The rest of this paper is as follows. In Section 2, we briefly review various parametric distributions. Then, we demonstrate the fitting results and determine the bestfit distributions via information criteria, goodnessoffit test, and graphical tail test for each stock index in Section 3. Section 4 gives the application results for the stock indices and ELS. Section 5 concludes the paper.
2. Distributions
As mentioned previously, we used 12 parametric distributions in this study: Cauchy, Laplace, normal, Student’s t, skew normal, skew Cauchy, skew Laplace, skew Student’s t, hyperbolic, normal inverse Gaussian, variance gamma, and general hyperbolic. The Cauchy distribution is named after the mathematician A.L. Cauchy (1789–1857). By using location μ and scale (σ > 0) as parameters, the probability density function (pdf) of the Cauchy distribution is
The random variable X has the Laplace distribution location μ and scale σ(>0) parameters if it has the following pdf:
A normal distribution is defined by two parameters: mean (or location) μ and variance (or scale) σ(>0). The pdf of the normal distribution is
The Student’s t distribution has heavier tails than a normal distribution and the degrees of freedom (ν ≥ 1) controls the heavytailedness. Student’s t distribution can be described by a normal and a Gamma distribution. Let Y be a random variable following the normal distribution with zeromean and the scale σ. Let U have a Gamma distribution in shape/rate parameterization (see [36]). Then, the random variableadmits the Student’s t distribution with mean μ, scale σ, and degrees of freedom parameter ν (see [37]). Its pdf iswhere Γ is the gamma function. Although the pdf of Student’s t can be generated by chisquare distribution, this distribution is the marginal posterior distribution for the normal mean with unknown variance and conjugate prior distribution (see [38]). Thus, threeparameter Student’s t distribution (5) is naturally used in many Bayesian inference problems.
Figure 1 illustrates the pdf of the Cauchy, Laplace, standard normal, and Student’s t distributions with different degrees of freedom. The pdf of Cauchy, Laplace, and Student’s t distribution has a fat tail compared to the normal distribution. Furthermore, the pdf of Student’s t converges to a normal distribution when the number of degrees of freedom ν tends to infinity. The case of ν = 1 is the Cauchy distribution.
The Cauchy, Laplace, normal, and Student’s t distribution cannot show the asymmetry of the pdf. Therefore, as an extension to the normal distribution, in order to accommodate asymmetry, the skew normal distribution was first suggested by Azzalini [8]. Let ϕ(⋅) and Φ(⋅) be the standard normal density function and the cumulative density function (cdf), respectively. A random variable X is said to have a skew normal distribution with the skew parameter . The pdf of the skew normal distribution is then given as
Location and scale parameters can be included via the linear transformation Y = μ + σX, which follows the skew normal distribution with location μ, scale σ, and shape parameter (denoted by SN(μ, σ^{2}, α)). The skew normal distribution has a kurtosis that is a little bit higher than the normal distribution, which is one of its limitations.
Similar to Azzalini’s study [8], Gupta and Huang [39] and Gupta et al. [40] introduced various skew distributions. One of them is the skew Cauchy distribution. The pdf of the skew Cauchy distribution is
Another technique of importing skewness into a symmetric distribution was developed by Fernández and Steel [41]. The idea was to convert a symmetric pdf into a skewed one by postulating inverse scale factors in the positive and negative orthants. Based on this technique, the pdf of an asymmetric Laplace distribution is given bywhere μ, σ(>0), and α(>0) are the location, scale, and skewness parameters, respectively. Figure 2 displays the pdf of skew Cauchy and skew Laplace distributions with varying α values. Through these figures, we can observe the asymmetric distributions.
(a)
(b)
According to Nurminen et al. [42], the univariate skew Student’s t distribution is parameterized by the location μ, scale σ, skew α, and degrees of freedom ν. The pdf of skew Student’s t distribution iswhere t(x; μ, σ, ν) is the pdf of Student’s t distribution given by (5) and . Also, T(·; 0, 1, ν) is the cdf of Student’s t distribution with degrees of freedom ν.
In Figure 3, the skew normal and skew Student’s t distributions with the skew parameter α = 0 are equivalent to the normal and Student’s t, respectively. Furthermore, according to the level of the skew parameter of the distribution, we can have positive or negative skewed pdf in the skew normal and skew Student’s t distributions.
(a)
(b)
The general hyperbolic distribution (ghyp) was introduced by BarndorffNielsen [43], and the ghyp can be parameterized in several ways. Following Prause et al. [44], the pdf of the ghyp iswhere K_{λ} denotes a modified Bessel function of the third kind with index λ (see [45]).
A remarkable aspect is that ghyp distributions envelop many special cases and limiting distributions, such as, for example, the Student’s t, normal, hyperbolic, and normal inverse Gaussian distributions (see [19]). We describe some of the special cases.
For λ = 1, we get the hyperbolic distribution. This type of distribution is also generated by the normal and generalized inverse Gaussian distributions. The hyperbolic distribution decreases exponentially, which is more slow than the normal distribution. The pdf of the hyperbolic distribution iswhere β is responsible for the steepness. The other μ, σ, and α are the location, scale, and skewness parameters, respectively.
We obtain the normal inverse Gaussian (NIG) distribution by letting in (11). The pdf of NIG iswhere σ > 0 and 0 < α < β.
On the contrary, the variance gamma distribution is also described by the normal and gamma distributions, and the tails of the variance gamma distribution decrease more slowly than the normal distribution. The pdf of variance gamma iswhere μ, σ, α, and λ(>0) are the location, scale, skewness, and shape parameters, respectively.
Figures 4–7 show the pdf of the hyp, NIG, ghyp, and variance gamma distributions, respectively. The plots can show a skewed and high kurtosis distribution by adjusting the skew parameter (α) and the kurtosis parameter (β, λ).
(a)
(b)
(a)
(b)
(a)
(b)
Table 1 summarizes the features of each distribution. These features are normality, kurtosis, and skewness, and they can be controlled by the skew Student’s t, hyperbolic, NIG, variance gamma, and generalized hyperbolic distributions.

3. Fitting and Selecting
3.1. Fitting Distribution
The four stock indices are considered for the period January 2013 to December 2017, amounting to 5 years of business days. Figure 8 illustrates the time series of price, log return, and histogram for each index. Furthermore, Table 2 describes basic statistics for the stock indices. By reading the skewness values presented in Table 2, we can observe that the return distributions for KOSPI 200 and S&P 500 show negative skewness. All return distributions are more fattailed relative to the normal distribution. In addition, Jarque–Bera statistics for all return series indicate that the four stock indices have nonnormal return distributions.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)

We can estimate the parameters of the distributions via maximum likelihood, and all models are implemented in the R packages fGarch, ghyp, sn, ald, VGAM, and MASS. In addition, the approach assumes that the time series is approximately independent and identically distributed. Hence, we plotted the autocorrelation function (ACF) of the log return series and squared the log return series for each stock index in Figure 9. Figure 9 indicates that while the log return series are serially uncorrelated, the squared log return series are serially correlated, which is in accordance with the stylized facts introduced by Cont [1]. In this study, we utilize the GARCH (1,1) model in order to remove the dependence in return series according to the arguments made by McNeil et al. [46]. We further present the definition of the GARCH (1,1) model.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
Definition 1. (GARCH (1,1)). The process (X_{t}) is GARCH (1,1) if it is covariance stationary and satisfies the following equations:where α_{0} > 0, α_{1} ≥ 0, β_{1} ≥ 0, and α_{1} + β_{1} < 1; z_{t} is a sequence of random variables following independent and identical distribution with zero mean and unit variance.
We create a filtered return series by subtracting the mean μ of the raw return series and then calibrating the GARCH (1,1) parameters (α_{0}, α_{1}, and β_{1}). We then obtain the filtered return series defined as follows: should be approximately independent and identically distributed (i.i.d).
The ACF of both filtered return series and squared filtered return series for each stock index is illustrated in Figure 10. Figure 10 shows that there is no dependence in both filtered return series and squared return series. Therefore, now we can apply the maximum likelihood method to the filtered return series of the four stock indices. The estimation results of the 12 distributions for each index are given in Tables 3–6 in Appendix.
The next step is to compare the distribution using several tools. In this study, we chose three tools to find the bestfit distribution for each stock index. The first such tool is the information criteria—Akaike information criterion and Bayesian information criterion. Following are the definitions of the information criteria.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
Definition 2. (Akaike information criterion [47]). Akaike information criterion (AIC) is a measure of goodness of fit defined aswhere k is the number of parameters to be estimated in the model, are the estimated parameters that maximize the likelihood (or loglikelihood), and is the maximum value of the loglikelihood.
Definition 3. (Bayesian information criterion [48]). Bayesian information criterion (BIC) or Schwarz criterion is a measure of goodness of fit defined aswhere n is the number of observations and are the estimated parameters that maximize the likelihood (or loglikelihood).
The AIC and BIC can be used to compare models based on different probability distributions. In the model selection application, the optimally fitted model is identified by the minimum value of AIC and BIC. We chose the fitted candidate distribution corresponding to the minimum value of AIC and BIC. This means that, for each index, we chose the bestfit distribution among the 12 distributions based on the AIC and BIC values given in Table 7.
The preferred distribution is the one with the lowest AIC and BIC values. Fortunately, for each stock index, the distribution has the lowest AIC and an equivalent BIC. Therefore, we chose the bestfit distribution for each stock index using the information criteria. Consequently, the bestfit distributions are Student’s t for HSCEI and variance gamma for KOSPI 200, S&P 500, and EURO STOXX 50.

3.2. GoodnessofFit Test
In order to enhance the robustness of the fitting results, we additionally performed the Kolmogorov–Smirnov test for the goodnessoffit test. According to the test results given in Table 8, we can determine which theoretical distribution differs significantly from the given return distribution for each stock index. Based on the Kolmogorov–Smirnov test, we determined that the distributions in bold are not able to describe the return distribution with the given significance level of 5%.

The Kolmogorov–Smirnov test uses the whole samples to calculate the statistics, which represent the maximum difference value between the empirical distribution function and the theoretical distribution function. However, in extreme cases, the left and the right tails of the return distribution are usually affected in terms of risk management. In other words, the tails of the return distribution and the risk management are interrelated. Therefore, we used a graphical left tail test for examining the fit in the tails. The graphical tests were performed as follows (see [49]):(i)Let denote the estimated cdf of the fitted distribution and (X_{(1)}, …, X_{(N)}) the order statistic of the historical data(ii)A plot of against X_{(t)} superimposed onto a plot of against X_{(t)} shows the left tail fit for the fitted distribution
In Figure 11, the circles correspond to the empirical data, the red line corresponds to the Cauchy distribution, the red dashed line to the Laplace, the red dashdotted line to the normal, the blue line to the Student’s t, the blue dashed line to the skew normal, the blue dashdotted line to the skew Cauchy, the black line to the skew Laplace, the black dashed line to the skew Student’s t, the black dashdotted line to the hyp, the green line to the NIG, the green dashed line to the variance gamma, and the green dashdotted line to the ghyp distributions.
(a)
(b)
(c)
(d)
Looking to each subfigure, we chose the bestfit candidate distributions for the left tail in the return distribution for each stock index. Therefore, skew Student’s t and ghyp distributions were chosen for HSCEI and the hyp, and variance gamma and ghyp distributions were chosen for KOSPI 200 and S&P 500. For EURO STOXX 50 and skew Student’s t, hyp and ghyp were chosen.
Based on the results from information criteria, Kolmogorov–Smirnov test, and the graphical tail test, we finally chose the bestfit distribution among the 12 parametric distributions for each stock index, which is shown in Table 9.

4. Risk Analysis
In order to measure the risk, we use the valueatrisk and the conditional tail expectation given in (19) and (20), respectively.
Definition 4. (valueatrisk (VaR)). VaR at a confidence level θ ∈ (0, 1) for loss L of a security or a portfolio is defined to bewhere F is the distribution function of loss.
Definition 5. (conditional tail expectation (CTE)). CTE at a confidence level θ ∈ (0, 1) for loss L of a security or a portfolio is defined to beAs a means of explanation, the VaR shows how much a portfolio loses within a certain time period. By contrast, the CTE indicates the expected loss whenever the occurred loss is greater than the VaR. In practice, the confidence level ranges from 95 to 99.5%, though the Basel committee recommends 99%.
4.1. Application to the Stock Index
In this section, we implement the backtest for VaR. The backtest period is 2013 to 2014. We compared the daily VaR estimates at the 99% confidence level. In order to forecast oneday ahead VaR, we used 250 business days, meaning that the window size is the last 250 observations. A rollingfixedwindow scheme is used to forecast oneday ahead VaR.
According to McNeil and Frey [50], the oneday ahead VaR under the GARCH model is given bywhere q_{θ} is the quantile at the given probability. and are the forecast value of standard deviation and mean, respectively, estimated by the GARCH model.
As discussed by Christoffersen [51], we will use an unconditional coverage test and a conditional convergence test to determine if the model is appropriate. The unconditional coverage test attempts to determine whether the observed ratio of exceptions is consistent with the ratio of expected exceptions according to the VaR model. On the contrary, the conditional coverage test investigates whether the total number of exceptions is equal to the expected one and the VaR exception process is independently distributed through time. A violation of the VaR model is defined bywhere L_{i} and VaR_{i}θ are the realized loss and VaR at ith day for given all N trading days.
We listed the results of the four indices from the VaR backtesting in Table 10. According to the backtest results, for the Student’s t, skew Laplace, hyp, NIG, variance gamma, and ghyp distributions, we cannot reject the null hypothesis for both unconditional and conditional coverage tests under confidence level 99%. In addition, the skew Laplace, hyp, NIG, variance gamma, and ghyp have actual violations smaller than the number of expected violations for the KOSPI 200 index. Most distributions have poor backtesting results for S&P 500. These results suggest that the optimal distribution for the risk measure of each stock index may vary.

4.2. Application to ELS
There are several kinds of ELS, for example, knockout, bullspread, digitalcall, hifive, and stepdown. In this study, we deal with twostock stepdown ELS because this type of ELS is the bestselling ELS product on the Korean derivative market. In order to calculate the VaR and CTE for the twostock stepdown ELS, we use the Monte Carlo simulation method. Since the twostock stepdown ELS incorporate twostock indices, we need to choose some pairs of stock indices. In this study, three pairs (HSCEI, KOSPI 200), (S&P 500, KOSPI 200), and (EURO STOXX 50, KOSPI 200) are selected for the underlying assets of the twostock stepdown ELS.
In this study, we consider a twostock stepdown ELS product with the following assumptions:(i)Maturity 3 years and the valuation period is 6 months(ii)Early redemption 90%, 90%, 90%, 85%, 85%, and 85%(iii)Knockin barrier 55%(iv)5.0% yield (annual)
To summarize, we calculate the VaR and CTE for the twostock stepdown ELS with three kinds of pairs following the above assumptions. For the Monte Carlo simulation method, the bestfit distribution is used to generate the return distribution for each stock index. The Monte Carlo simulation uses the following steps:(i)Generate the random variables under the bestfit distribution for each stock index(ii)Use the Cholesky decomposition and obtain the correlated two series of random variables(iii)According to filtering equation (16), obtain the defiltered returns generated by the bestfit distribution(iv)Calculate the VaR and CTE for the ELS
We defined a test period from 2013 to 2017. Because the ELS with a maturity of 3 years is taken for the investigation, we considered the number of periods of 3 years in terms of 750 business days from January 1, 2013, to December 31, 2017. When counting, we considered only the closing prices of the same business day for the four stock indices.
In order to use the Cholesky decomposition, we calculated the correlations between the filtered returns from the four stock indices. The correlations are given in Table 11. All the correlations are positive, which indicates that the world economy has a similar direction of trends, although there is a difference of magnitude within them. The correlation between S&P 500 and EURO STOXX 50 is the biggest, and the correlation between KOSPI 200 and S&P 500 is the smallest.

Figure 12 illustrates the relative prices of the four stock indices from 2005 to 2014. In Figure 12, all the relative prices seem to change similarly. Furthermore, sharp declines during the 2008 global financial crisis are evident.
In order to enhance the robustness of the simulation results, we calculated the VaR and CTE for direct investments using the historical data. The direct investment is used to make an equally weighted portfolio that consists of twostock indices. That is, we consider a portfolio of two stocks with weights 50% and 50% having an investment period of 3 years.
The VaR and CTE for the equally weighted portfolios are given in Table 12. Furthermore, we calculate the VaR and CTE for the twostock stepdown ELS using the historical data in Table 13. Tables 12 and 13 indicate that the indirect investment, the holding of the ELS, is generally riskier than the direct investment, especially for the pair (EURO STOXX 50, KOSPI 200), where the VaR and CTE for the ELS are higher than those of direct investment given by both confidence levels.


Furthermore, we investigated the number of autocallable cases using the historical data. The numbers and ratios of early repayment are illustrated in Table 14. From Table 14, we conclude that the first autocallable scored over 70% for the given observations. However, the ratio for loss depends heavily on the pair of stock indices.

The VaR and CTE calculated by the Monte Carlo simulation are given in Table 15. In addition, the investigation for the autocallable of ELS in the Monte Carlo simulation is given in Table 16. The Monte Carlo simulation results indicate that most autocallables occur first, and the ratio of loss is approximately 5%. Furthermore, from the simulation results, the inherited risk for ELS is larger than the risk calculated by the historical data.


According to the results above, the findings can be summarized as follows:(i)Direct investments are less risky than the indirect investments in cases of ELS(ii)Based on the historical data, the calculation of risk for ELS is inadequate to indicate the intrinsic risk of ELS to the investors(iii)The level of risk measures for ELS depends on a pair of stock indices
5. Summary
The motivation for conducting this study was to discover whether these distributions are appropriate for describing the return distribution. In order to find the bestfit distributions for the four stock indices, we implemented several steps. First, we calculated the information criteria, AIC and BIC. Second, we utilized Kolmogorov–Smirnov’s goodnessoffit test. Finally, we plotted the left tails of the fitted return distribution as a graphical test. Through these steps, we could choose the bestfit distribution for each stock index.
The empirical results provide a number of interesting conclusions, with useful practical implications. Our main findings can be summarized as follows.
First, we found the bestfit distribution for each stock index. The bestfit distributions are described by the general hyperbolic distribution, which can control both the kurtosis and the skewness. Therefore, we conclude that whether the distribution describes both kurtosis and skewness is crucial to find the bestfit return distribution for a stock index.
Second, we calculated the VaR for each stock index and implemented backtesting for the estimated VaR by each distribution. These test results indicate that the distribution that has both kurtosis and skewness is adequate for estimating the VaR.
Third, based on the backtest results, the long position for ELS is riskier than the direct investment for the portfolio consisting of the same stock indices. In other words, ELS hides the inherited risk of ELS’ payoff structure and attracts investors by giving the coupon higher than the return of money market.
The fourth conclusion is that the calculated measure for ELS depends on the pair of stock indices. In other words, the relationship between the twostock indices is the main factor for risk management of ELS.
Lastly, the Monte Carlo simulation results indicate that the ELS inherit a higher risk than the risk calculated from the historical data. In general, the backtest results are presented to the investors for the sales of ELS. Therefore, it is necessary to present the simulation results regarding the risk measures to assist the rational decision of the investors.
These results have at least two implementations. First, a timevarying correlation model can be used to calculate the VaR for ELS because we assume that the correlation between two returns is a constant for a given period. For example, dynamic conditional correlation (Engle [52]) can be considered. Second, the hedge performance for a portfolio with different distributions can be examined. The hedging is directly related with the business profit of portfolio management. Therefore, finding the proper hedge ratio for the portfolio depending on the return distribution should be considered.
Appendix
A. Parameter Estimation
The roles of the parameters are as follows:(i)μ: location parameter(ii)σ: scale parameter(iii)ν: degrees of freedom in Student’s t distribution(iv)α: skewness parameter(v)β: kurtosis parameter(vi)λ: shape parameter




Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The work of S.Y. Choi was supported by the Gachon University research fund of 2018 (GCU20180295) and by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (no. 2019R1G1A1010278). The research of J.H. Yoon was supported by the NRF of Korea grants NRF2017R1A5A1015722 and NRF2019R1A2C108931011.
References
 R. Cont, “Empirical properties of asset returns: stylized facts and statistical issues,” Quantitative Finance, vol. 1, no. 2, pp. 223–236, 2001. View at: Publisher Site  Google Scholar
 R. Huisman, R. Pownall, and K. Koedijk, “Varx: fat tails in financial risk management,” The Journal of Risk, vol. 1, no. 1, pp. 47–61, 1998. View at: Publisher Site  Google Scholar
 P. Verhoeven and M. McAleer, “Fat tails and asymmetry in financial volatility models,” Mathematics and Computers in Simulation, vol. 64, no. 34, pp. 351–361, 2004. View at: Publisher Site  Google Scholar
 Z. Rachev, B. RochevaIotovo, and S. Stoyanov, “Capturing fat tails,” Risk, vol. 23, no. 5, p. 72, 2010. View at: Google Scholar
 C.H. Lin and S.S. Shen, “Can the studentt distribution provide accurate value at risk?” The Journal of Risk Finance, vol. 7, no. 3, pp. 292–300, 2006. View at: Publisher Site  Google Scholar
 E. Platen and R. Rendek, “Empirical evidence on studenttLogreturns of diversified world stock indices,” Journal of Statistical Theory and Practice, vol. 2, no. 2, pp. 233–251, 2008. View at: Publisher Site  Google Scholar
 D. T. Cassidy, M. J. Hamp, and R. Ouyed, “Pricing European options with a log Student’s tdistribution: a Gosset formula,” Physica A: Statistical Mechanics and Its Applications, vol. 389, no. 24, pp. 5736–5748, 2010. View at: Publisher Site  Google Scholar
 A. Azzalini, “A class of distributions which includes the normal ones,” Scandinavian Journal of Statistics, vol. 12, no. 2, pp. 171–178, 1985. View at: Google Scholar
 T. Bodnar and A. K. Gupta, “Robustness of the inference procedures for the global minimum variance portfolio weights in a skewnormal model,” The European Journal of Finance, vol. 21, no. 1314, pp. 1176–1194, 2015. View at: Publisher Site  Google Scholar
 M. Taniguchi, A. Petkovic, T. Kase, T. DiCiccio, and A. C. Monti, “Robust portfolio estimation under skewnormal return processes,” The European Journal of Finance, vol. 21, no. 1314, pp. 1091–1112, 2015. View at: Publisher Site  Google Scholar
 M. C. Jones and M. J. Faddy, “A skew extension of the tdistribution, with applications,” Journal of the Royal Statistical Society: Series B (Statistical Methodology), vol. 65, no. 1, pp. 159–174, 2003. View at: Publisher Site  Google Scholar
 C. A. AbantoValle, V. H. Lachos, and D. K. Dey, “Bayesian estimation of a skewstudentt stochastic volatility model,” Methodology and Computing in Applied Probability, vol. 17, no. 3, pp. 721–738, 2015. View at: Publisher Site  Google Scholar
 C.T. Gao and X.H. Zhou, “Forecasting var and es using dynamic conditional score models and skew student distribution,” Economic Modelling, vol. 53, pp. 216–223, 2016. View at: Publisher Site  Google Scholar
 O. E. BarndorffNielsen, “Normal inverse Gaussian distributions and stochastic volatility modelling,” Scandinavian Journal of Statistics, vol. 24, no. 1, pp. 1–13, 1997. View at: Publisher Site  Google Scholar
 J. Fajardo and A. Farias, “Generalized hyperbolic distributions and brazilian data,” SSRN Electronic Journal, vol. 24, no. 2, 2002. View at: Google Scholar
 C. Necula, “Modeling heavytailed stock index returns using the generalized hyperbolic distribution,” Romanian Journal of Economic Forecasting, vol. 10, no. 2, pp. 118–131, 2009. View at: Google Scholar
 V. K. Socgnia and D. Wilcox, “A comparison of generalized hyperbolic distribution models for equity returns,” Journal of Applied Mathematics, vol. 2014, 2014. View at: Publisher Site  Google Scholar
 M. Takahashi, T. Watanabe, and Y. Omori, “Volatility and quantile forecasts by realized stochastic volatility models with generalized hyperbolic distribution,” International Journal of Forecasting, vol. 32, no. 2, pp. 437–457, 2016. View at: Publisher Site  Google Scholar
 E. Eberlein and U. Keller, “Hyperbolic distributions in finance,” Bernoulli, vol. 1, no. 3, pp. 281–299, 1995. View at: Publisher Site  Google Scholar
 E. Eberlein, “Application of generalized hyperbolic lévy motions to finance,” in Lévy Processes, pp. 319–336, Springer, Berlin, Germany, 2001. View at: Google Scholar
 D. B. Madan and E. Seneta, “The variance gamma (V.G.) model for share market returns,” The Journal of Business, vol. 63, no. 4, pp. 511–524, 1990. View at: Publisher Site  Google Scholar
 D. B. Madan, P. P. Carr, and E. C. Chang, “The variance gamma process and option pricing,” Review of Finance, vol. 2, no. 1, pp. 79–105, 1998. View at: Publisher Site  Google Scholar
 A. Tjetjep and E. Seneta, “Skewed normal variancemean models for asset pricing and the method of moments,” International Statistical Review, vol. 74, no. 1, pp. 109–126, 2006. View at: Publisher Site  Google Scholar
 T. Moosbrucker, “Pricing CDOS with correlated variance gamma distributions,” Journal of Fixed Income, vol. 12, pp. 1–30, 2006. View at: Google Scholar
 C. Adcock, M. Eling, and N. Loperfido, “Skewed distributions in finance and actuarial science: a review,” The European Journal of Finance, vol. 21, no. 1314, pp. 1253–1281, 2015. View at: Publisher Site  Google Scholar
 J. Nakajima, “Bayesian analysis of multivariate stochastic volatility with skew return distribution,” Econometric Reviews, vol. 36, no. 5, pp. 546–562, 2017. View at: Publisher Site  Google Scholar
 R. Vernic, “Multivariate skewnormal distributions with applications in insurance,” Insurance: Mathematics and Economics, vol. 38, no. 2, pp. 413–426, 2006. View at: Publisher Site  Google Scholar
 C. Bolance, M. Guillen, E. Pelican, and R. Vernic, “Skewed bivariate models and nonparametric estimation for the CTE risk measure,” Insurance: Mathematics and Economics, vol. 43, no. 3, pp. 386–393, 2008. View at: Publisher Site  Google Scholar
 M. Eling, “Fitting insurance claims to skewed distributions: are the skewnormal and skewstudent good models?” Insurance: Mathematics and Economics, vol. 51, no. 2, pp. 239–248, 2012. View at: Publisher Site  Google Scholar
 M. Eling, “Fitting asset returns to skewed distributions: are the skewnormal and skewstudent good models?” Insurance: Mathematics and Economics, vol. 59, pp. 45–56, 2014. View at: Publisher Site  Google Scholar
 K. Dowd, Beyond Value at Risk: The New Science of Risk Management, Wiley, Hoboken, NJ, USA, 1998.
 D. Duffie and J. Pan, “An overview of value at risk,” The Journal of Derivatives, vol. 4, no. 3, pp. 7–49, 1997. View at: Publisher Site  Google Scholar
 W. Hu and A. Kercheval, “Risk management with generalized hyperbolic distributions,” in Proceedings of the Fourth IASTED International Conference on Financial Engineering and Applications, pp. 19–24, ACTA Press, Berkeley, CA, USA, September 2007. View at: Google Scholar
 J. A. Lopez, “Methods for evaluating valueatrisk estimates,” Economic Policy Review, vol. 4, no. 3, 1998. View at: Google Scholar
 C. Acerbi and D. Tasche, “On the coherence of expected shortfall,” Journal of Banking & Finance, vol. 26, no. 7, pp. 1487–1503, 2002. View at: Publisher Site  Google Scholar
 M. H. DeGroot, Optimal Statistical Decisions, vol. 82, John Wiley & Sons, Hoboken, NJ, USA, 2005.
 M. Roth, On the Multivariate t Distribution, Linköping University Electronic Press, Linköping, Sweden, 2012.
 A. Gelman, J. B. Carlin, H. S. Stern, D. B. Dunson, A. Vehtari, and D. B. Rubin, Bayesian Data Analysis, Chapman and Hall/CRC, London, UK, 2013.
 A. K. Gupta and W.J. Huang, “Quadratic forms in skew normal variates,” Journal of Mathematical Analysis and Applications, vol. 273, no. 2, pp. 558–564, 2002. View at: Publisher Site  Google Scholar
 A. K. Gupta, F.C. Chang, and W.J. Huang, “Some skewsymmetric models,” Random Operators and Stochastic Equations, vol. 10, no. 2, pp. 133–140, 2002. View at: Publisher Site  Google Scholar
 C. Fernandez and M. F. J. Steel, “On bayesian modeling of fat tails and skewness,” Journal of the American Statistical Association, vol. 93, no. 441, pp. 359–371, 1998. View at: Publisher Site  Google Scholar
 H. Nurminen, T. Ardeshiri, R. Piche, and F. Gustafsson, “Robust inference for statespace models with skewed measurement noise,” IEEE Signal Processing Letters, vol. 22, no. 11, pp. 1898–1902, 2015. View at: Publisher Site  Google Scholar
 O. BarndorffNielsen, “Exponentially decreasing distributions for the logarithm of particle size,” Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences, vol. 353, no. 1674, pp. 401–419, 1977. View at: Publisher Site  Google Scholar
 K. Prause, “The generalized hyperbolic model: estimation, financial derivatives and risk measures,” University of Freiburg, Breisgau, Germany, 1999, Ph.D. thesis. View at: Google Scholar
 M. Abramowitz and I. A. Stegun, “Handbook of mathematical functions with formulas, graphs, and mathematical tables,” National Bureau of Standards Applied Mathematics Series 55, Dover, New York, NY, USA, 1972. View at: Google Scholar
 A. J. McNeil, R. Frey, and P. Embrechts, Quantitative Risk Management: Concepts, Techniques and Tools: Concepts, Techniques and Tools, Princeton University Press, Princeton, NJ, USA, 2015.
 H. Akaike, “Information theory and an extension of the maximum likelihood principle,” in Selected Papers of Hirotugu Akaike, pp. 199–213, Springer, Berlin, Germany, 1998. View at: Google Scholar
 G. Schwarz, “Estimating the dimension of a model,” The Annals of Statistics, vol. 6, no. 2, pp. 461–464, 1978. View at: Publisher Site  Google Scholar
 K. Aas and I. H. Haff, “The generalized hyperbolic skew student’s tdistribution,” Journal of Financial Econometrics, vol. 4, no. 2, pp. 275–309, 2006. View at: Publisher Site  Google Scholar
 A. J. McNeil and R. Frey, “Estimation of tailrelated risk measures for heteroscedastic financial time series: an extreme value approach,” Journal of Empirical Finance, vol. 7, no. 34, pp. 271–300, nov 2000. View at: Publisher Site  Google Scholar
 P. F. Christoffersen, “Evaluating interval forecasts,” International Economic Review, vol. 39, no. 4, pp. 841–862, 1998. View at: Publisher Site  Google Scholar
 R. Engle, “Dynamic conditional correlation,” Journal of Business & Economic Statistics, vol. 20, no. 3, pp. 339–350, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 SunYong Choi and JiHun Yoon. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.