In the renewal risk theory, the study of two sided jumps has been attracted by many researchers since its introduction. After the development of the distribution of modified inter time claim occurrence, the explicit expressions for ruin theory components in the literature under some assumptions, in this work, we examine probability density of the time of ruin, surplus immediately before ruin and deficit at ruin respectively under two sided risk process using some basic assumptions. Explicit expressions for distribution of interest are being derived.
The ruin problems and its analysis have been studied extensively by many authors ever since its inception. In the second half of 1900’s Dickson et.al, Gerber et al. etc. have made a remarkable contribution to risk process both renewal and classical as well. Recently more works have been carried down in somewhat new approach in risk renewal process, so-called two sided jump process. Gerber and Shiu [
The success of an insurance company is not only due to earnings in its conventional business but also due to intelligent investments of the money at its hand. Eventually risk invest can breed danger at the period when the market value of the assets is low and the company will not be able to recover from the loss because of market fluctuations. Therefore it is important to have a look on investments upon risk or risk free strategies. The investment of the money and settling down of the claim amount lead to a new branch of risk process called two sided risk process. More recently risk theory researches in two sided jumps field have got a remarkable attention. Substantial amount of works has been devoted to find the various ruin probability components under two sided risk model. For research on this kind of model, we refer to Perry, Stadje and Zacks [
The current paper aims at studying risk model with two sided jumps under some valid assumptions. In this paper we aim to derive a simple and unified expression for the densities of the time of ruin, the surplus immediately before ruin and deficit at ruin under two sided renewal risk process. The rest of the paper is organized as follows; In Section 2 we shall introduce some preliminary knowledge and focus on the main problem. In Section 3 we derive the densities under the condition suggested followed by the concluding remark.
Stochastic process with two sided jumps has captured a lot of attention in recent times. Various first crossing times, two sided passage problems, ascending/ descending and some other related quantities have been studied under different model. The classical or Sparre Anderson renewal model consider only one sided rather the claim amount along with number of claims under some basic assumptions. The premium amount assumed to be constant and the respective random variables such as number of claims, claim size are also assumed to be independent to each other. In two sided jump problems we let the surplus amount at any point of time to be a function of more random variables. In other words we deposit or invest the surplus amount in some other profit making portfolio which breeds one more random variable preferably in positive direction. Apart from negatively moving claim size variable we thus having a one more positively moving random gain variable. This leads to a new kind of risk model termed as risk model with two sided jumps. This model can be extended to risk model with multi dimensional moves by incorporating more significant variables.
Following the idea of Dong and Liu [
R ( t ) = u + p t + ∑ i = 1 M ( t ) X i − ∑ j = 1 N ( t ) Y j , t ≥ 0 (1)
where u ≥ 0,the initial surplus, p > 0, the constant premium rate, ∑ i = 1 M ( t ) X i Com-
pound Poisson process with intensity λ representing the total random gain (premium income or investment or annuity) until time t. Xi’s are independent and identically distributed random variables with common density ψ and mean μx. Here we assume that X follows an Erlangian (2, β) and Laplace trans-
forms of ψ be ψ ⌢ ( s ) = ∫ 0 ∞ e − s x ψ ( x ) d x . Therefore
ψ ( x ) = β 2 x e − β x , x > 0 , β > 0
and ψ ⌢ ( s ) = β 2 ( s + β ) 2 , whereas ∑ j = 1 N ( t ) Y j is accumulated claim process and {N(t)} is a Renewal process representing the number of claims up to time t with inter claim {Hi}. The Hi’s are i.i.d random variables with common density ξ and Laplace transform ξ ⌢ ( s ) = ∫ 0 ∞ e − s h ξ ( h ) d h . The claim sizes {Yj} are positive random variables with a common distribution function Q, density q and mean μy. Laplace transform q ^ ( s ) = ∫ 0 ∞ e − s y q ( y ) d y . Here we assume that Y follows Lindley distribution with parameter θ. Then, q ( y ) = θ 2 1 + θ ( 1 + y ) e − θ y , θ > 0 , [
Also we assume that { X i } , { Y j } , { M ( t ) } and { N ( t ) } are mutually independent. For ensuring a positive security loading condition, we take ( p + λ μ x ) ⋅ E ( H ) > μ y .
Let T = inf { t ≥ 0 , R ( t ) < 0 } be the ruin time, R(T_) be the surplus immediately before ruin, namely undershoot and | R ( T ) | be the deficit at ruin, termed as overshoot, again we define the probability of ruin
ψ ( u ) = P [ R ( T ) < 0 / R ( 0 ) = u ] , u ≥ 0.
For δ > 0 , due to Gerber and Shiu [
Φ ( u ) = E [ e − ( δ T ) w ( R ( T − ) , | R ( T ) | ) I ( T < ∞ ) / R ( 0 ) = u ] , u ≥ 0
where w ( x , y ) is a non-negative measurable function on [ 0 , ∞ ) × [ 0 , ∞ ) and I ( . ) is an indicator function. When w ( x , y ) = 1 and δ = 0 , Φ ( u ) becomes the probability of ruin,
ψ ( u ) = E [ I ( T < ∞ ) / R ( 0 ) = u ] = P [ T < ∞ / R ( 0 ) = u ] , u ≥ 0 .
Again, Let T n = ∑ i = 1 n H i be the time when nth claim occurs, T0 = 0. Since ruin only occurs at the epochs where claims occur, then we define the discrete time process R ˜ = { R ˜ n , n = 0 , 1 , 2 , ⋯ } and R ˜ 0 = 0
Again R ˜ n = R ( T n ) denotes the surplus immediately after the nth claim.
Now,
R ˜ n = u + p T n + ∑ i = 1 M ( T n ) X i − ∑ k = 1 n Y k = u + p T ˜ n − ∑ k = 1 n Y k
where T ˜ n = T n + ∑ i = 1 M ( T n ) X i / p with T ˜ 0 = 0 which corresponds to the Sparre Anderson risk model [
R ¯ ( t ) = u + p t − ∑ i = 1 N ¯ ( t ) Y i
where the initial surplus u and the claim size Yi are exactly the same as those in model (1). The counting number process N ¯ ( t ) denotes the number of claims up to time t with the modified inter claim times Z i = T ˜ i − T ˜ i − 1 . Clearly Zi are i.i.d random variables with a common density k. (say)
In this model we assume the initial inter claim time, H ~ exp(α), the number of claims M(T) follows Poisson(λ), the random gain X ~ Erlang(2, β) and finally we propose the claim amount Y follows Lindley distribution with parameter θ.
For u ≥ 0, let f ( x , y , t / u ) represents joint probability density function of R(T_), |R(T)| and T, then
∫ 0 ∞ ∫ 0 ∞ ∫ 0 ∞ f ( x , y , t / u ) d x d y d t = P ( T < ∞ / R ( 0 ) = u ) = ψ ( u )
taking T as time of ruin.
We let that x > u + p t , f ( x , y , t / u ) = 0 , or x ≤ u + p t , f ( x , y , t / u ) = 1 and
f ( u + p t , y , t / u ) d x d y d t = k ( t ) ⋅ 1 ⋅ f ( u + p t + y ) d y d t
where k(t) is the linear combination of exponential distribution (Rebello & Thampi, [
k ( t ) = a 1 e − R 1 t + a 2 e − R 2 t + a 3 e − R 3 t , t > 0
a1, a2, a3, R1, R2 and R3 are suitably chosen constants.
It is known, using P ( X ∪ Y ) = P ( X ) ⋅ P ( Y / X ) , that f ( x , y , t / u ) is the jpdf of R(T_) and T at the point (x, t) multiplied by the conditional probability density function of |R(T)| at y given that R(T_) = x and T = t.
Since the conditional probability |R(T_)| does not depend on t we have
γ ( x + y ) ∫ 0 ∞ γ ( x + y ) d x = γ ( x + y ) ϒ ¯ ( x ) (Gerber and Shiu, [
hence
f ( x , y , t / u ) = P ( R ( T _ ) ∩ T ) ⋅ P ( | R ( T ) | / R ( T _ ) ∩ T )
⇒ f ( x , y , t / u ) = ∫ 0 ∞ f ( x , z , t / u ) d z γ ( x + y ) ϒ ¯ ( x ) (2)
Define
p ( x , y / u ) = ∫ 0 ∞ e − δ t f ( x , y , t / u ) d t
where δ can be considered as a force of interest, as a dummy variable, and χ ( x / u ) = ∫ 0 ∞ p ( x , y / u ) d y . multiplying (2) by e − δ t and integrating with respect to t. i.e.,
∫ 0 ∞ e − δ t f ( x , y , t / u ) d t = ∫ 0 ∞ e − δ t ∫ o ∞ f ( x , z , t / u ) d z ⋅ γ ( x + y ) ϒ ¯ ( x ) d t
p ( x , y / u ) = χ ( x / u ) ⋅ γ ( x + y ) ϒ ¯ ( x )
Hence g ( y ) = ∫ 0 ∞ p ( x , y / u ) d x i.e., g ( y ) = ∫ 0 ∞ ∫ 0 ∞ e − δ t f ( x , y , t / 0 ) d t d x .
Using the definition of martingale; upon the stochastic process { R ( t ) } t ≥ 0 with stationary and independent increments, the process { e − δ t + ξ R ( t ) } t ≥ 0 is a martingale iff for each t > 0 its expectation at time t is equal to its initial value ( where ξ is a constant).
i.e., iff E ( e − δ t + ξ R ( t ) / R ( 0 ) = u ) = E ( e ξ u ) = e ξ u (3)
Define χ ⌢ be the Laplace transform of a function χ ( x ) , x > 0 ,
χ ⌢ ( ξ ) = ∫ 0 ∞ e − ξ x χ ( x ) d x
Then LHS of (3)
R ( e − δ t + ξ R ( t ) / R ( 0 ) = u ) = E ( e − δ t + ξ { U + p t − ∑ i = 1 N ¯ ( t ) X i } ) = E ( e − δ t ⋅ e ξ u ⋅ e ξ p t ⋅ e − ∑ i = 1 N ¯ ( t ) X i ) = exp { − δ t + ξ u + ξ p t + λ t ( χ ⌢ ( ξ ) − 1 ) }
(3) Implies − δ + ξ p + λ ( χ ⌢ ( ξ ) − 1 ) = 0 .
It tends to ∞ as ξ → ∞ and takes the value − δ at ξ = 0 . thus it has a unique non-negative root, ρ hence { e − δ t + ρ R ( t ) } t ≥ 0 is a martingale.
For x > u = R ( 0 ) , let T x = inf { t / R ( t ) = x } be the first time when the surplus reaches the level x.
Then R ( t ) ≤ x , 0 ≤ t ≤ T x and using optimal stopping theorem to { e − δ t + ρ R ( t ) } t ≥ 0 , we have
e ρ u = E ( e − δ T x + ρ R ( T x ) / R ( 0 ) = u ) = e ρ x ⋅ E ( e − δ T x / R ( 0 ) = u )
Or
E ( e − δ T x / R ( 0 ) = u ) = e − ρ ( x − u ) (4)
Using the idea of Gerber and Shiu [
For x > u, let π ( t ; u , x ) ; t > 0 denote the pdf of the random variable T.
Hence (4) implies ∫ 0 ∞ e − δ t π ( t ; u , x ) d t = e − ρ ( x − u ) . The differential π ( t ; u , x ) d t is the probability that the surplus process meets the level x first time between (t, t + dt) provides the surplus cannot reach x before t = x − u p . Hence for t < x − u p ; π ( t ; u , x ) = 0 and the density function π ( t ; u , x ) at t = x − u p equals a 1 R 1 ⋅ ( 1 − e − R 1 ( x − u ) / p ) + a 2 R 2 ⋅ ( 1 − e − R 2 ( x − u ) / p ) + a 3 R 3 ⋅ ( 1 − e − R 3 ( x − u ) / p ) or = ∑ i = 1 3 a i R i ⋅ ( 1 − e − R i ( x − u ) / p ) .
Using duality π ( t ; u , x ) = π ˜ ( t ; u , x ) , Now the differential f ( x , y , t / u ) d x d y d t can be interpreted as the probability of the event that ruin does not take place by the time t, that the surplus process meets the level x between (t, t + dt), but does not attain the level x + dx owing a claim within dx/p time units after Tx and that the size of this claim is x + y.
∴ f ( x , y , t / u ) d x d y d t = π ˜ ( t ; u , x ) d t ⋅ λ d x p ⋅ γ ( x + y ) d y
Or
f ( x , y , t / u ) = λ p ⋅ ∑ i = 1 3 a i R i ⋅ ( 1 − e − R i ( x − u ) / p ) ⋅ θ 2 1 + θ ( 1 + x + y ) e − θ ( x + y ) (5)
using (4), if we multiply (5) by e − δ t d t and integrate from t = 0 to t = ∞ and apply(4), we have
p ( x , y / u ) = λ p ⋅ e − ρ x ⋅ θ 2 1 + θ ( 1 + x + y ) e − θ ( x + y ) (6)
Or
χ ( x / 0 ) = λ p ⋅ e − θ x ( 1 + θ x θ + 1 ) ⋅ e − ρ x
and
g ( y ) = λ p ⋅ ∫ 0 ∞ θ 2 1 + θ e − θ ( x + y ) ( 1 + x + y ) ⋅ e − ρ x d x
CorollaryBy distinguishing whether or not ruin occurs at the first time when the surplus falls below the initial value u, we have
p ( x , y / u ) = ∫ 0 u p ( x , y / u − z ) ⋅ g ( z ) d z + p ( x − u , y + u / 0 ) , 0 < u < x (7)
Using (7)
p ( x , y / u ) = λ p e − ρ x γ ( x + y )
p ( x − u , y + u / 0 ) = λ p e − ρ ( x − u ) γ ( x + y ) = p ( x , y / 0 ) ⋅ e ρ u
∴ ( 7 ) ⇒ p ( x , y / u ) = ∫ 0 u p ( x , y / u − z ) ⋅ g ( z ) d z + p ( x , y / 0 ) ⋅ e ρ u ⋅ I ( u < x ) (8)
where I ( u < x ) is an indicator function. Integrating (8) with respect to y, we have χ ( x / u ) = ∫ 0 u χ ( x / u − z ) ⋅ g ( z ) d z + χ ( x / 0 ) ⋅ e ρ u ⋅ I ( u < x ) , defining
ς ( u ) = ∫ 0 u ς ( u − z ) g ( z ) d z + e ρ u I ( u < x )
then
χ ( x / u ) = χ ( x / 0 ) ς ( u )
This paper is prepared to learn the insurance process under two sided jump risk renewal processes. We investigated a reformulated insurance process in which random gain and claim amount are all uncertain renewal processes. This paper derives the distribution of time of ruin, deficit at ruin, surplus amount t before ruin under two sided risk theory setup. The explicit expressions for the same are being derived. The application of other feasible distributions may be considered as a scope of further study.
The authors would like to thank the anonymous referee for helpful comments and suggestions which improved the earlier version of the paper. This research was supported by University Grants Commission of India.
Rebello, J.J. and Thampi, K.K. (2017) The Distribution of the Time of Ruin, the Surplus Immediately before Ruin and Deficit at Ruin under Two Sided Risk Renewal Process. Journal of Ma- thematical Finance, 7, 624-632. https://doi.org/10.4236/jmf.2017.73032