The Compound Poisson risk model can be modeled by the quotation,
where and are independent. : amount of ith claim, and . u: initial surplus at time . c: rate of premium income per unit time. An additional assumption: and thus
where , loading coefficient and is the rate of compound Poisson Process.
In this paper, the theoretical assumptions are: 1) In a Compound Poisson model, the ruin probability with initial reserve satisfies,
and 2) In a compound Poisson model with exponential claim,
where and have estimated using simulation. The compound Poisson risk model has been estimated.
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 , loading coefficient , initial surplus 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.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.
It has repeated to construct more confidence interval for each .
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 and be defined by following model,
For simulation purpose set
where d is such that if and only if
Setting X and Y such that and . This setting is motivated from reinsurance problem. Let Z is the original risk and f, g are different reinsurance strategies.
Gatto and Mosimann  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 , 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 . 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 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 . The compound Poisson risk model modeled by (1), where is the amount of i-th claims with is the amount of claims generated by portfolio in time where is a Poisson process. We have and are independent. Again u is initial surplus at time . The ruin probability of compound Poisson surplus process is denoted by
We have an additional assumption that . Then c can be expressed as
where is a loading coefficient and is the rate of compound Poisson process. If , the ruin occurs with probability 1 which is generally not interested by insurer  . Intuitively, that means, averagely the premium should exceed the claim rate on one unit time period. We have from Wald’s Identity,
which can be proved by conditioning .
which implies that and consequently,
3. Theoretical Results
Theorem 3.1. In a compound Poisson model defined by (1), the ruin probability with no initial reserve satisfies
Theorem 3.2. In a compound Poisson model defined by (1) with
The purpose of simulation is to estimate infinite ruin probabilities at and compare to its theoretical value defined by (15). The steps are repeated for different u value to develop .
First, it has focused on estimating
using simulation. In each run, the process is simulated up to time
and the result of the i-th simulation run
is set to 1 if ruin occurs at or before time
and to 0 if this is not the case. After
The following algorithm has used to obtain of a single simulation run.
Step 1: Initialize the simulation: set .
Step 2: Draw an exponential inter-arrival time I with parameter and draw a claim size .
Step 3: If return . If not, set and if return . Set . Otherwise, return to step 2.
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: , with initial surplus and claim size are exponentially distributed with mean 1. An approximate 95% confidence interval of 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 . 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 defined by theorem 8.1 which has demonstrated on Figure 2.
Figure 1. Ruin probabilities plot.
5. Ruin Probabilities Comparison of Two Compound Poisson Model with Different Claims
Definition: Let X and Y be two random variable where X is smaller than Y in convex order denoted by such that
for all convex function , provided the expectation exists.
6. Theoretical Results
Theorem 6.1. For two random variable X and Y defined from (18)
provided that the expectation exist  .
Theorem 6.2. For two random variable X and Y defined from (18)
provided that expectation exist  .
The above result implies that means Y has more dispersion. The above result is important because this guarantees that the premium rate,
which makes easier to compare the ruin probabilities. A typical way to construct X and Y following order has seen from Ohlin’s Lemma.
7. Ohlin’s Lemma
Lemma 7.1. Let Z be a positive random variable and f and g be increasing functions. If g crosses f only once from below, i.e., there exist , such that
for all , and (22)
then  .
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 .
8. Compound Poisson Risk Models with Different Claims
It can be defined two compound Poisson risk model with different claims by (5) and by (6). Again and 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.
Theorem 8.1. Let and be defined in (4) and (5) respectively with generic claim size X and Y. If , then  . (27)
It has interested to see the above theoretical results by using simulation. The theoretical result has observed in simulation for different values. For simulation purpose we can set,
with where such that
Now compared with first simulation, X becomes and Y becomes . From Ohlin’s Lemma, we have .
The practical background of
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
is applicable on quota share (proportional insurance). Graph of has estimated as Figure 2.
Figure 2. Ruin probabilities plot for and .
The infinite ruin probabilities of has been estimated and observed. The infinite ruin probability at has been observed by using simulation. It has
found that is equal to its theoretical value, . Then the
comparison of and has estimated for two-compound Poisson risk models. It has been verified that for different u values. Eventually, the theorem defined by (27) has been verified from estimations.
 Gatto, R. and Mosimann, M. (2012) Four Approaches to Compute the Probability of Ruin in Compound. Poisson Risk Process with Diffusion. Mathematical and Computer Modelling, 55, 1169-1185.