Tutorials › Exam ASTAM - Advanced Short-Term Actuarial Mathematics › Simulation Counts for Stable Tail-Risk Estimates

Exam ASTAM · Simulation Precision and Design · Tutorial 669 of 979

Simulation Counts for Stable Tail-Risk Estimates

Plan simulation runs by translating a tail-risk precision requirement into a minimum number of observations beyond the target quantile.

Intermediate 12 min read

What You'll Learn

  • Relate a risk-measure confidence level to the probability of a simulated loss entering the tail.
  • Calculate the expected number of tail observations in an independent simulation.
  • Use exact binomial probabilities to select a simulation size for a target count and confidence.
  • Apply a conservative Chernoff bound when exact binomial calculations are inconvenient.
  • Recognize why a sufficient tail count is a planning diagnostic, not a guarantee of stable risk estimates.

Why Plan for Tail Observations?

A simulation may produce a precise-looking estimate of a high-quantile risk measure while containing very few losses beyond the estimated quantile. That is a warning: the estimate is being informed by a small number of observations in the part of the distribution that matters most. Before running a large model, it is useful to ask how many simulated losses are expected in that tail, and how likely the run is to contain at least a chosen minimum number.

This tutorial develops that count-based planning approach for value-at-risk (VaR) and tail value-at-risk (TVaR). It complements methods for quantile standard errors and confidence intervals: the goal here is to plan a simulation size around a required tail count, rather than to derive a full confidence interval for the risk measure.

Definition: For a loss \(X\) and confidence level \(q\), define \(\operatorname{VaR}_q(X)=\inf\{x:F_X(x)\ge q\}\). For a continuous loss distribution, the probability of exceeding this quantile is \(p=1-q\). A simulation of \(n\) independent losses has tail count \(K=\sum_{j=1}^n I(X_j>\operatorname{VaR}_q(X))\), where \(I(\cdot)\) is 1 when its condition holds and 0 otherwise. Under the stated continuity and independence assumptions, \(K\sim\operatorname{Bin}(n,p)\) and \(E[K]=np\).

For a continuous distribution, this count has a simple interpretation: each simulated loss independently has probability \(p\) of exceeding the true VaR. Thus \(np\) is the expected number of exceedances. It is not a guarantee that a particular simulation will contain that many. A planning requirement should therefore specify both a minimum count \(m\) and a probability \(1-\varepsilon\) of meeting it:

$$ P(K\ge m)\ge 1-\varepsilon. $$

The count requirement is a practical design criterion, not a universal definition of adequate precision. The appropriate \(m\) depends on the application, the severity of tail losses, and the desired precision. In particular, a minimum count alone does not determine the standard error of a tail average.

Exact Binomial Planning

When the exceedance probability \(p\) is known or specified by the model, the exact planning calculation is direct. For \(K\sim\operatorname{Bin}(n,p)\), choose the smallest integer \(n\) such that

$$ P(K\ge m) =1-\sum_{k=0}^{m-1}{n\choose k}p^k(1-p)^{n-k} \ge 1-\varepsilon. $$

This calculation accounts for the variability of the count. Merely choosing \(n\) so that \(np=m\) would make \(m\) the mean count, not a count reached with high probability. In practice, exact binomial probabilities can be evaluated with a statistical package or a reliable binomial probability table. For a small \(m\), the finite sum can also be calculated directly.

Worked Example: Exact Planning for a 90% VaR Tail

A hypothetical insurer uses independent simulated annual losses to estimate \(\operatorname{VaR}_{0.90}\). Assume the loss distribution is continuous. The simulation plan requires at least two losses above the true VaR with probability at least 95%. Find the smallest simulation size using the exact binomial calculation.

Method. The exceedance probability is \(p=1-0.90=0.10\). Therefore \(K\sim\operatorname{Bin}(n,0.10)\), and the requirement is \(P(K\ge2)\ge0.95\).

For any \(n\), the probability of fewer than two exceedances is the probability of zero plus the probability of one:

$$ P(K<2)=(0.90)^n+n(0.10)(0.90)^{n-1}. $$

Test \(n=45\). The two terms are

$$ P(K=0)=(0.90)^{45}=0.00872796, $$
$$ P(K=1)=45(0.10)(0.90)^{44} =4.5(0.00969774)=0.04363982. $$

Thus \(P(K<2)=0.05236778\), so \(P(K\ge2)=1-0.05236778=0.94763222\). This is below 0.95; 45 simulations do not meet the requirement.

Now test \(n=46\):

$$ P(K=0)=(0.90)^{46}=0.00785517, $$
$$ P(K=1)=46(0.10)(0.90)^{45} =4.6(0.00872796)=0.04014863. $$

Hence \(P(K<2)=0.04800380\), and \(P(K\ge2)=1-0.04800380=0.95199620\). The requirement is met at \(n=46\), while it fails at \(n=45\); because the probability of at least two exceedances cannot decrease as \(n\) increases, 46 is the smallest acceptable size. The common mistake avoided is treating the expected count \(np\) as if it were a probability guarantee.

A Conservative Bound for Larger Plans

For a large required count, repeatedly evaluating binomial sums can be inconvenient. A conservative alternative is the multiplicative Chernoff lower-tail bound. Let \(K\sim\operatorname{Bin}(n,p)\), let \(\mu=np\), and suppose \(0<m\le\mu\). Set \(\delta=1-m/\mu\), so \(0\le\delta<1\). Then

$$ P(K<m)\le \exp\left(-\frac{\mu\delta^2}{2}\right) = \exp\left(-\frac{(\mu-m)^2}{2\mu}\right). $$

The last expression bounds the chance of falling short of the target count. It follows that a sufficient condition for \(P(K\ge m)\ge1-\varepsilon\) is

$$ \exp\left(-\frac{(\mu-m)^2}{2\mu}\right)\le\varepsilon. $$

This bound may require more simulations than the exact binomial calculation. It is a guaranteed sufficient condition under the binomial model, not an estimate of the exact shortfall probability.

Derivation: Write \(K\) as a sum of \(n\) independent Bernoulli indicators, each with success probability \(p\). For \(t>0\), Markov's inequality applied to \(e^{-tK}\) gives \(P(K\le(1-\delta)\mu)\le e^{t(1-\delta)\mu}E[e^{-tK}]\). Independence gives \(E[e^{-tK}]=(1-p+pe^{-t})^n\). Optimizing this expression over \(t>0\) yields the binomial lower-tail Chernoff bound; the standard inequality for its exponent gives \(P(K\le(1-\delta)\mu)\le e^{-\mu\delta^2/2}\). Since \(K<m\) is contained in the event \(K\le m=(1-\delta)\mu\), this also bounds the probability of fewer than \(m\) observations.

If it is convenient to plan first in terms of a required mean count \(\mu\), let \(c=\ln(1/\varepsilon)\). Solving the sufficient-condition inequality for \(\mu\ge m\) gives

$$ \mu\ge m+c+\sqrt{c^2+2mc}, \qquad n\ge \left\lceil\frac{m+c+\sqrt{c^2+2mc}}{p}\right\rceil. $$

The ceiling is needed because the simulation size must be an integer. This formula plans for the chosen Chernoff guarantee; it does not give the smallest size under the exact binomial calculation.

Worked Example: A 99% VaR Tail-Count Guarantee

A hypothetical risk team plans independent simulations for \(\operatorname{VaR}_{0.99}\) of annual portfolio loss. It wants at least 20 exceedances and a probability of success of at least 95%. Assess a proposed run of \(n=3{,}500\) simulations using the Chernoff bound.

Method. Continuity gives \(p=1-0.99=0.01\). The mean count is \(\mu=np=3{,}500(0.01)=35\). The target \(m=20\) is below the mean, as required for the lower-tail bound. Then \(\delta=1-20/35=3/7\).

$$ P(K<20)\le \exp\left(-\frac{35(3/7)^2}{2}\right) =\exp\left(-\frac{315}{98}\right) =\exp\left(-\frac{45}{14}\right) \approx0.0402. $$

Therefore

$$ P(K\ge20)\ge1-\exp(-315/98)\approx0.9598>0.95. $$

The proposed run meets the target according to this conservative bound. The bound on the probability of fewer than 20 observations is about 4.02%, which is below the allowed 5%. The common arithmetic trap is to round the exponential too early or incorrectly: \(\exp(-315/98)\) is approximately 0.040184, which rounds to 0.0402, not 0.0401. The corresponding success lower bound rounds to 0.9598.

Planning Tail Value-at-Risk Simulations

For a continuous loss distribution, TVaR at level \(q\) is the expected loss conditional on exceeding \(\operatorname{VaR}_q(X)\). A simulation-based tail average therefore relies on observations from a region with probability \(1-q\). A first count-based plan is to require at least \(m\) such observations with a specified probability. The same binomial calculations apply, using \(p=1-q\).

This count requirement is especially useful as a minimum design check for a tail average: a run with only a handful of tail losses is unlikely to provide a dependable description of their average. But meeting the count target does not, on its own, establish a desired standard error. Tail-loss variability matters too; a few very large values can dominate the average even when the number of exceedances appears adequate.

Worked Example: Planning a TVaR Tail Count

A hypothetical benefits pool is modeled using independent annual simulated losses. For a continuous model, the team wants at least 40 losses beyond \(\operatorname{VaR}_{0.975}\), with a probability of at least 90%. Find a sufficient simulation size using the Chernoff bound.

Method. The tail probability is \(p=1-0.975=0.025\). Choose a mean count \(\mu=60\), which corresponds to

$$ n=\frac{\mu}{p}=\frac{60}{0.025}=2{,}400. $$

For \(m=40\), \(\delta=1-40/60=1/3\). The lower-tail bound is

$$ P(K<40)\le \exp\left(-\frac{60(1/3)^2}{2}\right) =\exp(-10/3) \approx0.0357. $$

It follows that \(P(K\ge40)\ge1-0.0357=0.9643\), which exceeds the required 0.90. Thus 2,400 independent simulations are sufficient under this bound. This is a conservative plan, not a claim that 2,400 is the smallest acceptable size. The common mistake avoided is using the 2.5% tail probability as though it were the expected tail count; the expected count is \(2{,}400(0.025)=60\).

Assumptions and Practical Limits

Before applying these calculations, check what each simulated observation represents. The binomial model requires independent, identically distributed losses and a fixed exceedance probability \(p\). The simple relation \(p=1-q\) for exceedance of VaR assumes a continuous distribution. If the loss distribution has an atom at the quantile, the exceedance probability may differ from \(1-q\); use the model's actual probability of exceeding the chosen threshold or specify a quantile convention that handles ties.

The formulas also concern ordinary independent simulation counts. With importance sampling, observations have weights and the raw number of simulated points in a region need not represent the number of equally informative tail observations. With correlated simulation output, the binomial model does not apply directly. In either case, an appropriate design must account for the weighting or dependence rather than treating the raw count as an independent binomial count.

A sound planning sequence is to define the quantile and its tail probability, set a justified minimum tail count, choose a target confidence, and then calculate \(n\) using an exact binomial probability or a clearly identified conservative bound. After the run, report the observed tail count alongside the risk estimate. For TVaR, also assess the variability of the tail losses; for VaR, count alone does not replace quantile-uncertainty analysis.

Key takeaway: For independent continuous-loss simulations, the number of observations beyond a level-\(q\) quantile is binomial with mean \(n(1-q)\). Plan against the probability of falling below a chosen tail count, not just against the expected count, and treat that count as a minimum diagnostic rather than a complete precision guarantee.

Check Your Understanding

For each problem, assume independent simulation losses and a continuous distribution unless the problem states otherwise.

  1. A team simulates \(n=600\) losses for \(\operatorname{VaR}_{0.95}\). What is the expected number of exceedances? (Answer: 30)
  2. A model uses 2,500 independent simulations for \(\operatorname{VaR}_{0.99}\). Find the expected number of losses beyond the true quantile. (Answer: 25)
  3. For \(\operatorname{VaR}_{0.98}\), a team wants at least 30 exceedances with at least 95% probability. Using the Chernoff mean-count formula, find a sufficient integer \(n\). (Answer: 2,337 simulations)
  4. For a TVaR calculation at level 0.975, require at least 50 exceedances with at least 95% probability. Using the Chernoff mean-count formula, find a sufficient integer \(n\). (Answer: 2,823 simulations)
  5. A continuous loss model has \(P(X>\operatorname{VaR}_{0.80}(X))=0.20\). What is the smallest \(n\) for which the probability of at least one exceedance is at least 90%? Use \(1-(0.80)^n\ge0.90\). (Answer: 11 simulations)
Clinical Biostats