Ruin Theory
Ruin Theory
Ruin theory
4.1 INTRODUCTION
In this chapter we focus again on collective risk models, but now in the long
term. We consider the development in time of the capital of an insurer.
This is a stochastic process which increases continuously because of the earned
premiums, and decreases stepwise because of the payment of claims. When the
capital becomes negative, we say that ruin occurs. Let denote the probability
that this ever happens, provided that the annual premium and the claims process
remain unchanged. This probability is a useful tool for the management since it
serves as an indication of the soundness of the insurer’s combined premiums and
claims process, given the available initial capital A high probability of
ruin indicates instability: measures such as reinsurance or raising some premiums
should be considered, or the insurer should attract extra working capital.
The probability of ruin enables one to compare portfolios with each other, but
we cannot attach any absolute meaning to the probability of ruin, as it doesn’t
actually represent the probability that the insurer will go bankrupt in the near
future. First of all, it might take centuries for ruin to actually happen. Moreover,
potential interventions in the process, for instance paying out dividends or raising
the premium for risks with an unfavorable claims performance, are ruled out in the
81
82 RUIN THEORY
difference between the earned premiums and the total payments up to any mo-
ment. Using this mgf, we will determine the value of the ruin probability in case
the claims are distributed according to variants of the exponential distribution.
Next, we will consider some approximations for the ruin probability.
where
with
A typical realization of the risk process is depicted in Figure 4.1. The random
variables denote the time points at which a claim occurs. The slope of
the process is if there are no claims; if, however, for some then the
capital drops by which is the size of the claim. Since in Figure 4.1, at time
the total of the incurred claims is larger than the initial
capital plus the earned premium the remaining surplus is less than
0. This state of the process is called ruin and the point in time at which this occurs
for the first time is denoted by T. So,
if for all
The random variable T is defective, as the probability of is positive. The
probability that ruin ever occurs, i.e., the probability that T is finite, is called the
84 RUIN THEORY
Before we turn to the claim process i.e., the total claims up to time we
first look at the process of the number of claims up to We will assume that
is a so-called Poisson process:
Definition 4.2.1 (Poisson process)
The process is a Poisson process if for some intensity the increments
of the process have the following property:
Next to this global definition of the claim number process, we can also consider
infinitesimal increments where the infinitesimal ‘number’
EXPONENTIAL UPPER BOUND 85
again is positive, but smaller than any real number larger than 0. For the Poisson
process we have:
Actually, these equalities are not really quite equalities: they are only valid if we
ignore terms of order
A third way to define such a process is by considering the waiting times
Because Poisson processes are memoryless, these waiting times are independent
exponential random variables, and they are also independent of the history of
the process. This can be shown as follows: if the history H represents an arbitrary
realization of the process up to time with the property that then
For most distributions, there is no explicit expression for the adjustment coeffi-
cient. To facilitate solving (4.10) by a spreadsheet or a computer program, one can
use the fact that see Exercise 4.3.2.
In the next theorem, we prove F. Lundberg’s famous exponential inequality for the
ruin probability. Surprisingly, the proof involves mathematical induction.
assume it has a size which has a probability Then the capital at that
moment equals Integrating over and yields
Note that if R is replaced by any other real number, the expression in square
brackets in (4.17) is unequal to 1, so in fact the adjustment coefficient is the unique
positive number R with the property that is a martingale.
In this section we give an expression for the ruin probability which involves the
mgf of U (T), i.e., the capital at the moment of ruin, conditionally given the event
RUIN PROBABILITY AND EXPONENTIAL CLAIMS 89
that ruin occurs in a finite time period. This expression enables us to give an exact
expression for the ruin probability in case of an exponential distribution.
Theorem 4.4.1 (Ruin probability)
The ruin probability for satisfies
From Remark 4.3.5, we know that the left-hand side equals For the first
conditional expectation in (4.19) we take and write, using
see also (4.17):
The second term vanishes if For the first term, note that
has an expected value and a variance
Because of Chebyshev’s inequality, it suffices to choose the function such
that We can for instance take
1. If then the chord in Figure 4.2 tends to a tangent line and, because of
Theorem 4.4.1, if then the ruin probability equals 1, see
Exercise 4.4.1.
4. It is quite plausible that the denominator of (4.18) has a finite limit for
say Then, of course, This yields the following asymptotic
approximation for for large we have
Notice that Lundberg’s exponential upper bound boils down to an equality here
except for the constant In this case, the denominator of (4.18) does not
depend on In general, however, it will depend on
In the discrete time model, we consider more general risk processes than the
compound Poisson process from the previous sections, but now only on the time
points 0,1,2,… Instead of we write Let denote the
profit between the time points and therefore
The last equation has a unique solution. This can be seen as follows: since E [G] > 0
and Pr[G < 0] > 0, we have and for while
so is a convex function.
Example 4.5.1 (Compound Poisson distributed annual claims)
In the special case that is a compound Poisson process, we have
where denotes the compound Poisson distributed total claims in year From
(4.11), we know that R satisfies the equation Hence,
Example 4.5.2 (Normally distributed annual claims)
If with then follows from:
92 RUIN THEORY
Combining this result with the previous example, we observe the following. If we
consider a compound Poisson process with a large Poisson parameter, i.e., with
many claims between the time points 0 , 1 , 2 , . . . , then will approximately follow
a normal distribution. Consequently, the adjustment coefficients will be close to
each other, so On the other hand, if we take
and in Exercise 4.3.2, then it turns out that is an upper bound for
R.
Analogously to Theorem 4.4.1, one can prove the following equality:
So in the discrete time model one can give an exponential upper bound for the ruin
probability, too, which is
Now assume that we take out a stop-loss reinsurance with For a reinsurance
with payment Y, the reinsurer asks a premium with If
the reinsurance premium amounts to
To determine the adjustment coefficient, we calculate the distribution of the profit
in one year which consists of the premium income minus the reinsurance
premium minus the retained loss. Hence,
where denotes the reinsurance premium. The reinsurer uses a loading factor
on the net premium. Assume that and for and
Furthermore, let and consider two values and
In case of proportional reinsurance the premium equals
In the table below, we give the results for different values of compared
with the same results in case of proportional reinsurance with the same expected
payment by the reinsurer:
For the loading factors of the reinsurer and the insurer are equal,
and the more reinsurance we take, the larger the adjustment coefficient is. If
the reinsurer’s loading factor equals then for the expected retained loss
is not less than the retained premium
Consequently, the resulting retained loading factor is not positive, and eventual ruin
is a certainty. The same phenomenon occurs in case of excess of loss reinsurance
with In the table below, this situation is denoted by the symbol *.
From the table we see that all adjustment coefficients for excess of loss coverage
(XL) are at least as large as those for proportional reinsurance (Prop.) with the
same expected payment. This is not a coincidence: by using the theory on ordering
of risks, it can be shown that XL coverage always yields the best R-value as well
as the smallest ruin probability among all reinsurance contracts with the same
expected value of the payment, see Example 10.4.4.
In this section we show that the non-ruin probability can be written as a com-
pound geometric distribution function. For this purpose, we consider the maximal
aggregate loss, i.e., the maximal difference between the payments and the earned
premium up to time
Since S(0) = 0, we have The event occurs if, and only if, a finite
point in time exists for which In other words, the inequalities
96 RUIN THEORY
Next, we consider the points where the surplus process reaches a new record low.
This happens necessarily at points in time when a claim is paid. Let the random
variables denote the amounts by which the record low is
less than the – 1st one, see Figure 4.3 where there are three new record lows,
assuming that the process drifts away to in the time period not shown. Let M
be the random number of new records. We have
From the fact that a Poisson process is memoryless, it follows that the probability
that a particular record low is the last one is the same every time. Hence, M follows
a geometric distribution. For the same reason, the amounts of the improvements
are independent and identically distributed. The parameter of M, i.e.,
the probability that the previous record is the last one, equals the probability to
avoid ruin starting with initial capital 0, hence it equals
So L has a compound geometric distribution. Both the value of the geometric
parameter and the distribution of conditionally given follow
from the following theorem:
BEEKMAN'S CONVOLUTION FORMULA 97
we can write
Substitute (4.43) into (4.42), subtract from both sides and divide by
Then we get
The double integrals in (4.45) can be reduced to single integrals as follows. For
the first double integral, exchange the order of integration, substitute
and again exchange the integration order. This leads to
Hence,
2. Assuming that there is at least one new record low, has the same distri-
bution as the amount with which ruin occurs starting from (if ruin
BEEKMAN’S CONVOLUTION FORMULA 99
occurs). So we have the following expression for the density function of the
record improvements:
where
and
4. The mgf of the maximal aggregate loss L which because of (4.38) is also
the mgf of the non-ruin probability is given by
To prove this, note that implies that the surplus eventually will drop below
the initial level, so
Two situations exist for which we can give expressions for the ruin probabilities.
In case of exponential distributions, and mixtures or combinations of these, an
analytical expression arises. For discrete distributions, we can derive an algorithm.
In the previous section, we derived the mgf with the non-ruin probability
In some cases, it is possible to identify this mgf, and thus give an
expression for the ruin probability. We will describe how this works for mix-
tures and combinations of two exponential distributions, see Section 3.7. Since
and
EXPLICIT EXPRESSIONS FOR RUIN PROBABILITIES 101
Note that, except for a constant, is a density function, see Exercise 4.8.1.
Now, if X is a combination or a mixture of two exponential distributions as in
(3.27), i.e., for some and it has density function
then the right-hand side of (4.61), after multiplying both the numerator and the
denominator by can be written as the ratio of two polynomials
in By using partial fractions, this can be written as a sum of terms of the form
corresponding to times an exponential distribution. We give two
examples to clarify this method.
Except for the constant this is the mgf of an exponential distribution. We con-
clude from (4.61) that – is equal to the density function of this distribution.
By using the boundary condition we see that for the exponential(1)
distribution
for
Notice that the right-hand side of this equation corresponds to the mgf of the
claims only if is to the left of the asymptotes, i.e., if If is larger, then
this mgf is hence we write instead of for these branches in
Figure 4.4. From this figure, one sees immediately that the positive roots and
are real numbers that satisfy
For other distributions than the ones above, it is difficult to calculate the exact
value of the ruin probability Furthermore, one may argue that this exact
value is not very important, since in case of doubt, other factors will be decisive.
So there is a need for good and simple approximations for the ruin probability.
First of all, we give some global properties of the ruin probability that should
preferably be satisfied by the approximations. Equation (4.50) yields
Since we have
hence
To fit the first two moments, the parameters and of the gamma cdf
must meet the following conditions:
4. Just as in the first approximation, one can replace the claims distribution
by another with a few moments in common, for which the corresponding
ruin probability can be easily calculated. A suitable candidate for such a
replacement is a mixture or combination of exponential distributions.
5. Another possible replacement is a discrete distribution. The ruin probabili-
ties can easily be computed from (4.70). For each claims distribution, one
can find a two-point distribution with the same first three moments. This
is not always possible in case of a mixture/combination of two exponential
distributions. Both methods yield good approximations.
6. From the theory of ordering of risks, it follows that one gets a lower bound for
the ruin probability if one replaces the claims distribution with expectation
by a one-point distribution on A simple upper bound can be obtained
if one knows the maximum value of the claims. If one takes a claims
distribution with probability for and probability for 0, then a
Poisson process arises which is equivalent to a Poisson process with claims
always equal to and claim number parameter instead of So, both
the lower bound and the upper bound can be calculated by using (4.70) with
4.10 EXERCISES
Section 4.2
1. Assume that the waiting times are independent and identically distributed random
variables with cdf and density function Given and for
some what is the conditional probability of a claim occurring between points in time
and (This generalization of a Poisson process is called a renewal process.)
2. Let be a Poisson process with parameter and let and
Show that and interpret these
formulas by comparing with
Section 4.3
Section 4.4
1. From Theorem 4.4.1, we know that if Why does this imply that if
3. For a compound Poisson process, it is known that the continuous ruin probability depends on
the initial capital in the following way: Determine the adjustment
coefficient for this process. Can anything be said about the Poisson parameter in this risk process?
What is
4. Assume that By looking at the event "non-ruin & no claim before with
denoting the premium income per unit of time, show that must hold.
5. For a certain risk process, it is given that and Which of the
numbers 0, 1 and 6 are roots of the adjustment coefficient equation
Which one is the real adjustment coefficient?
One of the four expressions below is the ruin probability for this process; determine which
expression is the correct one, and argue why the other expressions can’t be the ruin probability.
Section 4.5
Section 4.6
1. The claim process on some insurance portfolio is compound Poisson with and
The loading factor is Calculate the adjustment coefficient in case one takes out a
proportional reinsurance with a loading factor Calculate the relative loading
factor after this reinsurance. Which restrictions apply to
108 RUIN THEORY
2. For the same situation as in the previous exercise, but now with excess of loss coverage
write down the adjustment coefficient equation, and determine the loading factor after
reinsurance.
3. Assume that the claims per year are R(5,1) distributed and that A
reinsurer covers a fraction of each risk, applying a premium loading factor Give the
adjustment coefficient for the reinsured portfolio, as a function of Which value optimizes
the security of the insurer?
4. A total claims process is compound Poisson with and The relative
loading factor is One takes out a proportional reinsurance The relative
loading factor of the reinsurer equals 1. Determine the adjustment coefficient For which
values of is ruin not a certainty?
Section 4.7
1. What is the mgf of if the claims (a) are equal to with probability 1, and (b) have an
exponential distribution?
2. Prove that
3. In Exercises 4.4.3 and 4.4.6, what is
Section 4.8
11. The ruin processes of company 1 and 2 are both compound Poisson with intensities
and claims distributions exponential(3) and exponential(6), and loading factors
and The claims process of company 1 is independent of the one of company 2.
These companies decide to merge, without changing their premiums. Determine the intensity,
claims distribution and loading factor of the ruin process for the merged company. Assume
that both company 1 and 2 have an initial capital equal to 0, then obviously so does the
merged company. Compare the probabilities of the following events (continuous infinite ruin
probabilities): “both companies never go bankrupt” with “the merged company never goes
bankrupt”. Argue that, regardless of the values of the initial capitals and for the separate
companies, and consequently of for the merged company, the following holds: the event
“both companies never go bankrupt” has a smaller probability than “the merged company never
goes bankrupt”.
Section 4.9
1. Verify (4.72), (4.73), (4.75)/(4.76), and (4.80). Solve and from (4.80).
2. Work out the details of the final approximation.