Compound Poisson risk model has been simulated. It has started with exponential claim sizes. The simulations have checked for infinite ruin probabilities. An appropriate time window has been chosen to estimate and compare ruin probabilities. The infinite ruin probabilities of two-compound Poisson risk process have estimated and compared them with standard theoretical results.

It has assumed that the probability of ruin at time window of T = 1000 after 100,000 independent runs which has set just to make algorithm convenience. For computational convenience, it has chosen the rate
λ
=
0.5
, loading coefficient
θ
=
0.1
, initial surplus
u
=
0
and claim size which are exponentially distributed with mean 1. After simulation, an approximate 95% confidence interval based on 1,000,000 independent runs, and the confidence interval of
Ψ
(
0
)
: [0.9083541, 0.9094819] has obtained.

It helped to claim that the simulation which has done so far has very high accuracy since it is approximately equal to theoretical value of 0.909 from above given condition.

The time window and independent runs are changed as required. It has also estimated infinite ruin probabilities of two-compound Poisson Process and compared them. Let
U
X
(
t
)
and
U
Y
(
t
)
be defined by following model,

U
X
(
t
)
=
u
+
c
X
t
−
∑
i
=
1
N
(
t
)
X
i
(5)

and

U
Y
(
t
)
=
u
+
c
Y
t
−
∑
i
=
1
N
(
t
)
Y
i
(6)

where

{
X
i
}
i
i
d
︸
X
and
{
Y
i
}
i
i
d
︸
Y
and
i
=
1
,
2
,
3
,
⋯

where

c
X
=
(
1
+
θ
)
λ
E
[
X
]
and
c
Y
=
(
1
+
θ
)
λ
E
[
Y
]
(7)

For simulation purpose set

X
=
min
(
Z
,
d
)
(8)

and

Y
=
1
2
Z
(9)

with

Z
~
exp
(1)

where d is such that
E
[
X
]
=
E
[
Y
]
if and only if

(
1
−
e
−
d
)
=
1
2
(10)

Then

d
=
ln
2
(11)

Setting X and Y such that
X
=
f
(
Z
)
and
Y
=
g
(
Z
)
. This setting is motivated from reinsurance problem. Let Z is the original risk and f, g are different reinsurance strategies.

Gatto and Mosimann [2] have solved for four different approaches to compute probability of ruin of insurer compound Poisson surplus process with an additive Wiener perturbation. There is no ideal comparison of ruin probabilities which have supported by theoretical results (15), (16) and (27). In this paper, infinite ruin probabilities of two-compound Poisson risk process have estimated and compared them with standard theoretical results.

2. Compound Poisson Risk Model

Suppose that at time
t
=
0
, the insurer has an amount of money set aside for the portfolio. This amount of money is called the initial surplus and is denoted by u. We can further assume that throughout this work 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. We are also ignoring the expenses. The insurer’s surplus at any 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 compound Poisson risk model modeled by (1), where
{
X
i
}
i
i
d
︸
X
is the amount of i-th claims with
i
=
1
,
2
,
3
,
⋯
,
N
(
t
)
is the amount of claims generated by portfolio in time
[
0
,
t
]
where
N
(
t
)
is a Poisson process. We have
N
(
t
)
and
X
i
are independent. Again u is initial surplus at time
t
=
0
. The ruin probability of compound Poisson surplus process is denoted by

Ψ
(
u
)
=
P
[
inf
t
>
0
U
(
t
)
<
0
]
(12)

We have an additional assumption that
c
>
λ
E
[
X
]
. Then c can be expressed as

c
=
(
1
+
θ
)
λ
E
[
X
]
(13)

where
θ
>
0
is a loading coefficient and
λ
is the rate of compound Poisson process. If
c
<
λ
E
[
X
]
, the ruin occurs with probability 1 which is generally not interested by insurer [1] . Intuitively, that means, averagely the premium should exceed the claim rate on one unit time period. We have from Wald’s Identity,

The purpose of simulation is to estimate infinite ruin probabilities at
Ψ
(
0
)
and compare
Ψ
(
0
)
to its theoretical value defined by (15). The steps are repeated for different u value to develop
Ψ
(
u
)
.

First, it has focused on estimating
Ψ
(
u
,
τ
)
using simulation. In each run, the process is simulated up to time
τ
and the result of the i-th simulation run
Z
i
is set to 1 if ruin occurs at or before time
τ
and to 0 if this is not the case. After
runs, the Monte Carlo estimator is given by,

Z
=
1
/
n
∑
i
=
1
n
X
i
(17)

The following algorithm has used to obtain of a single simulation run.

Step 1: Initialize the simulation: set
U
(
0
)
=
u
.

Step 2: Draw an exponential inter-arrival time I with parameter
λ
and draw a claim size
X
~
F
.

For computational convenience, it has assumed that the probability of ruin time at time window of T = 1000 after 1,000,000 independent runs which are assigned just to make algorithm convenience. The parameters which are defined in (1) are chosen:
λ
=
0.5
,
θ
=
0.1
with initial surplus
u
=
0
and claim size are exponentially distributed with mean 1. An approximate 95% confidence interval of
Ψ
(
0
)
based on 1,000,000 independent runs, is given by [0.9083541, 0.9094819].

We expect to see that the theoretical curve lies in between the confidence intervals. The plot from above procedure has expected to see a good estimate. However, it is not the case always since the time window chosen may not be the sufficient to provide an accurate estimate. The simulation has repeated to see different confidence interval for each
Ψ
(
u
)
=
[
Ψ
L
(
u
)
,
Ψ
U
(
u
)
]
. Then four different ruin probabilities plots are estimated from simulation (Figure 1).

One main issue is that the Monte Carlo simulation used doesn’t work for small ruin probabilities. It could be the reason the graph doesn’t display the probabilities of ruin smaller than 0.15. We need to test for more surplus with sufficient time window to address this issue. For this, important sampling or other method is required. But the scope of this paper is to extend this estimation to compare ruin probabilities
Ψ
X
(
u
)
≤
Ψ
Y
(
u
)
defined by theorem 8.1 which has demonstrated on Figure 2.

The above setting is motivated from reinsurance problem. We can think Z as the original risk and f, g are different insurance strategies. In this context, besides increasing, f and g are assumed to be continuous and satisfy
0
≤
f
(
x
)
≤
x
.

8. Compound Poisson Risk Models with Different Claims

It can be defined two compound Poisson risk model with different claims
U
X
(
t
)
by (5) and
U
Y
(
t
)
by (6). Again
c
X
and
c
Y
which are defined in (7) are respective rate of premium income per unit time of the model defined by (5) and (6). We are interested in two quantities: the ruin time and infinite ruin probability of two different models defined by (5) and (6) respectively.

is the excess of loss insurance form. For example, customer’s loss in a car accident with an auto insurance where only the part in excess of certain level is covered by insurance company. The maximum loss for the customer has capped at d. Again the setting