CT 6-Ruin Theory
CT 6-Ruin Theory
Chapter 9
Ruin theory
Syllabus objectives
(iv) Explain the concept of ruin for a risk model. Calculate the adjustment coefficient
and state Lundberg’s inequality. Describe the effect on the probability of ruin of
changing parameter values and of simple reinsurance arrangements.
1. Explain what is meant by the aggregate claim process and the cash-flow
process for a risk.
3. Define a compound Poisson process and derive the moments and moment
generating function for such a process.
7. Describe the effect on the probability of ruin, in both finite and infinite
time, of changing parameter values.
0 Introduction
In the previous two chapters we used the collective risk model to look at the aggregate
claims S arising during a fixed period of time. S was given by the equation
S = X1 + X 2 + + X N , where N denotes the number of claims arising during the period.
In this chapter we will extend this model by treating S (t ) as a function of time. This
gives us the equation S (t ) = X1 + X 2 + + X N (t ) , where N (t ) denotes the number of
claims occurring before time t. N t is called a Poisson process and S t is called a
compound Poisson process. We can use S t to model claims received by an insurance
company and hence consider the probability that this insurance company is ruined. The
notation and the other basic concepts are covered in Section 1.
In Section 2 we will give formal definitions of both the Poisson process and the
compound Poisson process. You may have already met these in Subject CT4. We will
also introduce the concept of a premium security loading. Briefly, this is an additional
amount charged on an insurance premium to reduce the likelihood of an insurance
company becoming ruined.
Section 4 considers the effect of changing parameter values on the probability of ruin
for an insurance company. Finally, in Section 5 we will consider the impact on the
probability of ruin for an insurance company when reinsurance is introduced into the
equation.
1 Basic concepts
1.1 Notation
f (x)
lim =0
x 0 x
You can use this notation to simplify your working. For example, the function
2 3
g ( x) = 3 x + 0.5 x + 0.004 x can be rewritten as g ( x ) = 3 x + o( x ) , since
0.5 x 2 + 0.004 x3
0 as x 0 . Note that o( x) does not represent an actual number so
x
that co( x) ( c is a constant), -o( x) and o( x) are all equivalent.
Question 9.1
For the purposes of this course, you just need to be able to understand the notation.
However, if you wish to know more, then it’s covered in more detail in the Foundation
ActEd Course (FAC).
S (t ) the aggregate claims in the time interval [0, t], for all t 0.
You can think of a stochastic process as being a whole family of different random
variables. Consider a time line. On the line there are an infinite number of different
time intervals. For each interval of time, there is a random variable that corresponds to
the aggregate claim amount arising in that time interval. This is what we mean by a
stochastic process.
Of course, if you have previously studied Subject CT4, you should be familiar with
these ideas already.
N (t )
S (t ) = Â Xi
i =1
So we have just taken the idea of a compound distribution and generalised it to cover
different time periods.
The insurer of this portfolio will receive premiums from the policyholders. It is
convenient at this stage to assume, as will be assumed throughout this chapter,
that the premium income is received continuously and at a constant rate. Here is
some more notation:
so that the total premium income received in the time interval [0, t] is ct. It will
also be assumed that c is strictly positive.
Suppose that at time 0 the insurer has an amount of money set aside for this
portfolio. This amount of money is called the initial surplus and is denoted by U.
It will always be assumed that U 0. The insurer needs this initial surplus
because the future premium income on its own may not be sufficient to cover the
future claims. Here we are ignoring expenses. The insurer’s surplus at any
future time t (> 0) is a random variable since its value depends on the claims
experience up to time t. The insurer’s surplus at time t is denoted by U(t). The
following formula for U(t) can be written:
U (t ) = U + ct - S (t )
In words this formula says that the insurer’s surplus at time t is the initial
surplus plus the premium income up to time t minus the aggregate claims up to
time t. Notice that the initial surplus and the premium income are not random
variables since they are determined before the risk process starts. The above
formula is valid for t 0 with the understanding that U(0) is equal to U. For a
given value of t, U(t) is a random variable because S(t) is a random variable.
Hence U (t ) t 0 is a stochastic process, which is known as the cash flow
process or surplus process.
Figure 1
Figure 1 shows one possible outcome of the surplus process. Claims occur at
times T1, T2, T3, T4 and T5 and at these times the surplus immediately falls by
the amount of the claim. Between claims the surplus increases at constant rate
c per unit time. The model being used for the insurer’s surplus incorporates
many simplifications, as will any model of a complex real-life operation. Some
important simplifications are that it is assumed that claims are settled as soon as
they occur and that no interest is earned on the insurer’s surplus. Despite its
simplicity this model can give an interesting insight into the mathematics of an
insurance operation.
We are also assuming that there are no expenses associated with the process (or,
equivalently, that S (t ) makes allowance for expense amounts as well as claim
amounts), and that the insurer cannot vary the premium rate c .
We are also ignoring the possibility of reinsurance. Simple forms of reinsurance will be
incorporated into the model later in this chapter.
It can be seen from Figure 1 that the insurer’s surplus falls below zero as a result
of the claim at time T3. Speaking loosely for the moment, when the surplus falls
below zero the insurer has run out of money and it is said that ruin has occurred.
In this simplified model, the insurer will want to keep the probability of this
event, that is, the probability of ruin, as small as possible, or at least below a
predetermined bound. Still speaking loosely, ruin can be thought of as meaning
insolvency, although determining whether or not an insurance company is
insolvent is, in practice, a very complex problem. Another way of looking at the
probability of ruin is to think of it as the probability that, at some future time, the
insurance company will need to provide more capital to finance this particular
portfolio.
(U ) = P [U (t ) 0, for some t , 0 t ]
(U , t ) = P [U ( ) 0, for some , 0 t]
(U 2 , t ) (U1 , t ) (1.1)
(U 2 ) (U1 ) (1.2)
(U , t 1 ) (U , t 2 ) (U ) (1.3)
lim (U , t ) = (U ) (1.4)
t
The larger the initial surplus, the less likely it is that ruin will occur either in a
finite time period, hence (1.1), or an unlimited time period, hence (1.2).
For a given initial surplus U, the longer the period considered when checking for
ruin, the more likely it is that ruin will occur, hence (1.3).
Question 9.2
You may be wondering whether it is possible to find numerical values for these ruin
probabilities. In some very simple cases it is. However, for most practical situations,
finding an exact value for the probability of ruin is impossible. In some cases there are
useful approximations to (u) , even if calculation of an exact value is not possible.
The two probabilities of ruin considered so far have been continuous time
probabilities of ruin, so-called because they check for ruin in continuous time. In
practice it may be possible (or even desirable) to check for ruin only at discrete
intervals of time.
For a given interval of time, denoted h, the following two discrete time
probabilities of ruin are defined:
Figure 2
It can be seen from Figure 2 that in discrete time with h = 1, ruin does not occur
for this realisation of the surplus process before time 5, but ruin does occur (at
time 2½) in discrete time with h = ½.
Listed below are five relationships between different discrete time probabilities
of ruin for 0 U1 U2 and for 0 t1 t2 < . Formulae (1.5), (1.6), (1.7) and (1.8)
are the discrete time versions of formulae (1.1), (1.2), (1.3) and (1.4) above and
their intuitive explanations are similar. The intuitive explanation of (1.9) comes
from Figure 2.
h (U 2 , t ) h (U1 , t ) (1.5)
h (U 2 ) h (U1 ) (1.6)
h (U , t 1 ) h (U , t 2 ) h (U ) (1.7)
lim h (U , t ) = h (U ) (1.8)
t
h (U , t ) (U , t ) (1.9)
Question 9.3
Intuitively, it is expected that the following two relationships are true since the
probability of ruin in continuous time could be approximated by the probability
of ruin in discrete time, with the same initial surplus, U, and time horizon, t,
provided ruin is checked for sufficiently often, ie provided h is sufficiently small.
lim h (U , t ) = (U , t ) (1.10)
h 0+
lim h (U ) = (U ) (1.11)
h 0+
Formulae (1.10) and (1.11) are true but the proofs are rather messy and will not
be given here.
2.1 Introduction
In this section some assumptions will be made about the claim number process,
lN(t )qt 0
, and the claim amounts, lX i q i 1
. The claim number process will be
assumed to be a Poisson process, leading to a compound Poisson process
l q
S(t ) t 0 for aggregate claims. The assumptions made in this section will hold
for the remainder of this chapter.
If you have previously studied Subject CT3, the material in this section should be
familiar to you.
Recall from Subject CT3 that we use the term “Poisson process” to describe the number
of claims arising from a time period of length t.
If the number of claims N arising from a single time period has a Poisson distribution
with parameter l then the number of claims N (t ) which arise over a time period of
length t is a Poisson process, ie N (t ) has a Poisson distribution with parameter l t .
(ii) P N (t + h ) = r | N (t ) = r = 1 - l h + o(h )
P N (t + h ) = r + 1| N (t ) = r = l h + o(h ) (2.1)
P N (t + h) r + 1| N (t ) = r = o(h )
Condition (ii) implies that there can be a maximum of one claim in a very short time
interval h. It also implies that the number of claims in a time interval of length h does
not depend on when that time interval starts.
Example
Solution
The events in this case are occurrences of claim events (ie accidents, fires, thefts etc) or
claims reported to the insurer. The parameter represents the average rate of
occurrence of claims (eg 50 per day), which we are assuming remains constant
throughout the year and at different times of day. The assumption that, in a sufficiently
short time interval, there can be at most one claim is satisfied if we assume that claim
events cannot lead to multiple claims (ie no motorway pile-ups etc).
When studying a Poisson process the distribution of the time to the first claim
and the times between claims is often of particular interest.
This section will show that the time to the first claim has an exponential distribution
with parameter l .
Let the random variable T1 denote the time of the first claim. Then, for a fixed
value of t , if no claims have occurred by time t , T1 t . Hence:
P (T1 t ) = P (N (t ) = 0) = exp{ - l t }
This last step follows from the formula for the probability function of a Poisson( t)
distribution with x = 0.
And:
P (T1 t) = 1 exp{ t}
Important information
The time to the first claim in a Poisson process has an exponential distribution with
parameter l .
This section will show that the time between claims has an exponential distribution with
parameter l .
For i = 2,3, , let the random variable Ti denote the time between the (i - 1) th
and the i th claims. Then:
n n +1 n
P (Tn +1 t| Â Ti = r ) = P ( Â Ti t +r | Â Ti = r )
i =1 i =1 i =1
= P (N (t + r ) = n | N (r ) = n )
= P (N (t + r ) - N (r ) = 0 | N (r ) = n )
By condition (2.2) (ie using the independence of claim numbers in different time
periods):
P (N (t + r ) - N (r ) = 0 N (r ) = n ) = P (N (t + r ) - N (r ) = 0)
Finally:
P (N (t + r ) - N (r ) = 0) = P (N (t ) = 0) = exp{ - l t }
since the number of claims in a time interval of length r does not depend on
when that time interval starts (condition (2.1)). Thus inter-event times also have
an exponential distribution with parameter .
Important information
The time between claims in a Poisson process has an exponential distribution with
parameter l .
Note that the inter-event time is independent of the absolute time. In other words the
time until the next event has the same distribution, irrespective of the time since the last
event or the number of events that have already occurred. This is referred to as the
memoryless property of the exponential distribution.
Question 9.4
If reported claims follow a Poisson process with rate 5 per day (and the insurer has a
24 hour hotline), calculate:
(i) the probability that there will be fewer than 2 claims reported on a given day
(ii) the probability that another claim will be reported during the next hour.
In this section the Poisson process for the number of claims will be combined
with a claim amount distribution to give a compound Poisson process for the
aggregate claims.
This last assumption means that for any t 0 , the random variable N (t ) has a
Poisson distribution with parameter t , so that:
( t )k
P [N (t ) k] exp{ t} for k 0,1, 2,
k!
x
F ( x) = Ú f (t ) dt
-
mk = E [ X ik ] for k = 1, 2, 3,
Important information
mk = E[ X ik ] for k = 1, 2, 3,
Whenever the common moment generating function of the X i s exists, its value
at the point r will be denoted by M X (r ) .
In case you have forgotten, the definition of a moment generating function is:
M X (r ) = E[e rX ]
Note that we are using r for the dummy variable to avoid confusion with time t .
MS (r )= exp {l t (M X (r ) - 1)}
S = X1 + X 2 + + XN
E[ S ] = l E[ X ] = l m1 var[ S ] = l E[ X 2 ] = l m2
M S (r ) = M N [ln M X (r )]
c > m1 (2.3)
so that the insurer’s premium income (per unit time) is greater than the expected
claims outgo (per unit time).
Question 9.5
If we know the distribution of the aggregate claims S (t ) , we can often determine the
probability of ruin for the discrete model over a finite time horizon directly (without
reference to the models), by looking at the cashflows involved.
Example
The aggregate claims arising during each year from a particular type of annual insurance
policy are assumed to follow a normal distribution with mean 0.7 P and standard
deviation 2.0 P , where P is the annual premium. Claims are assumed to arise
independently. Insurers assess their solvency position at the end of each year.
A small insurer with an initial surplus of £0.1m expects to sell 100 policies at the
beginning of the coming year in respect of identical risks for an annual premium of
£5,000. The insurer incurs expenses of 0.2 P at the time of writing each policy.
Calculate the probability that the insurer will prove to be insolvent at the end of the
coming year. Ignore interest.
Solution
Using the information given, the insurer’s surplus at the end of the coming year will be:
= P( N [0.35m,(0.1m) 2 ] 0.5m)
Ê 0.5m - 0.35m ˆ
= 1- ÁË ˜¯ = 1 - 0.93319 = 0.067
0.1m
Question 9.6
If the insurer expects to sell 200 policies during the second year for the same premium
and expects to incur expenses at the same rate, calculate the probability that the insurer
will prove to be insolvent at the end of the second year.
In fact, the normal distribution is probably not a very realistic distribution to use for the
claim amount distribution in most portfolios, as it is symmetrical, whereas many claim
amount portfolios will have skewed underlying distributions. However, it is commonly
used in CT6 exam questions.
Example
The number of claims from a portfolio of policies has a Poisson distribution with
parameter 30 per year. The individual claim amount distribution is lognormal with
parameters m = 3 and s 2 = 1.1 . The rate of premium income from the portfolio is 1,200
per year.
If the insurer has an initial surplus of 1,000, estimate the probability that the insurer’s
surplus at time 2 will be negative, by assuming that the aggregate claims distribution is
approximately normal.
Solution
First we need the mean and variance of the aggregate claims in a two-year period. The
expected number of claims will be 60. So the mean and variance are (using the formulae
for the first two moments of a lognormal distribution):
Ruin will occur if S (2) is greater than the initial surplus plus premiums received. So we
want:
3, 400 - 2, 088.80
P S (2) 1, 000 + 2 1, 200 P N (0,1) = 1 - (2.8053) = 0.0025
218, 457
So far, we have used c to denote the rate of premium income per unit time, independent
of the claims outgo. In some circumstances it is more useful to think of the rate of
premium income as being related to the rate of claims outgo.
For the insurer to survive, the rate at which premium income comes in needs to be
greater than the rate at which claims are paid out. If this is not true, the insurer is
certain to be ruined at some point.
c (1 ) m1
The security loading is the percentage by which the rate of premium income exceeds the
rate of claims outgo. So, for the Poisson process outlined above, we have:
c = (1 + ) E ( S ) = (1 + )l m1
where is the security loading. is also sometimes called the “relative security
loading”. It might typically be a figure such as 0.2, ie 20%.
The insurer will need to adopt a positive security loading when pricing policies, in order
to cover expenses, profit, contingency margins and so on.
Note that this does not mean that ruin is impossible. It is quite possible for the actual
claims outgo to exceed substantially its expected value. So even in this situation the
insurer’s probability of ruin is non-zero.
Question 9.7
In summary:
For a compound Poisson process S (t ) , the mean and variance of the total claim amount
are given by:
E[ S (t )] = l t E ( X ) var[ S (t )] = l t E ( X 2 )
M S (t ) (r ) = exp l t ( M X (r ) - 1
2.6 A technicality
lim M X r = (2.4)
-
r
(For example, if the X i s have a range bounded above by some finite number,
then will be ; if the X i s have an exponential distribution with parameter ,
then will be equal to .)
Suppose for example that claim amounts have a continuous uniform distribution on the
interval (0,10) , so that they are bounded above by 10. Then the moment generating
function of the claim distribution is (from the Tables):
e10r 1
M X (r )
10r
This is defined for all positive values of r , and so in this case . We can see that
as r , the limit of the MGF is infinite. If the claim distribution is Exp( ) , the
MGF (as stated in the Tables) is:
1
M X (r ) (1 r / )
r
lim (l M X (r ) - cr ) = (2.5)
-
r
This is because , c and r would all have finite values in the limit.
Now it will be shown that (2.5) holds when is infinite. This requires a little
more care. First note that there is a positive number, say, such that:
P[ X i ] 0
The reason for this is that all claim amounts are positive.
So, if we pick a small enough number ( 0.01 maybe), we’re bound to get a
proportion of claims whose amount exceeds this.
M X (r ) er
M X (r ) = E (e rX ) = E (e rX | X ) P( X ) + E (e rX | X ) P( X )
0 + er p
Hence:
lim (l M X (r ) - cr ) lim (l e r p - cr ) =
r r
Important information
lim (l M X (r ) - cr ) =
r
(U ) exp{ RU }
The proof of Lundberg’s inequality is not on the syllabus for Subject CT6.
Don’t worry at this stage about what R actually represents. It will be defined shortly.
Until then just think of it as a parameter associated with the surplus process.
Note that if we can find a value for R , then Lundberg’s inequality tells us that we can
find an upper bound for the probability of ruin. This is a very useful result.
Figure 1 shows a graph of both exp{ -RU } and (U ) against U when claim
amounts are exponentially distributed with mean 1, and when the premium
loading factor is 10%. (The solution for R will be found in Section 3.2. The
formula for (U ) is given in Section 4.2.)
10
In fact it can be shown that the value of R in this case is .
11
It can be seen that, for large values of U , (U ) is very close to the upper bound,
so that (U ) exp{ RU } .
Figure 1
R can be interpreted as measuring risk. The larger the value of R , the smaller
the upper bound for (U ) will be. Hence, (U ) would be expected to decrease
as R increases. R is a function of the parameters that affect the probability of
ruin, and R ’s behaviour as a function of these parameters can be observed.
Note that R is an inverse measure of risk. Larger values of R imply smaller ruin
probabilities, and vice versa.
(i) the claim amount distribution is exponential with mean 10, and
Figure 2
Question 9.8
Note also that the value of R when claim amounts are exponentially distributed
is less than the value of R when all claim amounts are 10. Again, this result is
not surprising. Both claim amount distributions have the same mean, but the
exponential distribution has greater variability. Greater variability is associated
with greater risk, and hence a larger value of (U ) would be expected for the
exponential distribution, and a lower value of R . This example illustrates that R
is affected by the premium loading factor and by the characteristics of the
individual claim amount distribution. R is now defined and shown, in general, to
encapsulate all the factors affecting a surplus process.
The surplus process depends on the initial surplus, on the aggregate claims
process and on the rate of premium income. The adjustment coefficient is a
parameter associated with a surplus process which takes account of two of
these factors: aggregate claims and premium income. The adjustment
coefficient gives a measure of risk for a surplus process. When aggregate
claims are a compound Poisson process, the adjustment coefficient is defined in
terms of the Poisson parameter, the moment generating function of individual
claim amounts and the premium income per unit time.
l M X (r ) - l - cr = 0 (3.1)
l M X (R ) = l + cR (3.2)
Note that, although R relates to the aggregate claims, the MGF used in the definition is
for the individual claim amount.
It is probably not at all obvious to you at this stage why R is defined in this way. The
reason is bound up with the proof of Lundberg’s inequality, which you are not required
to know. Please accept the definition, so that you can find the value of R in simple
cases.
Note that equation (3.1) implies that the value of the adjustment coefficient
depends on the Poisson parameter, the individual claim amount distribution and
the rate of premium income. However, writing c (1 ) m1 gives:
M X (r ) 1 (1 )m1r
You can see from this equation that all the ’s have cancelled.
Important information
M X (r ) cr or M X (r ) = 1 + (1 + )m1r
Since the rate of claims outgo per unit time is m1 (the mean of a compound Poisson
distribution), then if a loading factor of is used, the rate of premium income will be
c (1 ) m1 .
Example
An insurer knows from past experience that the number of claims received per month
has a Poisson distribution with mean 15, and that claim amounts have an exponential
distribution with mean 500. The insurer uses a security loading of 30%. Calculate the
insurer’s adjustment coefficient and give an upper bound for the insurer’s probability of
ruin, if the insurer sets aside an initial surplus of 1,000.
Solution
M X (r ) = 1 + (1 + )m1r
Ê 1 ˆ
We have X ~ Exp Á ˜ , so that M X (r ) = (1 - 500r ) -1 (this comes from the Tables),
Ë 500 ¯
= 0.3 , and m1 = E[ X ] = 500 . Substituting these into the equation:
Rearranging:
r = 0.000462
(U ) e -0.000462 1,000
= 0.630