Next Article in Journal
Premiums for Long-Term Care Insurance Packages: Sensitivity with Respect to Biometric Assumptions
Previous Article in Journal
Acknowledgement to Reviewers of Risks in 2015
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ruin Analysis of a Discrete-Time Dependent Sparre Andersen Model with External Financial Activities and Randomized Dividends

Department of Statistics & Actuarial Science, University ofWaterloo, 200 University Avenue West, Waterloo, ON N2L 3G1, Canada
*
Author to whom correspondence should be addressed.
Submission received: 18 September 2015 / Accepted: 20 January 2016 / Published: 3 February 2016

Abstract

:
We consider a discrete-time dependent Sparre Andersen risk model which incorporates multiple threshold levels characterizing an insurer’s minimal capital requirement, dividend paying situations, and external financial activities. We focus on the development of a recursive computational procedure to calculate the finite-time ruin probabilities and expected total discounted dividends paid prior to ruin associated with this model. We investigate several numerical examples and make some observations concerning the impact our threshold levels have on the finite-time ruin probabilities and expected total discounted dividends paid prior to ruin.

1. Introduction

The classical Cramér–Lundberg model is a foundational mathematical representation of an insurer’s surplus process in risk theory. Despite its tractability, however, the model has limitations in terms of applications. Certainly, more complex models are desirable in modern industrial settings. This paper strives to contribute to the ever-growing literature on insurance risk models that reflect more general and realistic modelling approaches to ruin theory.
Bruno de Finetti [1] first introduced the notion of a dividend strategy and the idea of finding an optimal dividend payment strategy for the insurance risk model. This was followed by numerous other researchers who further explored the problem in a variety of contexts (for reviews of the area, the interested reader is directed to Albrecher and Thonhauser [2] and Avanzi [3]). In particular, Drekic and Mera [4] considered the ruin analysis of a threshold-based dividend payment strategy in a discrete-time Sparre Andersen model. Their analysis was an extension of the work by Alfa and Drekic [5], in which a Sparre Andersen insurance risk model in discrete time was analyzed as a doubly-infinite Markov chain to establish a computational procedure for calculating the joint probability distribution of the time of ruin, the surplus immediately prior to ruin, and the deficit at ruin.
In this paper, we focus on the development of a recursive computational procedure to calculate the finite-time ruin probabilities and expected total discounted dividends paid prior to ruin associated with a model which generalizes the single threshold-based risk model introduced by Drekic and Mera [4]. In actual fact, three additional threshold levels are introduced to depict a minimum surplus level control strategy and external financial activities related to both investment and loan undertakings. Readers are referred to, for example, Cai and Dickson [6] and Li [7] for other general investment strategies found in insurance risk models, where the latter paper examined an insurance risk model with risky investments under the assumption that the risky assets follow a Wiener process, and the former paper considered a Markov chain based interest rate model. Korn and Wiese [8] studied optimal investment strategies in an insurance risk model where they also assumed that the risky assets follow a Wiener process.
The remainder of the paper is organized as follows. In Section 2, we introduce notation and specify the fundamental components underlying our threshold-based risk model. Section 3 details the derivation of a recursive formula (namely, Equation (9)) for the finite-time ruin probability associated with our proposed risk model and demonstrates the simplification of the result to that of Drekic and Mera [4]. Section 4 presents the derivation of a similar recursive formula (namely, Equation (13)) to compute the expected total discounted dividends paid prior to ruin and likewise demonstrates the simplification of the result to that of Drekic and Mera [4]. Finally, Section 5 discusses some numerical examples and related findings.

2. Model Description and Assumptions

For t N (with N being the set of non-negative integers), we define U t as the insurer’s amount of surplus at time t (measured in discrete monetary units) and F t as the amount of funds present at time t in the external fund of the insurer, a separate monetary account the insurer holds to better manage its reserve through both investment activities and loan undertakings. In actual fact, U t represents the surplus level at the end of the time interval ( t - 1 , t ] , t Z + (with Z + being the set of positive integers), at which point any premiums, deposits, claims, or withdrawals corresponding to this time interval have been received/paid out. We adopt the convention that premiums are received at ( t - 1 ) + and any claims are applied at t - .
In what follows, we assume that it is the insurer’s policy to pay out all the outstanding debt before resuming investment activities, and that the insurer first utilizes its investment assets to make any adjustments to its surplus level before engaging in loan activities. To differentiate between investment activities and loan activities with respect to the external fund, we split the support set of { F t , t N } into two disjoint sets, namely F t ( 0 , ) and F t [ β , 0 ] , where β is a non-positive integer representing the minimum support value of { F t , t N } . In fact, β is one of four threshold levels we feature in our model with the understanding that - β represents the borrowing limit of the insurer. When F t > 0 , F t represents the insurer’s investment activities in which interest is assumed to be earned at a constant rate of κ 1 > 0 per period. Conversely, when β F t 0 , F t represents the insurer’s loan activities and interest expense accumulates at a constant rate of κ 2 > 0 per period.
We next introduce the remaining three threshold levels (to be denoted by 1 , 2 , and 3 ) and explain how they, along with β, define our risk model. To aid in the understanding, let U t - and F t - represent the surplus and external fund levels, respectively, immediately after a claim instance but before a withdrawal instance. Herein, withdrawal refers to any cash inflow from the external fund to the surplus process (whereas deposit refers to any cash outflow from the surplus process to the external fund). Firstly, the threshold 1 represents the insurer’s minimum acceptable surplus level, and if U t - (corresponding to the time interval ( t - 1 , t ] ) is below 1 due to a claim, we withdraw or borrow from F t - to bring U t - up to level 1 . However, if F t - = β , then we can neither withdraw nor borrow more from F t - even if U t - is below 1 . In addition, if F t - (corresponding to the time interval ( t - 1 , t ] ) drifts below β due to interest expense accumulation, we use U t - to pay back the difference at t - as a form of deposit so that F t is at least kept at its minimum support value of β.
On the other hand, 2 is a trigger point for investment activities. If U t 2 , a constant deposit of size d is paid to the external fund at t + . Note that the deposit and withdrawal amounts are also stochastic in the sense that they are dependent on the surplus process, and that a deposit can be made at both the left and right limits of a time interval. We denote the left limit deposit amount corresponding to the time interval ( t , t + 1 ] to be D 1 , t , the right limit deposit amount corresponding to the time interval ( t - 1 , t ] to be D 2 , t , and the withdrawal amount corresponding to the time interval ( t - 1 , t ] to be W t .
Lastly, as in Drekic and Mera [4], if U t 3 , a random dividend is paid out to shareholders at t + . We denote the random dividend paid at time t by D t and assume that { D t , t N } is a conditionally independent and identically distributed (iid) sequence of random variables given { U t , t N } . As a final requirement, we assume that 1 2 3 .
To sum up, premiums and left limit deposits corresponding to the time interval ( t , t + 1 ] , t N , are collected and paid out at t + according to the following respective (random) rates:
P t = c if  U t < 3 , X t if  U t 3 ,
and
D 1 , t = 0 if  U t < 2 , d if  U t 2 ,
where X t = c - D t and b i = P r { X t = i } , i = c 1 , c 1 + 1 , , c 2 , denotes the probability mass function (pmf) of X t . We refer to c Z + as the pure (constant) premium and assume that c 1 , c 2 { d , d + 1 , , c } where d c , c 1 c 2 , and i = c 1 c 2 b i = 1 . Clearly, c 1 and c 2 are the respective lower and upper support values of the distribution of the random premium amount X t . Note that, by assumption, the probability distribution of X t is identical for all values of t N . Let μ = E { X 0 } denote the common mean.
Withdrawals and right limit deposits corresponding to the time interval ( t - 1 , t ] , t Z + , are made at t - according to the following respective (random) rates:
W t = 0 if  U t - 1 , min { 1 - U t - , max { 0 , F t - - β } } if  U t - < 1 ,
and
D 2 , t = max { 0 , β - F t - } .
For illustrative purposes, Figure 1 and Figure 2 depict an example of the simultaneous evolution of both the surplus process and that of the external fund.
Figure 1. Sample evolution of the surplus process { U t , t N } .
Figure 1. Sample evolution of the surplus process { U t , t N } .
Risks 04 00002 g001
Figure 2. Sample evolution of the external fund process { F t , t N } .
Figure 2. Sample evolution of the external fund process { F t , t N } .
Risks 04 00002 g002
Beginning at time 0 with an initial surplus level of v { 1 , 1 + 1 , } and an initial external fund amount of g N , the insurer’s amount of surplus at time t is expressible as
U t = v + i = 0 t - 1 P i - i = 0 t - 1 D 1 , i + i = 1 t W i - i = 1 t D 2 , i - i = 1 N t Y i , t N ,
where N t is the number of claims occurring by time t and individual claim sizes { Y i , i Z + } are assumed to form an iid sequence of positive, integer-valued random variables. We assume that the number of claims process is a discrete-time renewal process with independent, positive, integer-valued interclaim times { W i , i Z + } , where W i is the time between the ( i - 1 ) -th and i-th claims (with the understanding that the 0-th claim occurs at time 0). In particular, { W i , i Z + } forms an iid sequence of positive random variables with common pmf a k = P r { W i = k } , k = 1 , 2 , , n a where n a Z + , and corresponding survival function A k = P r { W i > k } = 1 - j = 1 k a j . Furthermore, we assume that the pairs { ( W i , Y i ) , i Z + } are iid, so that the joint pmf of ( W i , Y i ) is given by
P r { W i = k , Y i = j } = a k α j ( k ) ,
where α j ( k ) = P r { Y i = j | W i = k } denotes the conditional pmf of Y i given W i . Such a structurization allows for possible dependence between interclaim times and claim sizes in Sparre Andersen models (e.g., see Cheung et al. [9]).

3. Calculation of Finite-Time Ruin Probabilities

We begin by examining the finite-time ruin probabilities associated with the risk model defined by Equation (1). To start with, ruin occurs if and only if U t < 0 for some t Z + and we denote T to be the time of ruin. In other words, T = min { t Z + | U t < 0 } with T = if U t 0 t Z + .
In what follows, we are interested in computing the quantity
Ψ ( v , g , n ) = P r { T n | U 0 = v , F 0 = g } = 1 - P r { T > n | U 0 = v , F 0 = g } , n N ,
which we refer to as the finite-time ruin probability. To aid in the computation of this quantity, we introduce the following related function:
σ ( u , f , n , m ) = P r { T > n | U 0 = u , F 0 = f , M 0 = m } , n N , u Z , f { β , β + 1 , } ,
where Z denotes the set of all integers and M t , referred to as the elapsed waiting time counter, represents the elapsed time at time t since the most recent claim occurrence and its values lie in the set { 0 , 1 , , n a - 1 } . With the introduction of this function, we remark that Ψ ( v , g , n ) = 1 - σ ( v , g , n , 0 ) .
First of all, assuming the occurrence of no claims and no right limit deposits, we need to identify when U t 2 and U t 3 for the first time. We introduce two functions to denote these time points, namely:
z t , u = 0 if  u 2 , min { 2 - u - 1 c + 1 , t } if  u < 2 ,
and
z t , u = 0 if  u 3 , min { 3 - u - c z t , u - 1 c - d + z t , u + 1 , t } if  u < 3 ,
where x , referred to as the floor function of x, yields the largest integer less than or equal to x.
To aid us in obtaining a mathematical expression for σ ( u , f , n , m ) , we have to examine how the process { F t , t N } evolves over time. Let us first assume that there are no claims or withdrawals to consider. Clearly, F t is a non-decreasing function of t if f N . On the other hand, if f N , F t could either be a strictly decreasing function of t or perhaps a convex function depending on the values of f, κ 2 and d. Consequently, if F t drifts below β, a deposit is forced to be made and this may cause ruin. Thus, in this model, ruin can occur due to either a claim or a deposit. This certainly adds more complexity in deriving a formula for σ ( u , f , n , m ) , and as a result, we have to introduce a few more functions. One such function is denoted by o t , u , f , representing the time point s { 1 , 2 , , t } at which F s is set to become greater than or equal to 0 for the first time. Obtaining this value is not difficult since, assuming the occurrence of no claims and that F s β s o t , u , f , F s becomes non-stochastic, the form of which we denote by:
F ˜ s , u , f = f ( 1 + κ 2 ) s + d z s , u , s κ 2 if  f < 0 , f ( 1 + κ 1 ) s + d z s , u , s κ 1 if  f 0 ,
where d k , l κ i , k , l N , represents the future value of deposits made at times k + , ( k + 1 ) + , , ( l - 1 ) + with respect to interest rate κ i per period, i = 1 , 2 . Clearly, we have d k , l κ i = 0 for l k and
d k , l κ i = d ( 1 + κ i ) l - k + d ( 1 + κ i ) l - k - 1 + + d ( 1 + κ i ) = d ( 1 + κ i ) [ ( 1 + κ i ) l - k - 1 ] κ i , l > k .
It subsequently follows that
o t , u , f = t if  F ˜ i + 1 , u , f < 0 i { 0 , 1 , , t } , min { i { 0 , 1 , , t } | F ˜ i + 1 , u , f 0 } otherwise .
In defining o t , u , f above, we utilize the value of F ˜ i + 1 , u , f instead of F ˜ i , u , f . This is because for f < 0 , the function F ˜ t , u , f , t N , up-crosses level 0 only if a positive amount of deposit is made to the external fund. We stated earlier that left deposits are made at the left limit point of a discrete-time interval. Thus, F ˜ i + 1 , u , f 0 for the first time implies that at time i + , there was a left deposit made to the external fund.
With the introduction of o t , u , f , we henceforth express the non-stochastic form of F t as
F ^ t , u , f = ( f ( 1 + κ 2 ) o t , u , f + d z o t , u , f , u , o t , u , f κ 2 ) ( 1 + κ 1 ) t - o t , u , f + d max { z t , u , o t , u , f } , t κ 1 , t N .
Note that Equation (2) involves the use of the floor function to calculate the (non-stochastic) value of the external fund at time t. Such an assumption can be viewed as conservative in nature, since any non-integer value of the external fund (which can arise due to interest accumulation) is essentially rounded down.
There is another important function we introduce next, as it represents the earliest time point when F t falls below β (again assuming the occurrence of no claims) due to interest expense accumulation. We denote this time point by c t , m , u , f and refer to it as a calling point. It is given by
c t , m , u , f = min { n a - m , t } if  F ^ i , u , f β i { 1 , 2 , , min { n a - m , t } } , min { i { 1 , 2 , , min { n a - m , t } } | F ^ i , u , f < β } otherwise .
Note that the above function depends on both t and m. As introduced earlier in this section, m represents the elapsed time at time 0 since the most recent claim occurrence.
With these preliminaries in place, we adopt the principle of conditioning on the first claim time as in Cossette et al. [10] or Drekic and Mera [4]. Measured from our initial time point (which we label as time 0), the lower limit of the time until the first claim occurs is 1, but its pmf is now conditional on the value of m. Morever, in evaluating σ ( u , f , n , m ) , we condition on first claim times ranging from 1 up to c n , m , u , f , and on the event that the time until the first claim occurs is greater than c n , m , u , f . For first claim time instances which take place at or before the calling point, the recursive process used is very similar to that of Drekic and Mera [4]. However, in the event that the time until the first claim occurs is greater than c n , m , u , f , the recursive process is performed differently. By doing so, we are essentially denoting c n , m , u , f to be the “new” initial time point, updating the parameters of the function σ, and proceeding with the recursive process. We further explain this situation after introducing some necessary boundary conditions for σ ( u , f , n , m ) , namely:
σ ( u , f , n , m ) = 0 if  u Z -   or   m = n a , 1 if  u N , n = 0 ,   and   m = 0 , 1 , , n a - 1 ,
where Z - in Equation (3) denotes the set of negative integers. By conditioning on the events outlined above, we get
σ ( u , f , n , m ) = k = 1 c n , m , u , f a k + m A m P r { T > n | U 0 = u , F 0 = f , M 0 = m , W 1 ( m ) = k } + A c n , m , u , f + m A m P r { T > n | U 0 = u , F 0 = f , M 0 = m , W 1 ( m ) > c n , m , u , f } ,
where W 1 ( m ) is the duration from our initial time point until the first claim occurs given that the elapsed waiting time at time 0 since the most recent claim is m.
At time k { 1 , 2 , , c n , m , u , f } , the elapsed waiting time counter is reset to 0 for the next recursion, n is reduced by k, and the “new” initial surplus and external fund amounts are determined by the size of the incurred claim and the premiums received up to time k. Specifically, we obtain
P r { T > n | U 0 = u , F 0 = f , M 0 = m , W 1 ( m ) = k } = l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u - d ( k - z k , u ) + l + F ^ k , u , f - β α j ( k + m ) σ ( u * , f * , n - k , 0 ) ,
where
u * = u + c z k , u - d ( k - z k , u ) + l - j + min { F ^ k , u , f - β , max { 0 , j - [ u + c z k , u - d ( k - z k , u ) + l - 1 ] } } ,
f * = max { β , min { F ^ k , u , f , F ^ k , u , f - j + [ u + c z k , u - d ( k - z k , u ) + l - 1 ] } } ,
and l denotes the value of the sum of the random premiums received up to time k, with corresponding pmf b l , k - z k , u representing the ( k - z k , u ) -fold convolution of b l with itself. To evaluate the pmf b l , r , r N , we define b l , 0 = δ l , 0 (where δ i , j , in general, denotes the Kronecker delta function of i and j),
b l , 1 = b l if  l = c 1 , c 1 + 1 , , c 2 , 0 otherwise ,
and for r = 2 , 3 , ,
b l , r = j = c 1 c 2 b j , 1 b l - j , r - 1 if  l = r c 1 , r c 1 + 1 , , r c 2 , 0 otherwise .
The reasoning behind the definitions of the above parameters is that we first consider whether the claim size is substantial enough for the surplus process to fall below its minimum support level 1 . If so, then j - [ u + c z k , u - d ( k - z k , u ) + l - 1 ] is a positive quantity and we consider whether the external fund is able to support the surplus process. We do this by comparing F ^ k , u , f - β and j - [ u + c z k , u - d ( k - z k , u ) + l - 1 ] , and choosing the minimum of the two quantities to ensure that the external fund does not fall below the maximum level of external funding allowed, β. If j - [ u + c z k , u - d ( k - z k , u ) + l - 1 ] is a non-positive quantity, then the surplus process is greater than or equal to 1 after the claim, in which case, we only need to consider whether F ^ k , u , f is below β. If so, F ^ k , u , f - β is less than 0, and we would subtract | F ^ k , u , f - β | from the surplus process and add it to the external fund to bring it up to β.
In situations when W 1 ( m ) > c n , m , u , f , we perform a recursion at c n , m , u , f to similarly acquire
P r { T > n | U 0 = u , F 0 = f , M 0 = m , W 1 ( m ) > c n , m , u , f } = l = ( c n , m , u , f - z c n , m , u , f , u ) c 1 ( c n , m , u , f - z c n , m , u , f , u , f ) c 2 b l , c n , m , u , f - z c n , m , u , f , u , f σ ( u , f , n - c n , m , u , f , c n , m , u , f + m ) ,
where
u = u + c z c n , m , u , f , u - d ( c n , m , u , f - z c n , m , u , f , u ) + l + min { 0 , F ^ c n , m , u , f , u , f - β }
and
f = max { F ^ c n , m , u , f , u , f , β } .
We remark that when W 1 ( m ) > c n , m , u , f , there is no claim size to consider at time c n , m , u , f . Thus, all we need to account for is whether F ^ c n , m , u , f , u , f falls below β. In this case, just enough funds would be withdrawn from the surplus process and added to the external fund to bring it up to β. However, note that F ^ c n , m , u , f , u , f may not necessarily be below β. If F ^ c n , m , u , f , u , f β , then c n , m , u , f = min { n a - m , n } and this implies that either n - c n , m , u , f = 0 or c n , m , u , f + m = n a in Equation (6). This yields an interesting outcome. Given that W 1 ( m ) c n , m , u , f and F ^ c n , m , u , f , u , f β , it must be that u 0 at time c n , m , u , f . However, if c n , m , u , f + m = n a , then σ ( u , f , n - c n , m , u , f , n a ) is set equal to 0 via Equation (3). On the other hand, if n < n a - m so that c n , m , u , f = n , then σ ( u , f , 0 , c n , m , u , f + m ) = 1 from Equation (3). Putting it altogether, we establish the following final formula for σ ( u , f , n , m ) :
σ ( u , f , n , m ) = k = 1 c n , m , u , f a k + m A m l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u - d ( k - z k , u ) + l + F ^ k , u , f - β α j ( k + m ) σ ( u * , f * , n - k , 0 ) + A c n , m , u , f + m A m l = ( c n , m , u , f - z c n , m , u , f , u ) c 1 ( c n , m , u , f - z c n , m , u , f , u ) c 2 b l , c n , m , u , f - z c n , m , u , f , u σ ( u , f , n - c n , m , u , f , c n , m , u , f + m ) .
Note that the determination of σ ( u , f , n , m ) via Equation (9) requires a double recursion in both n and m, with boundary conditions given by Equation (3) serving as the starting point. Moreover, if we assume that f = 0 , d = 0 , 1 = 0 , β = 0 , m = 0 , and α j ( k ) = P r { Y i = j | W i = k } = P r { Y i = j } = α j (i.e., Y i is independent of W i i Z + ), then the model under consideration is equivalent to the independent Sparre Andersen model studied by Drekic and Mera [4]. To verify this, we first observe that c n , 0 , u , 0 = min { n , n a } n Z + . If n < n a , then σ ( u , f , 0 , n ) = 1 and
A n A 0 P r ( T > n | U 0 = u , F 0 = 0 , M 0 = 0 , W 1 ( 0 ) > n ) = A n l = ( n - z n , u ) c 1 ( n - z n , u ) c 2 b l , n - z n , u = A n .
Conversely, if n n a , then σ ( u , f , n - n a , n a ) = 0 and
A n a A 0 P r ( T > n | U 0 = u , F 0 = 0 , M 0 = 0 , W 1 ( 0 ) > n a ) = A n a × 0 = 0 .
Thus, Equation (9) simplifies to become
σ ( u , 0 , n , 0 ) = A n + k = 1 min { n , n a } a k l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u + l α j σ ( u + c z k , u + l - j , 0 , n - k , 0 ) ,
which is consistent with the result in Drekic and Mera [4], p. 744.

4. Calculation of Expected Total Discounted Dividend Payments

Our next objective is to derive a corresponding recursive formula to compute the expected total discounted dividend payments made prior to ruin. The approach we employ essentially borrows from that of Dickson and Waters [11], Section 5. Let E { D v , g } denote the expected total discounted (i.e., to time 0 according to discount factor ν ( 0 , 1 ) per unit of time) dividends paid prior to ruin, where the random variable D v , g represents the total discounted dividends paid before ruin starting from an initial surplus of v and an initial level of g in the external fund. Moreover, we also introduce the analogous quantity E { D v , g , n } as the expected total discounted dividends paid before ruin occurs or strictly before time n Z + , whichever happens first.
In order to calculate E { D v , g } , we develop a computational procedure for calculating E { D v , g , n } and then use the fact that E { D v , g , n } E { D v , g } as n . To aid in the computation of E { D v , g , n } , we introduce a function (similar in nature to σ from the previous section) defined by
V ( u , f , n , m ) = E { D u , f , n | M 0 = m } , n Z + , u Z , f { β , β + 1 , } , m { 0 , 1 , , n a - 1 } .
Clearly, E { D v , g , n } = V ( v , g , n , 0 ) . As with the function σ in the previous section, the function V ( u , f , n , m ) also has its own set of boundary conditions, namely:
V ( u , f , n , m ) = 0 if  u Z -   or   m = n a , 0 if  0 u < 3   and   n = 1 , c - μ if  u 3   and   n = 1 .
We employ a similar approach as in the previous section by conditioning on values of W 1 ( m ) ranging from 1 up to c n - 1 , m , u , f , as well as the case when W 1 ( m ) > c n - 1 , m , u , f . By conditioning on these events, we immediately obtain
V ( u , f , n , m ) = k = 1 c n - 1 , m , u , f a k + m A m E { D u , f , n , m | W 1 ( m ) = k } + A c n - 1 , m , u , f + m A m E { D u , f , n , m | W 1 ( m ) > c n - 1 , m , u , f } .
For k { 1 , 2 , , c n - 1 , m , u , f } , an expected dividend payment of amount c - μ would occur at times z k , u + , ( z k , u + 1 ) + , , ( k - 1 ) + , followed by possible future dividend payments (starting from time k) once the initial claim is applied. Applying the appropriate conditioning arguments ultimately leads to
E { D u , f , n , m | W 1 ( m ) = k } = i = z k , u k - 1 ν i ( c - μ ) + ν k l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u - d ( k - z k , u ) + l + F ^ k , u , f - β α j ( k + m ) V ( u * , f * , n - k , 0 ) ,
where u * and f * are given by Equations (4) and (5), respectively. For W 1 ( m ) > c n - 1 , m , u , f , however, we need to reset the parameters of the function V as we did in our treatment of the finite-time ruin probabilities and base the recursion at c n - 1 , m , u , f . We also need to account for the expected dividend payments received at times z c n - 1 , m , u , f , u + , ( z c n - 1 , m , u , f , u + 1 ) + , , ( c n - 1 , m , u , f - 1 ) + . Thus, we eventually obtain
E { D u , f , n , m | W 1 ( m ) > c n - 1 , m , u , f } = i = z c n - 1 , m , u , f c n - 1 , m , u , f - 1 ν i ( c - μ ) + ν c n - 1 , m , u , f × l = ( c n - 1 , m , u , f - z c n - 1 , m , u , f , u ) c 1 ( c n - 1 , m , u , f - z c n - 1 , m , u , f , u ) c 2 b l , c n - 1 , m , u , f - z c n - 1 , m , u , f , u V ( u ^ , f ^ , n - c n - 1 , m , u , f , c n - 1 , m , u , f + m ) ,
where
u ^ = u + c z c n - 1 , m , u , f , u - d ( c n - 1 , m , u , f - z c n - 1 , m , u , f , u ) + l + min { 0 , F ^ c n - 1 , m , u , f , u , f - β } ,
and
f ^ = max { F ^ c n - 1 , m , u , f , u , f , β } .
Note that u ^ and f ^ are identical in form to Equations (7) and (8), respectively, with n simply replaced by n - 1 . Combining Equations (11) and (12) eventually yields the final overall formula
V ( u , f , n , m ) = k = 1 c n - 1 , m , u , f a k + m A m [ ( c - μ ) ν z k , u - ν k 1 - ν + ν k l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u - d ( k - z k , u ) + l + F ^ k , u - β α j ( k + m ) V ( u * , f * , n - k , 0 ) ] + A c n - 1 , m , u , f + m A m [ ( c - μ ) ν z c n - 1 , m , u , f , u - ν c n - 1 , m , u , f 1 - ν + ν c n - 1 , m , u , f × l = ( c n - 1 , m , u , f - z c n - 1 , m , u , f , u ) c 1 ( c n - 1 , m , u , f - z c n - 1 , m , u , f , u ) c 2 b l , c n - 1 , m , u , f - z c n - 1 , m , u , f , u V ( u ^ , f ^ , n - c n - 1 , m , u , f , c n - 1 , m , u , f + m ) ] .
Similar to the function σ of the previous section, we note that the use of Equation (13) requires a double recursion in both n and m, with boundary conditions given by Equation (10) serving as the starting point. As before, we consider the situation when f = 0 , d = 0 , 1 = 0 , β = 0 , m = 0 , and α j ( k ) = α j , so that the model under consideration is equivalent to the one analyzed by Drekic and Mera [4]. In this case, Equation (13) reduces to
V ( u , 0 , n , 0 ) = k = 1 min { n - 1 , n a } a k [ ( c - μ ) ν z k , u - ν k 1 - ν + ν k l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u + l α j V ( u + c z k , u + l - j , 0 , n - k , 0 ) ] + A min { n - 1 , n a } [ ( c - μ ) ν z n - 1 , u - ν n - 1 1 - ν + ν n - 1 l = ( n - 1 - z n - 1 , u ) c 1 ( n - 1 - z n - 1 , u ) c 2 b l , n - 1 - z n - 1 , u V ( u + c z n - 1 , u + l , 0 , 1 , n - 1 ) ] .
We remark that since c n - 1 , 0 , u , 0 = min { n - 1 , n a } , the square-bracketed term in Equation (14) that is pre-multiplied by A min { n - 1 , n a } matters only if c n - 1 , 0 , u , 0 = n - 1 since A n a = 0 . Thus, for convenience, we set c n - 1 , 0 , u , 0 = n - 1 inside this square-bracketed term in Equation (14). In addition, Equation (14) simplifies further to become:
V ( u , 0 , n , 0 ) = k = 1 min { n - 1 , n a } a k [ ( c - μ ) ν z k , u - ν k 1 - ν + ν k l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u + l α j V ( u + c z k , u + l - j , 0 , n - k , 0 ) ] if  z n - 1 , u < z n , u , k = 1 min { n - 1 , n a } a k [ ( c - μ ) ν z k , u - ν k 1 - ν + ν k l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u + l α j V ( u + c z k , u + l - j , 0 , n - k , 0 ) ] + A min { n - 1 , n a } i = z n - 1 , u n - 1 ν i ( c - μ ) if  z n - 1 , u = z n , u .
The logic behind Equation (15) is as follows. If z n - 1 , u < z n , u , then z n , u = n . This implies that there are no dividend payments before time n, and hence, Equation (12) becomes equal to 0. On the other hand, if z n - 1 , u = z n , u , then there is a guaranteed dividend payment at time n - 1 and
l = ( n - 1 - z n - 1 , u ) c 1 ( n - 1 - z n - 1 , u ) c 2 b l , n - 1 - z n - 1 , u V ( u + c z n - 1 , u + l , 0 , 1 , n - 1 ) = c - μ .
Thus, the quantity
A min { n - 1 , n a } i = z n - 1 , u n - 1 ν i ( c - μ )
in the second component of Equation (15) can be rewritten as
A n - 1 i = z n , u n - 1 ν i ( c - μ ) .
In addition, consider the expression
k = n n a a k i = z k , u n - 1 ν i ( c - μ ) .
If z n , u = z n + 1 , u = = z n a , u , then certainly we can replace z n + 1 , u , z n + 2 , u , , z n a , u with z n , u in Equation (16). Otherwise, z n , u = n , and since min { z n , u , z n + 1 , u , , z n a , u } = z n , u , Equation (16) evaluates to 0. In either case, Equation (16) becomes
k = n n a a k i = z k , u n - 1 ν i ( c - μ ) = k = n n a a k i = z n , u n - 1 ν i ( c - μ ) = A n - 1 i = z n , u n - 1 ν i ( c - μ ) .
Substituting these results into Equation (15) and simplifying ultimately yields
V ( u , 0 , n , 0 ) = ( c - μ ) k = 1 n a a k ν min { z k , u , n } - ν min { k , n } 1 - ν + k = 1 min { n - 1 , n a } a k ν k l = ( k - z k , u ) c 1 ( k - z k , u ) c 2 b l , k - z k , u j = 1 u + c z k , u + l α j V ( u + c z k , u + l - j , 0 , n - k , 0 ) .
Once again, the above formula coincides with that of Drekic and Mera ([4], p. 745).

5. Numerical Results

In this section, we implement our key recursive formulas (namely, Equations (9) and (13)) and investigate the behaviour of our proposed risk model through some numerical examples. To the best of our knowledge, we are unaware of any alternative computational methods to calculate these two ruin-related quantities of interest for the kind of discrete-time risk model considered in this paper. In fact, most methods in the literature have employed conditioning arguments similar to ours to analyze discrete-time risk models and ultimately formulate a recursive procedure to calculate such quantities.
In what follows, we focus specifically on the independent Sparre Andersen model (i.e., α j ( k ) = α j ) and introduce a set of interclaim time distributions to study, namely:
(a) a j = ( 2 / 11 ) ( 9 / 11 ) j - 1 if  j = 1 , 2 , , 24 , ( 9 / 11 ) 24 if  j = 25 , (b) a j = 1 / 10 , j = 1 , 2 , , 10 , (c) a j = 1 1 - ( 39 / 50 ) 25 25 j ( 11 / 50 ) j ( 39 / 50 ) 25 - j , j = 1 , 2 , , 25 , (d) a j = ( 0 . 645 ) ( 1 / 2 ) j + ( 0 . 355 ) ( 1 / 12 ) ( 11 / 12 ) j - 1 if  j = 1 , 2 , , 14 , ( 0 . 645 ) ( 1 / 2 ) 14 + ( 0 . 355 ) ( 1 / 12 ) ( 11 / 12 ) 14 if  j = 15 , ( 0 . 355 ) ( 1 / 12 ) ( 11 / 12 ) j - 1 if  j = 16 , 17 , , 49 , ( 0 . 355 ) ( 11 / 12 ) 49 if  j = 50 .
We observe that (a) is the pmf of a truncated geometric distribution with n a = 25 ; (b) is the pmf of a uniform distribution on { 1 , 2 , , 10 } ; (c) is the pmf of a zero-truncated binomial distribution with n a = 25 ; and (d) is the pmf of a mixture of two truncated geometric distributions with n a = 50 . We note that the means are essentially equal to 5.5 for all four interclaim time distributions, but their variability differs with (c) being the least variable and (d) being the most variable.
As for the random premium distribution, in effect, we consider a degenerate distribution with all the probability mass on 2 (i.e., c 1 = c 2 = 2 , so that b 2 = 1 ). In terms of the claim size distribution, we consider a discretized version of the Pareto distribution with mean 10 and pmf of the form
α j = 1 + j - 1 30 - 4 - 1 + j 30 - 4 , j Z + .
Incidentally, the above claim size pmf results from the application of the “lower bound” discretization method of Dickson ([12], p. 79) to approximate a Pareto distribution which is continuous. In fact, similar discretization ideas (e.g., see Dickson et al. [13] or Alfa and Drekic [5], pp. 306–307) can be employed to construct a discrete-time approximation to a continuous-time analogue of our proposed risk model. We also remark that the above set of interclaim time, random premium, and claim size distributions are all taken from Section 4 of Drekic and Mera [4]. In all our examples, we set v = 10 , g = 0 , c = 5 , κ 1 = 0 . 01 , d = 1 , and ν = 0 . 75 .
We make the following observations concerning the results in Table 1 through Table 3, in which interclaim time distribution (a) was used throughout:
(i)
In Table 1, we assumed that κ 2 = 0 . 02 , 1 = 0 , 2 = 20 , and 3 = 50 . Under these circumstances, changing the maximal level of external funding allowed resulted in a monotone behaviour in our two performance measures. As we decreased the value of β, the finite-time ruin probabilities decreased monotonically for all n 150 , whereas the expected total discounted dividends paid before ruin increased monotonically. It seems that the benefit of having more funds available outweighs the borrowing costs under the specific setting considered here. We point out that in Table 1 (as well as Table 2 through Table 5), the minimum time point n required to achieve convergence (to six significant digits) of V ( 10 , 0 , n , 0 ) to E { D 10 , 0 } is italicized and appears in parentheses next to its corresponding value.
(ii)
In Table 2, we assumed that κ 2 = 0 . 02 , β = - 10 , 2 = 25 , and 3 = 50 . Changing the minimal capital requirement level 1 resulted in a negative effect on the finite-time ruin probabilities. As we increased 1 , the finite-time ruin probabilities monotonically increased for all n 150 . In relation, the expected total discounted dividends paid before ruin increased as the value of 1 rose. Artificially requiring the level of the surplus process to be at a certain positive level prompts more borrowing and this generates higher interest expense. Thus, ruin is more likely to occur. The expected total discounted dividends paid prior to ruin increased since the surplus process is now more likely to reach the dividend payment trigger level 3 , as the surplus level is maintained at a higher level more often.
(iii)
In Table 3, we assumed that κ 2 = 0 . 02 , 1 = 0 , β = - 10 , and 3 = 50 . Raising the investment trigger level 2 resulted in increasing both the finite-time ruin probabilities and the expected total discounted dividend payments prior to ruin. As we increase 2 , we are delaying investments and this leads to a negative effect on the finite-time ruin probabilities since the external fund earns interest while the surplus process does not. Nevertheless, the surplus process is kept at a higher level as we increase 2 , and thus, the surplus process is more likely to reach 3 , resulting in higher expected total discounted dividend payments prior to ruin.
The following remarks are made concerning the results in Table 4, in which interclaim time distributions (a) to (d) were each studied:
(iv)
We assumed that κ 2 = 0 . 02 , 1 = 0 , 2 = 20 , 3 = 50 , and β = - 10 . In an effort to investigate the effects of variability in the choice of interclaim time distribution, we observed that the finite-time ruin probabilities were highest for interclaim time distribution (d) and lowest for interclaim time distribution (c) for all n 150 . The expected total discounted dividends paid prior to ruin ended up being highest for (a) and lowest for (d).
Finally, we make the following observations concerning the results in Table 5, in which interclaim time distribution (b) was used throughout:
(v)
In (i), we observed that the ability to borrow more money from the external fund had positive effects on both the finite-time ruin probabilities and expected total discounted dividends paid prior to ruin. This begs the question as to whether an insurer can continue to borrow more and more money and still produce a positive impact on the business. To investigate this matter further, we increased κ 2 to 0.30 and varied β from - 5 to - 30 in increments of size 5. Under this revised setting, not only has the monotone behaviour of the finite-time ruin probabilities changed, but the effects of decreasing the value of β have also changed. From Table 5, note that for n = 25 and onwards, β = - 10 produced the lowest ruin probabilities whereas β = - 25 yielded the highest ruin probabilities. On the other hand, the expected total discounted dividends paid before ruin were still highest for β = - 30 and lowest for β = - 5 , although the difference was rather minimal.
Table 1. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (a) with 1 = 0 , 2 = 20 , 3 = 50 , and κ 2 = 0 . 02
Table 1. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (a) with 1 = 0 , 2 = 20 , 3 = 50 , and κ 2 = 0 . 02
Ψ ( 10 , 0 , 25 ) Ψ ( 10 , 0 , 50 ) Ψ ( 10 , 0 , 75 ) Ψ ( 10 , 0 , 100 ) Ψ ( 10 , 0 , 150 ) E { D 10 , 0 }
β = 0 0.1748300.1966140.2046720.2078230.209558 0 . 248444 ( 60 )
β = - 4 0.1440860.1646620.1724980.1756380.177413 0 . 250317 ( 56 )
β = - 8 0.1199480.1393030.1468910.1500070.151815 0 . 251686 ( 59 )
β = - 12 0.1007260.1188990.1262300.1293160.131157 0 . 252692 ( 62 )
β = - 16 0.08522870.1022740.1093530.1124080.114284 0 . 253445 ( 67 )
β = - 20 0.07263600.08862590.09546790.09849810.100416 0 . 254016 ( 67 )
Table 2. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (a) with 2 = 25 , 3 = 50 , β = - 10 , and κ 2 = 0 . 02
Table 2. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (a) with 2 = 25 , 3 = 50 , β = - 10 , and κ 2 = 0 . 02
Ψ ( 10 , 0 , 25 ) Ψ ( 10 , 0 , 50 ) Ψ ( 10 , 0 , 75 ) Ψ ( 10 , 0 , 100 ) Ψ ( 10 , 0 , 150 ) E { D 10 , 0 }
1 = 0 0.1103500.1300560.1380290.1413930.143408 0 . 272456 ( 64 )
1 = 5 0.1115210.1318930.1402980.1439280.146174 0 . 274604 ( 56 )
1 = 10 0.1128730.1340820.1430260.1469880.149529 0 . 279519 ( 56 )
1 = 15 0.1145780.1368890.1465270.1509180.153847 0 . 291398 ( 59 )
1 = 20 0.1161510.1398170.1503250.1552670.158714 0 . 315121 ( 56 )
Table 3. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (a) with 1 = 0 , 3 = 50 , β = - 10 , and κ 2 = 0 . 02
Table 3. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (a) with 1 = 0 , 3 = 50 , β = - 10 , and κ 2 = 0 . 02
Ψ ( 10 , 0 , 25 ) Ψ ( 10 , 0 , 50 ) Ψ ( 10 , 0 , 75 ) Ψ ( 10 , 0 , 100 ) Ψ ( 10 , 0 , 150 ) E { D 10 , 0 }
2 = 5 0.1085460.1253090.1317340.1343180.135778 0 . 208521 ( 60 )
2 = 10 0.1087250.1257720.1323660.1350420.136572 0 . 210858 ( 57 )
2 = 15 0.1092650.1271350.1341370.1370110.138677 0 . 231600 ( 59 )
2 = 30 0.1113360.1322430.1408580.1445560.146820 0 . 331537 ( 56 )
2 = 45 0.1134480.1384420.1497530.1550120.158566 0 . 417485 ( 56 )
Table 4. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin when β = - 10 , 1 = 0 , 2 = 20 , 3 = 50 , and κ 2 = 0 . 02
Table 4. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin when β = - 10 , 1 = 0 , 2 = 20 , 3 = 50 , and κ 2 = 0 . 02
Ψ ( 10 , 0 , 25 ) Ψ ( 10 , 0 , 50 ) Ψ ( 10 , 0 , 75 ) Ψ ( 10 , 0 , 100 ) Ψ ( 10 , 0 , 150 ) E { D 10 , 0 }
(a)0.1098110.1285690.1360290.1391310.140955 0 . 252225 ( 59 )
(b)0.07397370.08937410.09551450.09805380.0995368 0 . 249026 ( 69 )
(c)0.05778120.07192400.07759490.07994440.0813180 0 . 247518 ( 56 )
(d)0.1985210.2255330.2365860.2413740.244328 0 . 227710 ( 58 )
Table 5. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (b) with 1 = 0 , 2 = 20 , 3 = 50 , and κ 2 = 0 . 30
Table 5. Finite-time ruin probabilities and expected total discounted dividends paid prior to ruin corresponding to interclaim time distribution (b) with 1 = 0 , 2 = 20 , 3 = 50 , and κ 2 = 0 . 30
Ψ ( 10 , 0 , 25 ) Ψ ( 10 , 0 , 50 ) Ψ ( 10 , 0 , 75 ) Ψ ( 10 , 0 , 100 ) Ψ ( 10 , 0 , 150 ) E { D 10 , 0 }
β = - 5 0.09733310.1167640.1254170.1299950.134825 0 . 247088 ( 58 )
β = - 10 0.08440380.1075210.1189640.1257260.133511 0 . 247093 ( 56 )
β = - 15 0.09310860.1284450.1400550.1440020.145997 0 . 247098 ( 56 )
β = - 20 0.09839120.1313070.1407460.1441960.146018 0 . 247106 ( 56 )
β = - 25 0.1036490.1322890.1410040.1442800.146028 0 . 247121 ( 59 )
β = - 30 0.1036290.1321960.1409670.1442670.146026 0 . 247143 ( 51 )

6. Conclusions

In this paper, we considered a discrete-time dependent Sparre Andersen risk model featuring multiple threshold levels in an effort to characterize an insurer’s minimal capital requirement, dividend paying scenarios, and external financial activities. In analyzing this model, we developed recursive computational procedures to calculate two particular performance measures of interest, namely finite-time ruin probabilities and expected total discounted dividends paid prior to ruin. Through a variety of numerical experiments performed, we were able to make some observations concerning the impact our threshold levels have on both of these performance measures.

Acknowledgments

Steve Drekic acknowledges the financial support from the Natural Sciences and Engineering Research Council of Canada through its Discovery Grants program (#238675-2010-RGPIN). The authors are also very grateful to the anonymous referees for their careful reading of the paper and their useful suggestions that have improved the overall presentation of the paper.

Author Contributions

Both authors contributed to all aspects of this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. B. De Finetti. “Su un’ impostazione alternativa dell teoria collettiva del rischio.” Trans. XVth Int. Congr. Actuar. 2 (1957): 433–443. [Google Scholar]
  2. H. Albrecher, and S. Thonhauser. “Optimality results for dividend problems in insurance.” RACSAM 103: 295–320. [CrossRef]
  3. B. Avanzi. “Strategies for dividend distribution: A review.” North Am. Actuar. J. 13 (2009): 217–251. [Google Scholar] [CrossRef]
  4. S. Drekic, and A.M. Mera. “Ruin analysis of a threshold strategy in a discrete-time Sparre Andersen model.” Methodol. Comput. Appl. Probab. 13 (2011): 723–747. [Google Scholar] [CrossRef]
  5. A.S. Alfa, and S. Drekic. “Algorithmic analysis of the Sparre Andersen model in discrete time.” ASTIN Bull. 37 (2007): 293–317. [Google Scholar] [CrossRef]
  6. J. Cai, and D.C.M. Dickson. “Ruin probabilities with a Markov chain interest model.” Insur. Math. Econ. 35 (2004): 513–525. [Google Scholar] [CrossRef]
  7. W. Li. “Ruin probability of the renewal model with risky investment and large claims.” Sci. China Ser. A Math. 52 (2009): 1539–1545. [Google Scholar]
  8. R. Korn, and A. Wiese. “Optimal investment and bounded ruin probability: Constant portfolio strategies and mean-variance analysis.” ASTIN Bull. 38 (2008): 423–440. [Google Scholar] [CrossRef]
  9. E.C.K. Cheung, D. Landriault, G.E. Willmot, and J.-K. Woo. “Structural properties of Gerber-Shiu functions in dependent Sparre Andersen models.” Insur. Math. Econ. 46 (2010): 117–126. [Google Scholar] [CrossRef]
  10. H. Cossette, D. Landriault, and E. Marceau. “Ruin probabilities in the discrete time renewal risk model.” Insur. Mathe. Econ. 38 (2006): 309–323. [Google Scholar] [CrossRef]
  11. D.C.M. Dickson, and H.R. Waters. “Some optimal dividends problems.” ASTIN Bull. 34 (2004): 49–74. [Google Scholar] [CrossRef]
  12. D.C.M. Dickson. Insurance Risk and Ruin. Cambridge, UK: Cambridge University Press, 2005. [Google Scholar]
  13. D.C.M. Dickson, A.D. Egidio dos Reis, and H.R. Waters. “Some stable algorithms in ruin theory and their applications.” ASTIN Bull. 25 (1995): 153–175. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Kim, S.S.; Drekic, S. Ruin Analysis of a Discrete-Time Dependent Sparre Andersen Model with External Financial Activities and Randomized Dividends. Risks 2016, 4, 2. https://doi.org/10.3390/risks4010002

AMA Style

Kim SS, Drekic S. Ruin Analysis of a Discrete-Time Dependent Sparre Andersen Model with External Financial Activities and Randomized Dividends. Risks. 2016; 4(1):2. https://doi.org/10.3390/risks4010002

Chicago/Turabian Style

Kim, Sung Soo, and Steve Drekic. 2016. "Ruin Analysis of a Discrete-Time Dependent Sparre Andersen Model with External Financial Activities and Randomized Dividends" Risks 4, no. 1: 2. https://doi.org/10.3390/risks4010002

APA Style

Kim, S. S., & Drekic, S. (2016). Ruin Analysis of a Discrete-Time Dependent Sparre Andersen Model with External Financial Activities and Randomized Dividends. Risks, 4(1), 2. https://doi.org/10.3390/risks4010002

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