Moments of the Discounted Aggregate Claims with Delay Inter-Occurrence Distribution and Dependence Introduced by a FGM Copula

Franck Adékambi

doi:10.5772/intechopen.88699

Abstract

In this chapter, with renewal argument, we derive higher simple moments of the Discounted Compound Delay Renewal Risk Process (DCDRRP) when introducing dependence between the inter-occurrence time and the subsequent claim size. To illustrate our results, we assume that the inter-occurrence time is following a delay-Poisson process and the claim amounts is following a mixture of Exponential distribution, we then provide numerical results for the first two moments. The dependence structure between the inter-occurrence time and the subsequent claim size is defined by a Farlie-Gumbel-Morgenstern copula. Assuming that the claim distribution has finite moments, we obtain a general formula for all the moments of the DCDRRP process.

Keywords

compound delay-Poisson process
discounted aggregate claims
moments
FGM copula
constant interest rate

chapter and author info

Author

Franck Adékambi*
- School of Economics, University of Johannesburg, South Africa

*Address all correspondence to: fadekambi@uj.ac.za

DOI: 10.5772/intechopen.88699

From the Edited Volume

Probability, Combinatorics and ControlEdited by Andrey Kostogryzov

Probability, Combinatorics and Control

Edited by Andrey Kostogryzov and Victor Korolev

Show +

1. Introduction

The classical Poisson model is attractive in the sense that the memoryless property of the exponential distribution makes calculations easy. Then the research was extended to ordinary Sparre-Andersen renewal risk models where the inter-claim times have other distributions than the exponential distribution. Dickson and Hipp [1, 2] considered the Erlang-2 distribution, Li and Garrido [3] the Erlang-n distribution, Gerber and Shiu [4] the generalized Erlang-n distribution (a sum of n independent exponential distributions with different scale parameters) and Li and Garrido [5] looked into the Coxian class distributions. One difficulty with these models is that we have to assume that a claim occurs at time 0, which is not the case in usual setting.

Albrecher and Teugels [6] considered modeling dependence with the use of an arbitrary copula. In a similar dependence model to Albrecher and Teugels as well, Asimit and Badescu [7] considered a constant force of interest and heavy tailed claim amounts.

Barges et al. [8] followed the idea of Albrecher and Teugels [6] and supposed that the dependence is introduced by a copula, the Farlie-Gumbel-Morgenstern (GGM) copula, between a claim inter-arrival time and its subsequent claim amount.

Adékambi and Dziwa [9] and Adékambi [10] provide a direct point of extension but assuming that the claim counting process to follow an unknown general distribution in a framework of dependence with random force of interest to calculate the first two moments of the present value of aggregate random cash flows or random dividends.

The discounted aggregate sum has also been applied in many other fields. For example, it can be used in health cost modeling, see Govorun and Latouche [11], Adékambi [12], or in reliability, in civil engineering, see Van Noortwijk and Frangopol [13].

The delayed or modified renewal risk model solves this problem by assuming that the time until the first claim has a different distribution than the rest of the inter-claim times. Not much research has been done for this model at this stage. Among the first works was Willmot [14] where a mixture of a “generalized equilibrium” distribution and an exponential distribution is considered for the distribution of the time until the first claim. Special cases of the model include the stationary renewal risk model and the delayed renewal risk model with the time until the first claim exponentially distributed. Our focus is to extend the work of Bargès et al. [8], Adékambi and Dziwa [9] and Adékambi [10] by allowing the counting process to follow a delay renewal risk process and thus derive a recursive formula of the moments of this subsequent Discounted Compound Delay Poisson Risk Value (DCDPRV).

For example, young performer companies typically have a high growth rate at the beginning, but as they mature their growth rate may decrease with the increasing scarcity of investment opportunities. That makes dividends dependent on the economic climate at the dividend occurrence time. Obviously the distribution of inter-dividends time in times of economic expansion and in times of economic crisis cannot be identically distributed. So it would be appropriate to use a delayed renewal model to model the distribution of the inter-dividend time. A delayed renewal process is just like an ordinary renewal process, except that the first arrival time is allowed to have a different distribution than the other inter-dividends times.

The chapter is organized as follows: In the second section, we present the model of the continuous time discounted compound delay-Poisson risk process that we use and give some notation. In Section 3, we present a general formula for all the moments of the DCDPRV process. A numerical example of the first two moments will then follow in Section 4.

2. The model

We use the same model as the one in Bargès et al. [8], where the instantaneous interest rate δis constant.

Define our risk model as follows:

The number of claims Ntt≥0and Ndtt≥0form, respectively, an ordinary and a delayed renewal process and, for k∈N=123…:

the positive claim occurrence times are given by Tk,
the positive claim inter-arrival times are given by τk=Tk−Tk−1,k∈N, and T0=0.
τkk≥2∼τ2are independent and identically distributed (i.i.d),

The kthrandom claim is given by Xk, and

Xkk∈Nare independent and identically distributed (i.i.d),
Xkτkk∈Nare mutually independent; and the higher moments, μk=EX1kof X1exist.

The discounted aggregate value at time t=0of the claims recorded over the period 0tyields, respectively, for the ordinary and the delayed renewal case:

Z0t=∑k=1Note−δTkXk,Zdt=∑k=1Ndte−δTkXk,E1

where Z0t=Zdt=0if N0t=Ndt=0.

2.1 The dependence

We introduce a specific structure of dependence based on the Farlie-Gumbel-Morgenstern (FGM) copula. The advantage of using the FGM copula and its generalizations lies in its mathematical manageability. The joint cumulative distribution function (c.d.f.) of Xiτi, the ith claim and its occurrence time is

FXi,τixv=CFXixFτiv=FXixFτiv+θFXixFτiv1‐FXix1‐Fτiv,E2

for xv∈R+∗R+and where FXixand Fτivare the marginals of Xiand τirespectively. Recall that the density of the FGM copula is

cθFGMuv=1+θ1‐2u1‐2v,E3

for uv∈01∗01so that the joint probability density function (p.d.f.) of Xiτiis

fXi,Tixv=cθFGMFXixFTivfXixfTiv=fXixfTiv+θfXixfTiv1‐2FXix1‐2FTiv,E4

where fXiand fτiare the p.d.f.’s of Xiand τirespectively.

With these hypotheses, we present in Section 3 recursive formula of the higher moments of this present value risk process, for a constant instantaneous interest rate.

3. Recursive expression for higher moments

It is often easier to calculate the moments of the random variable Zdtt≥0than finding its distribution. If the probability generation function of Zdtt≥0or its moment generating function (mgf) exists, it is possible to obtain the corresponding distribution by inversion of its mgf. Since, there is relatively little research devoted to the study of the distribution of the discounted compound renewal sums. We could then think about another technique other than the one proposed by the above authors by studying the moments of Zdtt≥0.

3.1 Delay renewal case

The mathematical expectation of total claims plays an important role in the determination of the pure premium, in addition to giving a measure of the central tendency of its distribution. The moments centered at the average of order 2, 3 and 4 are the other moments usually considered because they usually give a good indication of the pace of distribution, and these give us respectively a measure of the dispersion of the distribution around its mean, a measure of the asymmetry and flattening of the distribution considered.

Moments, whether simple, joined or conditional, may eventually be used to construct approximations of the distribution of the DCDPRV.

Theorem 3.1

The Laplace transform of the mth moment of Zdtt≥0is given by:

π˜Zdmr=1+λ2r+mδ+λ1−λ2r+mδ+λ1u˜mr=λ11+λ2r+mδ+λ1−λ2r+mδ+λ1×∑j=0m−1mjμm−j−θμm−j′−μm−jλ1+mδ+r+2θμm−j′−μm−j2λ1+mδ+rπ˜ZojrE5

where

π˜Zdmr=u˜mr+λ2mδu˜mr×Lτ1mδr+λ1−λ2mδ+λ1u˜m×Lτ1mδ+λ1r.E6

Proof

Conditioning on the arrival of the first claim leads to

πZdmt=EZmt=EEe−δsX1+e−δsZot−smτ1=s=∑j=0m−1mj∫0tfτ1se−mδsEXm−jτ1=sπZojt−sds+∫0tfτ1se−mδsπZomt−sds.E7

We have

EXm−jτ1=s=∫0∞xm−jfXτ1=sxdx=∫0∞xm−j1+θ1−2FXx1−2Fτ1sfXxdx=EXm−j+θ∫0∞xm−j1−2FXx1−2Fτ1sfXxdx=EXm−j+θ∫0∞xm−j2−2FXx1−2Fτ1sfXxdx−θ∫0∞xm−j1−2Fτ1sfXxdx=EXm−j1−θ1−2Fτ1s+θ1−2Fτ1s∫0∞m−jxm−j1−Fτ1sdx.E8

We let,

μm−j′=EX′m−j=∫0∞m−jxm−j−11−FXx2dx<∫0∞m−jxm−j−11−FXxdx=EXm−j<∞E9

such that the above equation becomes

EXm−jτ1=s=μm−j+θ1−2Fτ1sμm−j′−μm−j.E10

πZdmt=EZmt=EEe−δsX1+e−δsZot−smτ1=s=∑j=0m−1mj∫0tfτ1se−mδsμm−j+θ1−2Fτ1sμm−j′−μm−jπZojt−sds+∫0tfτ1se−mδsπZomt−sds.

Let us ∫0tfτ1se−mδsds=Hδt, ∫0tfτ2se−mδsds=Iδtthen

πZdmt=∑j=0m−1mj∫0tfτ1se−mδsμm−j+θ1−2Fτ1sμm−j′−μm−jπZojt−sds+Hmδ∗πZo.m=um+Hmδ∗um+Imδ∗πZo.m=um+Hmδ∗um+Hmδ∗Imδ∗πZo.m=um+Hmδ∗um+Hmδ∗Imδ∗um+Imδ∗πZo.m=um+Hmδ∗um+um∗∑k=1∞Hmδ∗Imδ∗kt=um+um∗∑k=0∞Hmδ∗Imδ∗kt=um+∫0tumt−se−mδsdmds,E11

where umt=∑j=0m−1mj∫0tfτ1se−mδsμj+θ1−2Fτ1sμj′−μjπZom−jt−sds.

We consider the case where the canonical random variable τ2has an Exponential distribution with parameter λ2>0and τ1has an Exponential distribution with parameter λ1>0.

That is, we have:

fτ1t=λ1e−λ1t,fτ2t=λ2e−λ2t,Lτ1λ1s=∫0∞e−svfτ1vdv=λ1λ1+s,Lτ2λ2s=λ2λ2+s.

mdt=λ2t+λ1−λ2λ11−eλ1tE12

The mth moment of Zdtis then given by,

πZdmt=um+∫0tumt−se−mδsdmds=um+λ2∫0tumt−se−mδsds+λ1−λ2∫0tumt−se−mδ+λ1sds=um+λ2mδ∫0tumt−smδe−mδsds+λ1−λ2mδ+λ1∫0tumt−smδ+λ1e−mδ+λ1sdsE13

Taking the Laplace transform of the above equation, we get:

π˜Zdmr=u˜mr+λ2mδu˜mr×Lτ1mδr+λ1−λ2mδ+λ1u˜m×Lτ1mδ+λ1rE14

But,

umt=∑j=0m−1mj∫0tfτ1se−mδsμj+θ1−2Fτ1sμj′−μjπZom−jt−sds=∑j=0m−1mj∫0tλ1e−λ1se−mδsμj+θ2e−λ1s−1μj′−μjπZom−jt−sds=λ1μj−θμj′−μjλ1+mδ∑j=0m−1mj∫0tλ1+mδe−λ1+mδsπZom−jt−sds+2θλ1μj′−μj2λ1+mδ∑j=0m−1mj∫0t2λ1+mδe−2λ1+mδsπZom−jt−sdsE15

Then the Laplace transform of umt, at r, will give:

u˜mr=λ1∑j=0m−1mjμj−θμj′−μjλ1+mδ+r+2θμj′−μj2λ1+mδ+rπ˜Zom−jrE16

Substituting Eq. (14) into Eq. (13), we have:

π˜Zdmr=1+λ2r+mδ+λ1−λ2r+mδ+λ1u˜mr=λ11+λ2r+mδ+λ1−λ2r+mδ+λ1∑j=0m−1mjμj−θμj′−μjλ1+mδ+r+2θμj′−μj2λ1+mδ+rπ˜Zom−jrE17

Solving the above equation for the ordinary case, where τ2k≥2∼τ2, we have:

π˜Zomr=λ2μmrr+δm+λ2+λ2r+δm+λ2∑k=1m−1Cmkμkπ˜Zom−kr+θμm′−μmλ2r+δmrr+δm+λ2r+δm+2λ2+θλ2r+δmr+δm+λ2r+δm+2λ2∑k=1m−1Cmkμk′−μkπ˜Zom−kr+λ2r+δm+λ2π˜Z0mrE18

Rearranging the above equation, we will get

π˜Zomr=λ2μmrr+δm+λ2r+δm∑k=1m−1Cmkμkπ˜Zom−kr+θλ2μm′−μmrr+δm+2λ2+θλ2r+δm+2λ2∑k=1m−1Cmkμk′−μkπ˜Zom−krE19

Corollary 3.1

The first moment of Zdtt≥0is given by:

πZdt=θλ1μ1′−μ1λ2+δδ+λ1δ+2λ2+λ1λ2+δδδ+λ1μ1+θλ1μ1′−μ1λ1−λ2λ1−2λ2−λ1−λ2μ11δ+λ1e−δ+λ1t−θλ11δ+2λ2λ2λ1−2λ2μ1′−μ1e−δ+2λ2t−2θλ11δ+2λ1μ1′−μ1e−δ+2λ1t−λ2δμ1e−δtE20

Proof:

From Theorem 3.1, we have:

π˜Zdr=λ1μ1rr+δ+λ1+λ1r+δ+λ1π˜Z0r+θμ1′−μ1λ1r+δrr+δ+λ1r+δ+2λ1E21

From Bargès et al. [8], we have

π˜Zor=λ2μ1rr+δ+θλ2μ1′−μ1rr+δ+2λ2E22

Substituting Eq. (22) into Eq. (21), yields

π˜Zdr=λ1μ1rr+δ+λ1+λ1r+δ+λ1λ2μ1rr+δ+θλ2μ1′−μ1rr+δ+2λ2+θλ1μ1′−μ1r+δrr+δ+λ1r+δ+2λ1=λ1λ2rr+δr+δ+λ1+λ1rr+δ+λ1μ1+θλ1μ1′−μ1λ2rr+δ+λ1r+δ+2λ2+r+δrr+δ+λ1r+δ+2λ1E23

with

λ1rr+δ+λ1=λ1δ+λ1.1r−λ1δ+λ1.1r+δ+λ1E24

λ1λ2rr+δr+δ+λ1=λ1λ2δδ+λ1.1r+λ2δ+λ1.1r+δ+λ1−λ2δ.1r+δE25

λ2rr+δ+λ1r+δ+2λ2=λ2δ+λ1δ+2λ21r+λ2λ1−2λ21δ+λ11r+δ+λ1−1δ+2λ21r+δ+2λ2E26

r+δrr+δ+λ1r+δ+2λ1=δδ+λ1δ+2λ11r+1δ+λ11r+δ+λ1−2δ+2λ11r+δ+2λ1E27

Substituting Eqs. (24), (25), (26) and (27) into Eq. (23), yields:

π˜Zdr=μ1λ1δ+λ1.1r−λ1δ+λ1.1r+δ+λ1+θμ1′−μ1λ2δ+λ1δ+2λ21r+λ2λ1−2λ21δ+λ11r+δ+λ1−1δ+2λ21r+δ+2λ2+θλ1μ1′−μ1δδ+λ1δ+2λ11r+1δ+λ11r+δ+λ1−2δ+2λ11r+δ+2λ1+μ1λ1λ2δδ+λ1.1r+λ2δ+λ1.1r+δ+λ1−λ2δ.1r+δE28

Rearranging the above equation, will give

πZdt=θλ1μ1′−μ1λ2+δδ+λ1δ+2λ2+λ1λ2+δδδ+λ1μ11r+θλ1μ1′−μ1λ1−λ2λ1−2λ2−λ1−λ2μ11δ+λ11r+δ+λ1−θλ11δ+2λ2λ2λ1−2λ2μ1′−μ11r+δ+2λ2−2θλ11δ+2λ1μ1′−μ11r+δ+2λ1−λ2δμ11r+δE29

Remark 1

If λ1=λ2then Eq. (29) becomes

π˜Zdr=θλμ1′−μ11δ+2λ+λδμ11r−1δ+2λθλμ1′−μ11r+δ+2λ−λδμ11r+δ=λδμ11r−1r+δ+θλμ1′−μ11δ+2λ1r−1r+δ+2λ=λμ1rr+δ+λθμ1′−μ1rr+δ+2λ,E30

which is exactly the result of Bargès et al. [8].

The inverse of the Laplace transform in Eq. (29) will give

πZdt=θλ1μ1′−μ1λ2+δδ+λ1δ+2λ2+λ1λ2+λ1δδδ+λ1μ1+θλ1μ1′−μ1λ2λ1−2λ2+1+λ2−λ1μ11δ+λ1e−δ+λ1t−θλ11δ+2λ2λ2λ1−2λ2μ1′−μ1e−δ+2λ2t−2θλ11δ+2λ1μ1′−μ1e−δ+2λ1t−λ2δμ1e−δtE31

Remarks 2

If θ=0and λ1≠λ2then

πZdt=λ1δλ2+δδ+λ1−λ2δe−δtμ1+λ2−λ1δ+λ1μ1e−δ+λ1t=λ21−e−δtδμ1+λ1−λ21−e−δ+λ1tδ+λ1μ1=λ2a¯tδ+λ1−λ2a¯tδ+λ1μ1E32

which is exactly the result of Léveillé et al. [15].

If λ1=λ2and θ≠0then

πZot=θλμ1′−μ11δ+2λ2+λδμ1−θλ1δ+2λμ1′−μ1e−δ+2λt−λδμ1e−δt=λδ1−e−δtμ1+θλ1−e−δ+2λtδ+2λμ1′−μ1E33

which is exactly the result of Bargès et al. [8].

If λ1=λ2and θ=0then

πZot=λδ1−e−δtμ1=λa¯tδμ1,E34

which is exactly the result of Léveillé et al. [15].

Corollary 3.2

The second moment of Zdtt≥0is given by the following development:

The result in Theorem 3.1 when n=2gives:

π˜Zd2r=2λ1λ1−λ2rr+2δ+λ1r+2δ+2λ11rμ2+2μ1π˜Zor+λ1r+2δ+λ1r+2δr+2δ+2λ2λ2r+2δ+2λ1+1π˜Z02r,E35

From Bargès et al. [8], we have.

π˜Zor=λ2μ1rr+δ+θλ2μ1′−μ1rr+δ+2λ2E36

and

π˜Zo2r=λ2μ2rr+2δ+θλ2μ2′−μ2rr+2δ+2λ2+2λ22μ12rr+δr+2δ+2θλ22μ1μ1′−μ1rr+δ+2λ2r+2δ+2θλ22μ1μ1′−μ1rr+2δ+2λ2r+δ+2θ2λ22μ1′−μ12rr+δ+2λ2r+2δ+2λ2E37

Substituting Eqs. (39) and (40) into Eq. (38), yields:

π˜Zd2r=2λ1λ1−λ2rr+2δ+λ1r+2δ+2λ11rμ2+2μ1λ2μ1rr+δ+θλ2μ1′−μ1rr+δ+2λ2+λ1r+2δ+λ1r+2δr+2δ+2λ2λ2r+2δ+2λ1+1×λ2μ2rr+2δ+θλ2μ2′−μ2rr+2δ+2λ2+2λ22μ12rr+δr+2δ+2θλ22μ1μ1′−μ1rr+δ+2λ2r+2δ+2θλ22μ1μ1′−μ1rr+2δ+2λ2r+δ+2θ2λ22μ1′−μ12rr+δ+2λ2r+2δ+2λ2,E38

and rearranging Eq. (38), will give:

π˜Zd2r=λ1μ2rr+2δ+λ1+2λ1λ2μ12rr+δr+2δ+λ1+2θλ1λ2μ1μ1′−μ1rr+δ+2λ2r+2δ+λ1+θλ1μ2′−μ2r+2δrr+2δ+λ1r+2δ+2λ1+2θλ1λ2μ1μ1′−μ1r+2δrr+δr+2δ+λ1r+2δ+2λ1+2θ2λ1λ2μ1′−μ12r+2δrr+δ+2λ2r+2δ+λ1r+2δ+2λ1+λ1λ2μ21rr+2δr+2δ+λ1+θλ1λ2μ2′−μ21rr+2δ+2λ2r+2δ+λ1+2λ1λ22μ121rr+δr+2δr+2δ+λ1+2θλ1λ22μ1μ1′−μ11rr+δ+2λ2r+2δr+2δ+λ1+2θλ22λ1μ1μ1′−μ11rr+2δ+2λ2r+δr+2δ+λ1+2θ2λ22λ1μ1′−μ121rr+δ+2λ2r+2δ+2λ2r+2δ+λ1,E39

which can be simplified to

π˜Zd2r=γ0r+γ1r+δ+γ2r+2δ+γ3r+2δ+λ1+γ4r+δ+2λ2+γ5r+2δ+2λ1+γ6r+2δ+2λ2,E40

with,

γ0=λ1μ22δ+λ1+2λ1λ2μ12δ2δ+λ1+2θλ1λ2μ1μ1′−μ12δ+λ1δ+2λ2+θδλ1μ2′−μ22δ+λ1δ+λ1+2θλ1λ2μ1μ1′−μ12δ+λ1δ+λ1+2δθ2λ1λ2μ1′−μ122δ+λ1δ+λ1δ+2λ2+λ1λ2μ22δ2δ+λ1+θλ1λ2μ2′−μ22δ+λ22δ+λ1+λ1λ22μ12δ22δ+λ1+θλ1λ22μ1μ1′−μ1δδ+2λ22δ+λ1+θλ1λ22μ1μ1′−μ1δδ+λ22δ+λ1+θ2λ1λ22μ1′−μ122δ+2λ2δ+λ22δ+λ1E41

γ1=−2λ1λ2μ12δδ+λ1−2θλ1λ2μ1μ1′−μ1δ+λ1δ+2λ1−2λ1λ22μ12δ2δ+λ1−2θλ1λ22μ1μ1′−μ1δδ+2λ2δ+λ1E42

γ2=−λ2μ22δ+λ22μ12δ2+θλ22μ1μ1′−μ1δδ−2λ2E43

γ3=−λ1μ22δ+λ1+2λ1λ2μ122δ+λ1δ+λ1+2θλ1λ2μ1μ1′−μ12δ+λ1λ1+δ−2λ2+θλ1μ2′−μ22δ+λ1−2θλ1λ2μ1μ1′−μ12δ+λ1δ+λ1−2θ2λ1λ2μ1′−μ122δ+λ1λ1+δ−2λ2+λ2μ22δ+λ1+θλ1λ2μ2′−μ22δ+λ1λ1−2λ2−2λ22μ122δ+λ1δ+λ1−2θλ22μ1μ1′−μ12δ+λ1λ1+δ−2λ2+2θλ1λ22μ1μ1′−μ12λ2−λ1δ+λ12δ+λ1+2θ2λ1λ22μ1′−μ122δ+λ1λ1+δ−2λ22λ2−λ1E44

γ4=−2θλ1λ2μ1μ1′−μ1δ+2λ2λ1+δ−2λ2+2θ2λ1λ2μ1′−μ122λ2−δδ+2λ2λ1+δ−2λ22λ1−2λ2+δ−2θλ1λ22μ1μ1′−μ1δ+2λ2δ−2λ2λ1+δ−2λ2−2θ2λ1λ22μ1′−μ12δδ+2λ2λ1+δ−2λ2E45

γ5=−θλ1μ2′−μ2δ+λ1+2θλ1λ2μ1μ1′−μ1δ+λ1δ+2λ1−2θ2λ1λ2μ1′−μ12δ+λ12λ1−δ−2λ2E46

γ6=−θλ1λ2μ2′−μ22δ+λ2λ1−2λ2+θλ1λ22μ1μ1′−μ1δ+λ2δ+2λ2λ1−2λ2+θ2λ1λ22μ1′−μ12δδ+λ2λ1−2λ2E47

Remark 2

When

λ1=λ2

π˜Zd2r=2λ2μ2rr+δr+2δ+2θλ2μ1μ1′−μ1rr+δ+2λr+2δ+λμ2rr+2δ+θλμ2′−μ2rr+2δ+2λ2θλ2μ1μ1′−μ1rr+δr+2δ+2λ+2θ2λ2μ1′−μ12rr+δ+2λr+2δ+2λ,E48

which is exactly the result of Bargès et al. [8].

The Laplace transform in Eq. (49), is a combination of terms of the form:

g˜r=1rα1+rα2+r…αn+r,E49

with ga function defined for all non-negative real numbers. As described in the proof of Theorem 1.1 in Baeumer [16], each of these terms can be expressed as a combinations of partial fraction such as g˜r=γ01r+γ11α1+r+…+γn1αn+rwhere.

γ0=1α1…αn, for i=1,…,n, γi=−1αi∏j=1,j≠i1αj−αi.

Since the inverse Laplace transform of 1αi+ris e−αit, it is easy to invert g˜and obtain

gt=γ0+γ1e−α1t+γ2e−α2t+…+γne−αnt.E50

Using Eq. (49) in Eq. (53), it results that

πZd2t=γ0+γ1e−δt+γ2e−2δt+γ3e−2δ+λ1t+γ4e−δ+2λ2t+γ5e−2δ+2λ1t+γ6e−2δ+2λ2t,t≥0,E51

where γii∈012…6are given by equation Eq. (50).

Remarks

If θ=0then

γ0θ=0=λ1μ22δ+λ1+2λ1λ2μ12δ2δ+λ1+λ1λ2μ22δ2δ+λ1+λ1λ22μ12δ22δ+λ1=λ12δ+λ12δ+λ22δμ2+λ1λ2δ2δ+λ12δ+λ2δμ12=λ12δ+λ22δ2δ+λ1μ2+2λ2μ12δ,E52

γ1θ=0=−2λ1λ2μ12δδ+λ1−2λ1λ22μ12δ2δ+λ1=−2λ1λ2μ12δδ+λ11+λ2δ,γ2θ=0=λ22μ12δ2−λ2μ22δE53

γ3θ=0=−λ1μ22δ+λ1+2λ1λ2μ122δ+λ1δ+λ1+λ2μ22δ+λ1−2λ22μ122δ+λ1δ+λ1=λ2−λ1μ22δ+λ1+2λ2λ1−λ2μ122δ+λ1δ+λ1,E54

γ4θ=0=γ5θ=0=γ6θ=0.E55

Then,

πZd2t=λ12δ+λ22δ2δ+λ1μ2+2λ2μ12δ−2λ1λ2μ12δδ+λ11+λ2δe−δt+λ22μ12δ2−λ2μ22δe−2δt+12δ+λ1λ2−λ1μ2+2λ2λ1−λ2μ122δ+λ1δ+λ1e−2δ+λ1t=λ12δ+λ22δ2δ+λ1−λ22δe−2δt+λ2−λ12δ+λ1e−2δ+λ1tμ2+λ2λ12δ+λ2δ22δ+λ1−2λ1λ2δδ+λ11+λ2δe−δt+λ22δ2e−2δt+2λ2λ1−λ22δ+λ1δ+λ1e−2δ+λ1tμ12=λ12δ+λ22δ2δ+λ1−λ22δ1−2δa¯t2δ+λ2−λ12δ+λ11−2δ+λ1a¯t2δ+λ1μ2+λ2λ12δ+λ2δ22δ+λ1−2λ1λ2δδ+λ11+λ2δe−δt+2λ22δ2e−2δt−λ22δ2e−2δtμ12+2λ2λ1−λ2δδ+λ1−2λ2λ1−λ2δ2δ+λ1μ12e−2δ+λ1t=λ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ2+λ2λ12δ+λ2δ22δ+λ1−2λ1λ2δδ+λ11+λ2δe−δt+2λ22δ2e−2δt−λ22δ21−2δa¯t2δμ12+2λ2λ1−λ2δδ+λ1e−2δ+λ1tμ12−2λ2λ1−λ21−2δ+λ1a¯t2δ+λ1μ12δ2δ+λ1=λ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ2+2λ2δλ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ12−eδt2λ1λ2δδ+λ11+λ2δ−2λ22δ2e−δt+2λ2λ1−λ2δδ+λ1e−δ+λ1tE56

πZd2t=λ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ2+2λ2δλ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ12−eδt2λ1λ2δδ+λ11+λ2δ−2λ22δ2e−δt+2λ2λ1−λ2δδ+λ1e−δ+λ1t=λ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ2+2λ2δλ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ12−eδt2λ1λ2δδ+λ11+λ2δ−2λ22δ21−δa¯t2δ+2λ2λ1−λ2δδ+λ11−δ+λ1a¯tδ+λ1=λ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ2+2λ2δλ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ12−2λ2δeδtλ2a¯tδ+λ1−λ2a¯tδ+λ1=λ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ2+2λ2δλ2a¯t2δ+λ1−λ2a¯t2δ+λ1−eδtλ2a¯tδ+λ1−λ2a¯tδ+λ1μ12.E57

To finally have:

πZd2t=λ2a¯t2δ+λ1−λ2a¯t2δ+λ1μ2+2λ2δλ2a¯t2δ+λ1−λ2a¯t2δ+λ1−eδtλ2a¯tδ+λ1−λ2a¯tδ+λ1μ12,E58

which is exactly the result of Léveillé et al.[15].

If λ1=λ2then

γ0λ1=λ2=λμ22δ+λ2μ12δ2+θλ2δδ+λ+θλ2δδ+2λμ1μ1′−μ1+θλμ2′−μ22δ+λ+θ2λ2μ1′−μ12δ+λδ+2λ,E59

γ1λ1=λ2=−2λ2μ12δ2−2θλ2μ1μ1′−μ1δδ+2λ, γ2λ1=λ2=−λμ22δ+λ2μ12δ2+θλ2μ1μ1′−μ1δδ−2λ,γ3λ1=λ2=0(60)

γ4λ1=λ2=−2θλ2μ1μ1′−μ1δ+2λδ−2λ−2θ2λ2μ1′−μ12δδ+2λ,E61

γ5λ1=λ2=−θλμ2′−μ2δ+λ+2θλ2μ1μ1′−μ1δ+λδ+2λ+2θ2λ2μ1′−μ12δδ+λ,E62

γ6λ1=λ2=θλμ2′−μ22δ+λ−θλ2μ1μ1′−μ1δ+λδ+2λ−θ2λ2μ1′−μ12δδ+λ.E63

Then,

πZd2t=λμ21δ−e−δt2δ+θλμ2′−μ212λ+2δ−e−2λ+2δt2λ+2δ+2λ2μ1212δ2−e−δtδ2+e−2δt2δ2+2θλ2μ1μ1′−μ112δ2λ+δ−e−2λ+δt2λ+δ−2λ+δ+e−2δt2δ−2λ+δ+2θλ2μ1μ1′−μ11δ2λ+2δ−e−δtδ2λ+δ+e−2λ+2δt2λ+2δ2λ+δ+2θ2λ2μ1′−μ1212λ+δ2λ+2δ−e−2λ+δtδ2λ+δ+e−2λ+2δtδ2λ+2δE64

which is exactly the result of Bargès et al. [8].

If λ1=λ2and θ=0then

πZo2t=λa¯t2δμ2+λa¯t2δμ12,E65

which is exactly the result of Léveillé et al. [15].

Remark 3.1

Noting for i=1,2,…,m,j=1,2,…,m, p=0,1and k∈N−0

ηkijp=ijθpλkEXj1−pEX′j−EXjpr+p×2λ+iδk.E66

We can rewrite π˜Zorand π˜Zo2ras

π˜Zor=1rηk110+ηk111,E67

π˜Zo2r=1rηk220+ηk221+ηk210+ηk211ηk110+ηk111=1rηk220+ηk221+ηk210ηk110+ηk210ηk111+ηk211ηk110+ηk211ηk111E68

The term π˜Zomrcan also be expressed using

π˜Zomr=1r∑i=1mi1j1p1…injnpn∈℘m,nηkinjnpn×…×ηki1j1p1,E69

where ℘m,n=i1j1p1,…,injnpn:i1=m,i1+…+in=m−1+n,i1>…>in,j1=m−1+n,j1+…+jn=m,j1>…>jn,p∈01, (70)

and

π˜Zdmr=2λ1λ1−λ2rr+mδ+λ1r+mδ+2λ11rμm+∑k=1m−1CmkEX1kπ˜Zom−kr+λ1r+mδ+λ1r+mδr+mδ+2λ2λ2r+mδ+2λ1+1π˜Z0mr,E71

4. Application

4.1 First two moments

For the numerical illustration, suppose that X∼pExpβ1=180+1−pExpβ2=1200, the inter-claim time distribution parameters λ1=2;4and λ2=1, the interest rate δ=3%(Tables 1–4). We use three different values for the copula parameter θ=−1;0;1, p=13and fix the time t=1;10;100. The mth moment of Xis

θ	t=1	t=10	t=100
−1	482.3375	4450	16,428
0	438.1057	4407.1	16,385
1	393.874	4364.2	16,342

Table 1.

EZdtfor λ1=1,λ2=10,δ=3%.

θ	t=1	t=10	t=100
−1	597.1633	4578.7	16,557
0	554.5237	4535.5	16,513
1	511.8841	4492.3	16,470

Table 2.

EZdtfor λ1=5,λ2=10,δ=3%.

θ	t=1	t=10	t=100
−1	33900.85	29646.71	7928.481
0	23972.65	20976.12	5711.548
1	359.8795	1036.223	1542.412

Table 3.

StdZdtfor λ1=1,λ2=10,δ=3%.

θ	t=1	t=10	t=100
−1	34027.1	29756.44	7954.606
0	24061.85	21053.72	5731.564
1	366.1484	1035.491	1545.671

Table 4.

StdZdtfor λ1=5,λ2=10,δ=3%.

μm=pm!β1m+1−pm!β2mand μm′=∫0∞mxm−11−FXx2dx=μm=pm!2β1m+1−pm!2β2m.(72)

4.2 Premium calculation

From the results in Section 4.1, we can compute the premium related to the risk of an insurance portfolio represented by Gt, depending on the premium calculation principles adopted by the insurance company. The loaded premium Zdtconsists in the sum of the pure premium EZdt, the expected value of the costs related to the portfolio, and a loading for the risk Mtas

Gt=EZdt+MtE73

The loading for the risk differs according to the premium calculation principles.

4.2.1 The expected value principle

Denote by θ>0the safety loading. The expected value principle defines the loaded premium as:

Gt=EZdt+θEZdt,E74

where Mt=θEZdt.

4.2.2 The variance principle

Denote by θ>0the safety loading. The variance principle defines the loaded premium as:

Gt=EZdt+θVarZdt,E75

where Mt=θVarZdt.

4.2.3 The standard deviation principle

Denote by θ>0the safety loading. The standard deviation principle defines the loaded premium as:

Gt=EZdt+θVarZdt,E76

where Mt=θVarZdt.

4.2.4 The quantile principle

The standard deviation principle defines the loaded premium as:

Gt=FZdt−11−ε,E77

where εis smallest (for example: ε=0.5%,1%,2.5%,5%).

In this case, the safety loading Mtis given by

Mt=FZdt−11−ε−EZdtE78

The principles of standard deviation and variance only require partial information on the distribution of the random variable, Zdt, i.e., its expectation and its variance.

Often, the actuary only has this information for different reasons (time constraints, information …).

If the actuary has more information about the random variable, Zdti.e., he knows the form of FZdt, then he can apply the quantile principle.

But he does not know much about FZdt, then he can approximate the distribution of Zdtusing the matching moments technique.

5. Conclusion

We have derived exact expressions for all the moments of the DCDPRV process using renewal arguments, again disproving the popular belief that renewal techniques cannot be applied in the presence of economic factors. Our results, for the DCDPRV process, are consistent: (i) with the results of Léveillé et al. [15] for θ=0,λ1≠λ2and for θ=0,λ1=λ2, (ii) with the results of Bargès et al. [8] for θ≠0,λ1=λ2.

Within this framework, further research is needed to get exact expressions (or approximations) of certain functional of the Zdtt≥0process, as stop-loss premiums and ruin probabilities.

Our models have applications in reinsurance, house insurance and car insurance. They can also be used in evaluation of health programs, finance, and other areas.

For example, consider the case of a male currently aged 25 who is starting a defined contribution (DC) pension plan and is planning to retire in, say, 40 years at the age of 65. He anticipates that when he reaches that age he will convert his accumulated pension fund into a life annuity in order to hedge his own longevity risk and avoid outliving his own financial resources. The value of his retirement income will depend not only on the value of his pension fund, but also on the price of annuities at the time. Other things being equal, this means that his retirement income prospects will be affected by the distribution on future annuity value: the greater the dispersion of that distribution, the riskier his retirement income will be. For the assessment of the accumulated pension fund and its variability our models can be used. We can suppose that this man makes a deposit to a bank account, and that the time between successive deposits follows a renewal process and the force of interest is stochastic. Our model allows us to calculate the accumulated pension fund and its variability at the age of 65.

Another possible application is in reliability, to model the net present value of aggregate equipment failures costs until its total breakdown. A piece of equipment is deemed to be beyond repair when the repair time exceeds a predetermined gap. Of course, another possible definition of total breakdown is when the cost of repair exceeds a predetermined gap. But, since the cost of repair is defined per unit time, the two definitions are somewhat equivalent.

Moments of the Discounted Aggregate Claims with Delay Inter-Occurrence Distribution and Dependence Introduced by a FGM Copula

Abstract

Keywords

chapter and author info

Author

Franck Adékambi*

From the Edited Volume

1. Introduction

2. The model

2.1 The dependence

3. Recursive expression for higher moments

3.1 Delay renewal case

4. Application

4.1 First two moments

Table 1.

Table 2.

Table 3.

Table 4.

4.2 Premium calculation

4.2.1 The expected value principle

4.2.2 The variance principle

4.2.3 The standard deviation principle

4.2.4 The quantile principle

5. Conclusion

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Over 21,000 IntechOpen readers like this topic

chapter statistics

More statistics for editors and authors

Related Content

This Book

Probability, Combinatorics and Control

Modelling the Information-Psychological Impact in Social Networks

Related Book

Probabilistic Modeling in System Engineering

Probabilistic Modelling in Solving Analytical Problems of System Engineering