Next Article in Journal
Bibliometric Analysis of the Literature on Measuring Techniques for Manipulating Financial Statements
Previous Article in Journal
Evaluation of Changes on World Stock Exchanges in Connection with the SARS-CoV-2 Pandemic. Survival Analysis Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Tail Probability of the Discounted Aggregate Claims under Homogeneous, Non-Homogeneous and Mixed Poisson Risk Model

Schoool of Economics, University of Johannesburg, Johannesburg 2006, South Africa
*
Author to whom correspondence should be addressed.
Risks 2021, 9(7), 122; https://doi.org/10.3390/risks9070122
Submission received: 28 April 2021 / Revised: 31 May 2021 / Accepted: 2 June 2021 / Published: 30 June 2021

Abstract

:
In this paper, we derive a closed-form expression of the tail probability of the aggregate discounted claims under homogeneous, non-homogeneous and mixed Poisson risk models with constant force of interest by using a general dependence structure between the inter-occurrence time and the claim sizes. This dependence structure is relevant since it is well known that under catastrophic or extreme events the inter-occurrence time and the claim severities are dependent.

1. Introduction

In actuarial science, the discounted sum of losses over time can be considered as a randomly weighted sum of a sequence of independent and identically distributed (i.i.d.) random variables. Those independent and identically distributed random variables represent the claim amounts in successive development years, while the weights represent the discounted factors. Similar studies that have assumed independence between the inter-occurrence time and the forthcoming claim amounts include Taylor (1979), Delbaen and Haezendonck (1987), Waters (1983), Sundt and Teugels (1995), Cai (2002) and Yuen et al. (2006) who not only assumed independence between inter-occurrence time but have also derived the asymptotic tail probabilities associated with free interest risk model. In addition, Boogaert et al. (1988), Willmot (1989), Leveillé and Garrido (2001a, 2001b), Leveillé and Adékambi (2011, 2012) and Léveillé et al. (2010) have not only assumed the independence between the inter-occurrence time and the claim sizes but have also generalized the moments of the aggregate renewal sums. Leveillé and Hamel (2018) have used non-homogeneous and doubly compound Poisson process with stochastic force of interest and have derived the moments generating functions of the discounted aggregate claims. Moreover, they have found the inverse of the moments generating functions numerically and analytically. Geluk and Tang (2009) have studied the asymptotic behavior of the tail probabilities of a sum of real-value random variables. They have assumed that the distributions of the random variables belong to the subexponential class. Using asymptotically independence assumptions, they have showed that the asymptotic behavior of the tail probabilities is the same as that of the independence case. Although the independence assumption appears to be custom in risk theory, it has recently been subject to severe criticism due to the fact that it does not provide reliable estimates under extreme events, for example, see Boudreault et al. (2006) and Zhang and Yang (2011). Thus, the use of dependent risk models has exhibited some positive interests in the literature. Yang et al. (2012) have also studied the asymptotic behavior of the tail probabilities of a randomly weighted sum of subexponential random variables under dependence structure. In their paper, they have assumed that the random weights and the corresponding summand are dependent and the sequence of the pairs of random variables are i.i.d. Liu and Gao (2015) have considered the problem of the tail behavior of the discounted aggregate sum in dependent risk models with constant force of interest. They have assumed that the claim sizes are of an upper tail asymptotic independence structure. Moreover, they have assumed that the claim sizes and the inter-occurrence time have a specific dependence structure and they have derived an asymptotic formula in the case where the claim size distribution belongs to the intersection of long-tailed distribution class and dominant variation class. Albrecher and Boxma (2004, 2005) have developed a dependent risk model in which they have assumed that the inter-claim time and the forthcoming claim amounts are governed by a semi-Markov process. Further to that, Albrecher and Teugels (2006) have used a random walk process to model the dependency between the inter-claim time and the forthcoming claim amounts. Their results are more elaborate as they have modeled this dependency using a copula and have derived exponential estimates of the ruin probabilities in finite and infinite time horizons. Yang et al. (2014) have considered a nonstandard risk model in which they have assumed a Lévy process for the interest rate. Assuming a specific dependency among the claims, they have derived the tail probability of the discounted aggregate claim. Lu et al. (2018) have extended the work of Yang et al. (2014) by considering the multi-risk model. They incorporated the dependence structure among the claim severities and by assuming a multidimensional Lévy process for the force of interest, they have derived the asymptotic tail distribution of the discounted aggregate claim. This paper attempts to assess the impact of the dependency between the inter-occurrence time and the forthcoming claim severities through a copula approach and derives the asymptotic tail probabilities. The closest study to ours is that of Asimit and Badescu (2010) who have considered in a constant force of interest environment a heavy tailed distribution of the discounted aggregate claim in which the dependence structure is modeled via copulas. They have derived numerically the asymptotic tail probabilities, the asymptotic finite time ruin probabilities and asymptotic approximations for some common risk measures associated with the discounted aggregate claims distribution, without taking into account any effect of the interest rate. Asimit and Badescu (2010) have also derived a convolution form of the asymptotic tail probability of the aggregate claim under homogeneous Poisson risk model when the force of interest rate δ 0 and a closed-form expression has only been derived for the special case of δ = 0 . In the present paper, we have extended their result by finding a closed-form expression of the tail probability under homogeneous, mixed and non-homogeneous Poisson processes (Theorems 2–7), assuming that the counting process following Poisson distribution may be restrictive.
The model inputs are set as follows:
  • The counting processes { N ( t ) , t 0 } , { N n h ( t ) , t 0 } and { N m h ( t ) , t 0 } form respectively homogeneous, non-homogeneous and mixed Poisson processes.
  • Non-negative random variables { V k , k N } that represent the inter-occurrence times.
  • Non-negative random variables { T k , k N } such that T k = j = 1 k V j , represent the claim occurrence times.
  • Positive i.i.d random variables { B k , k N } that represent the claim severities. It is assumed that { B k , V k , k N 1 } are independent of { B 1 , V 1 } .
The discounted aggregate claim process for homogeneous, non-homogeneous or mixed Poisson models over a finite time horizon ( 0 ; t ] is given by
Z δ ( t ) = k = 1 N ( t ) B k e δ T k , Z δ n h ( t ) = k = 1 N n h ( t ) B k e δ T k , Z δ m h ( t ) = k = 1 N m h ( t ) B k e δ T k ,
where Z δ ( t ) = Z δ n h ( t ) = Z δ m h ( t ) = 0 if N ( t ) = N n h ( t ) = N m h ( t ) = 0 and δ is the force of interest.
The remaining is structured as follows. Section 2 reviews the type of dependence structure used in the study. Section 3 gives a closed-form expression of the tail probability for the homogeneous, non-homogeneous and mixed Poisson risk models and finally in Section 4 we apply these results to derive some risk measures such as the asymptotic finite time ruin probability and the value-at-risk and also provide some numerical illustrations for these risk measures.

2. Preliminaries

This paper considers the same dependence structure for the k-th inter-occurrence time V k and its corresponding claim size B k as in Asimit and Badescu (2010) and derived the tail probability for the aggregate claim.
The following assumption (from Asimit and Badescu (2010)) characterizes the type of dependence structure under study.
Assumption A1.
The pairs of random variables V i , B i , i = 1 , 2 , are mutually i.i.d. Moreover, there exists a positive and locally bounded function l ( ) such that the relation
Pr B 1 > b | V 1 = v 1 Pr B 1 > b l ( v 1 ) ,
holds uniformly for all v 1 ( 0 , T ] as b + where T is the time horizon.
Here uniformly is understood as
lim b sup v ( 0 , t ] Pr B 1 > b | V 1 = v 1 Pr B 1 > b l ( v 1 ) 1 = 0 .
As in Asimit and Badescu (2010), Assumption 1 is motivated by the fact that under its premises, one can study a wide class of dependence structure given in terms of a copula.
From Sklar’s theorem (see Sklar 1959), if H is a joint cumulative distribution function with continuous marginals H 1 and H 2 respectively, there exists a unique copula C such that
H ( b , v ) = C ( H 1 ( b ) , H 2 ( v ) ) , ( b , v ) D o m ( H ) .
So, for example, the survival copula is defined as the copula relative to the joint survival function and is given by
C ^ ( u 1 , u 2 ) = u + v 1 + C ( 1 u 1 , 1 u 2 ) , ( u 1 , u 2 ) [ 0 , 1 ] 2 .
For more details on the theory of copulas, the reader is referred to the book of Nelsen (1999).
It follows (see Asimit and Badescu 2010) that if C ^ 2 ( u 1 , u 2 ) = C ^ ( u 1 , u 2 ) u 2 exists, then Assumption 1 can be rewritten as
lim u 1 0 sup u 2 [ e λ T , 1 ) C ^ 2 ( u 1 , u 2 ) u 1 h ( u 2 ) 1 = 0 , where h ( v ) : = l log v λ .
In this paper, we consider only the case of heavy-tailed claim size distributions. Let us denote by F the cumulative distribution function of a non-negative random variable B, F ¯ its tailed distribution and F * n ¯ the n-fold convolution of F ¯
(1)
The distribution of B is said to belong to the subexponential class ( S ) if and only if
lim b F * 2 ¯ ( b ) F ¯ ( b ) = 2 .
(2)
The distribution of B is said to belong to regularly varying class ( R V α ) if and only if
lim b F ¯ ( b y ) F ¯ ( b ) = y α , α > 0 .
For more information on heavy-tailed distributions and its applications, the reader is referred to Su et al. (2007) and references therein.

3. Asymptotic Expression of the Cumulative Distribution of Z δ ( T )

Recall from Section 1 that Z δ ( T ) is the random variable that represents the discounted aggregate claim process. If the actuary knows the form of Z δ ( T ) , then he can apply the quantile principle to calculate the premium or risk measures such as the ruin probability or the value-at-risk.

3.1. Homogeneous Poisson Risk Process

Let us recall the following theorem from Asimit and Badescu (2010) for the tail probability of Z δ ( t ) .
Theorem 1.
Consider a homogeneous compound Poison risk process with constant force of interest δ > 0 such that the cumulative function ( F ) of B belongs to R V α . If Assumption 1 is satisfied for any t ( 0 ; T ] , T is time horizon, then
P Z δ ( T ) > b K δ ( T ) P ( B 1 > b ) , b + ,
where
K δ ( T ) = e λ T n = 1 λ n 0 T 0 v n 0 v 2 i = 1 n l ( v i ) e α δ j = 1 i v j d v 1 d v n .
Theorem 1 gives the convolution form of the tail probability of the aggregate claim. Although this result is interesting, it has a drawback as a closed analytical expression is still unknown, thus the multiple integrals need to be evaluated numerically and may cost in time. To address the issue, the following theorems derive and solve the differential equation of K δ ( T ) given above.
Theorem 2.
Under the assumption of Theorem 1, the deterministic function K δ ( T ) satisfied the following differential equation
d 2 K δ ( T ) d T 2 = λ e α δ T α δ λ d K δ ( T ) d T + λ l ( T ) e ( α δ + λ ) T ,
where K δ ( 0 ) = 0 and K δ ( T ) T | T = 0 = λ l ( 0 ) .
Remark 1.
If δ = 0 then
d 2 K 0 ( T ) d T 2 = λ l ( T ) e λ T .
The solution of this differential equation is consistent with the result of Asimit and Badescu (2010).
Proof. 
From Equation (3) of Theorem 1 we have for any T 0 ,
K δ ( T ) = e λ T n = 1 λ n D n i = 1 n l ( v i ) e α δ j = 1 i v j d v = e λ T n = 1 λ n 0 T 0 v n 0 v 2 i = 1 n l ( v i ) e α δ j = 1 i v j d v 1 d v n = e λ T n = 1 λ n 0 T 0 v n 0 v 2 i = 1 n l ( v i ) j = 1 i e α δ v j d v 1 d v n = e λ T 0 T λ l ( v ) e α δ v d v + n = 2 λ n 0 T 0 v n 0 v 2 i = 1 n l ( v i ) e α δ j = 1 i v j d v 1 d v n ,
where, at n = 1 , the preceding notations mean that we integrate only between 0 and T.
By taking the first derivative of the preceding expression with respect to T we have,
d K δ ( T ) d T = λ K δ ( T ) + e λ T λ l ( T ) e α δ T + e λ T λ n = 2 λ n 1 0 T 0 v n 1 0 v 2 i = 1 n 1 l ( v ) j = 1 i e α δ v j + l ( v ) e α δ j = 1 n 1 v j e α δ T d v 1 d v n 1 = λ K δ ( T ) + λ l ( T ) e ( α δ + λ ) T + e λ T λ n = 2 λ n 1 0 T 0 v n 1 0 v 2 i = 1 n 1 l ( v ) j = 1 i e α δ w j d v 1 d v n 1 + e ( α δ + λ ) T λ n = 2 λ n 1 0 T 0 v n 1 0 v 2 l ( v ) e α δ j = 1 n 1 v j d v 1 d v n 1 .
Rearranging this expression gives,
d K δ ( T ) d T = λ K δ ( T ) + λ l ( T ) e ( α δ + λ ) T + e λ T λ n = 1 λ n 0 T 0 v n 0 v 2 i = 1 n l ( v ) j = 1 i e α δ v j d v 1 d v n + e ( α δ + λ ) T λ n = 1 λ n 0 T 0 v n 0 v 2 l ( v ) e α δ j = 1 n v j d v 1 d v n = λ K δ ( T ) + λ l ( T ) e ( α δ + λ ) T + λ K δ ( T ) + e ( α δ + λ ) T λ n = 1 λ n 0 T 0 v n 0 v 2 l ( v ) e α δ j = 1 n v j d v 1 d v n = λ e ( α δ + λ ) T l ( T ) + n = 1 λ n 0 T 0 v n 0 v 2 l ( v ) e α δ j = 1 n v j d v 1 d v n .
As previously, the derivative of the preceding expression with respect to T yields,
d 2 K δ ( T ) d T 2 = α δ + λ d K δ ( T ) d T + λ e ( α δ + λ ) T { l ( T ) + λ l ( T ) e α δ T + λ n = 2 λ n 1 0 T 0 v n 1 0 v 2 l ( v ) e α δ j = 1 n 1 v j e α δ T d v 1 d v n 1 } = α δ + λ d K δ ( T ) d T + λ e ( α δ + λ ) T { l ( T ) + λ l ( T ) e α δ T + λ e α δ T n = 1 λ n 0 T 0 v n 0 v 2 l ( v ) e α δ j = 1 n v j d v 1 d v n } .
From Equation (5) we have,
λ e α δ T n = 1 λ n 0 T 0 v n 0 v 2 l ( v ) e α δ j = 1 n v j d v 1 d v n = e λ T d K δ ( T ) d T λ l ( T ) e α δ T .
Then, substituting (7) into (6) leads to,
d 2 K δ ( T ) d T 2 = α δ + λ d K δ ( T ) d T + λ e ( α δ + λ ) T l ( T ) + λ l ( T ) e α δ T + e λ T d K δ ( T ) d T λ l ( T ) e α δ T = α δ + λ d K δ ( T ) d T + λ e ( α δ + λ ) T l ( T ) + e λ T d K δ ( T ) d T = α δ + λ K δ ( T ) T + λ l ( T ) e ( α δ + λ ) T + λ e α δ T d K δ ( T ) d T = λ e α δ T α δ λ d K δ ( T ) d T + λ l ( T ) e ( α δ + λ ) T ,
which proves the statement. □
The next theorem gives an analytical expression of K δ ( T ) in a homogeneous Poisson risk model.
Theorem 3.
Under the assumption of Theorem 2, if we assume further that K δ ( t ) > 0 for every t [ 0 , T ] , then the solution of the differential equation given by expression (4) of Theorem 2 is
K δ ( T ) = 0 T e λ α δ e α δ v + ( α δ + λ ) v d v + λ 0 T l ( v ) e ( α δ + λ ) v 1 + α δ + λ T v d v ,
where
lim δ 0 λ α δ e α δ v = + .
Proof. 
(1)
Firstly let us solve Equation (4) of Theorem 2 without the second member.
From d 2 K δ ( T ) d T 2 = λ e α δ T α δ λ d K δ ( T ) d T ,
we have,
d d T ln d L δ ( T ) d T = λ e α δ T α δ λ ln d L δ ( T ) d T = λ α δ e α δ T + ( α δ + λ ) T + μ 0 d L δ ( T ) d T = e λ α δ e α δ T + ( α δ + λ ) T + μ 1 L δ ( T ) = 0 T e λ α δ e α δ v + ( α δ + λ ) v d v + μ .
By letting μ = 0 we have,
L δ ( T ) = 0 T e λ α δ e α δ v + ( α δ + λ ) v d v .
(2)
Secondly let us find the solution for the second member of Equation (4).
Let for any T 0 ,
H ( T ) = λ 0 T l ( v ) e ( α δ + λ ) v 1 + α δ + λ T v d v .
One can easily check that
d 2 H ( T ) d T 2 = λ l ( T ) e ( α δ + λ ) T .
Therefore, the final solution is
K δ ( T ) = L δ ( T ) + H ( T ) ,
which proves the statement. □
Remark 2.
In the case of a mixed Poisson compound model, the differential equation of K δ ( T ) given Λ = λ is again given by Equation (4)
We can now state the theorem that gives the closed-form expression of K δ ( T ) in a mixed Poisson compound risk model.
Theorem 4.
Consider a mixed Poisson compound risk model. Assume that the intensity rate Λ is a continuous positive random variable such that its moment generating function M Λ ( s ) exists then, the closed-form expression of K δ ( T ) is given by
K δ ( T ) = 0 T M Λ v e α δ v α δ e α δ v + E l ( v ) e ( α δ + Λ ) v Λ 1 + α δ + Λ T v d v .
Proof. 
From Theorem 3, we have
K δ ( T ) | Λ = λ = 0 T e λ α δ e α δ v + ( α δ + λ ) v + λ l ( v ) e ( α δ + λ ) v 1 + α δ + λ T v d v ,
therefore,
K δ ( T ) = E 0 T e Λ α δ e α δ v + ( α δ + Λ ) v + Λ l ( v ) e ( α δ + Λ ) v 1 + α δ + Λ T v d v ,
and using Fubini’s theorem, we can interchange the expectation and the integral.
Hence,
K δ ( T ) = 0 T E e Λ α δ e α δ v + ( α δ + Λ ) v + E Λ l ( v ) e ( α δ + Λ ) v 1 + α δ + Λ T v d v = 0 T E e e α δ v α δ Λ + Λ v e α δ v + E Λ l ( v ) e ( α δ + Λ ) v 1 + α δ + Λ T v d v = 0 T E e e α δ v α δ + v Λ e α δ v + E l ( v ) e ( α δ + Λ ) v Λ 1 + α δ + Λ T v d v = 0 T M Λ v e α δ v α δ e α δ v + E l ( v ) e ( α δ + Λ ) v Λ 1 + α δ + Λ T v d v .
 □
Using specific copulas given in terms of a closed analytical expression of l ( ) one can reduce Equation (9).
Corollary 1.
Under assumption of Theorem 4. If the dependence structure is modeled by Gumbel–Barnett copula then,
K δ ( T ) = 0 T M Λ v e α δ v α δ e α δ v d v + ( 1 + θ ) n = 0 ( T ) n + 1 ( n + 1 ) ! × T k = 0 n + 1 n + 1 k α δ n + 1 k M Λ ( k + 1 ) ( 0 ) + ( 1 n ) k = 0 n n k α δ n k M Λ ( k + 1 ) ( 0 )
The expression of K δ ( T ) is similar for the other copulas.
Proof. 
Using the result in Section 3.2, we have for these copulas.
Λ 0 T l ( v ) e ( α δ + Λ ) v 1 + α δ + Λ T v d w = Λ ( 1 + θ ) { 2 1 e ( α δ + Λ ) T α δ + Λ + ( α δ + Λ ) T 1 e ( α δ + Λ ) T α δ + Λ + T e ( α δ + Λ ) T } .
Moreover,
2 Λ ( 1 + θ ) 1 e ( α δ + Λ ) T α δ + Λ = 2 ( 1 + θ ) n = 0 ( T ) n + 1 ( n + 1 ) ! k = 0 n n k Λ k + 1 α δ n k Λ ( 1 + θ ) ( α δ + Λ ) T 1 e ( α δ + Λ ) T α δ + Λ = ( 1 + θ ) T n = 0 ( T ) n + 1 ( n + 1 ) ! k = 0 n + 1 n + 1 k Λ k + 1 α δ n + 1 k Λ ( 1 + θ ) T e ( α δ + Λ ) T = ( 1 + θ ) n = 0 ( n + 1 ) ( T ) n + 1 ( n + 1 ) ! k = 0 n n k Λ k + 1 α δ n k .
The result follows by taking the expectation and interchanging the sum and the expectation. □
In the next subsection we give a closed-form expression of the constant K δ ( T ) where the dependence structure is modeled by some specific well-known copulas.

3.2. Examples with Some Copulas

Recall that e x = n = 0 x n n ! and ( a + b ) n = k = 0 n n k a k b n k . Let
H = 0 T e λ α δ e α δ v + ( α δ + λ ) v d v
therefore,
H = 0 T n = 0 λ α δ e α δ v + ( α δ + λ ) v n n ! d v = 0 T n = 0 ( 1 ) n n ! λ α δ e α δ v + ( α δ + λ ) v n d v = 0 T n = 0 ( 1 ) n n ! k = 0 n n k λ α δ e α δ v k ( α δ + λ ) v n k d v = 0 T n = 0 ( 1 ) n n ! k = 0 n n k α δ + λ n k λ α δ e α δ v k v n k d v = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k 0 T e k α δ v v n k d v = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k α δ k k 1 n 0 α δ k T e x x n k d x = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k α δ k k 1 n γ ( n + 1 k , α δ T ) ,
where
γ ( α , x ) = 0 x t α 1 e t d t is the incomplete gamma function .
then,
H = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k α δ k k 1 n γ ( n + 1 k , α δ T )
Example 1.
Ali–Mikhail–Haq copula:
C ( u 1 , 2 ) = u 1 u 2 1 θ ( 1 u 1 ) ( 1 u 2 ) , θ [ 1 , 1 ) l ( v ) = 1 θ 1 2 e λ v .
If the dependence structure is modeled by the Ali–Mikhail–Haq copula then,
K δ ( T ) = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k α δ k k 1 n γ ( n + 1 k , α δ T ) + λ T 1 + θ 2 θ α δ + λ α δ + 2 λ + 2 λ ( 1 + θ ) a ¯ ( α δ + λ ) T 2 λ 2 θ α δ + 2 λ a ¯ ( α δ + 2 λ ) T .
This expression can be easily evaluated numerically.
Proof. 
λ 0 T l ( v ) e ( α δ + λ ) v ( 1 + ( α δ + λ ) ( T v ) d v = λ 0 T l ( v ) e ( α δ + λ ) v d v + λ ( α δ + λ ) 0 T l ( v ) e ( α δ + λ ) v ( T v ) d v ,
hence,
λ 0 T l ( v ) e ( α δ + λ ) v d v = ( 1 + θ ) λ 0 T e ( α δ + λ ) v d v 2 θ λ 0 T e ( α δ + 2 λ ) v d v = ( 1 + θ ) λ 1 e ( α δ + λ ) T α δ + λ 2 λ θ 1 e ( α δ + 2 λ ) T α δ + 2 λ = ( 1 + θ ) λ a ¯ ( ( α δ + λ ) T ) 2 λ θ a ¯ ( ( α δ + 2 λ ) T ) .
Let x = T v then,
0 T l ( v ) e ( α δ + λ ) T ( T v ) d v = 0 T x l ( T x ) e ( α δ + λ ) T x d v = e ( α δ + λ ) T 0 T x l ( T x ) e ( α δ + λ ) x d x = e ( α δ + λ ) T 0 T x 1 + θ 1 2 e λ ( T x ) e ( α δ + λ ) x d x = e ( α δ + λ ) T ( 1 + θ ) 0 T x e ( α δ + λ ) x d x 2 θ e λ T 0 T x e ( α δ + 2 λ ) x d x .
Using integration by parts yields
λ ( α δ + λ ) 0 T l ( v ) e ( α δ + λ ) T ( T v ) d v = λ T 1 + θ 2 θ α δ + λ α δ + 2 λ + λ ( 1 + θ ) a ¯ ( α δ + λ ) T + 2 θ α δ + λ α δ + 2 λ a ¯ ( α δ + 2 λ ) T .
The result follows by combining (11)–(13). □
Remark 3.
For the FGM copula, l ( . ) is the same as for the Ali–Mikhail–Haq copula; therefore, K δ ( T ) will be the same.
Example 2.
Gumbel–Barnett copula:
C ( u 1 , u 2 ) = u 1 u 2 exp ( θ ln u 1 ln u 2 ) , 0 θ 1 , l ( v ) = 1 θ θ ln 1 e λ v .
If the dependence structure is modeled by the Gumbel–Barnett copula then,
K δ ( T ) = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k α δ k k 1 n γ ( n + 1 k , α δ T ) + λ ( 1 + θ ) 2 + ( α δ + λ ) T a ¯ ( α δ + λ ) T + T e ( α δ + λ ) T + λ θ k = 2 1 k 1 α δ + λ α δ + λ k T e ( α δ + λ k ) T + a ¯ ( α δ + λ k ) T + 1 + ( α δ + λ ) T a ¯ ( α δ + λ k ) T .
Proof. 
Let
κ ( v ) = e ( α δ + λ ) v ( 1 + ( α δ + λ ) ( T v ) ) ,
then using (14),
λ 0 T l ( v ) κ ( w ) d v = λ ( 1 + θ ) 0 T e ( α δ + λ ) v ( 1 + ( α δ + λ ) ( T v ) ) d v = λ 1 + θ ) ( 1 + ( α δ + λ ) T 0 T e ( α δ + λ ) v d v λ ( 1 + θ ) ( α δ + λ ) 0 T v e ( α δ + λ ) v d v ,
using integration by parts yields
λ 0 T l ( v ) κ ( v ) d v = λ ( 1 + θ ) 2 + ( α δ + λ ) T a ¯ ( α δ + λ ) T + T e ( α δ + λ ) T .
Moreover,
0 T ln 1 e λ v κ ( v ) d v = 0 T n = 1 e n λ v n e ( α δ + λ ) v 1 + ( α δ + λ ) ( T v ) d v = n = 1 [ 1 + ( α δ + λ ) T 0 T e ( α δ + ( n + 1 ) λ ) v ( α δ + λ ) 0 T v e ( α δ + ( n + 1 ) λ ) v ] = n = 1 1 n [ 1 + ( α δ + λ ) T a ¯ ( α δ + ( n + 1 ) λ ) T + α δ + λ α δ + ( n + 1 ) λ ) T e ( α δ + ( n + 1 ) λ ) T + a ¯ ( α δ + ( n + 1 ) λ ) T ] .
By taking k = n + 1 in the above result and combining it with (11) and (15) yields the result. □
Example 3.
Marshall–Olkin copula:
C ( u 1 , u 2 ) = min u 1 1 θ 1 u 2 , u 1 u 2 1 θ 2 0 θ 1 , θ 1 1 with l ( v ) = 1 θ 1 .
K δ ( T ) is given by
K δ ( T ) = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k α δ k k 1 n γ ( n + 1 k , α δ T ) + λ ( 1 + θ 1 ) 2 + ( α δ + λ ) T a ¯ ( α δ + λ ) T + T e ( α δ + λ ) T .
Example 4.
Frank copula:
C ( u 1 , u 2 ) = 1 θ ln 1 + e θ u 1 1 e θ u 2 1 e θ 1 , θ R * , l ( v ) = θ e θ 1 e θ e θ 1 e λ v .
If the dependence is modeled by the Frank copula then,
K δ ( T ) = n = 0 ( 1 ) n n ! k = 0 n n k λ α δ k α δ + λ n k α δ k k 1 n γ ( n + 1 k , α δ T ) + λ θ e θ 1 e θ n = 0 θ n n ! k = 1 n + 1 n k 1 ( 1 ) k 1 { α δ + λ α δ + λ k a ¯ ( α δ + λ k ) T + T e ( α δ + λ k ) T + ( 1 + ( α δ + λ ) T ) a ¯ ( α δ + λ k ) T } .
Proof. 
Recall that
l ( v ) = θ e θ 1 e θ e θ 1 e λ v = θ e θ 1 e θ n = 0 θ n n ! j = 0 n n j ( 1 ) j e j λ v .
Then using (14)
0 T g ( v ) κ ( v ) d v = θ e θ 1 e θ n = 0 θ n n ! j = 0 n n j ( 1 ) j { ( 1 + ( α δ + λ ) T ) × 0 T e ( α δ + ( j + 1 ) λ ) v d v ( α δ + λ ) 0 T v e ( α δ + ( j + 1 ) λ ) v d v } , = θ e θ 1 e θ n = 0 θ n n ! j = 0 n n j ( 1 ) j { 1 + ( α δ + λ ) T a ¯ ( α δ + ( j + 1 ) λ ) T + α δ + λ α δ + ( j + 1 ) λ a ¯ ( α δ + ( j + 1 ) λ ) T + T e ( α δ + ( j + 1 ) λ ) T } .
By letting k = j + 1 we have,
λ 0 T l ( v ) κ ( v ) d v = λ θ e θ 1 e θ n = 0 θ n n ! k = 1 n + 1 n k 1 ( 1 ) k 1 { 1 + ( α δ + λ ) T × a ¯ ( α δ + k λ ) T + α δ + λ α δ + k λ a ¯ ( α δ + k λ ) T + T e ( α δ + k λ ) T } ,
The result follows by combining (11) and (16). □

3.3. Non-Homogeneous Poisson Risk Model

In this subsection, we extend our result by considering a non-homogeneous compound Poisson process where the intensity rate is a deterministic function λ ( t ) , which we assume to be piece-wise continuous and bounded on finite intervals.
Theorem 5.
Consider a non-homogeneous Poisson compound renewal model with constant force of interest ( δ > 0 ) such that the cumulative distribution function F R V α . If Assumption 1 is fulfilled for any time t ( 0 , T ] then,
P Z δ n h ( T ) > b K δ ( T ) P ( B 1 > b ) , b
where
K δ ( T ) = e m ( T ) n = 1 0 T λ ( v n ) 0 v 2 λ ( v 1 ) i = 1 n l ( v i ) e α δ j = 1 i v j d v 1 d v n ,
and m ( T ) = 0 T λ ( s ) d s is mean function (or cumulative intensity function).
Before proving Theorem 5, let us recall the following useful Lemmas.
Lemma 1
(From Ruwet (2007)). Consider a non-homogeneous Poisson process N ( t ) with an intensity rate λ ( t ) as a function of t. Then the conditional distribution function of the first n occurrence times of the claims, assuming that N ( t ) = n , has the same distribution as the order statistics of the n i.i.d. uniform random variables U 1 , U 2 , , U n with common distribution
F U i ( u ) = m ( u ) m ( t ) , 0 < u < t , where m ( t ) = 0 t λ ( s ) d s
i.e.,
f T 1 , T 2 , , T n | N ( t ) = n ( t 1 , t 2 , , t n | n ) = n ! m n ( t ) i = 1 n λ ( t i ) , t 1 < t 2 < < t n .
Lemma 2
(From Asimit and Badescu (2010)). Let B, X 1 , X 2 , be a sequence of independent non-negative random variables with F the cumulative distribution of B such that F S . Assume further that there exists a constant M such that P ( X i > x ) M P ( B > b ) holds for all b > 0 and any i = 1 , 2 , . Then, for any ϵ > 0 there exists A < such that
s n : = P i = 1 n X i > x P ( B > b ) A 1 + ϵ n ,
holds for any integer n.
Proof of Theorem 5.
By conditioning on the number of claims and inter-occurrence time by time T, the cumulative distribution functions of Z δ ( T ) follow the integral equation
P ( Z δ n h ( T ) > b ) = n = 1 P ( N ( T ) = n ) D n P k = 1 n B k e δ i = 1 k V i > b | V = v , N ( T ) × P ( V = v | N ( T ) = n ) d v = n = 1 P ( N ( T ) = n ) D n P k = 1 n B k e δ i = 1 k V i > b | V = v × f V 1 , , V n | N ( T ) ( t 1 , t n | n ) d v .
Due to the fact that
P k = 1 n B k e δ i = 1 k V i > b | V = v P k = 1 n B k > b | V = v ,
by the fact that B k | V k = v k are i.i.d. and F R V α therefore belongs to S , we can apply Lemma 2. This is true since Assumption 1 implies that there exists a constant M > 0 such that for all b > 0 and v ( 0 , T ] , we have
P B 1 e δ V 1 > b | V 1 = v 1 P B 1 > b | V 1 = v 1 M P ( B 1 > b ) .
From Lemma 2, it follows that for any ϵ > 0 there exists A > 0 such that,
P k = 1 n B k e δ i = 1 k V i > b | V = v P k = 1 n B k > b | V = v A 1 + ϵ n P ( B 1 > b )
holds uniformly for all b > 0 , v ( 0 , T ] n and n = 1 , 2 , .
Then
n = 1 P ( N ( T ) = n ) D n P k = 1 n B k e δ i = 1 k V i > b | V = v P ( V = v | N ( T ) = n )
n = 1 P ( N ( T ) = n ) D n P k = 1 n B k > b | V = v P ( V = v | N ( T ) = n )
n = 1 P ( N ( T ) = n ) D n A 1 + ϵ n P ( B 1 > b ) P ( V = v | N ( T ) = n )
A P ( B 1 > b ) n = 1 P ( N ( T ) = n ) = A P ( B 1 > b ) P N ( T ) ( 1 + ϵ ) < .
Since P N ( T ) exists for any t > 0 .
Using Lemma 1 together with the fact that P ( N ( T ) = n ) = m ( T ) n n ! e m ( T ) we have
P ( Z δ n h ( T ) > b ) = n = 1 e m ( T ) D n i = 0 n λ ( v i ) P k = 1 n B k e δ i = 1 k V i > b | V = v d v .
The result follows by applying the Dominated Convergence Theorem in (18). □
The following theorem gives the expression of the differential equation of the tail probability of the discounted aggregate claims under a non-homogeneous Poisson risk model.
Theorem 6.
Under the assumption of Theorem 5, if the the intensity function λ ( t ) 0 for all t [ 0 , T ] , then the deterministic function K δ ( T ) given in Theorem 5 satisfied the following differential equation
d 2 K δ ( T ) d T 2 = λ ( T ) λ ( T ) + λ ( T ) e α δ T ( λ ( T ) + α δ ) d K δ ( T ) d T + λ ( T ) l ( T ) e m ( T ) + α δ T ,
where K δ ( 0 ) = 0 , K δ ( T ) T | T = 0 = λ ( 0 ) l ( 0 ) and m ( T ) = 0 T λ ( s ) d s .
Remark 4.
If λ ( t ) = λ for all t [ 0 , T ] , then λ ( t ) = 0 and m ( T ) = λ T ; therefore, the result in Theorem 6 is consistent with the result of Theorem 2
Proof. 
From Equation (17) of Theorem 5 we have for any T 0 ,
K δ ( T ) = e m ( T ) n = 1 D n k = 1 n λ ( v k ) i = 1 n l ( v i ) e α δ j = 1 i v j d v = e m ( T ) n = 1 0 T λ ( v n ) 0 v n λ ( v n 1 ) 0 v 2 λ ( v 1 ) i = 1 n l ( v i ) e α δ j = 1 i v j d v 1 d v n = e m ( T ) n = 1 0 T λ ( v n ) 0 v n λ ( v n 1 ) 0 v 2 λ ( v 1 ) i = 1 n l ( v i ) j = 1 i e α δ v j d v 1 d v n = e m ( T ) { 0 T λ ( v ) l ( v ) e α δ v d v + n = 2 0 T 0 v n λ ( v n 1 ) 0 v 2 λ ( v 1 ) i = 1 n l ( v i ) e α δ j = 1 i v j d v 1 d v n } .
As in Theorem 2, the derivative of the preceding expression with respect to T yields.
d K δ ( T ) d T = m ( T ) K δ ( T ) + e m ( T ) λ ( T ) l ( T ) e α δ T + e m ( T ) λ ( T ) n = 2 0 T λ ( v n 1 ) 0 v 2 λ ( v 1 ) i = 1 n 1 l ( v ) j = 1 i e α δ v j + l ( v ) e α δ j = 1 n 1 v j e α δ T d v 1 d v n 1 = λ ( T ) K δ ( T ) + λ ( T ) l ( T ) e ( m ( T ) + α δ T ) + e m ( T ) λ ( T ) { n = 2 0 T λ ( v n 1 ) 0 v 2 λ ( v 1 ) i = 1 n 1 l ( v ) j = 1 i e α δ v j d v 1 d v n 1 + e α δ T n = 2 0 T λ ( v n 1 ) 0 v 2 λ ( v 1 ) l ( v ) e α δ j = 1 n 1 v j d v 1 d v n 1 } = λ ( T ) K δ ( T ) + λ ( T ) l ( T ) e ( m ( T ) + α δ T ) + e m ( T ) λ ( T ) n = 1 0 T λ ( v n ) 0 v 2 λ ( v 1 ) i = 1 n l ( v ) j = 1 i e α δ v j d v 1 d v n + e ( m ( T ) + α δ T ) λ ( T ) n = 1 0 T λ ( v n ) 0 v 2 λ ( v 1 ) l ( v ) e α δ j = 1 n v j d v 1 d v n = λ ( T ) K δ ( T ) + λ ( T ) l ( T ) e ( m ( T ) + α δ T ) + λ ( T ) K δ ( T + e ( m ( T ) + α δ T ) λ ( T ) ) × n = 1 0 T λ ( v n ) 0 v n λ ( v n 1 ) 0 v 2 λ ( v 1 ) l ( v ) e α δ j = 1 n v j d v 1 d v n = λ ( T ) e ( m ( T ) + α δ T ) l ( T ) + + λ ( T ) e ( m ( T ) + α δ T ) n = 1 0 T λ ( v n ) 0 v 2 λ ( v 1 ) l ( v ) e α δ j = 1 n v j d v 1 d v n .
By taking the derivative of the preceding expression with respect to T we have,
d 2 K δ ( T ) d T 2 = λ ( T ) λ ( T ) α δ + λ ( T ) d K δ ( T ) d T + λ ( T ) e ( m ( T ) + α δ T ) { l ( T ) + λ ( T ) l ( T ) e α δ T + λ ( T ) n = 2 0 T λ ( v n 1 ) 0 v 2 λ ( v 1 ) l ( v ) e α δ j = 1 n 1 v j e α δ T d v 1 d v n 1 } = λ ( T ) λ ( T ) α δ + λ ( T ) d K δ ( T ) d T + λ ( T ) e ( m ( T ) + α δ T ) { l ( T ) + λ ( T ) l ( T ) e α δ T + λ ( T ) e α δ T n = 1 0 T λ ( v n ) 0 v 2 λ ( v 1 ) l ( v ) e α δ j = 1 n v j d v 1 d v n } .
From Equation (20) we have,
λ ( T ) e α δ T n = 1 0 T λ ( v n ) 0 v 2 λ ( v 1 ) l ( v ) e α δ j = 1 n v j i = 1 n d v i = e m ( T ) d K δ ( T ) d T λ ( T ) l ( T ) e α δ T .
Then, substituting (22) into (21) leads to,
d 2 K δ ( T ) d T 2 = λ ( T ) λ ( T ) α δ + λ ( T ) d K δ ( T ) d T + λ ( T ) e ( m ( T ) + α δ T ) { l ( T ) + λ ( T ) l ( T ) e α δ T + e m ( T ) d K δ ( T ) d T λ ( T ) l ( T ) e α δ T } = λ ( T ) λ ( T ) α δ + λ ( T ) d K δ ( T ) d T + λ ( T ) e ( m ( T ) + α δ T ) l ( T ) + e m ( T ) d K δ ( T ) d T = λ ( T ) λ ( T ) + λ ( T ) e α δ T α δ + λ ( T ) d K δ ( T ) d T + λ ( T ) l ( T ) e ( m ( T ) + α δ T ) .
 □
The following theorem gives the solution of the differential equation of K δ ( T ) where the intensity function is not a constant.
Theorem 7.
Under the assumption of Theorem 6. If we assume further that K δ ( t ) > 0 for every t [ 0 , T ] , then
K δ ( T ) = 0 T λ ( v ) e m ( v ) + α δ v e 0 v λ ( s ) e α δ s d s d v + 0 T 0 v λ ( s ) l ( s ) e m ( s ) + α δ s d s d v .
Proof. 
(1)
Firstly let us solve Equation (19) of Theorem 6 without the second member.
From d 2 K δ ( T ) d T 2 = λ ( T ) λ ( T ) + λ ( T ) e α δ T α δ λ ( T ) d K δ ( T ) d T ,
we have,
d d T ln d L δ ( T ) d T = λ ( T ) λ ( T ) + λ ( T ) e α δ T α δ λ ( T ) ln d L δ ( T ) d T = ln λ ( T ) m ( T ) + α δ T + 0 T λ ( s ) e α δ s d s + μ 0 d L δ ( T ) d T = λ ( T ) e m ( T ) + α δ T e 0 T λ ( s ) e α δ s d s + μ 1 L δ ( T ) = 0 T λ ( v ) e m ( v ) + α δ v e 0 v λ ( s ) e α δ s d s d v + μ .
By letting μ = 0 we have,
L δ ( T ) = 0 T λ ( v ) e m ( v ) + α δ v e 0 v λ ( s ) e α δ s d s d v .
(2)
Secondly let us find the solution of Equation (19) of Theorem (6) for the second member.
Let for any T 0 ,
H ( T ) = 0 T 0 v λ ( s ) l ( s ) e m ( s ) + α δ s d s d v .
One can easily check that
d 2 H ( T ) d T 2 = λ ( T ) l ( T ) e m ( T ) + α δ T .
Therefore, the final solution is
K δ ( T ) = L δ ( T ) + H ( T ) ,
which proves the statement. □

4. Applications

In this section we use the expression of K δ ( T ) to derive an upper bound of the ruin probability and also to evaluate some risk measure such as the value-at-risk.

4.1. Ruin Probability

Finding the exact analytical expression of the ruin probability that can be computationally tractable is an extremely complex task. For simple cases of the collective risk model, exact formulas and approximations have been found for the ruin probability. However, when the effect of the force of interest is incorporated into the collective risk model, the calculations become more difficult. In the case of the discounted compound Poisson risk model, with a constant force of interest, a differential equation has been obtained for the ruin probability and a solution has been found when the amount of the claims follows an exponential distribution. For more details, see Sundt and Teugels (1995). The calculation of the ruin probability is also possible if we know the distribution of our risk process. The expressions obtained by Leveillé and Hamel (2018); Léveillé et al. (2010) for the distribution of the discounted aggregated renewal sum confirm that it will be very difficult to find an explicit expression of the probability of ruin in the framework defined by the previous authors. Hence the use of bounds for an estimation of the probability of ruin. Several bounds have been proposed in the framework of Sparre Andersen’s model, see Cai and Dickson (2003), when the force of interest is constant and Cai (2002) when the force of interest follows a time series. Other authors such as Albrecher and Teugels (2006), Tang (2005) and Wei and Yang (2004) have incorporated the dependency into the risk model and have derived asymptotic formulas for finite and infinite time ruin probabilities. In this subsection, we present some asymptotic results of the ruin probability for homogeneous and non-homogeneous compound Poisson risk model with a heavy tailed distribution for the claim severities and we also assume that the dependence structure is given in Section 3.2. This subsection provides numerical illustration with the use of copulas.
Let
U δ ( t ) = b e δ t + C δ ( t ) e δ t Z δ ( t ) ,
where b is the surplus at time t = 0 and C δ ( t ) = 0 t e δ ( t s ) d C ( s ) represents the accumulated amount of premiums received at time t. The non-decreasing and right continuous stochastic process, { C ( s ) } s 0 with C ( 0 ) = 0 represents the total amount of premiums accumulated to time s. We then define the time to ruin as
T δ , b = inf t > 0 : U δ ( t ) < 0 | U δ ( 0 ) = b ,
therefore, the finite-time ruin probability is given by
ψ δ ( b , T ) = P ( T δ , b T ) .
Clearly
P Z δ ( T ) b + e δ T C δ ( T ) ψ δ ( b , T ) P Z δ ( T ) b , holds for δ 0 .
Using the long-tailed property of subexponential distributions, see Embrechts et al. (1997), in (26) leads to the following Corollary of Theorems 1 and 5.
Corollary 2.
Consider a homogeneous and a non-homogeneous compound renewal risk process with constant force of interest δ > 0 such that Assumption 1 is fulfilled for any t ( 0 ; T ] , assume further that F R V α . In addition, if C δ ( T ) < then
ψ δ ( b , T ) K δ ( T ) P ( B 1 b ) , b + .
Proof. 

Numerical Illustration

In this subsection, we provide numerical illustrations for the asymptotic finite-time ruin probability. For this purpose we assume that the counting process is a homogeneous, a non-homogeneous and a mixed Poisson process. Furthermore, we assume that the dependence structure is given by some copulas defined in Section 3.2. However, the forthcoming claim amounts are assumed to have a Weibull distribution with cumulative distribution given by:
F B 1 ( b ) = 1 exp b 1 τ
Insurance claims are known to be right-skewed. Moreover, under extreme events scenarios the frequency of the claim occurrence is relatively low; therefore, we choose our parameters accordingly.
In Table 1, we consider a homogeneous compound Poisson process where λ = 0.5 , the constant force of interest rate is assumed to be δ = 20 % , τ = 2 and the time horizon T 20 , 40 .
In Table 2, we consider a mixed homogeneous compound Poisson process where Λ has gamma distribution with density
f Λ ( λ ) = β a Γ ( a ) λ a 1 exp ( β λ ) , for λ 0 ,
where a = 15 , β = 20 .
In Table 3, we consider a non-homogeneous Poisson risk model where the intensity function is assumed to be
λ ( t ) = 0.5 + 0.2 exp ( 0.2 t ) ,
and the other parameters were chosen as in Table 1.
Table 1, Table 2 and Table 3 give numerical values for the ruin probability respectively for the Ali–Mikhail–Haq and Frank copulas.
Our numerical illustrations show (see Table 1 and Table 2) that when the time to maturity increases, the ruin probability for the homogeneous and mixed Poisson processes increases. However, the inverse effect is observed for an increase of the copula’s parameter for the Frank copula. For non-homogeneous risk model (see Table 3), increasing the parameter α decreases the ruin probability in the case of the Ali–Mikhail–Haq copula.

4.2. Risk Management

In the movement to redefine the solvency margin according to risks for all financial risks, after the bank industry and the implementation of the Basel II agreement, it was the turn of insurance industry to see its regulations adapted to integrate risk. After Solvency I, which provided for a solvency margin determined as a function of percentages on premiums and claims, insurance regulations moved to more complex rules integrating risk, either by applying a standard formula or by using an internal model. Note that the subprime crisis, far from having put this reform on hold, convinced the European Commission to legislate to protect itself against so-called systemic risks (chain failures of players—banks, insurers, etc.—in the financial world).
The second pillar aims to set qualitative standards for monitoring risks internally in companies and how the supervisory authority should exercise its supervisory powers in this context. The identification of the “riskiest” companies is an objective, and the supervisory authorities will have in their power the possibility of requiring these companies to hold a capital higher than the amount suggested by the calculation of the so-called solvency capital, requirement (SCR) and/or reduce their exposure to risks. See Hoiring (2013), Gatzert and Wesker (2012) and Boonen (2017). This capital requirement covers all the risks that an insurer faces and is defined as the value-at-risk (VaR) of the basic own funds subject to a confidence level of 99.5 % on a 1-year period. In this subsection we provide an asymptotic expression for the value-at-risk.

4.2.1. Value-at-Risk

Considered as the most popular risk measure in financial industry, the value-at-risk (VaR) at level p for a loss variable L represents the p-quantile, defined as in Jorion (2001):
VaR p ( L ) = inf b R : P ( L > b ) 1 p .
An alternative definition of the value-at-risk (VaR) is the α -quantile of the distribution of the discounted aggregate sum of claims. Using Theorems 1 and 5, this quantile can be expressed as:
VaR 1 p Z δ ( T ) VaR 1 p K δ ( B 1 ) , for p 0 ,
provided that the density function of Z δ ( T ) is continuously close enough to the right tail.

4.2.2. Numerical Illustration

In this subsection, we provide numerical illustrations for the value-at-risk where Weibull distribution is considered for the claim severities using some copulas defined in Section 3.2. With the same parameters as before, we compute the value-at-risk at level p = 0.5 % .
In Table 4 and Table 5, we note that the value-at-risk is higher in the case of a negative value of θ for the Frank copula and the inverse effect is observed for the Ali–Mikhail–Haq copula. For the non-homogeneous Poisson process (see Table 6), the value-at risk is positively correlated to the parameter α in the case of the Frank copula and negatively correlated to the parameter α for the Ali–Mikhail–Haq copula.

5. Discussion and Conclusions

Under the model assumption, the numerical simulations show that the ruin is high for the Frank copula in comparison with the Ali–Mikhail–Haq copula in both homogeneous, mixed homogeneous and non-homogeneous Poisson processes. Moreover, the ruin is higher when the claim number process is modeled by a non-homogeneous process than the other process understudies. This finding was further confirmed by the simulation of the value-at-risk as shown in Table 6. One of the most important advantages of simulation in social science is the fact that simulation allows to approach the conditions of the experiment. It is very generally accepted that only experimentation can ensure the existence of cause-and-effect. Indeed, it is by manipulating the cause under experimental conditions (i.e., controlling other variables) that one can ensure a variation of the effect. Another major advantage of the simulation approach is that it helps with the discovery. Therefore, although the findings in this study are interesting and important from a theoretical point of view, further research must be done in terms of empirical studies with real data to confirm our findings.

Author Contributions

This doctoral research paper was done by K.E. under the supervision of F.A. Both authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Global Excellence and Stature (GES) 4.0 scholarship of the University of Johannesburg (UJ).

Acknowledgments

The authors would like to thank the editor and four anonymous referees for many constructive comments which greatly helped to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Albrecher, Hansjörg, and Jef L. Teugels. 2006. Exponential Behavior in the Presence of Dependence in Risk Theory. Journal of Applied Probability 43: 257–73. [Google Scholar] [CrossRef] [Green Version]
  2. Albrecher, Hansjörg, and Onno J. Boxma. 2004. A ruin model with dependence between claim sizes and claim intervals. Insurance: Mathematics and Economics 35: 245–54. [Google Scholar] [CrossRef] [Green Version]
  3. Albrecher, Hansjörg, and Onno J. Boxma. 2005. On the discounted penalty function in a Markov dependent risk model. Insurance: Mathematics and Economics 37: 650–72. [Google Scholar] [CrossRef] [Green Version]
  4. Asimit, Alexandru V., and Andrei L. Badescu. 2010. Extremes on the discounted aggregate claims in a time dependent risk model. Scandinavian Actuarial Journal, 93–104. [Google Scholar] [CrossRef]
  5. Boogaert, P., J. Haezendonck, and F. Delbaen. 1988. Limit theorems for the present value of the surplus of an insurance portfolio. Insurance: Mathematics and Economics 72: 131–38. [Google Scholar] [CrossRef]
  6. Boonen, Tim J. 2017. Solvency II solvency capital requirement for life insurance companies based on expected shortfall. European Actuarial Journal. [Google Scholar] [CrossRef] [Green Version]
  7. Boudreault, Mathieu, Helene Cossette, David Landriault, and Etienne Marceau. 2006. On a Risk Model with Dependence between inter-claim Arrivals and Claim Sizes. Scandinavian Actuarial Journal 5: 265–85. [Google Scholar] [CrossRef]
  8. Cai, Jun. 2002. Ruin probabilities with dependent rates of interest. Journal of Applied Probability 39: 312–23. [Google Scholar] [CrossRef]
  9. Cai, Jun, and David Dickson. 2003. Upper bounds for ultimate ruin probabilities in the Sparre Andersen model with interest. In Insurance: Mathematics and Economics. Amsterdam: Elsevier, vol. 32, pp. 61–71. [Google Scholar]
  10. Delbaen, Freddy, and Jean Haezendonck. 1987. Classical risk theory in an economic environment. Insurance: Mathematics and Economics 6: 85–116. [Google Scholar] [CrossRef]
  11. Embrechts, Paul, Claudia Klüppelberg, and Thomas Mikosch. 1997. Modelling of Extremal Events for Insurance and Finance. Berlin: Springer. [Google Scholar]
  12. Gatzert, Nadine, and Hannah Wesker. 2012. A Comparative Assessment of Basel II/III and Solvency II. The Geneva Papers 37: 539–70. [Google Scholar] [CrossRef] [Green Version]
  13. Geluk, Jaap, and Qihe Tang. 2009. Asymptotic Tail Probabilities of Sums of Dependent Subexponential Random Variables. Journal of Theoretical Probability 22: 871–82. [Google Scholar] [CrossRef] [Green Version]
  14. Hoiring, Dirk. 2013. Will Solvency II Market Risk Requirements Bite? The Impact of Solvency II on Insurers’ Asset Allocation. The Geneva Papers on Risk and Insurance Issues and Practice 38: 250–73. [Google Scholar] [CrossRef]
  15. Jorion, Philippe. 2001. Value at Risk, 2nd ed. New York: McGraw-Hill. [Google Scholar]
  16. Léveillé, Ghislain, and Emmanuel Hamel. 2018. Conditional, Non-Homogeneous and Doubly Stochastic Compound Poisson Processes with Stochastic Discounted Claims. Methodology and Computing in Applied Probability. [Google Scholar] [CrossRef]
  17. Léveillé, Ghislain, and Franck Adékambi. 2011. Covariance of discounted compound renewal sums with a stochastic interest rate. Scandinavian Actuarial Journal 2: 138–53. [Google Scholar] [CrossRef]
  18. Leveillé, Ghislain, and Franck Adékambi. 2012. Joint moments of discounted compound renewal sums. Scandinavian Actuarial Journal 1: 40–55. [Google Scholar] [CrossRef]
  19. Léveillé, Ghislain, and José Garrido. 2001a. Moments of compound renewal sums wit discounted claims. Insurance: Mathematics and Economics 28: 217–31. [Google Scholar]
  20. Léveillé, Ghislain, and José Garrido. 2001b. Recursive moments of compound renewal sums with discounted claims. Scandinavian Actuarial Journal 2: 98–110. [Google Scholar] [CrossRef]
  21. Léveillé, Ghislain, José Garrido, and Ya Fang Wang. 2010. Moment generating functions of compound renewal sums with discounted claims. Scandinavian Actuarial Journal 3: 165–84. [Google Scholar] [CrossRef]
  22. Liu, Xijun, and Qingwu Gao. 2015. Uniformly asymptotic behavior for the tail probability of discounted aggregate claims in the time-dependent risk model with upper tail asymptotically independent claims. Communications in Statistics Theory and Methods. [Google Scholar] [CrossRef]
  23. Lu, Dawei, Lixin Song, and Fuqi Li. 2018. Uniform asymptotics for discounted aggregate claims in dependent multi-risk model. In Communications in Statistics-Theory and Methods. Oxford: Taylor & Francis, vol. 48, pp. 781–93. [Google Scholar]
  24. Nelsen, Roger B. 1999. An Introduction to Copulas. New York: Springer. [Google Scholar]
  25. Ruwet, Christel. 2007. Processus de Poisson. Belgique: Université de Liège. [Google Scholar]
  26. Sklar, A. 1959. Fonctions de Répartition à n Dimensions et Leurs Marges. Paris: Publications de l’Institut de Statistique de l’ Université de Paris, vol. 8, pp. 229–31. [Google Scholar]
  27. Su, Chun, Zhishui Hu, Yu Chen, and Hanying Liang. 2007. A wide class of heavy-tailed distributions and its applications. Frontiers of Mathematics in China 2: 257–86. [Google Scholar] [CrossRef]
  28. Sundt, Bjørn, and Jozef L. Teugels. 1995. Ruin estimate under interest force. Insurance: Mathematics and Economics 16: 7–22. [Google Scholar] [CrossRef]
  29. Tang, Qithe. 2005. The Finite-time Ruin Probability of the Compound Poisson Model with Constant Interest Force. Journal of Applied Probability 42: 608–19. [Google Scholar] [CrossRef]
  30. Taylor, Gregory Clive. 1979. Probability of ruin under inflationary conditions or under experience rating. Astin Bulletin 10: 149–62. [Google Scholar] [CrossRef] [Green Version]
  31. Waters, Howard R. 1983. Probability of ruin for a risk process with claims cost inflation. Scandinavian Actuarial Journal 3: 148–64. [Google Scholar] [CrossRef]
  32. Wei, Li, and Hai-liang Yang. 2004. Explicit Expressions for the Ruin Probabilities of Erlang Risk Processes with Pareto Individual Claim. Acta Mathematicae Applicatae Sinica, English Series 20: 495–506. [Google Scholar] [CrossRef]
  33. Willmot, Gordon E. 1989. The total claims distribution under inflationary conditions. Scandinavian Actuarial Journal 10: 1–12. [Google Scholar] [CrossRef]
  34. Yang, Yang, Kaiyong Wang, and Dimitrios G. Konstantinides. 2014. Uniform asymptotics for discounted aggregate claims in dependent risk models. Journal of Applied Probability 51: 669–84. [Google Scholar] [CrossRef]
  35. Yang, Yang, Remigijus Leipus, and Jonas Šiaulys. 2012. Tail probability of randomly weighted sums of subexponential random variables under dependence structure. Statistics and Probability Letters 82: 1727–36. [Google Scholar] [CrossRef]
  36. Yuen, Kam C., Guojing Wang, and Rong Wu. 2006. On the renewal risk process with stochastic interest. Stochastic Processes and their Applications 116: 1496–510. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, Zhi-min, and Hu Yang. 2011. The compound Poisson risk model with dependence under a multi-layer dividend strategy. Applied Mathematics Journal, Chinese University 26: 1–13. [Google Scholar] [CrossRef]
Table 1. Asymptotic finite-time ruin probability for homogeneous Poisson risk model with interest rate.
Table 1. Asymptotic finite-time ruin probability for homogeneous Poisson risk model with interest rate.
α α = 2.5
θ θ = 0.5 θ = 0.75
TT = 20T = 40T = 20T = 40
b50100501005010050100
ψ δ A l i ( b , T ) 0.00740.00040.01450.00080.01120.00060.02180.0012
ψ δ F r a n k ( b , T ) 0.01850.00100.03650.00200.01710.00090.03390.0018
Table 2. Asymptotic ruin probability for mixed homogeneous Poisson risk model with interest rate.
Table 2. Asymptotic ruin probability for mixed homogeneous Poisson risk model with interest rate.
α α = 2.5
θ θ = 0.5 θ = 0.75
TT = 20T = 40T = 20T = 40
b50100501005010050100
ψ δ A l i ( b , T ) 0.01130.00060.02250.00120.01550.00080.03060.0016
ψ δ F r a n k ( b , T ) 0.02710.00150.05400.00290.02590.00140.05170.0028
Table 3. Asymptotic ruin probability for non-homogeneous Poisson risk model with interest rate.
Table 3. Asymptotic ruin probability for non-homogeneous Poisson risk model with interest rate.
θ θ = 0.5
α α = 2.5 α = 5
TT = 20T = 40T = 20T = 40
b50100501005010050100
ψ δ A l i ( b , T ) 0.03100.00170.06220.00330.02490.00130.04980.0027
ψ δ F r a n k ( b , T ) 0.03650.00190.07430.00400.03650.00200.07440.0040
Table 4. Asymptotic value-at-risk for homogeneous Poisson risk model with interest rate.
Table 4. Asymptotic value-at-risk for homogeneous Poisson risk model with interest rate.
θ θ = 0.5 θ = 0.75
Time horizon T = 10 T = 20 T = 40 T = 10 T = 20 T = 40
Ali VaR 99.5 % 46.524555.762866.2020752.364862.065573.0105
Frank VaR 99.5 % 59.493270.230782.088858.213568.935680.7418
Table 5. Asymptotic value-at-risk for mixed homogeneous Poisson risk model with interest rate.
Table 5. Asymptotic value-at-risk for mixed homogeneous Poisson risk model with interest rate.
θ θ = 0.5 θ = 0.75
Time horizon T = 10 T = 20 T = 40 T = 10 T = 20 T = 40
Ali VaR 99.5 % 51.969462.228173.526056.862667.329778.9272
Frank VaR 99.5 % 65.308876.787589.304464.464175.976988.4878
Table 6. Asymptotic value-at-risk for non-homogeneous Poisson risk model with interest rate.
Table 6. Asymptotic value-at-risk for non-homogeneous Poisson risk model with interest rate.
α α = 2.5 α = 5
Time horizon T = 10 T = 20 T = 40 T = 10 T = 20 T = 40
Ali VaR 99.5 % 67.152379.131992.018763.835675.310287.7768
Frank VaR 99.5 % 69.305282.040795.454769.358282.068595.4694
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adékambi, F.; Essiomle, K. Asymptotic Tail Probability of the Discounted Aggregate Claims under Homogeneous, Non-Homogeneous and Mixed Poisson Risk Model. Risks 2021, 9, 122. https://doi.org/10.3390/risks9070122

AMA Style

Adékambi F, Essiomle K. Asymptotic Tail Probability of the Discounted Aggregate Claims under Homogeneous, Non-Homogeneous and Mixed Poisson Risk Model. Risks. 2021; 9(7):122. https://doi.org/10.3390/risks9070122

Chicago/Turabian Style

Adékambi, Franck, and Kokou Essiomle. 2021. "Asymptotic Tail Probability of the Discounted Aggregate Claims under Homogeneous, Non-Homogeneous and Mixed Poisson Risk Model" Risks 9, no. 7: 122. https://doi.org/10.3390/risks9070122

APA Style

Adékambi, F., & Essiomle, K. (2021). Asymptotic Tail Probability of the Discounted Aggregate Claims under Homogeneous, Non-Homogeneous and Mixed Poisson Risk Model. Risks, 9(7), 122. https://doi.org/10.3390/risks9070122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop