# Interval Estimation of Value-at-Risk Based on Nonparametric Models

Hussein Khraibani, Bilal Nehme, Olivier Strauss

► To cite this version:

Hussein Khraibani, Bilal Nehme, Olivier Strauss. Interval Estimation of Value-at-Risk Based on Nonparametric Models. Econometrics, 2018, 6 (4), pp.47-77. 10.3390/econometrics6040047. lirmm-01951982

HAL Id: lirmm-01951982

https://hal-lirmm.ccsd.cnrs.fr/lirmm-01951982

Submitted on 11 Dec 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Interval Estimation of Value-at-Risk Based on Nonparametric Models

Hussein Khraibani¹*, Bilal Nehme² and Olivier Strauss³

¹ Department of Applied Mathematics, Faculty of Sciences, Lebanese University, Beirut 2038 1003, Lebanon

² Department of Economics, Faculty of Economic Sciences & Business Administration, Lebanese University, Beirut 2038 1003, Lebanon; Bilal.Nehme@ul.edu.lb

³ Department of Robotics, LIRMM University of Montpellier II, 61 rue Ada, 34392 Montpellier CEDEX 5, France; strauss@lirmm.fr

* Correspondence: hussein.khraibani@gmail.com; Tel.: +961-70453332

Received: 6 December 2018; Accepted: 6 December 2018; Published: 10 December 2018

<figure><img src="image_2.png" /></figure>

**Abstract:** Value-at-Risk (VaR) has become the most important benchmark for measuring risk in portfolios of different types of financial instruments. However, as reported by many authors, estimating VaR is subject to a high level of uncertainty. One of the sources of uncertainty stems from the dependence of the VaR estimation on the choice of the computation method. As we show in our experiment, the lower the number of samples, the higher this dependence. In this paper, we propose a new nonparametric approach called *maxitive kernel estimation* of the VaR. This estimation is based on a coherent extension of the kernel-based estimation of the cumulative distribution function to convex sets of kernel. We thus obtain a convex set of VaR estimates gathering all the conventional estimates based on a kernel belonging to the above considered convex set. We illustrate this method in an empirical application to daily stock returns. We compare the approach we propose to other parametric and nonparametric approaches. In our experiment, we show that the interval-valued estimate of the VaR we obtain is likely to lead to more careful decision, i.e., decisions that cannot be biased by an arbitrary choice of the computation method. In fact, the imprecision of the obtained interval-valued estimate is likely to be representative of the uncertainty in VaR estimate.

**Keywords:** risk measures; quantile estimation; financial time series; Value-at-Risk; choquet integral; possibility theory; maxitive kernel; kernel estimation; parametric models

## 1. Introduction

Controlling financial risk is an important issue for financial institutions. For the necessity of risk management, the first task is to measure risk. Value-at-Risk (VaR) is probably the most widely used risk measurement in financial institutions. It has made its way into the Basel II Capital-Adequacy framework. VaR is an estimate of the largest loss over a specified time horizon at a particular probability level. For example, if the daily VaR of an investment portfolio is $10 million with a 95% confidence level, this means that we can be 95% confident that the portfolio will not lose more than $10 million over the next day. In a more formal way, the VaR of a portfolio at a probability level $\alpha$ can be defined as McNeil et al. (2005, p. 38).

$$VaR_{\alpha} = \inf\{x, F_X(x) = P(X \le x) \ge \alpha\} = F_X^{-1}(\alpha), \quad (1)$$

where $X$ is a random variable representing daily percent return for a total stock index with Cumulative Distribution Function (CDF) $F_X$. As McNeil et al. (2005) point out, VaR is simply a quantile of the corresponding loss distribution, which makes its computation easy. Various methods have been proposed in the relevant literature to compute the VaR. Each method differs in the methodology used,
the hypotheses, and the way the models are implemented. A highly questionable fact is that each method leads to different results. Thus the choice of computational method is likely to have a large impact on the way a financial institution manages its credit portfolio.

The most commonly used approach is the nonparametric Historical Simulation (HS) method described in Linsmeier and Pearson (1997). It is based on the empirical CDF of the historically simulated returns by attributing an equal probability weight to each day's return. The HS approach is easy to implement but suffers from two major drawbacks: First, the success of the approach depends on the ability to collect a large series of data; second, this is an unconditional model and thus we need a number of extreme scenarios in the historical record to provide more informative estimates of the tail of the loss distribution McNeil et al. (2005). Moreover, estimates of extreme quantiles are inefficient since extrapolation beyond past observations is impossible. To avoid these drawbacks, several parametric and nonparametric approaches have been developed. For example, Butler and Schachter (1998) propose the kernel estimators in conjunction with the historical simulation method. Also, Charpentier and Oulidi (2010) calculate the VaR by using several nonparametric estimators based on Beta kernel. This study shows that those estimators improve the efficiency of traditional ones, not only for light tailed, but also heavy tailed distributions.

Numerous recent VaR models are referred to as parametric. From the point of view of the risk manager, estimate VaR assuming normal distribution of asset returns is inappropriate and lead to underestimate the left tail at low probability level. For this reason, most recent research papers deal with going beyond the normal model and attempt to capture the related phenomena of heavy and long tails and asymmetric form of the returns series. EVT (Extreme Value Theory) and the so-called GHYP (Generalized HYPerbolic) distributions are among the most widely used. The main advantage of hypothesizing a GHYP distribution is its ability to account for the statistical properties of financial market data such as volatility clustering, asymmetry and heavy-tail phenomena (see McNeil et al. (2005) for an introduction and Paolella and Polak (2015) for a recent application). Kuester et al. (2006) use an EVT-based approach and focuses on the long tails of the return distribution. Braione and Scholtes (2016) study the performance of forecasting VaR under different parametric distributional assumptions and show the predominance and the predictive ability for the skewed and heavy-tailed distributions in the univariate case. The main drawback of using a parametric model is the high dependency of the obtained method to the hypothesized distribution model. To lower this dependency, some authors have proposed a group of semiparametric methods that are based on extreme value theory. These methods have been successfully used by financial analysts in estimating VaR Danielsson and De Vries (2000). However, all the above mentioned attempts have the common weakness of a low robustness to modelling leading to variations in the VaR estimation. To enable this method to be used in a prudent manner, it would be of instrumental interest to estimate the variation of the computed VaR w.r.t. the modelling. Some attempts have been proposed to achieve such a goal. For example, Butler and Schachter (1998) propose a method to measure the precision of a VaR estimate. Jorion (1996) suggests that VaR always be reported with confidence intervals and shows that it is possible to improve the efficiency of VaR estimates using their standard errors. Kupiec (1995) proposes a method for quantifying the uncertainty, in the estimated VaR, induced by the fact that the return distribution is unknown. Pritsker (1997) proposes to estimate a standard error by using a Monte Carlo VaR analysis.

In this paper, we propose a completely new approach to estimate the variations in the VaR induced by choosing a particular model in a kernel-based approach. The idea is that an estimate of the CDF of the daily percent return can be used to compute the VaR as suggested by Equation (1). Such an estimate can be obtained by using a kernel-based approach Silvermann (1986). In Kernel Density Estimation (KDE), the role of the kernel is to achieve an interpolation to lower the impact of sampling on the obtained estimation. Roughly speaking, it smoothes the empirical CDF. However, as in the above mentioned methods, there is a systematic bias induced by the choice of the kernel used to estimate the CDF. Our proposal is to focus on the new nonparametric approach developed by
Loquin and Strauss (2008a) to estimate the CDF. This approach makes use of the ability of a new kind of interpolating kernel, called **maxitive kernel**, to represent a convex set of conventional kernels—that we call **summative kernels** Loquin and Strauss (2008b). Within this approach, the estimate of the CDF is interval-valued. Such an interval-valued estimate of the CDF, also called a p-box Destercke et al. (2007), is the convex set of all the CDF that would have been computed by using the convex set of conventional summative kernel represented by the maxitive kernel. We show that this approach can advantageously be used to compute with accuracy the corresponding convex set of kernel-based VaR estimates. This approach has advantages over Monte Carlo based approaches since its complexity is comparable to that of classical kernel estimation. Moreover, within this approach the bounds are exact while Monte Carlo approaches provide an inner estimation of those bounds.

This paper is structured as follows. Section 2 reviews some nonparametric and parametric approaches and bootstrap methods used to compute the confidence interval of VaR. In Section 3, we introduce the empirical distribution function and the kernel cumulative estimator based on summative kernel. We define in Section 4 the maxitive kernel, which forms the basis of our approach, and it is shown how an interval-valued estimation of VaR, based on maxitive kernel, can be derived that has relevant properties in this context. Section 5 presents and discusses our empirical findings. Firstly, we show how the choice of kernel function affects the VaR estimates. Secondly, we investigate the performance of the interval-valued proposed in this paper by comparing it to three very competitive approaches: The simple Normal VaR, the HS VaR, and the GHYP VaR. Finally, Section 6 concludes this paper with some further remarks.

Throughout this paper, we consider that the observations belong to a convex and compact subset (universe) $\Omega$ of $\mathbb{R}$, called *the reference set*. $P(\Omega)$ is the collection of all measurable subsets of $\Omega$. It naturally contains the empty set and is closed under complementation and countable unions. Then $(\Omega, P(\Omega))$ is a measurable space. Let $L(\Omega)$ be the set of functions defined on $\Omega$ with values in $\mathbb{R}$.

## 2. A Review of Some Statistical Approaches and Nonparametric Bootstrap Methods

### 2.1. Historical Simulation

The HS method calculates the Value-at-Risk using real historical data of asset returns and captures the non-normal distribution of the returns. HS is a nonparametric method because it doesn't make a specific assumption about distribution of returns. However HS method assumes that the distribution of past returns is a good and complete representation of expected future returns.

The VaR with $\alpha$% of probability level is calculated as the $\alpha$% percentile of the sorted data return values. For example, with a returns data with 1000 observations, the 1% VaR estimate is simply the negative of the 10th sample order statistic. This HS VaR can be defined as follows:

$$ \text{VaR}_{\alpha} = \text{percentile} \left( \{r_t\}_{t=1}^T, \alpha\% \right), $$

where $r_t$ is the asset return at time t.

### 2.2. GHYP Parametric Distribution

The generalized hyperbolic distribution (GHYP) was introduced by Barndorff-Nielsen (1978) to fit financial returns. The GHYP is an asymmetric heavy-tailed distribution that can account for the extreme events and cater for skewness embedded in the data. It has since been applied in diverse disciplines such as physics, biology, financial mathematics (see Eberlein and Keller (1995); Sørensen and Bibby (1997); Paolella (2007, chp. 9)).

The probability density function (pdf) of univariate GHYP distribution with the parameterization of Eberlein et al. (1998) is given, for $x \in \mathbb{R}$, by:
$$f_{\text{GHD}}(x; \lambda, \chi, \psi, \mu, \sigma^2, \gamma) = a K_{\lambda-1/2} \left( \sqrt{\left(\chi + \frac{(x-\mu)^2}{\sigma^2}\right) \left(\psi + \frac{\gamma^2}{\sigma^2}\right)} \right) \exp\left(\gamma \frac{(x-\mu)}{\sigma^2}\right) \left(\chi + \frac{(x-\mu)^2}{\sigma^2}\right) \left(\psi + \frac{\gamma^2}{\sigma^2}\right)^{\lambda/2-1/4},$$

with the normalizing constant

$$a = \frac{(\chi\psi)^{-\lambda/2}\psi^{\lambda}(\psi + \frac{\gamma^2}{\sigma^2})^{1/2-\lambda}}{(2\pi)^{1/2}\sigma K_{\lambda}(\sqrt{\chi\psi})},$$

where

(a) $K_\lambda$ denotes a modified Bessel function of the third kind with index $\lambda$.

(b) $\lambda$ defines the subclasses of GHYP and is related to the fail flatness.

(c) $\chi$ and $\psi$ determine the distribution shape; in general, the larger those parameters are, the closer the distribution is to the normal distribution.

(d) $\mu$ is the location parameter.

(e) $\sigma$ is the dispersion parameter (standard deviation)

(f) $\gamma$ is the skewness parameter (if $\gamma = 0$, the distribution reduces to the symmetric generalized hyperbolic distribution).

The GHYP family contains many special cases known under special names, listed as follows:

* **Hyperbolic Distribution (HYP)**

If $\lambda = 1$, we get the hyperbolic distribution. However, HYP distribution is characterized by having a hyperbolic log-density function whereas the logdensity for the normal distribution is a parabola. Thus, one may expect the HYP distribution to be coherent alternatives for heavy tailed data.

* **Normal Inverse Gaussian (NIG) Distribution**

If $\lambda = -\frac{1}{2}$, then the distribution is known as normal inverse gaussian (NIG). NIG distribution is also widely used in financial modeling.

* **Variance Gamma (VG) Distribution**

If $\lambda > 0$ and $\chi = 0$, then we obtain the limiting case which is known as variance gamma (VG) distribution.

* **The Skewed Student's -Distribution (St)**

When $\lambda < 0$ and $\chi = 0$ we get another limiting case called the generalized hyperbolic skew student's t distribution because it generalizes the usual Student t distribution, obtained from the skewed-t by setting the skewness parameter $\gamma = 0$.

In order to estimate the unknown parameters $(\lambda, \chi, \psi, \mu, \sigma^2, \gamma)$, one can use the maximum likelihood (ML) with numerical optimization method (EM-based algorithm).

## 2.3. Nonparametric Bootstrap Approach

The bootstrap method is used in an important number of statistics topics based on building a sampling distribution for a statistic by resampling from the data at hand. The term bootstrap was coined by Efron (1979), and is an allusion to the expression *pulling oneself up by one's bootstraps*—in this case, using the sample data as a population from which repeated samples are drawn Fox (2002). The nonparametric bootstrap allows us to empirically estimate the sampling distribution of a statistic $\theta$—such as a mean, median, standard deviation, or quantile (VaR)—or an estimator $\hat{\theta}$ without making assumptions about the form of the population, and without deriving the sampling distribution explicitly. The basic idea of the nonparametric bootstrap is as follows: Assuming a data set $S = (x_1, \dots, x_n)$ is available. We draw a sample of size $n$ from among the elements of $S$, sampling with
replacement. This result is called bootstrap sample $S_1^* = (x_{11}^*,...,x_{1n}^*)$. We repeat this procedure a large number of times, B, selecting many bootstrap samples; the bth such bootstrap sample is denoted $S_b^* = \{x_{b1}^*,...,x_{bn}^*\}$. Next, we compute the statistic $\theta$ for each of the bootstrap samples; that is $\theta_b^* = t(x_b^*)$, with t denoting some function, that we can use to estimate from data. Then the distribution of $\theta_b^*$ around the original estimate $\theta$ is analogous to the sampling distribution of the estimator $\theta$ around the population parameter. The variance is estimated by the sample variance (but for the bootstrap sample of $\hat{\theta}$) and the bias is estimated by the difference between the average of the bootstrap sample and the original $\hat{\theta}$. Then, the bootstrap estimates of bias and variance are given approximately by:

$$ \text{Bias}(\hat{\theta}^*) \approx \frac{1}{B} \sum_{b=1}^{B} \hat{\theta}_i^* - \hat{\theta} = \hat{\theta}^* - \hat{\theta} \quad \text{and} \quad \text{var}(\hat{\theta}^*) \approx \frac{1}{B-1} \sum_{b=1}^{B} (\hat{\theta}_i^* - \hat{\theta}^*)^2. $$

There are various methods for constructing bootstrap confidence intervals. The normal-theory interval assumes that the distribution of $\hat{\theta}$ is normally distributed (which is often approximately the case for statistics in sufficiently large samples), and uses the bootstrap estimate of sampling variance, and perhaps of bias, to construct a 100(1 - $\alpha$)-percent confidence interval of the form:

$$ [(\hat{\theta} - \hat{\theta}^*) - z_{1-\frac{\alpha}{2}} se(\theta^*), (\hat{\theta} - \hat{\theta}^*) + z_{1-\frac{\alpha}{2}} se(\theta^*)], $$

where $se(\hat{\theta}^*) = \sqrt{var(\hat{\theta}^*)}$ is the bootstrap estimate of the standard error of $\hat{\theta}$, and $z_{1-\frac{\alpha}{2}}$ is the $1-\frac{\alpha}{2}$ quantile of the standard-normal distribution (e.g., 1.96 for a 95-percent confidence interval, where $\alpha = 0.05$).

We can also estimate confidence intervals using percentiles of the sample distribution: The upper and lower bounds of the confidence interval would be given by the percentile points (or quantiles) of the sample distribution of parameter estimates. This alternative approach, called the bootstrap percentile interval, is to use the empirical quantiles of $\theta_b^*$ to form a confidence interval $[\theta_{(lo)}^*, \theta_{(up)}^*]$, where $\theta_{(1)}^*,...,\theta_{(r)}^*$ are the ordered bootstrap replicates of the statistic; $lo = [(r+1)\alpha/2]$; up = $[(r+1)(1-\alpha/2)]$; and $[\cdot]$ is the nearest integer function. This basic percentile interval approach is limited itself, particularly if parameter estimators are biased Dowd (2005). It is therefore often better to use the bias-corrected, accelerated (or BC$_a$) percentile intervals. In order to find the BC$_a$, we firstly calculate $z = \Phi^{-1}\rho$, where $\Phi^{-1}(\cdot)$ is the standard-normal quantile function, and $\rho = \frac{\sum_{b=1}^{n} (\theta_b^* \le \theta)}{r+1}$ is the adjusted proportion of bootstrap replicates at or below the original-sample estimate $\theta$. If the bootstrap sampling distribution is symmetric, and if $\theta$ is unbiased, then this proportion will be close to 0.5, and $z$ will be close to 0. Now, let $\theta_{(-i)}$ be the value of $\theta$ produced when the ith observation is deleted from the sample there are n of these quantities. Let $\bar{\theta} = \frac{1}{n} \sum_{i=1}^{n} \theta_{(-i)}$ be the average of the $\theta_{(-i)}$. We calculate

$$ \beta = \frac{\sum_{i=1}^{n} (\theta_{(-i)} - \bar{\theta})^3}{6 \left[ \sum_{i=1}^{n} (\theta_{(-i)} - \bar{\theta}^2) \right]^{\frac{3}{2}}}. $$

We calculate

$$ \beta_1 = \Phi \left[ z + \frac{z - z_{1-\frac{\alpha}{2}}}{1 - \beta(z - z_{1-\frac{\alpha}{2}})} \right] \quad \text{and} \quad \beta_2 = \Phi \left[ z + \frac{z + z_{1-\frac{\alpha}{2}}}{1 - \beta(z + z_{1-\frac{\alpha}{2}})} \right], $$

where $\Phi(\cdot)$ is the standard-normal cumulative distribution function. The values $\beta_1$ and $\beta_2$ are used to locate the endpoints of the corrected percentile confidence interval: $[\theta_{(lo^*)}^*, \theta_{(up^*)}^*]$, where $lo^* = [r\beta_1]$ and $up^* = [r\beta_2]$. When $\beta$ and $z$ are both 0, $\beta_1 = \Phi(-z_{1-\frac{\alpha}{2}}) = \Phi(z_{\frac{\alpha}{2}}) = \frac{\alpha}{2}$ and $\beta_2 = \Phi(z_{1-\frac{\alpha}{2}}) = 1 - \frac{\alpha}{2}$, which corresponds to the (uncorrected) percentile interval.

For more background on the bootstrap approach and a broader array of applications see Efron and Tibshirani (1993); Davison and Hinkley (1997).
The bootstrap approach can be used to estimate the Value-at-Risk (VaR) and a confidence intervals for VaR. A resampling method based on the bootstrap and a bias-correction to improving the Value-at-Risk (VaR) forecasting ability of the normal-GARCH model has been developed in Hartz et al. (2006). As mentioned in Dowd (2005), if we have a data set of n observations, we create a new data set by taking n drawings, each taken from the whole of the original data set. Each new data set created in this way gives us a new VaR estimate. We then create a large number of such data sets and estimate the VaR of each. The resulting VaR distribution function enables us to obtain estimates of the confidence interval for our VaR. For estimate confidence intervals using a bootstrap approach, we produce a bootstrapped histogram of resample-based VaR estimates, and then read the confidence interval from the quantiles of this histogram. A very good discussion on this and other improvements can be seen in Dowd (2005, chp. 4). In Section 5.3, we use the bias corrected method to estimate the confidence interval of VaR obtained from HS, Normal and GHYP distribution. The purpose of this simulation exercise is to compare these methodologies with our maxitive kernel approach, in order to identify that our interval-valued estimation of VaR perform properly.

# 3. Summative Kernels

## 3.1. Summative Kernels and Probability Distributions

Kernels are used in signal processing and nonparametric statistics to define a weighted neighborhood around a location $u \in \Omega$. In kernel regression, it is used to estimate the conditional expectation of a random variable (see e.g., Silvermann 1986; Wand and Jones 1995).

**Definition 1.** A summative kernel—or conventional kernel—is a positive function $\kappa: \Omega \to \mathbb{R}^+$ that verifies the summativity property

$$\int_{\Omega} \kappa(u)du = 1.$$

From a basic summative kernel $\kappa$, we can define a summative kernel $\kappa_{\Delta}^x$ translated in $x \in \Omega$ and dilated with a bandwidth $\Delta > 0$ by

$$\forall u \in \Omega, \kappa_{\Delta}^{x}(u) = \frac{1}{\Delta}\kappa(\frac{u-x}{\Delta}).$$

By convention, $\kappa(u) = \kappa_1^0(u)$.

The most used summative kernels in functional estimation are usually unimodal, symmetric, and centered (i.e., defining a weighted neighborhood around the origin). When $\kappa(u)$ is symmetric and unimodal function, then the following conditions are fulfilled

$$\int_{\Omega} u\kappa(u)du = 0 \text{ and } \int_{\Omega} u^2\kappa(u)du = c > 0.$$

For example, the Epanechnikov kernel proposed by Epanechnikov (1969) illustrated in Figure 1 is defined by

$$\kappa(u) = \begin{cases} \frac{3}{4}(1-u^2) & \text{if } |u| \le 1, \\ 0 & \text{otherwise.} \end{cases} \qquad (2)$$

In the following, $\mathcal{K}(\Omega)$ represents the set of unimodal, symmetric, and centered kernels. A summative kernel $\kappa$ can be seen as a probability distribution inducing a probability measure $P_{\kappa}: \mathcal{P}(\Omega) \to [0, 1]$ on the measurable space $(\Omega, \mathcal{P}(\Omega))$. $P_{\kappa}$ is defined by

$$\forall A \subseteq \Omega, P_{\kappa}(A) = \int_{A} \kappa(u)du.$$
<figure><img src="image_3.png" /><figcaption>Figure 1. Epanechnikov kernel with a bandwidth $\Delta = 1$.</figcaption></figure>

Considering the summative kernel $\kappa$, we can define what we call the *summative expectation* as follows:

**Definition 2.** Let $s$ be a function of $L(\Omega)$ and let $\kappa$ be a summative kernel of $K(\Omega)$. The summative expectation of $s$, in the neighborhood defined by $\kappa$, is the classical expectation of $s$ w.r.t. $P_\kappa$:

$$ \mathbb{E}_{\kappa}(s) = \int_{\Omega} sdP_{\kappa} = \int_{\Omega} s(u)\kappa(u)du, \qquad (3) $$

since $dP_\kappa = \kappa(u)du$.

## 3.2. Empirical Distribution Function

Let $(x_1, \dots, x_n)$ be a finite set of observations of $n$ random variables $(X_1, \dots, X_n)$ i.i.d with unknown pdf $f$ and CDF $F$. One can estimate $F$ by the empirical CDF $F_n$ defined by

$$ F_n(x) = \frac{1}{n} \sum_{i=1}^{n} 1(x_i \le x), \qquad (4) $$

where $1$ is the indicator function, namely $1(x_i \le x) = 1$ if $x_i \le x$ and zero otherwise. Obviously $1(x_i \le x) = 1$ with probability $P(x_i \le x) = F(x)$, and $1(x_i \le x) = 0$ with probability $P(x_i > x) = 1 - F(x)$. Then $1(x_i \le x)$ is a Bernoulli random variable with success probability $F(x)$. Since the $x_i$ are independent, so are the $1(x_i \le x)$. Thus $Y_n = \sum_{i=1}^{n} 1(x_i \le x)$ (a sum of independent $Ber\{F(x)\}$ random variables) is a $Bin\{n, F(x)\}$ random variable, and

$$ F_n(x) \sim \frac{1}{n} Bin\{n, F(x)\}. $$

Then $\forall x \in \Omega$,

$$ \mathbb{E}(F_n(x)) = F(x) \text{ and } \mathrm{Var}(F_n(x)) = \frac{F(x)(1-F(x))}{n}. $$

This implies that the empirical CDF is convergent in probability to the true CDF:

$$ \forall x \in \Omega, F_n(x) \xrightarrow{P} F(x). $$

Thus the empirical estimate of Value-at-Risk
$$VaR_{\alpha,n} = \inf \{x, F_n(x) = P(X \le x) \ge \alpha\}, \qquad (5)$$

converges towards the true VaR$_{α}$.

### 3.3. Summative Kernel Cumulative Estimator

The empirical distribution function (Equation (4)) is not smooth as it jumps by $\frac{1}{n}$ at each point $x_i$ ($i = 1...n$). A kernel estimator of the CDF, introduced by authors such as Nadaraya (1964); and Watson and Leadbetter (1964), is a smoothed version of the empirical distribution estimator. Such an estimator arises as an integral of the Parzen-Rosenblatt kernel density estimator (see Rosenblatt 1956; Parzen 1962).

**Definition 3.** The summative-kernel cumulative estimator of the CDF (also called the Parzen-Rosenblatt kernel cumulative estimator), based on the summative kernel $\kappa$ is given, in each point $x \in \Omega$, by

$$\hat{F}_{\kappa}^{n}(x) = \frac{1}{n} \sum_{i=1}^{n} \Gamma(x - x_i), \qquad (6)$$

where $\Gamma(x) = \int_{-\infty}^{x} \kappa(u)du$.

Let $\kappa \in \mathcal{K}(\Omega)$ be a summative kernel function of the second order with support $[-1, 1]$, the function $\Gamma(u)$ verifies the following properties

$$\Gamma(u) = \begin{cases} 0 & \text{if } u \in [-\infty, -1], \\ 1 & \text{if } u \in [1, \infty]. \end{cases}$$

$$\begin{cases} \int_{-1}^{1} \Gamma^{2}(u)du \le \int_{-1}^{1} \Gamma(u)du = 1, \\ \int_{-1}^{1} \Gamma(u)\kappa(u)du = \frac{1}{2}, \\ \int_{-1}^{1} u\Gamma(u)\kappa(u)du = \frac{1}{2}(1 - \int_{-1}^{1} \Gamma^{2}(u)du). \end{cases}$$

Note that if $\kappa$ is the pdf of a probability measure $P_\kappa$, then $\Gamma$ is the cumulative distribution function of $P_\kappa$. For example, $\kappa$ being the Epanechnikov kernel defined by (2), the function $\Gamma$, depicted in Figure 2, has the following expression:

$$\Gamma(u) = \begin{cases} 0 & \text{if } u \le -1, \\ -\frac{1}{4}u^3 + \frac{3}{4}u + \frac{1}{2} & \text{if } |u| \le 1, \\ 1 & \text{if } u \ge 1. \end{cases} \qquad (7)$$

Let $\kappa \in \mathcal{K}(\Omega)$ be a summative kernel with support $[-1, 1]$, for a fixed $x$, the bias and the variance of $\hat{F}_{\kappa}^{n}(x)$ are given by

$$\mathbb{E}(\hat{F}_{\kappa}^{n}(x)) - F(x) = \frac{\Delta^2}{2} f'(x) \int_{\Omega} u^2 \kappa(u) du + o(\Delta^2),$$

$$\mathrm{var}(\hat{F}_{\kappa}^{n}(x)) = \frac{1}{n} \left\{ F(x)(1-F(x)) + \Delta f(x) \left( \int_{-1}^{+1} \kappa^2(u) du - 1 \right) + o(\Delta) \right\}.$$

Some theoretical properties of the estimator $\hat{F}_{\kappa}^{n}$ have been investigated, among others, by Winter (1973); Winter (1979); Sarda (1993); Yamato (1973); and Jones (1990). Some properties have long been known, e.g., the uniform convergence of $\hat{F}_{\kappa}^{n}$ towards $F$ when the pdf $f = dF$ is
continuous Nadaraya (1964); and Yamato (1973) or without conditions on f Singh et al. (1983). The asymptotic expression of the Mean Integrated Squared Error (MISE) (i.e., $\int_{\Omega}(\hat{F}_{\kappa}^{n}(x) - F(x))dx)^2$ has been studied in Swanepoel (1988). For a continuous pdf f, it has been proved that the best kernel is the uniform kernel of bandwidth $\Delta > 0$ defined by $\kappa(u) = \frac{1}{2\Delta}1_{[-\Delta,\Delta]}(u)$ and for a discontinuous function f, in a finite number of points, the best kernel is the exponential kernel of bandwidth $\Delta > 0$ defined by $\kappa(u) = \frac{\Delta}{2}e^{-\Delta|u|}$.

<figure><img src="image_4.png" /><figcaption>Figure 2. Integrated Epanechnikov $\kappa$ kernel with a bandwidth $\Delta = 1$.</figcaption></figure>

The method of choosing the optimal value of bandwidth in kernel estimation of the cumulative distribution function is of crucial interest. Many procedures—such as plug-in and cross validation—have been proposed in the relevant literature (see e.g., Sarda 1993; Polansky and Baker 2000; Altman and Leger 1995) to choose (estimate) the optimal bandwidth for estimating the CDF of the random process underlying a sample. As mentioned in Quintela Del Rio and Estevez-Perez (2013), the Polansky and Baker plug-in bandwidth is given by

$$ \Delta_{\text{PB}} = \left( \frac{\zeta(\kappa)}{-n\xi_2^2(\kappa)\hat{V}_2(g_2)} \right)^{1/3}, $$

where

$$ \xi_2(\kappa) = \int_{\Omega} u^2 \kappa(u) du \quad ; \quad \zeta(\kappa) = 2 \int_{\Omega} u \kappa(u) \Gamma(u) du, $$

$$ \hat{V}_r(g) = \frac{1}{n^2 g^{r+1}} \sum_{i=1}^{n} \sum_{j=1}^{n} H^{(r)} \left( \frac{x_i - x_j}{g} \right) $$

estimates

$$ V_r = \int_{\Omega} f^{(r)}(u) f(u) du, $$

where $r \ge 2$ is an even integer and

$$ g_2 = \left( \frac{2H^{(2)}(0)}{-n\xi_2^2(H)V_4} \right)^{1/5}. $$

The kernel function H is not necessarily equal to $\kappa$. An iterative method for calculating the plug-in bandwidth has been proposed by Polansky and Baker (2000). As mentioned in Quintela Del Rio and
Estevez-Perez (2013), the plug-in bandwidth is calculated as follows: let $b > 0$ be an integer, firstly, we calculate $\hat{V}_{2b+2}$ using

$$ \hat{V}_r = \frac{(-1)^{r/2}r!}{(2\hat{\sigma}(x_i))^{r+1}(r/2)!(3.14)^{1/2}}, $$

where $\sigma$ is the standard deviation of the data which can be estimated by $\hat{\sigma}(x_i) = \min\{\hat{s}, \frac{Q_3-Q_1}{1.349}\}$, with $\hat{s}$ is the sample standard deviation, and $Q_1, Q_3$ denoting the first and third quartile, respectively. Secondly, begin from $j = b$ to $j = 1$, we calculate $\hat{V}_{2j}(\hat{s}_{2j})$, where

$$ \hat{s}_{2j} = \left( \frac{2H^{(2j)}(0)}{-n\hat{\zeta}_2(H)\hat{V}_{2j+2}} \right)^{1/(2j+3)}, $$

with

$$ \hat{V}_{2j+2} = \begin{cases} \hat{V}_{2b+2}, & \text{if } j = b, \\ \hat{V}_{2j+2}(\hat{s}_{2j+2}), & \text{if } j < b. \end{cases} $$

Thirdly, the plug-in bandwidth is

$$ \hat{\Delta}_{\text{PB}} = \left( \frac{\zeta(\kappa)}{-n\hat{\zeta}_2^2(\kappa)\hat{V}_{2j}(\hat{s}_{2j})} \right)^{1/3}. $$

In practice, for most applications, we consider $b = 2$ Quintela Del Rio and Estevez-Perez (2013). Sarda (1993) introduced a cross-validation method to estimate the optimal bandwidth which minimize the MISE:

$$ CV(\Delta) = \sum_{i=1}^{n} \left( F_n(x_i) - \hat{F}_{\kappa;-i}^{n}(x_i) \right)^2, $$

where $F_n$ is the empirical distribution function (Equation (4)) and $\hat{F}_{\kappa;-i}^n$ is the kernel cumulative distribution estimator computed by leaving out $x_i$:

$$ \hat{F}_{\kappa;-i}^{n}(u) = \frac{1}{n-1} \sum_{j \neq i} \Gamma(u - x_j). $$

Adrian et al. (1998) proposed a modified cross-validation method which minimizes the function

$$ CV(\Delta) = \frac{1}{n} \sum_{i=1}^{n} \int_{\Omega} \left( \mathbb{1}(u - x_i) - \hat{F}_{\kappa;-i}^{n}(u) \right)^2 du. $$

For more details about bandwidth selection of summative kernel CDF estimation, we refer our reader to Quintela Del Rio and Estevez-Perez (2013).

Now, if we suppose that all parameters have been chosen appropriately for $\hat{F}_{\kappa}^n$ to be a consistent estimate of the CDF, then a kernel estimator of VaR, denoted $\hat{VaR}_{\kappa,n}$, can be easily obtained by inverting the equation $\hat{F}_{\kappa}^n(x) = \alpha$. In that case, $\hat{VaR}_{\kappa,n}$ satisfies:

$$ \frac{1}{n} \sum_{i=1}^{n} \Gamma(\hat{VaR}_{\kappa,n} - x_i) = \alpha. \qquad (8) $$

The kernel VaR estimator $\hat{VaR}_{\kappa,n}$ can be seen as a smoothed version of the empirical distribution function of VaR ($VaR_{\kappa,n}$) defined by (Equation (5)). For further details about the properties of the kernel VaR estimator (see e.g., Gourieroux et al. 2000; Chen and Tang 2005; Sheather and Marron 1990). However, the validity of such an estimate highly depends on the appropriateness of the chosen kernel and bandwidth. In the next Section, we propose another approach that allows to estimate the
dependance of the obtained estimate to the parameters of the method, and therefore to increase the robustness of the decision process based on the data.

# 4. Maxitive Kernels

## 4.1. Maxitive Kernels and Possibility Distributions

In many applications the summative kernels and their bandwidth are chosen in a very empirical way. As proposed in Loquin and Strauss (2008b), the empirical character of choosing a kernel could be taken in consideration by taking a family of kernels rather than one kernel. This family of kernels can be represented by using a maxitive kernel.

**Definition 4.** A maxitive kernel is a positive function $\pi : \Omega \to [0, 1]$ that verifies the following maxitivity property

$$ \sup_{u \in \Omega} \pi(u) = 1. $$

From a basic maxitive kernel $\pi$, we can define a maxitive kernel $\pi_{\Delta}^{x}$ translated in $x \in \Omega$ and dilated with a bandwidth $\Delta > 0$ by

$$ \forall u \in \Omega, \pi_{\Delta}^{x}(u) = \pi(\frac{u-x}{\Delta}). $$

By convention, $\pi(u) = \pi_{1}^{0}(u)$.

For example, the triangular maxitive kernel proposed in Loquin and Strauss (2008b) illustrated in Figure 3 is defined by

$$ \pi(u) = \begin{cases} 1 - |u| & \text{if } |u| \le 1, \\ 0 & \text{otherwise.} \end{cases} \qquad (9) $$

<figure><img src="image_5.png" /><figcaption>Figure 3. Triangular maxitive kernel with a bandwidth $\Delta = 1$.</figcaption></figure>

A maxitive kernel can be seen as a possibility distribution Loquin and Strauss (2008b), inducing two dual non-additive confidence measures, a possibility measure $\Pi_{\pi}$, and a necessity measure $N_{\pi}$ Dubois and Prade (1988); and De Cooman (1997) defined, $\forall A \subseteq \Omega$, by
$$ \Pi_{\pi}(A) = \sup_{x \in A} \pi(x) \text{ (possibility),} $$

$$ N_{\pi}(A) = 1 - \Pi_{\pi}(A^{c}) \text{ (necessity),} $$

with $A^c$ being the complementary set of A in $\Omega$.

A maxitive kernel $\pi$ is said to dominate a summative kernel $\kappa$ Loquin and Strauss (2008b), if the possibility measure $\Pi_{\pi}$ dominates the probability measure $P_{\kappa}$, i.e.,

$$ \forall A \subseteq \Omega, P_{\kappa}(A) \leq \Pi_{\pi}(A). $$

In that sense, a maxitive kernel defines the convex set of summative kernels it dominates. This set is denoted $M(\pi)$ and defined by

$$ M(\pi) = \{ \kappa \in K(\Omega) / \forall A \subseteq \Omega, N_{\pi}(A) \leq P_{\kappa}(A) \leq \Pi_{\pi}(A) \}. \quad (10) $$

The specificity of a maxitive kernel $\pi$, defined by its integral, i.e., $Sp(\pi) = \int_{\Omega} \pi(u)du$, is a measure of the information contained by the possibility measure associated to $\pi$ Loquin and Strauss (2008b). A maxitive kernel $\pi_1$ is at least as informative as another one $\pi_2$ if $\forall u \in \Omega, \pi_1(u) \leq \pi_2(u)$ Dubois (2006). In that case $\pi_1$ is at least as specific as $\pi_2$. It also characterizes the amount of summative kernels dominated by $\pi$ in the sense that, if $\forall u \in \Omega, \pi_1(u) \leq \pi_2(u)$, then $Sp(\pi_1) \leq Sp(\pi_2)$ and $M(\pi_1) \subseteq M(\pi_2)$. Moreover, if $\exists u_0 \in \Omega$, such that $\pi_1(u_0) < \pi_2(u_0)$, then $Sp(\pi_1) < Sp(\pi_2)$ and $M(\pi_1) \subset M(\pi_2)$.

In this context, Dubois et al. (2004) proved that the triangular maxitive kernel, defined by (Equation (9)), with a support $[-\Delta, +\Delta]$ is the most specific maxitive kernel that dominates all summative symmetric and unimodal kernels whose support belongs to $[-\Delta, +\Delta]$.

## 4.2. Choquet Integrals and Maxitive Expectation

We start with the definition of a capacity or **non-additive measure** which generalize the notion of additive measure, i.e., probability. The notion of capacity was introduced by Gustave Choquet in 1953 and has played an important role in game theory, fuzzy set theory, Dempster-Shafer theory and many others (see Denneberg 1994; Shafer 1976; Choquet 1953).

**Definition 5.** Let $(\Omega, \mathcal{P}(\Omega))$ be the measurable space. A capacity $v$ on $(\Omega, \mathcal{P}(\Omega))$ is a set function $v : \mathcal{P}(\Omega) \to [0, 1]$ satisfying:

* $v$ is normalized (i.e., $v(\emptyset) = 0$ and $v(\Omega) = 1$).

* $v$ is monotone (i.e., $\forall A, B \in \mathcal{P}(\Omega), A \subseteq B \Rightarrow v(A) \leq v(B)$).

One of the most important concepts closely related to additive measures is integration. It has a natural generalization to non-additive measure theory. Historically the first applied integral with respect to non-additive measures is the Choquet Integral formalized by Gustave Choquet.

**Definition 6.** Let $(\Omega, \mathcal{P}(\Omega))$ be the measurable space and $v : \mathcal{P}(\Omega) \to [0, 1]$ a capacity. Let $s$ be a bounded function of $L(\Omega)$. The continuous Choquet integral of $s$ w.r.t. $v$ is defined by

$$ C_v(s) = \int_{\Omega} sdv = \int_{-\infty}^{0} \left( v(\{u \in \Omega : s(u) \geq z\}) - 1 \right) dz + \int_{0}^{\infty} v(\{u \in \Omega : s(u) \geq z\}) dz, \quad (11) $$

where the integral on the right hand side is a Riemann integral.
**Definition 7.** Let $v$ be a capacity on $(\Omega, \mathcal{P}(\Omega))$ and $x = (x_1, \dots, x_n) \in \Omega^n$ be a discrete function $[1, \dots, n] \to \Omega$ of $n$ samples. The discrete Choquet integral of $x$ w.r.t. $v$ is defined by:

$$ \mathbb{C}_v(x) = \int_{\Omega} x d\nu = \sum_{i=1}^{n} (x_{\tau(i)} - x_{\tau(i-1)}) v(A_{\tau(i)}), \quad (12) $$

where $\tau$ is the permutation on $[1, \dots, n]$, such that $x_{\tau(1)} \le \dots \le x_{\tau(n)}$, $A_{\tau(i)} := \{\tau(i), \dots, \tau(n)\}$ and $x_{\tau(0)} = 0$ by convention.

Note that, a probability $P$ being a special case of (additive) capacity, the Choquet integral w.r.t. $P$ coincides with the classical expected value w.r.t. $P$.

Recall that, if $X$ is a random variable defined on a probability space $(\Omega, \mathcal{P}(\Omega), P)$, then the classical expected value of $X$ is given by

$$ \mathbb{E}_P(X) = \int_{\Omega} X dP = \int_{\Omega} X(u) dP(u), \quad (13) $$

where $P: \mathcal{P}(\Omega) \to [0, 1]$ is a probability measure on the measurable space $(\Omega, \mathcal{P}(\Omega))$ and the integral is a Lebesgue-Stieltjes integral.

The notion of expectation w.r.t. a summative kernel can be extended to a maxitive kernel, as defined in Loquin and Strauss (2008b).

**Definition 8.** Let $\pi$ be a maxitive kernel and let $s$ be a bounded function of $L(\Omega)$. The expectation of $s$ w.r.t $\pi$ is defined by:

$$ \bar{\mathbb{E}}_{\pi}(s) = [\mathbb{E}_{\pi}(s), \bar{\mathbb{E}}_{\pi}(s)] = [\mathbb{C}_{N_{\pi}}(s), \mathbb{C}_{\Pi_{\pi}}(s)], \quad (14) $$

where $\Pi_{\pi}$ (resp. $N_{\pi}$) is the possibility (resp. necessity) measure induced by the maxitive kernel $\pi$ and $\mathbb{C}$ is the Choquet integral.

An important property, that will be used in the sequel, is that the maxitive interval-valued expectation of $s$ w.r.t. $\pi$ is the convex set of all the summative precise-valued expectations w.r.t. all the summative kernels dominated by $\pi$ Loquin and Strauss (2008b).

**Property 1.** Let $\pi$ be a maxitive kernel and let $\mathcal{M}(\pi)$ be the set of the summative kernels it dominates (see Equation (10)). Let $s$ be a bounded function of $L(\Omega)$. We have

$$ \forall g \in \bar{\mathbb{E}}_{\pi}(s), \exists \kappa \in \mathcal{M}(\pi): \mathbb{E}_{\kappa}(s) = g, \quad (15) $$

and

$$ \forall \kappa \in \mathcal{M}(\pi), \mathbb{E}_{\kappa}(s) \in \bar{\mathbb{E}}_{\pi}(s), \quad (16) $$

where $\mathbb{E}_{\kappa}(s)$ (resp. $\bar{\mathbb{E}}_{\pi}(s)$) is the summative (resp. maxitive) expectation of $s$ w.r.t. $\kappa$ (resp. $\pi$).

### 4.3. Maxitive Kernel Cumulative and VaR Estimator

As mentioned in Loquin and Strauss (2008a) (Theorem 4), the summative kernel cumulative estimator $\hat{F}_{\kappa}^n$, defined by (Equation (3)), in each point $x \in \Omega$, can be written as the summative expectation of the empirical CDF $F_n$, defined by (Equation (4)), in a probabilistic neighborhood defined by the summative kernel $\kappa$ translated in $x$

$$ \hat{F}_{\kappa}^{n}(x) = \mathbb{E}_{\kappa^{x}}(F_{n}). \quad (17) $$

Based on this reformulation, Loquin and Strauss (2008a) propose an extension of the summative kernel estimator of the CDF. This extension called **maxitive kernel estimator of the CDF** or **interval-valued estimation of the CDF**.
**Definition 9.** Let $F_n$ be the empirical CDF based on the sample $(x_1, \dots, x_n)$. Let $\pi$ be a maxitive kernel and let $\bar{\mathbb{E}}_{\pi}(\cdot)$ be the maxitive expectation w.r.t. $\pi$ (Expression (14)). The maxitive kernel cumulative estimator is defined, $\forall x \in \Omega$, by

$$ \bar{\mathbb{E}}_{\pi}^{n}(x) = [\underline{\mathbb{E}}_{\pi}^{n}(x), \bar{\mathbb{E}}_{\pi}^{n}(x)] = \bar{\mathbb{E}}_{\pi^{x}}(F_{n}) = [\underline{\mathbb{E}}_{\pi^{x}}(F_{n}), \bar{\mathbb{E}}_{\pi^{x}}(F_{n})]. \quad (18) $$

The computation of $\bar{\mathbb{E}}_{\pi}^{n}$ involves two Choquet integrals. This computation Loquin and Strauss (2008a) is given, for all $x \in \Omega$, by

$$ \underline{\mathbb{E}}_{\pi}^{n}(x) = \underline{\mathbb{E}}_{\pi^{x}}(F_{n}) = \mathbb{C}_{N_{\pi^{x}}}(F_{n}) = \frac{1}{n} \sum_{i=1}^{n} \left( (1 - \pi^{x}(x_i)) \mathbb{1}(x_i \le x) \right), \quad (19) $$

$$ \bar{\mathbb{E}}_{\pi}^{n}(x) = \bar{\mathbb{E}}_{\pi^{x}}(F_{n}) = \mathbb{C}_{\Pi_{\pi^{x}}}(F_{n}) = \frac{1}{n} \sum_{i=1}^{n} \left( \pi^{x}(x_i) \mathbb{1}(x \le x_i) + \mathbb{1}(x_i \le x) \right). \quad (20) $$

The estimation of the CDF is usually computed on $p$ regularly spaced points of $\Omega$. Let $\{y_j\}_{j \in \{1, \dots, p\}}$ be those $p$ points. The algorithm of compute of $\{\bar{\mathbb{E}}_{\pi}^{n}(y_j)\}_{j \in \{1, \dots, p\}}$, in each point $\{y_j\}_{j \in \{1, \dots, p\}}$ of $\Omega$, is given by Algorithm 1.

**Algorithm 1:** Computation of $\{[\underline{\mathbb{E}}_{\pi}^{n}(y_j), \bar{\mathbb{E}}_{\pi}^{n}(y_j)]\}_{j \in \{1, \dots, p\}}$.

**Data:** The considered location $\{y_j\}_{j \in \{1, \dots, p\}}$, the observations $\{x_i\}_{i \in \{1, \dots, n\}}$ and the maxitive kernel $\pi$.

**Result:** $\{[\underline{\mathbb{E}}_{\pi}^{n}(y_j), \bar{\mathbb{E}}_{\pi}^{n}(y_j)]\}_{j \in \{1, \dots, p\}}$.

1. begin

2. for $j = 1$ to $p$ do

3. $\bar{F}_{\pi}^{n}(y_j) \leftarrow 0, \underline{F}_{\pi}^{n}(y_j) \leftarrow 0$

4. for $i = 1$ to $n$ do

5. if $y_j \le x_i$ then

6. $\bar{F}_{\pi}^{n}(y_j) = \bar{F}_{\pi}^{n}(y_j) + 1$

7. else

8. $\pi_i = \pi^{y_j}(x_i)$

9. $\bar{F}_{\pi}^{n}(y_j) = \bar{F}_{\pi}^{n}(y_j) + \pi_i$

10. end

11. if $y_j > x_i$ then

12. $\underline{F}_{\pi}^{n}(y_j) = \underline{F}_{\pi}^{n}(y_j) + 1$

13. else

14. $\pi_i = \pi^{y_j}(x_i)$

15. $\underline{F}_{\pi}^{n}(y_j) = \underline{F}_{\pi}^{n}(y_j) + \pi_i$

16. end

17. end

18. $\bar{F}_{\pi}^{n}(y_j) = \frac{1}{n}\bar{F}_{\pi}^{n}(y_j)$

19. $\underline{F}_{\pi}^{n}(y_j) = 1 - \frac{1}{n}\underline{F}_{\pi}^{n}(y_j)$

20. end

21. end

The interval-values estimate $\bar{\mathbb{E}}_{\pi}^{n}$ based on a maxitive kernel $\pi$ (Expression (9)) is the most specific interval containing all summative estimate $\hat{F}_{\kappa}^{n}$ (Expression (6)) based on a summative kernel $\kappa$ dominated by $\pi$, i.e., such that $\kappa \in \mathcal{M}(\pi)$.
**Property 2.** Let $\pi$ be a maxitive kernel on $\Omega$ and $\mathcal{M}(\pi)$ the set of all the summative kernels dominated by $\pi$. We have

$$ \forall x \in \Omega, \forall \kappa \in \mathcal{M}(\pi), \hat{F}_{\kappa}^{n}(x) \in \bar{F}_{\pi}^{n}(x), $$

and

$$ \forall g \in \bar{F}_{\pi}^{n}(x), \exists \kappa \in \mathcal{M}(\pi), \hat{F}_{\kappa}^{n}(x) = g. $$

Notice that, due to Property 2, both $\bar{F}_{\pi}$ and $\hat{F}_{\pi}$ are also estimates of the sought after CDF. Clearly, $\bar{F}_{\pi}^{n}$ and $\hat{F}_{\pi}^{n}$ are two strictly increasing continuous functions on $\Omega$, and so they both have an inverse. Then, we can infer immediately an maxitive interval-valued estimation of the VaR $\overline{\text{VaR}}_{\alpha, \pi} = [\text{VaR}_{\alpha, \pi}, \overline{\text{VaR}}_{\alpha, \pi}]$:

$$ \overline{\text{VaR}}_{\alpha, \pi} = \inf \{ x, \bar{F}_{\pi}(x) = P(X \le x) \ge \alpha \} = \bar{F}_{\pi}^{-1}(\alpha), $$

$$ \overline{\text{VaR}}_{\alpha, \pi} = \inf \{ x, \bar{F}_{\pi}(x) = P(X \le x) \ge \alpha \} = \bar{F}_{\pi}^{-1}(\alpha), $$

that inherits from the properties of the CDF: $\overline{\text{VaR}}_{\alpha, \pi} = \{\overline{\text{VaR}}_{\alpha, \kappa} | \kappa \in \mathcal{M}\{\pi\}\}$.

The MATLAB program computing the interval-valued estimation of VaR is outlined in Appendix A to the present paper.

# 5. Experiment: Empirical Results

## 5.1. Data and Experimental Process

The data used in the present study consisted of daily closing prices collected between January 2010 and December 2016 for four stock indexes from developed countries: S&P500, DJI, Nikkei225, and CAC40. The daily close prices were converted to daily log-returns. We took the first differences of the natural logarithm of the daily prices. For an observed price $P_t$, the corresponding one-day log-return on day $t$ was defined by: $r_t = \ln(\frac{P_t}{P_{t-1}})$. Table 1 contains descriptive statistics for the sample return of the four considered stock indexes. We observed that the returns were negatively skewed and characterized by heavy tails since the kurtoses were significantly greater than 3. The series had a distribution with tails that were significantly fatter than those of normal distribution. The Jarque-Bera test also indicated that hypothesizing the four data to be normally distributed can be rejected with a very low significance level.

<table><caption>**Table 1.** Descriptive statistics of daily returns on 4 stock indexes of a financial institution observed between 1 January 2010 and 1 September 2016.</caption><thead><tr><th>Index</th><th>S&P500</th><th>Nikkei225</th><th>DJI</th><th>CAC40</th></tr></thead><tbody><tr><td>No. of observations</td><td>1702</td><td>1673</td><td>1702</td><td>1581</td></tr><tr><td>Min</td><td>-6.895</td><td>-11.15</td><td>-5.7</td><td>-5.634</td></tr><tr><td>Max</td><td>4.63</td><td>10.66</td><td>4.15</td><td>9.22</td></tr><tr><td>Mean</td><td>0.0379</td><td>0.027</td><td>0.032</td><td>0.0038</td></tr><tr><td>Standard deviation</td><td>0.991</td><td>1.467</td><td>0.92</td><td>1.379</td></tr><tr><td>Skewness</td><td>-0.441</td><td>-0.468</td><td>-0.384</td><td>-0.015</td></tr><tr><td>Kurtosis</td><td>7.08</td><td>9.77</td><td>6.416</td><td>5.85</td></tr><tr><td></td><td>1239.5</td><td>3262.7</td><td>869.91</td><td>535.17</td></tr><tr><td>Jarque-Bera</td><td>p-value</td><td>p-value</td><td>p-value</td><td>p-value</td></tr><tr><td></td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.000</td></tr></tbody></table>

Due to these properties, the normal distribution would be an inappropriate model for calculating the daily VaR—as mentioned in the introduction. We thus rather envisage a nonparametric asymmetric
approach based on kernel estimation of the CDF. However, as previously mentioned, estimating the VaR using a kernel based approach can be highly biased, since the choice of the kernel has an important impact on assessing the VaR estimate at a given probability level.

In this section, we propose to gauge this impact by estimating the VaR at a several probability levels with different types of kernels, whose bandwidths have been adapted to the available datasets. We confirmed that no kernel can be considered as optimal in this context and show that the new maxitive interval-valued nonparametric approach we propose leads to a more cautious behavior when applied to real data. In fact, the bounds of the interval-valued estimate of the VaR, for a given probability level, can be instrumental in applications such as reserve requirements for banks. We then show that this approach allows discarding some parametric approaches since they are not adapted to heavy tailed data like those we consider.

## 5.2. Nonparametric Estimate of the VaR: Comparing Summative and Maxitive Approaches

In this experiment, we have considered estimating the VaR for low probability levels—namely at levels lower than 0.1—for each of the four datasets. In the first experiment, we have focused on DJI and Nikkei225 indexes. The CDF(s) have been estimated by considering four of the most used summative kernels whose bandwidth have been adapted using the rule of thumb method: Epanechnikov, normal, biweight, and triweight (linear)—see Silvermann (1986) for their analytical expressions. Figure 4 illustrates the impact of the kernel choice in the VaR estimation. This impact is more marked with high volatility stocks (Nikkei225). For example, referring to Figure 4a, it can be seen that the VaR estimate based on the Epanechnikov kernel is absolutely greater than the one based on the normal kernel at the probability level of 1% while the opposite situation occurs at the probability level of 2.5%. Referring to each plots of Figure 4, it is obvious that no kernel can be seen as always providing the upper or lower valued CDF. Since the CDF curves intersect each other for different probability levels then the risk of bias in the VaR estimate is relatively great, especially when the number of observations in the data set is low—compare Figure 4b–d. In this context, no kernel function can be considered as more or less conservative than the others. Also, other indexes would have produced similar results.

This first experiment shows that the high dependence of the CDF estimate to the chosen kernel shape can highly impact and bias the VaR estimate. The goal of the second experiment aims to show that none of the four considered kernel can be considered as optimal in this context. This experiment needs a ground truth which is not available. A resampling methodology was carried out to estimate such a ground truth. This methodology consists of five steps:

* **Step 1:** Fit a GARCH model to the four stock returns data. The fitted model considered for the returns $r_t$ is an AR(1)$^1$-GARCH(1,1) given by:

$$r_t = \mu_t + \varepsilon_t, \\ \varepsilon_t = \sigma_t z_t, \quad \forall t = 1, \dots, n, \\ \sigma_t^2 = \eta_0 + \eta_1 \varepsilon_{t-1}^2 + \eta_2 \sigma_{t-1}^2,$$

where $r_t$ is the return value at time $t$, $z_t$ is a standard Normal random variable and $n$ is the sample size ($n = 200$ and $1500$). The conditional mean $\mu_t$, is assumed to follow an AR(1) model given by: $\mu_t = \phi_0 + \phi_1 r_{t-1}$. By definition $\varepsilon_t$ is serially uncorrelated with mean zero, but a time varying conditional variance equal to $\sigma_t^2$. The three positive parameters $\eta_0, \eta_1, \eta_2$ and the two parameters $\phi_0, \phi_1$ are, respectively, the parameters of the GARCH(1,1) and the AR(1) models. The ML method provides a systematic way to adjust the parameters of the model to give the best fit. Table 2 lists the fitted models for each daily stock returns data.

<footnote n="1">¹ Autoregressive of order 1.</footnote>
<table><caption>Table 2. Fitted GARCH models.</caption><thead><tr><th>Sample Size</th><th>Index</th><th>μ<sub>t</sub> ~ AR(1)</th><th>GARCH(1,1)</th></tr></thead><tbody><tr><td rowspan="4">200</td><td>S&P500</td><td>r<sub>t</sub> = 0.032 - 0.132r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.08 + 0.23ε<sup>2</sup><sub>t-1</sub> + 0.663σ<sup>2</sup><sub>t-1</sub></td></tr><tr><td>Nikkei225</td><td>r<sub>t</sub> = -0.05 - 0.08r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.179 + 0.12ε<sup>2</sup><sub>t-1</sub> + 0.827σ<sup>2</sup><sub>t-1</sub></td></tr><tr><td>DJI</td><td>r<sub>t</sub> = 0.033 - 0.116r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.08 + 0.25ε<sup>2</sup><sub>t-1</sub> + 0.63σ<sup>2</sup><sub>t-1</sub></td></tr><tr><td>CAC40</td><td>r<sub>t</sub> = -0.05 - 0.031r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.032 + 0.044ε<sup>2</sup><sub>t-1</sub> + 0.93σ<sup>2</sup><sub>t-1</sub></td></tr><tr><td rowspan="4">1500</td><td>S&P500</td><td>r<sub>t</sub> = 0.046 - 0.063r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.078 + 0.207ε<sup>2</sup><sub>t-1</sub> + 0.69σ<sup>2</sup><sub>t-1</sub></td></tr><tr><td>Nikkei225</td><td>r<sub>t</sub> = 0.068 - 0.0043r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.088 + 0.176ε<sup>2</sup><sub>t-1</sub> + 0.814σ<sup>2</sup><sub>t-1</sub></td></tr><tr><td>DJI</td><td>r<sub>t</sub> = 0.039 - 0.069r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.07 + 0.22ε<sup>2</sup><sub>t-1</sub> + 0.68σ<sup>2</sup><sub>t-1</sub></td></tr><tr><td>CAC40</td><td>r<sub>t</sub> = 0.025 - 0.027r<sub>t-1</sub> + ε<sub>t</sub></td><td>σ<sup>2</sup><sub>t</sub> = 0.03 + 0.076ε<sup>2</sup><sub>t-1</sub> + 0.906σ<sup>2</sup><sub>t-1</sub></td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Cumulative distribution function (CDF) Kernel smoothed estimators for the DJI daily returns with sample size 200 (a) and 650 (b); and for the Nikkei225 daily returns with sample size 200 (c) and 650 (d). The black curve is the Epanechnikov distribution; red curve is the Triweight distribution. The blue curve corresponds to the Biweight distribution and the green curve is the Normal distribution.</figcaption></figure>

* **Step 2:** Generate 1000 simulated samples from each returns data using the coefficients obtained from the above fitted models. The main reason behind proposing GARCH models to simulate
our real data is the dependence properties and the volatility phenomena of stock returns; see Engle (1982).

In a first step, we generate an i.i.d series, z_t, by the random generator in MATLAB (R2010a) and another series σ_t. The random numbers sampled were all assumed to be normally distributed with expectation zero and unit variance. Then, the innovations for the GARCH series have been obtained via the equation ε_t = σ_t z_t. Finally the generation of r_t from the AR(1) process is straightforward. Therefore we run the series for 1500 and 200 times in each of the 1000 simulated samples.

* **Step 3:** Calculate the VaR estimates for each sample data generated in step two using the 4 kernel functions. These VaRs are calculated at five chosen levels of probability (α = 1%, 2.5%, 5%, 7.5%, 10%) with 2 times horizons (n = 200 and n = 1500).

* **Step 4:** Assess the performance of each kernel function by comparing their VaRs with the empirical VaRs. The performance criterion we examine is the mean square error (MSE), i.e.:

$$ \text{MSE} = \mathbb{E} \left( \left( \hat{Q}_{(n)}(p) - Q(p) \right)^2 \right) \approx \frac{1}{m} \sum_{k=1}^{m} \left( \hat{Q}_{n}^{(k)}(p) - Q(p) \right)^2, \quad (21) $$

where $\hat{Q}_{(n)}(p)$ is a vector of quantiles obtained from the simulated samples.

* **Step 5:** Tabulate the results that lead us to some important conclusions.

The results of the MSE of each α (%) VaR estimates are reported in Tables 3 and 4. From Tables 3 and 4, for S&P500, the best kernel function that estimate the VaR at level of 10%, 7.5%, 5% and 2.5% is the Epanechnikov kernel, while the Normal kernel is the best at 1% probabilities level. For Nikkei225, it seems that the Normal perform better than the other kernel functions. For DJI, the Epanechnikov kernel is more accurate to estimate the VaR, while for CAC40 the Normal and is matched more accurately. From Table 4, when n = 200, it appears that the Biwieight and Epanechnikov kernel are more likely to estimate VaR more accurately at several probability levels. However, the performance of kernel functions rapidly declines as n and α get smaller.

<table><caption>Table 3. Mean square error results for several types of kernel functions and n = 1500.</caption><thead><tr><th>α</th><th>Index</th><th>Epanechnikov</th><th>Normal</th><th>Biweight</th><th>Triweight</th><th>Best</th></tr></thead><tbody><tr><td rowspan="4">10%</td><td>S&P500</td><td><strong>0.01285</strong></td><td>0.01287</td><td>0.01286</td><td>0.01286</td><td>Epanechnikov</td></tr><tr><td>Nikkei225</td><td>0.226</td><td><strong>0.223</strong></td><td>0.227</td><td>0.2265</td><td>Normal</td></tr><tr><td>DJI</td><td>0.00535</td><td><strong>0.00533</strong></td><td>0.00534</td><td>0.0535</td><td>Normal</td></tr><tr><td>CAC40</td><td>0.0256</td><td><strong>0.0255</strong></td><td>0.02551</td><td>0.02551</td><td>Normal</td></tr><tr><td rowspan="4">7.5%</td><td>S&P500</td><td><strong>0.03575</strong></td><td>0.03580</td><td>0.035764</td><td>0.035763</td><td>Epanechnikov</td></tr><tr><td>Nikkei225</td><td>0.0316</td><td><strong>0.0315</strong></td><td>0.0318</td><td>0.0318</td><td>Normal</td></tr><tr><td>DJI</td><td>0.01075</td><td>0.01075</td><td><strong>0.01074</strong></td><td>0.01075</td><td>Biweight</td></tr><tr><td>CAC40</td><td>0.0339</td><td><strong>0.034</strong></td><td>0.0341</td><td>0.0342</td><td>Normal</td></tr><tr><td rowspan="4">5%</td><td>S&P500</td><td><strong>0.07288</strong></td><td>0.07295</td><td>0.07289</td><td>0.07289</td><td>Epanechnikov</td></tr><tr><td>Nikkei225</td><td>0.0435</td><td><strong>0.042</strong></td><td>0.0435</td><td>0.0436</td><td>Normal</td></tr><tr><td>DJI</td><td><strong>0.03877</strong></td><td>0.03879</td><td>0.03878</td><td>0.03882</td><td>Epanechnikov</td></tr><tr><td>CAC40</td><td>0.078</td><td><strong>0.0784</strong></td><td>0.0782</td><td>0.0781</td><td>Normal</td></tr><tr><td rowspan="4">2.5%</td><td>S&P500</td><td><strong>0.01763</strong></td><td>0.01764</td><td>0.017635</td><td>0.017637</td><td>Epanechnikov</td></tr><tr><td>Nikkei225</td><td>0.0858</td><td><strong>0.0857</strong></td><td>0.0859</td><td>0.08578</td><td>Normal</td></tr><tr><td>DJI</td><td><strong>0.141</strong></td><td>0.1409</td><td>0.1408</td><td>0.1409</td><td>Epanechnikov</td></tr><tr><td>CAC40</td><td>0.132</td><td><strong>0.131</strong></td><td>0.133</td><td>0.134</td><td>Normal</td></tr><tr><td rowspan="4">1%</td><td>S&P500</td><td>0.05804</td><td><strong>0.05801</strong></td><td>0.05804</td><td>0.058017</td><td>Normal</td></tr><tr><td>Nikkei225</td><td><strong>0.246</strong></td><td>0.2467</td><td>0.2481</td><td>0.2482</td><td>Epanechnikov</td></tr><tr><td>DJI</td><td><strong>0.2404</strong></td><td>0.2405</td><td>0.2406</td><td>0.2407</td><td>Epanechnikov</td></tr><tr><td>CAC40</td><td>0.363</td><td>0.363</td><td>0.364</td><td><strong>0.362</strong></td><td>Triweight</td></tr></tbody></table>
<table><caption>**Table 4.** Mean square error results for several types of kernel functions and $n = 200$.</caption><thead><tr><th>α</th><th>Index</th><th>Epanechnikov</th><th>Normal</th><th>Biweight</th><th>Triweight</th><th>Best</th></tr></thead><tbody><tr><td rowspan="4">10%</td><td>S&P500</td><td>0.0337</td><td><strong>0.0335</strong></td><td>0.0336</td><td>0.0338</td><td>Normal</td></tr><tr><td>Nikkei225</td><td>0.766</td><td><strong>0.765</strong></td><td>0.767</td><td>0.766</td><td>Normal</td></tr><tr><td>DJI</td><td><strong>0.0441</strong></td><td>0.0443</td><td>0.0442</td><td>0.0442</td><td>Epanechnikov</td></tr><tr><td>CAC40</td><td>0.0347</td><td>0.0348</td><td><strong>0.0346</strong></td><td>0.0348</td><td>Biweight</td></tr><tr><td rowspan="4">7.5%</td><td>S&P500</td><td><strong>0.0632</strong></td><td>0.0634</td><td>0.0633</td><td>0.064</td><td>Epanechnikov</td></tr><tr><td>Nikkei225</td><td>0.938</td><td>0.937</td><td>0.939</td><td><strong>0.936</strong></td><td>Triweight</td></tr><tr><td>DJI</td><td>0.0614</td><td>0.0615</td><td><strong>0.0612</strong></td><td>0.0613</td><td>Biweight</td></tr><tr><td>CAC40</td><td><strong>0.0588</strong></td><td>0.0591</td><td>0.0592</td><td>0.0594</td><td>Epanechnikov</td></tr><tr><td rowspan="4">5%</td><td>S&P500</td><td>0.1169</td><td>0.1173</td><td><strong>0.1167</strong></td><td>0.1183</td><td>Biweight</td></tr><tr><td>Nikkei225</td><td>1.131</td><td>1.129</td><td>1.131</td><td><strong>1.127</strong></td><td>Triweight</td></tr><tr><td>DJI</td><td>0.0975</td><td>0.0977</td><td><strong>0.0974</strong></td><td>0.0977</td><td>Biweight</td></tr><tr><td>CAC40</td><td><strong>0.148</strong></td><td>0.149</td><td>0.147</td><td>0.0148</td><td>Epanechnikov</td></tr><tr><td rowspan="4">2.5%</td><td>S&P500</td><td>0.28</td><td>0.284</td><td><strong>0.281</strong></td><td>0.283</td><td>Biweight</td></tr><tr><td>Nikkei225</td><td>1.453</td><td>1.451</td><td>1.453</td><td><strong>1.448</strong></td><td>Triweight</td></tr><tr><td>DJI</td><td><strong>0.252</strong></td><td>0.255</td><td>0.253</td><td>0.254</td><td>Epanechnikov</td></tr><tr><td>CAC40</td><td>0.2905</td><td>0.2888</td><td><strong>0.287</strong></td><td>0.2889</td><td>Biweight</td></tr><tr><td rowspan="4">1%</td><td>S&P500</td><td><strong>0.869</strong></td><td>0.821</td><td>0.87</td><td>0.881</td><td>Epanechnikov</td></tr><tr><td>Nikkei225</td><td><strong>1.977</strong></td><td>1.995</td><td>1.981</td><td>1.980</td><td>Epanechnikov</td></tr><tr><td>DJI</td><td>0.702</td><td><strong>0.691</strong></td><td>0.696</td><td>0.693</td><td>Normal</td></tr><tr><td>CAC40</td><td>0.794</td><td>0.791</td><td><strong>0.788</strong></td><td>0.793</td><td>Biweight</td></tr></tbody></table>

The simulations results showed that, for large sample sizes (here $n = 1500$), all kernel functions are consistent estimators, i.e., MSE values are close to 0. On the other hand, MSE values are larger for the smaller sample size (here $n = 200$). This demonstrates that choosing a particular kernel, when the sample size is low, is risky. Therefore a coherent interval-valued of VaR as we propose is likely to provide a more careful decision in this context.

### 5.3. *Interval Estimation of Value-at-Risk and Some Numerical Comparisons*

Here we apply the maxitive kernel estimator presented in Section 4.3 to obtain the lower and upper bound for VaR corresponding to each four stock indexes for three probability levels ($\alpha = 10\%, 5\%, 1\%$). To obtain the optimal bounds for the VaR, the bandwidth has been, once again, chosen using the most popular methods such as biased cross-validation method and plug-in method presented in Section 3.3. As shown in Table 5, the maximum available bandwidth is taken to insure that all kernel estimators are inside the interval given by the maxitive kernel estimator. In order to examine the performance of the maxitive kernel estimation method, we divide the data into three samples: the first sample corresponds to data with 6 year time horizons, the second sample is chosen with time horizons of 3 years and the third correspond to one year time horizon. In a first step, after choosing the best-fit distribution from the family of GHYP distribution by using the AIC<fnref n="2" /> criteria, we compute the VaR across the candidate distribution. Next, we compute the VaR confidence interval based on these three distributions using the bias-corrected and accelerated (BCa) bootstrap method. Finally, we compare our proposed maxitive interval-valued with those based on the bootstrap technique under different samples sizes. So far, the bootstrap confidence interval of VaR for the three distributions: GHYP, Normal, and HS are presented together in Tables 6–8. With these punctual VaRs at hand for comparison purposes, we evaluate the performance of our approach. In this section, we have chosen to illustrate our results and to discuss the benefit of our approach in two ways. In Tables 6–8, we show the

<footnote n="2">² Akaike Information Criterion. AIC = $-2 \ln(L(GHYP)) + 2\zeta$, where $\zeta$ is the number of estimated parameters and $L(GHYP)$ is the likelihood of the GHYP model.</footnote>
explicit results for the minimum and maximum value, as well as the width of the interval estimation of VaR for several probability levels between 1% and 10%. Figure 5 shows the PP-plots, i.e., the $\hat{F}(x_i)$ of all models against the $F(x_i)$. Based on the numerical results we can formulate several conclusions: Evidently, the lower and upper bounds increase with the probability level $\alpha$. However, it is important to note that the intervals estimation of VaRs for the two indexes Nikkei225 and CAC40 are larger than those of DJI and S&P500 indexes. Note also that the influence of the volatility stock market is much more important than the influence of the sample size. Also, the widths of maxitive VaR intervals are rather tight than those derived from the GHYP, HS, and normal distributions. For example, for the short time horizon 1 year with the high volatility stock (Nikkei225) and at the 1% of probability level the width of the interval-valued of VaR is 1.073 while the widths of the BCa (99%) confidence interval based on HS, normal, and GHYP distributions are respectively 4.236, 1.925, and 4.246 respectively (See Table 6). Furthermore, we can remark that the maxitive interval-valued estimation of the VaR is on the left side of the normal VaR and this shows the ability of our approach to model very dangerous financial risks while the normal distribution is not consistent with tail-thickness and right tail risk. This results indicate that our approach is more accurate and informative especially for the smaller sample size.

<table><caption>Table 5. Optimal bandwidth estimation for the four stocks return data.</caption><thead><tr><th>Return Data</th><th>Time Horizon</th><th>Bandwidth Selection Method</th><th>Epanechnikov</th><th>Normal</th><th>Biweight</th><th>Triweight</th><th>Maxitive Kernel Bandwidth</th></tr></thead><tbody><tr><td rowspan="4">S&P500</td><td rowspan="2">6 years</td><td>Plug-in</td><td>0.189</td><td>0.084</td><td>0.245</td><td>0.250</td><td rowspan="2">0.250</td></tr><tr><td>Cross validation</td><td>0.174</td><td>0.060</td><td>0.190</td><td>0.211</td></tr><tr><td rowspan="2">3 years</td><td>Plug-in</td><td>0.228</td><td>0.101</td><td>0.269</td><td>0.294</td><td rowspan="2">0.294</td></tr><tr><td>Cross validation</td><td>0.197</td><td>0.118</td><td>0.211</td><td>0.267</td></tr><tr><td rowspan="4">Nikkei225</td><td rowspan="2">1 year</td><td>Plug-in</td><td>0.318</td><td>0.141</td><td>0.376</td><td>0.405</td><td rowspan="2">0.405</td></tr><tr><td>Cross validation</td><td>0.277</td><td>0.092</td><td>0.338</td><td>0.397</td></tr><tr><td rowspan="2">6 years</td><td>Plug-in</td><td>0.332</td><td>0.147</td><td>0.405</td><td>0.440</td><td rowspan="2">0.440</td></tr><tr><td>Cross validation</td><td>0.330</td><td>0.110</td><td>0.331</td><td>0.412</td></tr><tr><td rowspan="4">DJI</td><td rowspan="2">3 years</td><td>Plug-in</td><td>0.406</td><td>0.180</td><td>0.480</td><td>0.540</td><td rowspan="2">0.540</td></tr><tr><td>Cross validation</td><td>0.307</td><td>0.101</td><td>0.390</td><td>0.430</td></tr><tr><td rowspan="2">1 year</td><td>Plug-in</td><td>0.840</td><td>0.374</td><td>0.996</td><td>1.118</td><td rowspan="2">1.118</td></tr><tr><td>Cross validation</td><td>0.841</td><td>0.382</td><td>0.994</td><td>1.098</td></tr><tr><td rowspan="4">CAC40</td><td rowspan="2">6 years</td><td>Plug-in</td><td>0.173</td><td>0.077</td><td>0.205</td><td>0.224</td><td rowspan="2">0.224</td></tr><tr><td>Cross validation</td><td>0.150</td><td>0.050</td><td>0.180</td><td>0.208</td></tr><tr><td rowspan="2">3 years</td><td>Plug-in</td><td>0.228</td><td>0.101</td><td>0.270</td><td>0.294</td><td rowspan="2">0.294</td></tr><tr><td>Cross validation</td><td>0.189</td><td>0.113</td><td>0.210</td><td>0.223</td></tr><tr><td rowspan="4">CAC40</td><td rowspan="2">1 year</td><td>Plug-in</td><td>0.360</td><td>0.160</td><td>0.425</td><td>0.468</td><td rowspan="2">0.468</td></tr><tr><td>Cross validation</td><td>0.326</td><td>0.148</td><td>0.381</td><td>0.422</td></tr><tr><td rowspan="2">6 years</td><td>Plug-in</td><td>0.323</td><td>0.143</td><td>0.382</td><td>0.431</td><td rowspan="2">0.446</td></tr><tr><td>Cross validation</td><td>0.374</td><td>0.080</td><td>0.390</td><td>0.446</td></tr><tr><td rowspan="4">CAC40</td><td rowspan="2">3 years</td><td>Plug-in</td><td>0.423</td><td>0.188</td><td>0.498</td><td>0.564</td><td rowspan="2">0.564</td></tr><tr><td>Cross validation</td><td>0.433</td><td>0.144</td><td>0.501</td><td>0.544</td></tr><tr><td rowspan="2">1 year</td><td>Plug-in</td><td>0.93</td><td>0.416</td><td>1.107</td><td>1.250</td><td rowspan="2">1.250</td></tr><tr><td>Cross validation</td><td>0.934</td><td>0.411</td><td>1.150</td><td>1.230</td></tr></tbody></table>
<table><caption>**Table 6.** A comparison between the bound estimates of the daily VaR based on the maxitive kernel estimation method and bootstrap confidence intervals of Value-at-Risk (VaRs) based on three different distributions: Historical simulation method (HS), Normal, and Generalized HYPerbolic (GHYP) distributions. These VaR methods are applied to four daily stock returns at several probability levels a% and time horizon of 1 year.</caption><thead><tr><th rowspan="2">Returns</th><th rowspan="2">α</th><th colspan="3">Interval-Valued of VaR</th><th rowspan="2">pdf</th><th rowspan="2">Estimated VaR</th><th colspan="6">Bootstrap Confidence Intervals of VaRs</th></tr><tr><th>L.B.</th><th>U.B.</th><th>Width</th><th>BCa (90%)</th><th>BCa (95%)</th><th>BCa (99%)</th></tr></thead><tbody><tr><td rowspan="6">S&P500</td><td rowspan="3">10%</td><td rowspan="3">-1.092</td><td rowspan="3">-0.751</td><td rowspan="3">0.341</td><td>HS</td><td>-0.933</td><td>-1.089</td><td>-0.773</td><td>0.316</td><td>-1.126</td><td>-0.758</td><td>0.368</td><td>-1.162</td><td>-0.665</td><td>0.497</td></tr><tr><td>Normal</td><td>-1.117</td><td>-1.303</td><td>-0.933</td><td>0.370</td><td>-1.336</td><td>-0.893</td><td>0.443</td><td>-1.379</td><td>-0.825</td><td>0.554</td></tr><tr><td>GHYP</td><td>-0.839</td><td>-0.993</td><td>-0.724</td><td>0.269</td><td>-1.016</td><td>-0.694</td><td>0.322</td><td>-1.094</td><td>-0.673</td><td>0.421</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-1.605</td><td rowspan="3">-1.152</td><td rowspan="3">0.453</td><td>HS</td><td>-1.325</td><td>-1.552</td><td>-1.146</td><td>0.406</td><td>-1.581</td><td>-1.109</td><td>0.472</td><td>-1.632</td><td>-1.070</td><td>0.562</td></tr><tr><td>Normal</td><td>-1.440</td><td>-1.716</td><td>-1.256</td><td>0.460</td><td>-1.757</td><td>-1.234</td><td>0.523</td><td>-1.822</td><td>-1.131</td><td>0.691</td></tr><tr><td>GHYP</td><td>-1.264</td><td>-1.486</td><td>-1.053</td><td>0.433</td><td>-1.555</td><td>-1.015</td><td>0.540</td><td>-1.595</td><td>-0.900</td><td>0.695</td></tr><tr><td rowspan="6">Nikkei225</td><td rowspan="3">1%</td><td rowspan="3">-3.240</td><td rowspan="3">-2.117</td><td rowspan="3">1.123</td><td>HS</td><td>-2.196</td><td>-2.548</td><td>-1.938</td><td>0.610</td><td>-2.631</td><td>-1.876</td><td>0.755</td><td>-2.773</td><td>-1.736</td><td>1.037</td></tr><tr><td>Normal</td><td>-2.052</td><td>-2.399</td><td>-1.818</td><td>0.581</td><td>-2.456</td><td>-1.788</td><td>0.668</td><td>-2.469</td><td>-1.678</td><td>0.791</td></tr><tr><td>GHYP</td><td>-2.281</td><td>-2.589</td><td>-1.997</td><td>0.592</td><td>-2.683</td><td>-1.963</td><td>0.720</td><td>-2.701</td><td>-1.851</td><td>0.850</td></tr><tr><td rowspan="6">DJI</td><td rowspan="3">10%</td><td rowspan="3">-5.013</td><td rowspan="3">-3.236</td><td rowspan="3">1.777</td><td>HS</td><td>-4.331</td><td>-5.055</td><td>-3.713</td><td>1.342</td><td>-5.176</td><td>-3.590</td><td>1.586</td><td>-5.906</td><td>-3.419</td><td>2.487</td></tr><tr><td>Normal</td><td>-2.401</td><td>-2.805</td><td>-2.100</td><td>0.705</td><td>-2.831</td><td>-1.906</td><td>0.925</td><td>-2.970</td><td>-1.768</td><td>1.202</td></tr><tr><td>GHYP</td><td>-3.889</td><td>-4.588</td><td>-3.434</td><td>1.154</td><td>-4.772</td><td>-3.376</td><td>1.396</td><td>-4.875</td><td>-3.144</td><td>1.731</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-6.760</td><td rowspan="3">-5.010</td><td rowspan="3">1.750</td><td>HS</td><td>-5.797</td><td>-6.852</td><td>-5.196</td><td>1.656</td><td>-7.054</td><td>-5.121</td><td>1.933</td><td>-7.054</td><td>-4.789</td><td>2.265</td></tr><tr><td>Normal</td><td>-3.061</td><td>-3.625</td><td>-2.677</td><td>0.948</td><td>-3.687</td><td>-2.601</td><td>1.086</td><td>-3.806</td><td>-2.431</td><td>1.375</td></tr><tr><td>GHYP</td><td>-5.567</td><td>-6.458</td><td>-4.806</td><td>1.652</td><td>-6.886</td><td>-4.743</td><td>2.143</td><td>-7.009</td><td>-4.615</td><td>2.394</td></tr><tr><td rowspan="6">CAC40</td><td rowspan="3">1%</td><td rowspan="3">-10.170</td><td rowspan="3">-9.097</td><td rowspan="3">1.073</td><td>HS</td><td>-9.259</td><td>-10.31</td><td>-7.788</td><td>2.522</td><td>-10.710</td><td>-7.639</td><td>3.071</td><td>-11.320</td><td>-7.084</td><td>4.236</td></tr><tr><td>Normal</td><td>-4.301</td><td>-5.233</td><td>-3.774</td><td>1.459</td><td>-5.409</td><td>-3.743</td><td>1.666</td><td>-5.512</td><td>-3.587</td><td>1.925</td></tr><tr><td>GHYP</td><td>-9.907</td><td>-11.53</td><td>-8.766</td><td>2.764</td><td>-11.760</td><td>-8.549</td><td>3.211</td><td>-11.970</td><td>-7.724</td><td>4.246</td></tr><tr><td rowspan="6">10%</td><td rowspan="3">-1.099</td><td rowspan="3">-0.766</td><td rowspan="3">0.333</td><td>HS</td><td>-0.897</td><td>-1.081</td><td>-0.752</td><td>0.329</td><td>-1.157</td><td>-0.708</td><td>0.449</td><td>-1.184</td><td>-0.648</td><td>0.536</td></tr><tr><td>Normal</td><td>-1.081</td><td>-1.276</td><td>-0.938</td><td>0.338</td><td>-1.362</td><td>-0.915</td><td>0.447</td><td>-1.376</td><td>-0.858</td><td>0.518</td></tr><tr><td>GHYP</td><td>-0.838</td><td>-1.016</td><td>-0.694</td><td>0.322</td><td>-1.040</td><td>-0.673</td><td>0.367</td><td>-1.061</td><td>-0.623</td><td>0.438</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-1.540</td><td rowspan="3">-1.094</td><td rowspan="3">0.446</td><td>HS</td><td>-1.262</td><td>-1.550</td><td>-1.097</td><td>0.453</td><td>-1.576</td><td>-0.999</td><td>0.577</td><td>-1.138</td><td>-0.669</td><td>0.469</td></tr><tr><td>Normal</td><td>-1.394</td><td>-1.597</td><td>-1.181</td><td>0.416</td><td>-1.642</td><td>-1.142</td><td>0.500</td><td>-1.741</td><td>-1.078</td><td>0.663</td></tr><tr><td>GHYP</td><td>-1.237</td><td>-1.459</td><td>-1.085</td><td>0.374</td><td>-1.470</td><td>-1.039</td><td>0.431</td><td>-1.496</td><td>-0.989</td><td>0.507</td></tr><tr><td rowspan="6">1%</td><td rowspan="3">-2.949</td><td rowspan="3">-1.837</td><td rowspan="3">1.112</td><td>HS</td><td>-2.006</td><td>2.303</td><td>-1.727</td><td>-4.030</td><td>-2.349</td><td>-1.669</td><td>0.680</td><td>-2.439</td><td>-1.568</td><td>0.871</td></tr><tr><td>Normal</td><td>-1.982</td><td>-2.295</td><td>-1.760</td><td>0.535</td><td>-2.315</td><td>-1.690</td><td>0.625</td><td>2.396</td><td>-1.588</td><td>-3.984</td></tr><tr><td>GHYP</td><td>-2.166</td><td>-2.458</td><td>-1.919</td><td>0.539</td><td>-2.532</td><td>-1.863</td><td>0.669</td><td>-2.646</td><td>-1.715</td><td>0.931</td></tr><tr><td rowspan="6">10%</td><td rowspan="3">-3.555</td><td rowspan="3">-2.154</td><td rowspan="3">1.401</td><td>HS</td><td>-2.670</td><td>-3.075</td><td>-2.246</td><td>0.829</td><td>-3.160</td><td>-2.198</td><td>0.962</td><td>-3.184</td><td>-2.074</td><td>1.110</td></tr><tr><td>Normal</td><td>-1.920</td><td>-2.184</td><td>-1.616</td><td>0.568</td><td>-2.305</td><td>-1.561</td><td>0.744</td><td>-2.388</td><td>-1.515</td><td>0.873</td></tr><tr><td>GHYP</td><td>-2.702</td><td>-3.258</td><td>-2.099</td><td>1.159</td><td>-3.257</td><td>-2.234</td><td>1.023</td><td>-3.258</td><td>-2.099</td><td>1.159</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-4.202</td><td rowspan="3">-3.184</td><td rowspan="3">1.018</td><td>HS</td><td>-3.605</td><td>-4.142</td><td>-3.102</td><td>1.040</td><td>-4.156</td><td>-3.045</td><td>1.111</td><td>-4.415</td><td>-2.940</td><td>1.475</td></tr><tr><td>Normal</td><td>-2.446</td><td>-2.700</td><td>-1.991</td><td>0.709</td><td>-2.818</td><td>-1.961</td><td>0.857</td><td>-2.978</td><td>-1.961</td><td>1.017</td></tr><tr><td>GHYP</td><td>-3.576</td><td>-4.084</td><td>-3.106</td><td>0.978</td><td>-4.150</td><td>-2.936</td><td>1.214</td><td>-4.433</td><td>-2.608</td><td>1.825</td></tr><tr><td rowspan="6">1%</td><td rowspan="3">-6.051</td><td rowspan="3">-4.771</td><td rowspan="3">1.280</td><td>HS</td><td>-5.100</td><td>-5.759</td><td>-4.618</td><td>1.141</td><td>-6.129</td><td>-4.477</td><td>1.652</td><td>-6.125</td><td>-4.477</td><td>1.648</td></tr><tr><td>Normal</td><td>-3.433</td><td>-3.836</td><td>-2.979</td><td>0.857</td><td>-3.903</td><td>-2.863</td><td>1.040</td><td>-4.124</td><td>-2.757</td><td>1.367</td></tr><tr><td>GHYP</td><td>-5.463</td><td>-6.391</td><td>-4.860</td><td>1.531</td><td>-6.710</td><td>-4.821</td><td>1.889</td><td>-6.717</td><td>-4.411</td><td>2.306</td></tr></tbody></table>
<table><caption>**Table 7.** A comparison between the bound estimates of the daily VaR based on the maxitive kernel estimation method and bootstrap confidence intervals of VaRs based on three different distributions: Historical simulation method (HS), Normal, and GHYP distributions. These VaR methods are applied to four daily stock returns at several probability levels a% and time horizon of 3 years.</caption><thead><tr><th rowspan="3">Returns</th><th rowspan="3">α</th><th colspan="3">Interval-Valued of VaR</th><th rowspan="3">pdf</th><th rowspan="3">Estimated VaR</th><th colspan="9">Bootstrap Confidence Intervals of VaRs</th></tr><tr><th rowspan="2">L.B.</th><th rowspan="2">U.B.</th><th rowspan="2">Width</th><th colspan="3">BCa (90%)</th><th colspan="3">BCa (95%)</th><th colspan="3">BCa (99%)</th></tr><tr><th>L.B.</th><th>U.B.</th><th>Width</th><th>L.B.</th><th>U.B.</th><th>Width</th><th>L.B.</th><th>U.B.</th><th>Width</th></tr></thead><tbody><tr><td rowspan="6">S&P500</td><td rowspan="3">10%</td><td rowspan="3">-0.996</td><td rowspan="3">-0.732</td><td rowspan="3">0.264</td><td>HS</td><td>0.860</td><td>-0.941</td><td>-0.765</td><td>0.176</td><td>-0.958</td><td>-0.753</td><td>0.205</td><td>-1.001</td><td>-0.737</td><td>0.264</td></tr><tr><td>Normal</td><td>-1.088</td><td>-1.199</td><td>-1.002</td><td>0.197</td><td>-1.221</td><td>-0.986</td><td>0.235</td><td>-1.250</td><td>-0.956</td><td>0.294</td></tr><tr><td>GHYP</td><td>-0.834</td><td>-0.933</td><td>-0.763</td><td>0.170</td><td>-0.946</td><td>-0.746</td><td>0.200</td><td>-0.969</td><td>-0.709</td><td>0.260</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-1.394</td><td rowspan="3">-1.130</td><td rowspan="3">0.264</td><td>HS</td><td>-1.253</td><td>-1.382</td><td>-1.130</td><td>0.252</td><td>-1.363</td><td>-1.145</td><td>0.218</td><td>-1.429</td><td>-1.085</td><td>0.344</td></tr><tr><td>Normal</td><td>-1.402</td><td>-1.527</td><td>-1.301</td><td>0.226</td><td>-1.548</td><td>-1.286</td><td>0.262</td><td>-1.598</td><td>-1.240</td><td>0.358</td></tr><tr><td>GHYP</td><td>-1.207</td><td>-1.316</td><td>-1.112</td><td>0.204</td><td>-1.335</td><td>-1.090</td><td>0.245</td><td>-1.376</td><td>-1.040</td><td>0.336</td></tr><tr><td rowspan="6">Nikkei225</td><td rowspan="3">1%</td><td rowspan="3">-2.230</td><td rowspan="3">-1.979</td><td rowspan="3">0.251</td><td>HS</td><td>-2.059</td><td>-2.208</td><td>-1.922</td><td>0.286</td><td>-2.235</td><td>-1.897</td><td>0.338</td><td>-2.308</td><td>-1.851</td><td>0.457</td></tr><tr><td>Normal</td><td>-1.993</td><td>-2.163</td><td>-1.854</td><td>0.309</td><td>-2.192</td><td>-1.829</td><td>0.363</td><td>-2.268</td><td>-1.788</td><td>0.480</td></tr><tr><td>GHYP</td><td>-2.067</td><td>-2.246</td><td>-1.936</td><td>0.310</td><td>-2.275</td><td>-1.909</td><td>0.366</td><td>-2.313</td><td>-1.832</td><td>0.481</td></tr><tr><td rowspan="6">DJI</td><td rowspan="3">10%</td><td rowspan="3">-3.024</td><td rowspan="3">-2.283</td><td rowspan="3">0.741</td><td>HS</td><td>-2.602</td><td>-2.937</td><td>-2.386</td><td>0.551</td><td>-3.021</td><td>-2.341</td><td>0.680</td><td>-3.119</td><td>-2.231</td><td>0.888</td></tr><tr><td>Normal</td><td>-1.976</td><td>-2.176</td><td>-1.756</td><td>0.420</td><td>-2.235</td><td>-1.721</td><td>0.514</td><td>-2.350</td><td>-1.612</td><td>0.738</td></tr><tr><td>GHYP</td><td>-2.517</td><td>-2.773</td><td>-2.254</td><td>0.519</td><td>-2.816</td><td>-2.209</td><td>0.607</td><td>-2.929</td><td>-2.090</td><td>0.839</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-4.324</td><td rowspan="3">-3.616</td><td rowspan="3">0.708</td><td>HS</td><td>-3.911</td><td>-4.359</td><td>-3.566</td><td>0.793</td><td>-4.486</td><td>-3.503</td><td>0.983</td><td>-4.647</td><td>-3.346</td><td>1.301</td></tr><tr><td>Normal</td><td>-2.539</td><td>-2.851</td><td>-2.298</td><td>0.553</td><td>-2.885</td><td>-2.246</td><td>0.639</td><td>-3.028</td><td>-2.170</td><td>0.858</td></tr><tr><td>GHYP</td><td>-3.832</td><td>-4.338</td><td>-3.498</td><td>0.840</td><td>-4.384</td><td>-3.423</td><td>0.961</td><td>-4.441</td><td>-3.292</td><td>1.149</td></tr><tr><td rowspan="6">CAC40</td><td rowspan="3">1%</td><td rowspan="3">-6.602</td><td rowspan="3">-6.052</td><td rowspan="3">0.550</td><td>HS</td><td>-6.273</td><td>-6.980</td><td>-5.697</td><td>1.283</td><td>-7.053</td><td>-5.599</td><td>1.454</td><td>-7.308</td><td>-5.425</td><td>1.883</td></tr><tr><td>Normal</td><td>-3.595</td><td>-4.012</td><td>-3.263</td><td>0.749</td><td>-4.077</td><td>-3.182</td><td>0.895</td><td>-4.279</td><td>-3.081</td><td>1.198</td></tr><tr><td>GHYP</td><td>-7.421</td><td>-8.262</td><td>-6.738</td><td>1.524</td><td>-8.432</td><td>-6.632</td><td>1.800</td><td>-8.695</td><td>-6.273</td><td>2.422</td></tr><tr><td rowspan="6">10%</td><td rowspan="3">-0.976</td><td rowspan="3">-0.743</td><td rowspan="3">0.233</td><td>HS</td><td>-0.851</td><td>-0.927</td><td>-0.763</td><td>0.164</td><td>-0.945</td><td>-0.755</td><td>0.190</td><td>-0.986</td><td>-0.742</td><td>0.244</td></tr><tr><td>Normal</td><td>-1.072</td><td>-1.168</td><td>-0.968</td><td>0.200</td><td>-1.184</td><td>-0.953</td><td>0.231</td><td>-1.202</td><td>-0.926</td><td>0.276</td></tr><tr><td>GHYP</td><td>-0.830</td><td>-0.906</td><td>-0.758</td><td>0.148</td><td>-0.930</td><td>-0.748</td><td>0.182</td><td>-0.981</td><td>-0.722</td><td>0.259</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-1.373</td><td rowspan="3">-1.161</td><td rowspan="3">0.212</td><td>HS</td><td>-1.245</td><td>-1.349</td><td>-1.148</td><td>0.201</td><td>-1.377</td><td>-1.136</td><td>0.241</td><td>-1.410</td><td>-1.088</td><td>0.322</td></tr><tr><td>Normal</td><td>-1.379</td><td>-1.489</td><td>-1.286</td><td>0.203</td><td>-1.506</td><td>-1.263</td><td>0.243</td><td>-1.564</td><td>-1.230</td><td>0.334</td></tr><tr><td>GHYP</td><td>-1.193</td><td>-1.300</td><td>-1.102</td><td>0.198</td><td>-1.319</td><td>-1.086</td><td>0.233</td><td>-1.375</td><td>-1.062</td><td>0.313</td></tr><tr><td rowspan="6">1%</td><td rowspan="3">-1.975</td><td rowspan="3">-1.752</td><td rowspan="3">0.223</td><td>HS</td><td>-1.816</td><td>-1.959</td><td>-1.694</td><td>0.265</td><td>-1.972</td><td>-1.664</td><td>0.308</td><td>-2.021</td><td>-1.622</td><td>0.399</td></tr><tr><td>Normal</td><td>-1.956</td><td>-2.081</td><td>-1.815</td><td>0.266</td><td>-2.108</td><td>-1.790</td><td>0.318</td><td>-2.173</td><td>-1.726</td><td>0.447</td></tr><tr><td>GHYP</td><td>-2.025</td><td>-2.169</td><td>-1.891</td><td>0.278</td><td>-2.197</td><td>-1.856</td><td>0.341</td><td>-2.256</td><td>-1.830</td><td>0.426</td></tr><tr><td rowspan="6">10%</td><td rowspan="3">-2.311</td><td rowspan="3">-1.658</td><td rowspan="3">0.653</td><td>HS</td><td>-1.963</td><td>-2.152</td><td>-1.832</td><td>0.320</td><td>-2.148</td><td>-1.771</td><td>0.377</td><td>-2.230</td><td>-1.742</td><td>0.488</td></tr><tr><td>Normal</td><td>-1.630</td><td>-1.781</td><td>-1.507</td><td>0.274</td><td>-1.806</td><td>-1.482</td><td>0.324</td><td>-1.836</td><td>-1.428</td><td>0.408</td></tr><tr><td>GHYP</td><td>-1.968</td><td>-2.137</td><td>-1.837</td><td>0.300</td><td>-2.161</td><td>-1.804</td><td>0.357</td><td>-2.235</td><td>-1.738</td><td>0.497</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-3.095</td><td rowspan="3">-2.481</td><td rowspan="3">0.614</td><td>HS</td><td>-2.740</td><td>-2.952</td><td>-2.563</td><td>0.389</td><td>-2.960</td><td>-2.537</td><td>0.423</td><td>-3.062</td><td>-2.435</td><td>0.627</td></tr><tr><td>Normal</td><td>-2.093</td><td>-2.255</td><td>-1.941</td><td>0.314</td><td>-2.295</td><td>-1.910</td><td>0.385</td><td>-2.359</td><td>-1.880</td><td>0.479</td></tr><tr><td>GHYP</td><td>-2.756</td><td>-2.968</td><td>-2.588</td><td>0.380</td><td>-2.968</td><td>-2.518</td><td>0.450</td><td>-3.041</td><td>-2.448</td><td>0.593</td></tr><tr><td rowspan="6">1%</td><td rowspan="3">-4.855</td><td rowspan="3">-4.137</td><td rowspan="3">0.718</td><td>HS</td><td>-4.428</td><td>-4.770</td><td>-4.161</td><td>0.609</td><td>-4.806</td><td>-4.112</td><td>0.694</td><td>-4.903</td><td>-4.005</td><td>0.898</td></tr><tr><td>Normal</td><td>-2.961</td><td>-3.157</td><td>-2.745</td><td>0.412</td><td>-3.198</td><td>-2.237</td><td>0.961</td><td>-3.348</td><td>-2.677</td><td>0.671</td></tr><tr><td>GHYP</td><td>-4.519</td><td>-4.870</td><td>-4.268</td><td>0.602</td><td>-4.874</td><td>-4.178</td><td>0.696</td><td>-5.103</td><td>-4.034</td><td>1.069</td></tr></tbody></table>
<table><caption>**Table 8.** A comparison between the bound estimates of the daily VaR based on the maxitive kernel estimation method and bootstrap confidence intervals of VaRs based on three different distributions: Historical simulation method (HS), Normal, and GHYP distributions. These VaR methods are applied to four daily stock returns at several probability levels a% and time horizon of 6 years.</caption><thead><tr><th rowspan="3">Returns</th><th rowspan="3">α</th><th colspan="3">Interval-Valued of VaR</th><th rowspan="3">pdf</th><th rowspan="3">Estimated VaR</th><th colspan="9">Bootstrap Confidence Intervals of VaRs</th></tr><tr><th rowspan="2">L.B.</th><th rowspan="2">U.B.</th><th rowspan="2">Width</th><th colspan="3">BCa (90%)</th><th colspan="3">BCa (95%)</th><th colspan="3">BCa (99%)</th></tr><tr><th>L.B.</th><th>U.B.</th><th>Width</th><th>L.B.</th><th>U.B.</th><th>Width</th><th>L.B.</th><th>U.B.</th><th>Width</th></tr></thead><tbody><tr><td rowspan="12">S&P500</td><td rowspan="3">10%</td><td rowspan="3">-1.189</td><td rowspan="3">-0.941</td><td rowspan="3">0.248</td><td>HS</td><td>-1.070</td><td>-1.156</td><td>-1.007</td><td>0.149</td><td>-1.168</td><td>-0.994</td><td>0.174</td><td>-1.201</td><td>-0.968</td><td>0.233</td></tr><tr><td>Normal</td><td>-1.233</td><td>-1.314</td><td>-1.155</td><td>0.159</td><td>-1.331</td><td>-1.141</td><td>0.190</td><td>-1.377</td><td>-1.111</td><td>0.266</td></tr><tr><td>GHYP</td><td>-1.029</td><td>-1.103</td><td>-0.959</td><td>0.144</td><td>-1.121</td><td>-0.946</td><td>0.175</td><td>-1.152</td><td>-0.913</td><td>0.239</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-1.679</td><td rowspan="3">-1.443</td><td rowspan="3">0.236</td><td>HS</td><td>-1.574</td><td>-1.674</td><td>-1.481</td><td>0.193</td><td>-1.693</td><td>-1.463</td><td>0.230</td><td>-1.729</td><td>-1.432</td><td>0.297</td></tr><tr><td>Normal</td><td>-1.593</td><td>-1.699</td><td>-1.502</td><td>0.197</td><td>-1.720</td><td>-1.483</td><td>0.237</td><td>-1.754</td><td>-1.449</td><td>0.305</td></tr><tr><td>GHYP</td><td>-1.548</td><td>-1.655</td><td>-1.468</td><td>0.187</td><td>-1.678</td><td>-1.450</td><td>0.228</td><td>-1.699</td><td>-1.410</td><td>0.289</td></tr><tr><td rowspan="3">1%</td><td rowspan="3">-2.939</td><td rowspan="3">-2.772</td><td rowspan="3">0.167</td><td>HS</td><td>-2.824</td><td>-2.979</td><td>-2.679</td><td>0.300</td><td>-3.007</td><td>-2.656</td><td>0.351</td><td>-3.086</td><td>-2.609</td><td>0.477</td></tr><tr><td>Normal</td><td>-2.269</td><td>-2.407</td><td>-2.153</td><td>0.254</td><td>-2.430</td><td>-2.128</td><td>0.302</td><td>-2.475</td><td>-2.083</td><td>0.392</td></tr><tr><td>GHYP</td><td>-2.810</td><td>-2.974</td><td>-2.665</td><td>0.309</td><td>-3.007</td><td>-2.637</td><td>0.370</td><td>-3.074</td><td>-2.582</td><td>0.492</td></tr><tr><td rowspan="12">Nikkei225</td><td rowspan="3">10%</td><td rowspan="3">-2.717</td><td rowspan="3">-2.037</td><td rowspan="3">0.680</td><td>HS</td><td>-2.360</td><td>-2.532</td><td>-2.207</td><td>0.325</td><td>-2.561</td><td>-2.179</td><td>0.382</td><td>-2.629</td><td>-2.135</td><td>0.494</td></tr><tr><td>Normal</td><td>-1.853</td><td>-1.996</td><td>-1.728</td><td>0.268</td><td>-2.021</td><td>-1.702</td><td>0.319</td><td>-2.120</td><td>-1.658</td><td>0.462</td></tr><tr><td>GHYP</td><td>-2.331</td><td>-2.497</td><td>-2.184</td><td>0.313</td><td>-2.544</td><td>-2.159</td><td>0.385</td><td>-2.607</td><td>-2.119</td><td>0.488</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-3.698</td><td rowspan="3">-2.991</td><td rowspan="3">0.707</td><td>HS</td><td>-3.258</td><td>-3.486</td><td>-3.053</td><td>0.433</td><td>-3.535</td><td>-3.180</td><td>0.355</td><td>-3.627</td><td>-2.958</td><td>0.669</td></tr><tr><td>Normal</td><td>-2.386</td><td>-2.565</td><td>-2.240</td><td>0.325</td><td>-2.596</td><td>-2.207</td><td>0.389</td><td>-2.697</td><td>-2.154</td><td>0.543</td></tr><tr><td>GHYP</td><td>-3.307</td><td>-3.533</td><td>-3.111</td><td>0.422</td><td>-3.571</td><td>-3.073</td><td>0.498</td><td>-3.653</td><td>-3.009</td><td>0.644</td></tr><tr><td rowspan="3">1%</td><td rowspan="3">-5.964</td><td rowspan="3">-5.284</td><td rowspan="3">0.680</td><td>HS</td><td>-5.589</td><td>-5.970</td><td>-5.252</td><td>0.718</td><td>-6.057</td><td>-5.206</td><td>0.851</td><td>-6.190</td><td>-5.055</td><td>1.135</td></tr><tr><td>Normal</td><td>-3.386</td><td>-3.614</td><td>-3.171</td><td>0.443</td><td>-3.659</td><td>-3.127</td><td>0.532</td><td>-3.746</td><td>-3.035</td><td>0.711</td></tr><tr><td>GHYP</td><td>-5.959</td><td>-6.389</td><td>-5.643</td><td>0.746</td><td>-6.451</td><td>-5.593</td><td>0.858</td><td>-6.631</td><td>-5.483</td><td>1.148</td></tr><tr><td rowspan="12">DJI</td><td rowspan="3">10%</td><td rowspan="3">-1.023</td><td rowspan="3">-0.815</td><td rowspan="3">0.208</td><td>HS</td><td>-0.917</td><td>-0.985</td><td>-0.856</td><td>0.129</td><td>-0.995</td><td>-0.845</td><td>0.150</td><td>-1.017</td><td>-0.823</td><td>0.194</td></tr><tr><td>Normal</td><td>-1.147</td><td>-1.226</td><td>-1.080</td><td>0.146</td><td>-1.236</td><td>-1.067</td><td>0.169</td><td>-1.266</td><td>-1.038</td><td>0.228</td></tr><tr><td>GHYP</td><td>-0.909</td><td>-0.972</td><td>-0.853</td><td>0.119</td><td>-0.994</td><td>-0.842</td><td>0.152</td><td>-1.022</td><td>-0.822</td><td>0.200</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-1.459</td><td rowspan="3">-1.252</td><td rowspan="3">0.207</td><td>HS</td><td>-1.341</td><td>-1.432</td><td>-1.273</td><td>0.159</td><td>-1.451</td><td>-1.260</td><td>0.191</td><td>-1.467</td><td>-1.227</td><td>0.240</td></tr><tr><td>Normal</td><td>-1.482</td><td>-1.569</td><td>-1.400</td><td>0.169</td><td>-1.587</td><td>-1.386</td><td>0.201</td><td>-1.618</td><td>-1.357</td><td>0.261</td></tr><tr><td>GHYP</td><td>-1.336</td><td>-1.422</td><td>-1.261</td><td>0.161</td><td>-1.437</td><td>-1.249</td><td>0.188</td><td>-1.468</td><td>-1.223</td><td>0.245</td></tr><tr><td rowspan="3">1%</td><td rowspan="3">-2.436</td><td rowspan="3">-2.192</td><td rowspan="3">0.244</td><td>HS</td><td>-2.292</td><td>-2.426</td><td>-2.183</td><td>0.243</td><td>-2.451</td><td>-2.161</td><td>0.290</td><td>-2.505</td><td>-2.124</td><td>0.381</td></tr><tr><td>Normal</td><td>-2.109</td><td>-2.223</td><td>-2.005</td><td>0.218</td><td>-2.253</td><td>-1.984</td><td>0.269</td><td>-2.292</td><td>-1.938</td><td>0.354</td></tr><tr><td>GHYP</td><td>-2.326</td><td>-2.450</td><td>-2.214</td><td>0.236</td><td>-2.474</td><td>-2.186</td><td>0.288</td><td>-2.533</td><td>-2.144</td><td>0.389</td></tr><tr><td rowspan="12">CAC40</td><td rowspan="3">10%</td><td rowspan="3">-2.501</td><td rowspan="3">-1.888</td><td rowspan="3">0.613</td><td>HS</td><td>-2.197</td><td>-2.319</td><td>-2.090</td><td>0.229</td><td>-2.349</td><td>-2.066</td><td>0.283</td><td>-2.391</td><td>-2.024</td><td>0.367</td></tr><tr><td>Normal</td><td>-1.764</td><td>-1.863</td><td>-1.665</td><td>0.198</td><td>-1.880</td><td>-1.650</td><td>0.230</td><td>-1.919</td><td>-1.623</td><td>0.296</td></tr><tr><td>GHYP</td><td>-2.218</td><td>-2.351</td><td>-2.113</td><td>0.238</td><td>-2.372</td><td>-2.089</td><td>0.283</td><td>-2.416</td><td>-2.053</td><td>0.363</td></tr><tr><td rowspan="3">5%</td><td rowspan="3">-3.577</td><td rowspan="3">-2.976</td><td rowspan="3">0.601</td><td>HS</td><td>-3.235</td><td>-3.405</td><td>-3.091</td><td>0.314</td><td>-3.438</td><td>-3.062</td><td>0.376</td><td>-3.510</td><td>-3.005</td><td>0.505</td></tr><tr><td>Normal</td><td>-2.265</td><td>-2.384</td><td>-2.152</td><td>0.232</td><td>-2.410</td><td>-2.124</td><td>0.286</td><td>-2.456</td><td>-2.091</td><td>0.365</td></tr><tr><td>GHYP</td><td>-3.142</td><td>-3.309</td><td>-3.005</td><td>0.304</td><td>-3.339</td><td>-2.978</td><td>0.361</td><td>-3.406</td><td>-2.931</td><td>0.475</td></tr><tr><td rowspan="3">1%</td><td rowspan="3">-5.477</td><td rowspan="3">-5.023</td><td rowspan="3">0.454</td><td>HS</td><td>-5.110</td><td>-5.367</td><td>-4.880</td><td>0.487</td><td>-5.412</td><td>-4.825</td><td>0.587</td><td>-5.532</td><td>-4.729</td><td>0.803</td></tr><tr><td>Normal</td><td>-3.205</td><td>-3.370</td><td>-3.054</td><td>0.316</td><td>-3.407</td><td>-3.030</td><td>0.377</td><td>-3.470</td><td>-2.949</td><td>0.521</td></tr><tr><td>GHYP</td><td>-5.281</td><td>-5.542</td><td>-5.053</td><td>0.489</td><td>-5.578</td><td>-5.013</td><td>0.565</td><td>-5.677</td><td>-4.921</td><td>0.756</td></tr></tbody></table>
Next, in order to inspect the goodness-of-fit of the used models a graphical tool (PP-plot) is constructed (Figure 5) to compare the empirical cumulative function to the fitted cumulative functions. This plot confirms that the Epanechnikov kernel function and the GHYP distribution give a good global fit for the 4 returns data. The points of the PP-plot are close to the 45 degree line and they also lie within the interval estimation using maxitive kernel.

In contrast, from the same figure, the PP-plot of the normal cumulative function against the empirical cumulative function shows that the left end pattern is above the 45 degree line and the right end is below it. Thus, the normal distribution underestimates the VaR at low probability level. This is due to the fact that the normal distribution ignores the presence of fat tails in the actual distribution. Based on these results, we can conclude that the GHYP and the maxitive kernel method provide a better fit than a normal distribution to market return data. Thus, the maxitive kernel method seems to be a good choice to estimate the risk for VaR.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** A graphical tool (PP-plot) of the theoretical CDF of Epanechnikov (in black), Generalized HYPerbolic (GHYP, in green), and normal (in violet) versus the empirical CDF for the four returns data sets: (a) Daily CAC40 (b) Daily DJI (c) Daily Nikkei225 and (d) Daily S&P500. In each of the four plots, the empirical CDF on the horizontal axis and the theoretical CDF on the vertical axis. The highest (in red) and lowest (in blue) dotted lines correspond respectively to the upper and lower bounds of the interval CDF using maxitive kernel.</figcaption></figure>
## 6. Conclusions

Using an estimate of the Value-at-Risk (VaR) based on a small-sized sample may pose a risk to financial application due to the high dependence of this VaR estimate to the computation method. In fact, computing a VaR estimate can be performed in many ways including parametric and nonparametric approaches. We have shown, with experiments based on daily closing prices data of four stock indexes, that no method can be said to be optimal to achieve this estimate. Moreover, the bias induced by choosing a particular method is particularly sensitive for estimates based on small samples. In this paper, we have presented a new method for computing the VaR which highly differs from its competitors in the sense that it is interval-valued. This interval-valued VaR estimate is the set of all estimates that could have been obtained, with the same data sample, by using a set of kernel-based estimation methods. In our experiments, we noted that the output of parametric methods, like the GHYP VaR estimate, always belong to the interval-valued VaR we propose while others, like the normal VaR estimate, do not. It appears that VaR estimates belonging to the interval-valued VaR estimate we propose are likely to be less risky than those that do not. Moreover, a wide interval-valued VaR estimate is a marker of a high risk for a trader since it reflects thick tails, pronounced skewness, and excess kurtosis of financial asset price returns.

The interval-valued estimation, based on maxitive kernel, that we propose in this paper is a convex envelope of the kernel VaR estimation. Although the VaR is probably one of the most popular tool in risk management, an alternative measure for Value-at-Risk which satisfied the conditions for a coherent risk measure has been proposed (see e.g., Rockafellar and Uryasev 2000; Artzner et al. 1999). This risk measure is called Expected Shortfall (ES; also known as conditional VaR or average VaR). Indeed, the Basel Committee published, in January 2016 Basel Committee on Banking Supervision (2016), revised standards for minimum capital requirements for market risk which include a shift from Value-at-Risk to expected shortfall as the preferred risk measure. Then, the future research should be conducted into finding an interval-valued estimation of ES, based on maxitive kernel, and to compare it with the interval-valued estimation of VaR and the ES for many distributions. In this regard, Broda and Paolella (2011) presents easily-computed expressions for evaluating the expected shortfall for numerous distributions which are now commonly used for modelling asset returns. Another interesting avenue of research would be to construct a CoVaR interval combining GARCH model with maxitive kernel.

**Author Contributions:** All authors contributed equally to the paper.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Appendix A. Matlab Code for Maxitive Interval-Valued Estimation of VaR

### % Data definition

Data = load('S&P500_daily.txt');

NbData = length(Data);

### % Definition of visualization samples

SampleSize = length(Data);

kernelBandwidth = 0.25;

MIN = min(Data);

MAX = max(Data);

delta = (MAX - MIN) / SampleSize;

MIN = MIN - 50*delta;

MAX = MAX + 50*delta;

delta = (MAX - MIN) / SampleSize;

Absc = sort(Data);

SampleSize = length(Absc);

bandwidth = kernelBandwidth;

### % Summative kernel estimation of CDF
### % Epanechnikov kernel

```fortran
for n=1:SampleSize
    FonctionEscalier1 = KernelIntegralEpanechnikov(Absc(n)-Data, 0, 0.305);
    CDFEpanechnikov(n) = sum(FonctionEscalier1) / NbData;
end
```

### % Cosine kernel

```fortran
for n=1:SampleSize
    FonctionEscalier6 = KernelIntegralCosine(Absc(n)-Data, 0, 0.7);
    CDFBCosine(n) = sum(FonctionEscalier6) / NbData;
end
```

### % Biweight kernel

```fortran
for n=1:SampleSize
    FonctionEscalier4 = KernelIntegralBiweight(Absc(n)-Data, 0, 0.05);
    CDFBiweight(n) = sum(FonctionEscalier4) / NbData;
end
```

### % Triweight kernel

```fortran
for n=1:SampleSize
    FonctionEscalier5 = KernelIntegralTriweight(Absc(n)-Data, 0, 0.06);
    CDFTriweight(n) = sum(FonctionEscalier5) / NbData;
end
```

### % Interval-valued estimation of CDF

```fortran
cumul_upp_PLUS = zeros(SampleSize, 1);
cumul_low_PLUS = zeros(SampleSize, 1);
cumul_upp Moins = zeros(SampleSize, 1);
cumul_low Moins = zeros(SampleSize, 1);
```

```fortran
for n=1:SampleSize
```

```fortran
for k=1:NbData
    if Data(k) <= (Absc(n))
        cumul_upp_PLUS(n) = cumul_upp_PLUS(n) + 1;
    else
        Possibility = Max_kernel_Triangular(Data(k), Absc(n), bandwidth);
        cumul_upp_PLUS(n) = cumul_upp_PLUS(n) + Possibility;
end

if Data(k) >= (Absc(n))
    cumul_low_PLUS(n) = cumul_low_PLUS(n) + 1;
else
    Possibility = Max_kernel_Triangular(Data(k), Absc(n), bandwidth);
    cumul_low_PLUS(n) = cumul_low_PLUS(n) + Possibility;
end
```

```fortran
end
```

```fortran
cumul_low_PLUS = NbData - cumul_low_PLUS;
cumul_low_PLUS = cumul_low_PLUS / NbData;
cumul_upp_PLUS = cumul_upp_PLUS / NbData;
```
```matlab
y = zeros(size(x));
```

```matlab
xxx = x.*x.*x;
```

```matlab
xxxxx = x.*x.*x.*x.*x;
```

```matlab
y = ( ( x > -1 ) .* ( x <= 0 ) .* ((15/16)*x - (5/8)*xxx + (3/16)*xxxxx + 0.5) );
```

```matlab
y = y + ( ( x < 1 ) .* ( x >= 0 ) .* ((15/16)*x - (5/8)*xxx + (3/16)*xxxxx + 0.5) + (x>=1) );
```

## Triangular Maxitive Kernel

```matlab
% Triangulaire maxitive kernel
```

```matlab
function [y] = Max_kernel_Triangular(x, centre, delta)
```

```matlab
x = (x-centre)/delta; y = zeros(size(x));
```

```matlab
x = abs(x);
```

```matlab
Mask = x <= 1;
```

```matlab
y = ( 1.0 - x ) .* Mask;
```

```matlab
return;
```

## References

* Adrian, Bowman, Peter Hall, and Tania Prvan. 1998. Bandwidth selection for the smoothing of distribution functions. *Biometrika* 85: 799-808.

* Altman, Naomi, and Christian Leger. 1995. Bandwidth selection for kernel distribution function estimation. *Journal of Statistical Planning and Inference* 46: 195-214. [CrossRef]

* Artzner, Philippe, Freddy Delbaen, Jean-Marc Eber, and David Heath. 1999. Coherent measures of risk. *Mathematical Finance* 9: 203-28. [CrossRef]

* Barndorff-Nielsen, Ole. 1978. Hyperbolic distributions and distributions on hyperbolae. *Scandinavian Journal of Statistics* 5: 151-57.

* Basel Committee on Banking Supervision. 2016. Minimum Capital Requirements for Market Risk. Technical Report. Available online: http://www.bis.org (accessed on 3 November 2018).

* Braione, Manuela, and Nicolas Scholtes. 2016. Forecasting value-at-risk under different distributional assumptions. *Econometrics* 4: 1-27. [CrossRef]

* Broda, Simon A., and Marc S. Paolella. 2011. Expected shortfall for distributions in finance. In *Statistical Tools for Finance and Insurance*. Edited by Pavel Čížek, Wolfgang Karl Härdle and Rafał Weron. Berlin: Springer, pp. 57-99.

* Butler, J. S., and Barry Schachter. 1998. Estimating value-at-risk with a precision measure by combining kernel estimation with historical simulation. *Review of Derivatives Research* 1: 371-90.

* Charpentier, Arthur, and Abder Oulidi. 2010. Beta kernel quantile estimators of heavy-tailed loss distributions. *Statistics and Computing* 20: 35-55. [CrossRef]

* Chen, Song Xi, and Cheng Yong Tang. 2005. Nonparametric inference of value-at-risk for dependent financial returns. *Journal of Financial Econometrics* 3: 227-55. [CrossRef]

* Choquet, Gustave 1953. Théorie des capacités. *Annales de l'Institut Fourier* 5: 131-295. [CrossRef]

* Danielsson, Jon, and Casper De Vries. 2000. Value-at-risk and extreme returns. *Annales d'économie et de Statistique* 60: 236-69. [CrossRef]

* Davison, Anthony Christopher, and David Victor Hinkley. 1997. *Booststrap Methods and Their Applications*. New York: Cambridge University Press.

* De Cooman, Gert. 1997. Possibility theory, 1: The measure-and integral-theoretic groundwork. *International Journal of General Systems* 25: 291-323. [CrossRef]

* Denneberg, Dieter 1994. *Non Additive Measure and Integral*. Dordrecht: Kluwer Academic Publishers.

* Destercke, Sebastien, Didier Dubois, and Eric Chojnacki. 2007. On the relationships between random sets, possibility distributions, p-boxes and clouds. Paper presented at 28th Linz Seminar on Fuzzy Set Theory, Linz, Austria, February 6-10.

* Dowd, Kevin. 2005. *Measuring Market Risk*, 2nd ed. New York: John Wiley & Sons Inc.
Dubois, Didier. 2006. Possibility theory and statistical reasoning. *Computational Statistics and Data Analysis* 51: 47–69. [CrossRef]

Dubois, Didier, and Henri Prade. 1988. *Théorie des possibilités: Applications à la représentation des connaissances en informatique*. Paris: Masson.

Dubois, Didier, Henri Prade, Laurent Foulloy, and Gilles Mauris. 2004. Probability-possibility transformations, triangular fuzzy sets, and probabilistic inequalities. *Reliable Computing* 10: 273–97. [CrossRef]

Eberlein, Ernst, and Ulrich Keller. 1995. Hyperbolic distributions in finance. *Bernoulli* 1: 281–99. [CrossRef]

Eberlein, Ernst, Ulrich Keller, and Karsten Prause. 1998. New insights into smile, mispricing, and value at risk: The hyperbolic model. *The Journal of Business* 71: 371–405. [CrossRef]

Efron, Bradley. 1979. Bootstrap methods: Another look at the jackknife. *The Annals of Statistics* 7: 1–26. [CrossRef]

Efron, Bradley, and Robert J. Tibshirani. 1993. *An Introduction to the Bootstrap*. London: Chapman and Hall.

Engle, F. Robert. 1982. Autoregressive conditional heteroscedasticity with estimates of the variance of united kingdom inflation. *Econometrica* 50: 987–1007. [CrossRef]

Epanechnikov, V. A. 1969. Nonparametric estimation of a multidimensional probability density. *Theory of Probability and its Applications* 14: 153–58. [CrossRef]

Fox, John. 2002. Bootstrapping Regression Models Appendix to an R and S-Plus Companion to Applied Regression. Available online: http://statweb.stanford.edu/~owen/courses/305a/FoxOnBootingRegInR.pdf (accessed on 10 December 2018).

Gourieroux, Christian, Jean-Paul Laurent, and Olivier Scaillet. 2000. Sensitivity analysis of values at risk. *Journal of Empirical Finance* 7: 225–45. [CrossRef]

Hartz, Christoph, Stefan Mittnik, and Marc Paolella. 2006. Accurate value-at-risk forecasting based on the normal-garch model. *Computational Statistics and Data Analysis* 51: 2295–312. [CrossRef]

Jones, Michael Chris. 1990. The performance of kernel density functions in kernel distribution function estimation. *Statistics & Probability Letters* 37: 129–32.

Jorion, Philippe. 1996. Risk2: Measuring the risk in value at risk. *Financial Analysts Journal* 52: 47–56. [CrossRef]

Kuester, Keith, Stefan Mittnik, and Marc S. Paolella. 2006. Value-at-risk prediction: A comparison of alternative strategies. *Journal of Financial Econometrics* 4: 53–89. [CrossRef]

Kupiec, Paul. 1995. Techniques for verifying the accuracy of risk measurement models. *Journal of Derivatives* 3: 73–84. [CrossRef]

Linsmeier, Thomas J., and Neil D. Pearson. 1997. Quantitative disclosures of market risk in the sec release. *Pearson Accounting Horizons* 11: 107–35.

Loquin, Kevin, and Olivier Strauss. 2008a. Imprecise functional estimation: The cumulative distribution case. In *Soft Methods in Probability and Statistics*. Edited by Didier Dubois, Maria Asuncion Lubiano, Henri Prade, Przemysław Grzegorzewski and Olgierd Hryniewicz. Advanced in Soft Computing. Heidelberg: Springer. vol. 48, pp. 175–82.

Loquin, Kevin, and Olivier Strauss. 2008b. On the granularity of summative kernels. *Fuzzy Sets and Systems* 159: 1952–72. [CrossRef]

McNeil, Alexander J., Rudiger Frey, and Paul Embrechts. 2005. *Quantitative Risk Management*. Princeton: Princeton University Press.

Nadaraya, Elizbar. A. 1964. Some new estimates for distribution function. *Theory of Probability and its Application* 9: 497–500. [CrossRef]

Paolella, Marc S. 2007. *Intermediate Probability: A Computational Approach*. Chichester: John Wiley & Sons.

Paolella, Marc S., and Paweł Polak. 2015. Comfort: A common market factor non-gaussian returns model. *Journal of Econometrics* 187: 593–605. [CrossRef]

Parzen, Emanuel. 1962. On estimation of a probability density function and mode. *The Annals of Mathematical Statistics* 33: 1065–76. [CrossRef]

Polansky, Alan M., and Edsel R. Baker. 2000. Multistage plug-in bandwidth selection for kernel distribution function estimates. *Journal of Statistical Computation and Simulation* 65: 63–80. [CrossRef]

Pritsker, M. 1997. Evaluating value at risk methodologies: Accuracy versus computational time. *Journal of Financial Services Research* 12: 201. [CrossRef]

Quintela Del Rio, A., and G. Estevez-Perez. 2013. Nonparametric kernel distribution function estimation with kerdiest: An r package for bandwidth choice and applications. *Journal of Statistical Software* 50: 1–21.
Rockafellar, R. Tyrrell, and Stanislav Uryasev. 2000. Optimization of conditional value-at-risk. *Journal of Risk* 2: 21–41. [CrossRef]

Rosenblatt, Murray. 1956. Remarks on some nonparametric estimates of a density function. *The Annals of Mathematical Statistics* 27: 832–37. [CrossRef]

Sarda, Pascal. 1993. Smoothing parameter selection for smooth distribution functions. *Journal of Statistical Planning and Inference* 35: 65–75. [CrossRef]

Shafer, Glenn. 1976. *A Mathematical Theory of Evidence*. Princeton: Princeton University Press.

Sheather, Simon J., and J. Stephen Marron. 1990. Kernel quantile estimators. *Journal of the American Statistical Association* 85: 410–16. [CrossRef]

Silvermann, Bernard Walter. 1986. *Density Estimation for Statistics and Data Analysis*. London: Chapman and Hall.

Singh, Radhey S., Theo Gasser, and Bhola Prasad. 1983. Nonparametric estimates of distributions functions. *Communication in Statistics—Theory and Methods* 12: 2095–108. [CrossRef]

Sørensen, Michael, and Bo Martin Bibby. 1997. A hyperbolic diffusion model for stock prices. *Finance and Stochastics* 1: 25–41.

Swanepoel, Jan W. H. 1988. Mean integrated squared error properties and optimal kernels when estimating a distribution function. *Communications in Statistics—Theory and Methods* 17: 3785–99. [CrossRef]

Wand, M. P., and M. C. Jones. 1995. *Kernel Smoothing*. London: Chapman and Hall.

Watson, Geoffrey S., and M. R. Leadbetter. 1964. Hazard analysis II. *Sankhya Series A* 26: 101–16.

Winter, B. 1973. Strong uniform consistency of integrals of density estimators. *Canadian Journal of Statistics* 1: 247–53. [CrossRef]

Winter, B. 1979. Convergence rate of perturbed empirical distribution functions. *Journal of Applied Probability* 16: 163–73. [CrossRef]

Yamato, Hajime 1973. Uniform convergence of an estimator of a distribution function. *Bulletin on Mathematical Statistics* 15: 69–78.

<figure><img src="image_8.png" /></figure>

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).