Loading [MathJax]/jax/output/SVG/jax.js
Research article

Exact expression of ultimate time survival probability in homogeneous discrete-time risk model

  • Received: 04 September 2022 Revised: 15 November 2022 Accepted: 29 November 2022 Published: 13 December 2022
  • MSC : 60G50, 60J80, 91G05

  • In this work, we set up the generating function of the ultimate time survival probability φ(u+1), where

    φ(u)=P(supn1ni=1(Xiκ)<u),

    uN0,κN and the random walk {ni=1Xi,nN} consists of independent and identically distributed random variables Xi, which are non-negative and integer-valued. We also give expressions of φ(u) via the roots of certain polynomials. The probability φ(u) means that the stochastic process

    u+κnni=1Xi

    is positive for all nN, where a certain growth is illustrated by the deterministic part u+κn and decrease is given by the subtracted random part ni=1Xi. Based on the proven theoretical statements, we give several examples of φ(u) and its generating function expressions, when random variables Xi admit Bernoulli, geometric and some other distributions.

    Citation: Andrius Grigutis. Exact expression of ultimate time survival probability in homogeneous discrete-time risk model[J]. AIMS Mathematics, 2023, 8(3): 5181-5199. doi: 10.3934/math.2023260

    Related Papers:

    [1] Lihong Guan, Xiaohong Wang . A discrete-time dual risk model with dependence based on a Poisson INAR(1) process. AIMS Mathematics, 2022, 7(12): 20823-20837. doi: 10.3934/math.20221141
    [2] Huifang Yuan, Tao Jiang, Min Xiao . The ruin probability of a discrete risk model with unilateral linear dependent claims. AIMS Mathematics, 2024, 9(4): 9785-9807. doi: 10.3934/math.2024479
    [3] Carey Caginalp . A minimization approach to conservation laws with random initialconditions and non-smooth, non-strictly convex flux. AIMS Mathematics, 2018, 3(1): 148-182. doi: 10.3934/Math.2018.1.148
    [4] You Lv . Asymptotic behavior of survival probability for a branching random walk with a barrier. AIMS Mathematics, 2023, 8(2): 5049-5059. doi: 10.3934/math.2023253
    [5] Zhipeng Liu, Cailing Li, Zhenhua Bao . On a dependent risk model perturbed by mixed-exponential jump-diffusion processes. AIMS Mathematics, 2025, 10(4): 9882-9899. doi: 10.3934/math.2025452
    [6] J.-C. Cortés, A. Navarro-Quiles, J.-V. Romero, M.-D. Roselló . First-order linear differential equations whose data are complex random variables: Probabilistic solution and stability analysis via densities. AIMS Mathematics, 2022, 7(1): 1486-1506. doi: 10.3934/math.2022088
    [7] Xian Wen, Haifeng Huo, Jinhua Cui . The optimal probability of the risk for finite horizon partially observable Markov decision processes. AIMS Mathematics, 2023, 8(12): 28435-28449. doi: 10.3934/math.20231455
    [8] Clara Burgos, Juan Carlos Cortés, Elena López-Navarro, Rafael Jacinto Villanueva . Probabilistic analysis of linear-quadratic logistic-type models with hybrid uncertainties via probability density functions. AIMS Mathematics, 2021, 6(5): 4938-4957. doi: 10.3934/math.2021290
    [9] Mengrao Ma, Linmin Hu, Yuyu Wang, Fang Luo . Discrete-time stochastic modeling and optimization for reliability systems with retrial and cold standbys. AIMS Mathematics, 2024, 9(7): 19692-19717. doi: 10.3934/math.2024961
    [10] Khaled M. Alqahtani, Mahmoud El-Morshedy, Hend S. Shahen, Mohamed S. Eliwa . A discrete extension of the Burr-Hatke distribution: Generalized hypergeometric functions, different inference techniques, simulation ranking with modeling and analysis of sustainable count data. AIMS Mathematics, 2024, 9(4): 9394-9418. doi: 10.3934/math.2024458
  • In this work, we set up the generating function of the ultimate time survival probability φ(u+1), where

    φ(u)=P(supn1ni=1(Xiκ)<u),

    uN0,κN and the random walk {ni=1Xi,nN} consists of independent and identically distributed random variables Xi, which are non-negative and integer-valued. We also give expressions of φ(u) via the roots of certain polynomials. The probability φ(u) means that the stochastic process

    u+κnni=1Xi

    is positive for all nN, where a certain growth is illustrated by the deterministic part u+κn and decrease is given by the subtracted random part ni=1Xi. Based on the proven theoretical statements, we give several examples of φ(u) and its generating function expressions, when random variables Xi admit Bernoulli, geometric and some other distributions.



    The study of the sum of independent and identically distributed random variables ni=1Xi is hardly avoidable in probability theory and related fields. This sequence of sums {ni=1Xi,nN} is called the random walk. Let us define the stochastic process

    W(n):=u+κnni=1Xi,nN, (1.1)

    where uN0:=N{0}, κN and random variables Xi,iN are independent, identically distributed, non-negative and integer-valued. If κ=1, the defined process (1.1) is known as a discrete-time risk model; see [1]. Allowing κN, we call the process (1.1) a generalized premium discrete-time risk model; see [2]. Such types of processes appear in insurance mathematics (ruin theory), arguing that they describe an insurer's wealth in time moments nN, where u means the initial surplus (also called capital or reserve), κ denotes the premium rate (earnings per unit of time), i.e., (n+1)κnκ=κ, and the random walk {ni=1Xi,nN} represents the expenses (payoffs) caused by random size claims. Then, one can become curious to know whether the initial surplus u plus the gained premiums κn are sufficient to cover the incurred random expenses ni=1Xi. More precisely, one aims to know whether W(n)>0 for all n{1,2,,T} when T is some fixed natural number or T. The positivity of W(n) is of course associated with the probability. For the model given in (1.1), we define the finite time survival probability:

    φ(u,T):=P(Tn=1{W(n)>0})=P(sup1nTni=1(Xiκ)<u),TN,

    and the ultimate time survival probability:

    φ(u):=P(n=1{W(n)>0})=P(supn1ni=1(Xiκ)<u). (1.2)

    Both φ(u,T) and φ(u) are nothing but distribution functions of the provided integer-valued sequence of sums of random variables; these functions are left-continuous, non-decreasing and step functions if we allow uR. Also, φ()=1 if EX<κ; see Section 2. In particular, φ(0) is interpreted as the ultimate time survival probability when an insurer starts the activity with no initial surplus, i.e., when u=0. Then, the insurer maintains chances to "persist alive" if the payoff's size in the first moment of time n=1 is less than κ, i.e., if X1<κ.

    Calculation of φ(u,T) is simple; see, for instance, [2,Theorem 1]. Let us turn to the ultimate time survival probability φ(u). The law of total probability and rearrangements in (1.2) imply

    φ(u)=u+κi=1xu+κiφ(i); (1.3)

    see [2,page 3].

    By setting u=0 in (1.3), we get

    φ(0)=xκ1φ(1)+xκ2φ(2)++x0φ(κ); (1.4)

    to calculate the probability φ(κ) when x0>0, we must know the initial ones φ(0),φ(1),,φ(κ1). The calculation of φ(u), when u=κ,κ+1,, using the recurrence equation (1.3), requires φ(0),φ(1),,φ(κ1) too. The needed quantity of these initial values is X distribution-dependent, as some of the probabilities x0,x1,,xκ1 may equal to zero, cf. (1.4) when P(X>j)=1 for some j0. The paper [2] deals with finding the mentioned initial values φ(0),φ(1),,φ(κ1), and it is shown there that they can be found by calculating the limits of certain recurrent sequences. For instance, if κ=2 and x0>0, then it follows by (1.4) that

    φ(0)=x1φ(1)+x0φ(2),

    where (see [3,pages 2 and 3])

    φ(0)=φ()limnγn+1γn|βnγnβn+1γn+1|,φ(1)=φ()limnβnβn+1|βnγnβn+1γn+1|; (1.5)
    β0=1,β1=0,βn=1x0(βn2n1i=1xniβi),forn2,γ0=0,γ1=1,γn=1x0(γn2n1i=1xniγi),forn2,

    and φ()=1 if EX<2.

    Calculating the limits in (1.5) and aiming to prove that the provided determinant 2×2 never vanishes; in paper [3], it was proved the connection to the solutions of s2=GX(s), where sC,|s|1 and GX(s) is the probability-generating function of the random variable X. On top of that, it was realized in [3] that the values of φ(0) and φ(1) in (1.5) can be derived by using the classical stationarity property for the distribution of the maximum of a reflected random walk; see [4,Chapter Ⅵ,Section 9]. Using the mentioned stationarity property, the generating function of φ(u+1),uN0 for κ=2 was found in [3,Theorem 5]; however, this required the finiteness of the second moment of the random variable X, i.e., EX2<. In this article, we extend the work in [3] and find the generating function of φ(u+1),uN0 for an arbitrary κN. Moreover, we show that the requirement of EX2< is redundant and provide exact expressions of φ(u),uN0 via solutions of systems of linear equations which are based on the roots of sκ=GX(s) and Vandermonde-like matrices.

    For the short overview of the literature, we mention that the references [1,5,6,7,8,9,10,11,12,13] are known as the classical ones on the wide subject of renewal risk models, while [14,15,16] might be mentioned as the recent ones. The main reason for so much literature is that the ruin theory, being random walk-based, is heavily dependent on the random walk's structural assumptions, such as the independence of random variables, their distributions, etc. This work is also closely related to branching and Galton-Watson processes and queueing theory; see [17] and related papers. See also [18] or [19, Figure 1] on random walk occurrence in number theory. Last but not least, it is worth mentioning that Vandermonde matrices have a broad range of occurrences, from pure mathematics to many other applied sciences; see [20] and related works.

    Let

    M:=supn1(ni=1(Xiκ))+,

    where x+=max{0,x}, xR is the positive part function and the random variables Xi and κN are the same as in the model (1.1). Let us denote the probability mass function of the random variable M by

    πi:=P(M=i),iN0.

    Then, the ultimate time survival probability definition (1.2) implies that

    φ(u+1)=ui=0πi=P(Mu)foralluN0. (2.1)

    In general, the random variable M can be extended, i.e., P(M=)>0; however, the condition EX<κ ensures

    limuφ(u)=P(M<)=1;

    see [2,Lemma 1]. This condition EX<κ is called the net profit condition, and it is crucial because survival is impossible, i.e., φ(u)=0 for all uN0, if EXκ, except for a few trivial cases when P(X=κ)=1; see [2,Theorem 9]. Intuitively, it is clear that long-term survival by the model (1.1) is impossible if the threatening claim amount X on average is equal or greater to the collected premium κ per unit of time.

    For sC, let us denote the generating function of φ(1),φ(2), as follows:

    Ξ(s):=i=0φ(i+1)si,|s|<1

    and the probability-generating functions of the random variables X and M:

    GX(s):=i=0xisi,GM(s):=i=0πisi,|s|1.

    Then, Ξ(s) and GM(s), for |s|<1, satisfy the relation

    Ξ(s)=i=0φ(i+1)si=i=0ij=0πjsi=j=0πji=jsi=j=0πjsj1s=GM(s)1s. (2.2)

    In many examples, the radius of convergence of GX(s) or GM(s) is larger than one. See [3,Lemma 8] for more properties of the probability-generating function in |s|1.

    In this section, based on the previously introduced notations and relation (1s)Ξ(s)=GM(s) in (2.2), we formulate the main results of the work.

    Theorem 3.1. Let us consider the model defined in (1.1) and suppose that the net profit condition EX<κ holds. Then, the probability mass functions of the random variables M and X satisfy the following two equalities:

    GM(s)(sκGX(s))=κ1i=0πiκ1ij=0xj(sκsi+j),|s|1, (3.1)
    κEX=κ1i=0πiκ1ij=0xj(κij). (3.2)

    We prove Theorem 3.1 in Section 5.

    Equality (3.1) implies the following relation among the probabilities π0, π1, .

    Corollary 3.1. Let πi=P(M=i),iN0 and FX(u)=ui=0xi,uN0 be the distribution function of the random variable X. Then, for κN, the following equalities hold:

    πκx0=π0κ1i=0πiFX(κi),πnx0=πnκκ1i=0πixni,n=κ+1,κ+2,. (3.3)

    Proof of Corollary 3.1. The n-th derivative of both sides of the equality (3.1) and s0 gives

    or

    πnx0=πnκn1i=0πixni,n=κ+1,κ+2,.

    Let us turn to the survival probabilities φ(1),φ(2), generating function Ξ(s). It is easy to see that the equalities (2.2) and (3.1) imply

    Ξ(s)=κ1i=0πiκ1ij=0xj(sκsi+j)(1s)(sκGX(s)). (3.4)

    Therefore, in a similar way that the recurrence equation (1.3) requires the initial values of φ(0), φ(1), , φ(κ1), the generating function Ξ(s) in (3.4) (the equality (3.3) as well) requires π0,π1,,πκ1, κN. These probabilities can be solved by using the relations (3.1) and (3.2) and this is achievable as provided in Items (i)(iv) below:

    (i) We can choose |s|1 such that the left-hand side of (3.1) vanishes, i.e., the roots of sκ=GX(s).

    If the net profit condition GX(1)=EX<κ holds and the greatest common divisor of powers of s in sκ=GX(s) is one, there are exactly κ1 roots of sκ=GX(s) in |s|<1 when counted with their multiplicities. This fact is implied by Rouché's theorem and estimate |GX(s)|1<|λsκ| when λ>1 and |s|=1, which means that, because of the fundamental theorem of algebra, both functions λsκ and λsκGX(s) have κ zeros in |s|<1. When λ1+, there is always one root out of those κ in |s|<1 migrating to s=1 (s=1 is always the root of sκ=GX(s)), and some to other boundary points |s|=1 (roots of unity) if the greatest common divisor of powers of s in sκ=GX(s) is greater than one; see [21,Chapter 10], [22,Remark 10] and [3,Section 4,Lemmas 9 and 10 therein].

    (ii) Let α1 be a root of sκ=GX(s) in |s|1 and denote π:=(π0,π1,,πκ1)T as the column vector. Then, by (3.1) and

    (αj+αj+1++ακ1)(α1)=ακαj,j{0,1,,κ1},

    it holds that

    0=(κ1j=0xj(ακαj),κ2j=0xj(ακαj+1),,x0(ακακ1))π=(κ1j=0xjκ1i=jαi,κ2j=0xjκ1i=j+1αi,,x0ακ1)π=(κ1j=0αjFX(j),κ2j=0αj+1FX(j),,ακ1x0)π=κ1i=0πiκ1ij=0αj+iFX(j),

    where FX(u) is the distribution function of X.

    (iii) Let α1,,ακ11 be the roots of sκ=GX(s) in |s|1. Then, by (i), (ii) and (3.2),

    (κ1j=0αj1FX(j)κ2j=0αj+11FX(j)ακ11x0κ1j=0αjκ1FX(j)κ2j=0αj+1κ1FX(j)ακ1κ1x0κ1j=0xj(κj)κ2j=0xj(κj1)x0)(π0πκ2πκ1)=(00κEX). (3.5)

    If Aπ=B denotes the system (3.5), x0>0 and α1,α2,,ακ11 are the roots of multiplicity one, then, according to Lemma 4.2 proved in Section 4, the determinant |A|0 and, therefore, π=A1B.

    (iv) Suppose the root α1 of sκ=GX(s) in |s|1 is of multiplicity l{2,3,,κ1},κ3. Then, according to the equality (3.1) in Theorem 3.1 and (ii), the derivatives

    dmdsm(κ1i=0πiκ1ij=0sj+iFX(j))|s=α=0forallm{0,1,,l1}, (3.6)

    and, to avoid identical lines in matrix A, we can set up the modified system (3.5) by replacing its lines (except the last one) with the corresponding derivatives (3.6). If x0>0, such a modified main matrix A remains non-singular, as proved in Lemma 4.3 of Section 4.

    Note 1: The condition x0>0 does not lose generality. If P(X>j)=1 for some j{0,1,,κ2},κ2 and the net profit condition remains valid (note that P(X>κ1) implies EXκ), then there is a reduction the order of recurrence in (1.3) and, consequently, some terms in the sums of (3.1) and (3.2) equal zero, causing corresponding adjustments in the system (3.5) or its modified version described in (iv). We then end up dividing by some xj+1 instead of x0 where needed. For instance, if x0=0 and x1>0, we then can express φ(κ1) from (1.4) dividing by x1. See also [2,Theorem 7] and Corollary 3.2 when x0=0. Also, the both sides of sκ=GX(s) can be canceled by some power of s0 if P(X>j)=1 for some j{0,1,,κ2},κ2.

    We further denote by |A| the determinant of the matrix A where Mi,j, i,j{1,2,,κ}, κN are its minors and the matrix A is the main matrix in (3.5) or its modification replacing the coefficients by derivatives, as described in (iv).

    The equality (3.4) and thoughts listed in (i)(iv) allow us to formulate the following statement.

    Theorem 3.2. Let |s|<1 and sκGX(s)0. If the net profit condition EX<κ holds, then the survival probability-generating function is given by

    Ξ(s)=κEXGX(s)sκκ1i=0˜πiκ1ij=0sj+iFX(j), (3.7)

    where ˜πi=πi/(κEX),

    ˜π0=(1)κ+1Mκ,1|A|,˜π1=(1)κ+2Mκ,2|A|,,˜πκ1=Mκ,κ|A|

    and the matrix A is created as provided in (i)–(iv).

    Moreover, the initial values for the recurrence equation (1.3), including φ(κ), are

    φ(0)=κEX|A|κi=1(1)κ+iMκ,iFX(κi),φ(u)=κEX|A|ui=1(1)κ+iMκ,i,u=1,2,,κ.

    We prove Theorem 3.2 in Section 5.

    Note 2: We agree that, for κ=1, the matrix A=(x0), its determinant |A|=x0 and the minor M1,1=1. Recall that x0 gets replaced by some xj+1 if P(X>j)=1 for some j{0,1,,κ2},κ2 and the net profit condition holds; see Note 1.

    The next statement provides possible expressions of ˜π0,˜π1,,˜πκ1 and φ(0),φ(1),,φ(κ), κN.

    Theorem 3.3. Suppose that x0>0 and α1,α2,,ακ11 are the roots of multiplicity one of sκ=GX(s) in |s|1. Then, the values ~πi=πi/(κEX) for i=0,1,,κ1 admit the following representation:

    ˜π0=1x0κ1j=1αjαj1,˜π1=1j1<<jκ2κ1αj1αjκ2x0κ1j=1(αj1)FX(1)x0˜π0,˜π2=1j1<<jκ3κ1αj1αjκ3x0κ1j=1(αj1)FX(2)x0˜π0FX(1)x0˜π1,˜πκ1=(1)κ+1x0κ1j=11αj11x0κ2i=0˜πiFX(κ1i),κ2,

    and the initial values for the recurrence equation (1.3), including φ(κ), are

    ˜φ(0)=(1)κ+1κ1j=11αj1,˜φ(1)=1x0κ1j=1αjαj1,˜φ(2)=FX(1)x0˜φ(1)+κ1j=11/x0αj1(κ1j=1αj1j1<<jκ2κ1αj1αjκ2),˜φ(3)=FX(1)x0˜φ(2)FX(2)x0˜φ(1)+κ1j=11/x0αj1×(κ1j=1αj1j1<<jκ2κ1αj1αjκ2+1j1<<jκ3κ1αj1αjκ3),˜φ(κ)=1x0κ1i=1FX(κi)˜φ(i)+κ1j=11/x0αj1×(κ1j=1αj1j1<<jκ2κ1αj1αjκ2+1j1<<jκ3κ1αj1αjκ3++(1)κ+1),

    κ2, where

    ˜φ(i)=φ(i)(κEX),i{0,1,,κ}.

    Note that 0j=1()=1j1<j0()=1 in Theorem 3.3, and we prove this theorem in Section 5.

    In view of Theorems 3.2 and 3.3, we give several separate expressions on Ξ(s).

    Corollary 3.2. If κ=1, then

    Ξ(s)=1EXGX(s)s.

    If κ=2 and x0>0, then

    Ξ(s)=2EXα1αsGX(s)s2,

    where α[1,0) is the unique root of GX(s)=s2.

    If κ=2, x0=0 and x1>0, then

    Ξ(s)=2EX˜GX(s)s,

    where ˜GX(s)=i=0xi+1si,|s|1.

    Proof of Corollary 3.2. The provided Ξ(s) expressions are implied by Theorem 3.2. Recall that s2=GX(s),x0>0 has the unique real root α[1,0). In addition, when x0>0, then α=1 is the root of s2=GX(s) iff P(X2N0)=1; see [3,Section 4 and Corollary 15 therein] and the description (i) in Section 3.

    In this section, we formulate and prove several auxiliary statements needed to derive the main results stated in Section 3.

    Lemma 4.1. The random variable

    M=supn1(ni=1(Xiκ))+,

    where x+=max{0,x} is the positive part of xR, admits the following distribution property:

    (M+Xκ)+d=M.

    Proof. The proof is straightforward according to the definition of M and basic properties of the maximum. Indeed,

    (M+Xκ)+=max{0,max{0,supn1ni=1(Xiκ)}+Xκ}d=max{0,max{X1κ,supn2ni=1(Xiκ)}}d=max{0,supn1,ni=1(Xiκ)}=M.

    See also, [22,Lemma 5.2], [3,Lemma 25] and [4,page 198].

    Lemma 4.2. Let α1,,ακ11 be the roots of multiplicity one of sκ=GX(s) in the region |s|1, and suppose that the probability x0 is positive. Then, the determinant |A| of the main matrix in (3.5) is

    |A|=xκ0(1)κ+1κ1j=1(αj1)1i<jκ1(αjαi)0.

    Proof. Let us calculate the determinant

    |A|=|κ1j=0αj1FX(j)κ2j=0αj+11FX(j)ακ21x0+ακ11FX(1)ακ11x0κ1j=0αjκ1FX(j)κ2j=0αj+1κ1FX(j)ακ2κ1x0+ακ1κ1FX(1)ακ1κ1x0κ1j=0xj(κj)κ2j=0xj(κj1)2x0+x1x0|.

    We first put forward x0 from the last column. Then, multiplying the last column by FX(κ1),FX(κ2),,FX(1), respectively, and subtracting it from the first, the second, etc., columns, we obtain

    |A|=x0|κ2j=0αj1FX(j)κ3j=0αj+11FX(j)ακ21x0ακ11κ2j=0αjκ1FX(j)κ3j=0αj+1κ1FX(j)ακ2κ2x0ακ1κ1κ2j=0xj(κj1)κ3j=0xj(κj2)x01|.

    Proceeding the similar with the penultimate column of the last determinant (to put forward x0 and rearrange) and so on and applying the basic determinant properties, we obtain that

    |A|=xκ0|1α1ακ111ακ1ακ1κ1111|=xκ0(1)κ+1|α11α211ακ111α21α221ακ121ακ11α2κ11ακ1κ11|=xκ0(1)κ+1κ1j=1(αj1)|1α1ακ211α2ακ221ακ1ακ2κ1|.

    The last determinant is nothing but the well-known Vandermonde determinant; see for example [23,Section 6.1]. Thus,

    |A|=xκ0(1)κ+1κ1j=1(αj1)1i<jκ1(αjαi)0,

    because the roots α1,α2,,ακ1 are distinct and lie in the region |s|1, s1.

    Lemma 4.3. Let |s|1. Suppose some roots α1,,ακ11 of GX(s)=sκ are multiple, and assume that the probability x0 is positive. Then, the modified main matrix in (3.5), after replacing its lines (except the last one) by the derivatives (3.6), remains non-singular.

    Proof. In short, the statement follows because the derivative is a linear mapping. More precisely, let us assume that α1 is of multiplicity two. Then, there exists such sufficiently close to zero δR{0} that the matrix with the replaced second line

    (κ1j=0αj1FX(j)κ2j=0αj+11FX(j)ακ11x0κ1j=0(α1+δ)jFX(j)κ2j=0(α1+δ)j+1FX(j)(α1+δ)κ1x0κ1j=0αjκ1FX(j)κ2j=0αj+1κ1FX(j)ακ1κ1x0κ1j=0xj(κj)κ2j=0xj(κj1)x0) (4.1)

    is non-singular, see the expression of the determinant in Lemma 4.2. Then, subtracting the second line from the first in (4.1), dividing the first line by δ afterward and letting δ0, we get the desired line replacement using the derivative.

    The proof is analogous for higher derivatives and/or more multiple roots.

    In this section, we prove the statements formulated in Section 3. Let us start with the proof of Theorem 3.1.

    Proof of Theorem 3.1. By Lemma 4.1 and the rule of total expectation,

    GM(s)=Es(M+Xκ)+=E(E(s(M+Xκ)+|M))=κ1i=0πiEs(i+Xκ)++sκGX(s)i=κπisi=κ1i=0πi(Es(X+iκ)+siκGX(s))+sκGX(s)GM(s),

    which implies the equality (3.1):

    GM(s)(sκGX(s))=κ1i=0πi(Es(X+iκ)++κsiGX(s))=κ1i=0πiκ1ij=0xj(sκsi+j).

    To prove the second equality (3.2) in Theorem 3.1, we take the derivative of both sides of (3.1) with respect to s:

    S1+S2:=GM(s)(sκGX(s))+GM(s)(κsκ1GX(s))=κ1i=0πiκ1ij=0xj(κsκ1(i+j)si+j1)=:S3.

    We now let s1 in the last equality. It is easy to see that

    lims1S3=κ1i=0πiκ1ij=0xj(κij)

    and

    lims1S2=κEX,

    because the net profit condition EX<κ holds. Before calculating lims1S1, we observe that EX2=EM= and EX2<EM<; see [24,Theorems 5 and 6]. Therefore, the requirement EX2< implies lims1S1=0 immediately. However, lims1S1=0 despite EM=. Indeed, if GM(s) as s1, then

    lims1S1=lims1sκGX(s)1/GM(s)=lims1κsκ1GX(s)GM(s)/(GM(s))2,

    where

    lim sups1(GM(s))2GM(s)NN1i=Nπi

    for any N{2,3,}; see [22,Lemma 5.5]. Thus, the equality (3.2) follows and the theorem is proved.

    Proof of Theorem 3.2. For sκGX(s)0, the equality (3.4) and division by 1s (see (ii) in Section 3) imply

    Ξ(s)=κ1i=0πiκ1ij=0sj+iFX(j)GX(s)sκ=1GX(s)sκ(κ1j=0sjFX(j),κ2j=0sj+1FX(j),,sκ1x0)(π0π1πκ1).

    By the system (3.5), including its modified version described in (iv) in Section 3, and the recalled notations π=(π0,π1,,πκ1)T and ˜πi=πi/(κEX), we obtain

    π=1|A|(M1,1M1,2(1)1+κM1,κM2,1M2,2(1)2+κM2,κ(1)κ+1Mκ,1(1)κ+2Mκ,2Mκ,κ)T(00κEX)=κEX|A|((1)κ+1Mκ,1(1)κ+2Mκ,2Mκ,κ)=(κEX)(˜π0˜π1˜πκ1).

    Thus, the expression of Ξ(s) in (3.7) follows.

    The claimed equalities on φ(u) for u=1,,κ are evident due to the obtained expression of π and φ(u+1)=ui=0πi,uN0 provided in (2.1). It can be seen that the recurrence equation (1.3) yields

    φ(0)=κi=1xκiφ(i)=κ1i=0πiFX(κ1i).

    Proof of Theorem 3.3. We calculate the minors Mκ,1,Mκ,2,,Mκ,κ of the following matrix:

    A=(κ1j=0αj1FX(j)κ2j=0αj+11FX(j)ακ11x0κ1j=0αjκ1FX(j)κ2j=0αj+1κ1FX(j)ακ1κ1x0κ1j=0xj(κj)κ2j=0xj(κj1)x0).

    Following the calculation of determinant |A| in the proof of Lemma 4.2, we get

    Mκ,1=|κ2j=0αj+11FX(j)κ3j=0αj+21FX(j)ακ11x0κ2j=0αj+1κ1FX(j)κ3j=0αj+2κ1FX(j)ακ1κ1x0|=xκ10|α1α21ακ11ακ1α2κ1ακ1κ1|=xκ10κ1i=1αi1i<jκ1(αjαi).

    Note that Mκ,1 is defined for κ1 and M1,1=1 by agreement. The next one

    Mκ,2=|κ1j=0αj1FX(j)κ3j=0αj+21FX(j)ακ11x0κ1j=0αjκ1FX(j)κ3j=0αj+2κ1FX(j)ακ1κ1x0|=xκ20|x0+α1FX(1)α21ακ11x0+ακ1FX(1)α2κ1ακ1κ1|=xκ101i<jκ1(αjαi)1j1<<jκ2κ1αj1αjκ2+FX(1)x0Mκ,1.

    Similarly as before, Mκ,2 is defined for κ2 only, and M2,2=x0+FX(1)α, where α[1,0) is the unique root of s2=GX(s); see (i) in Section 3 and [3,Section 4 and Corollary 15 therein].

    Proceeding,

    Mκ,3=|κ1j=0αj1FX(j)κ2j=0αj+11FX(j)κ4j=0αj+31FX(j)ακ11x0κ1j=0αjκ1FX(j)κ2j=0αj+1κ1FX(j)κ4j=0αj+3κ1FX(j)ακ1κ1x0|=xκ20|x0+α21FX(2)α1α31ακ11x0+α2κ1FX(2)ακ1α3κ1ακ1κ1|+FX(1)x0Mκ,2=xκ101i<jκ1(αjαi)1j1<<jκ3κ1αj1αjκ3FX(2)x0Mκ,1+FX(1)x0Mκ,2,κ3,

    and so on until the last minor:

    Mκ,κ=|κ1j=0αj1FX(j)κ2j=0αj+11FX(j)x0ακ21+ακ11FX(1)κ1j=0αjκ1FX(j)κ2j=0αj+1κ1FX(j)x0ακ2κ1+ακ1κ1FX(1)|=|κ1j=0αj1FX(j)κ2j=0αj+11FX(j)x0ακ21κ1j=0αjκ1FX(j)κ2j=0αj+1κ1FX(j)x0ακ2κ1|+FX(1)x0Mκ,κ1=xκ101i<jκ1(αjαi)+(1)κFX(κ1)x0Mκ,1+(1)κ+1FX(κ2)x0Mκ,2++(1)2κ1FX(2)x0Mκ,κ2+FX(1)x0Mκ,κ1.

    The statement on expressions of ˜π0,˜π1,,˜πκ1 follows dividing the obtained minors (multiplied by 1 where needed) by the determinant |A|.

    We now prove the claimed formulas of φ(0),φ(1),,φ(κ), κN. By the recurrence equation (1.3) with u=0, φ(u+1)=ui=0πi, uN0 and the already proved expression of πκ1,κN in Theorem 3.3,

    φ(0)=κ1i=0φ(i+1)xκ1i=κ1i=0πiFX(κ1i)=κEX(1)κ+1κ1j=11αj1.

    The formula for φ(1) is evident because φ(1)=π0, where the expression of π0 is already proved in Theorem 3.3, too. The rest is clear by calculating the sum φ(u+1)=ui=0πi, uN0, where πi are given in the first part of Theorem 3.3.

    In this section, we give several examples illustrating the applicability of theoretical statements formulated in Section 3. The required numerical computations were performed by using Wolfram Mathematica [25]. As mentioned in Section 1, in [2], it has been proved that the required initial values for the recurrence equation (1.3) can be approximately found by calculating certain recurrent limits, while results of this work in Section 3, in many instances, provide exact closed-form expressions of the survival probabilities. Therefore, in some considered examples here, we check if the calculated exact value of φ matches the previously known approximate one.

    Example 6.1. Suppose the random claim amount X is Bernoulli-distributed, i.e., 1P(X=0)=p=P(X=1),0<p<1 and the premium κN. We find the ultimate time survival probability-generating function Ξ(s) and calculate φ(u),uN0.

    If κ=1, in view of the first part of Corollary 3.2 and the recurrence equation (1.3), it is trivial that Ξ(s)=1/(1s),|s|<1 and φ(0)=x0φ(1)=1p, φ(u)=1, uN. In other words, the ultimate time survival is guaranteed if the initial surplus uN and the maximal claim size is one in the model u+nni=1Xi.

    If κ2, it is easy to understand that u+κnni=1Xi>0 for all nN, uN0, regardless of the size of Xi; consequently, φ1, Ξ(s)=1/(1s),|s|<1.

    Example 6.2. Suppose that the random claim amount X is distributed geometrically with the parameter p(0,1), i.e., P(X=i)=p(1p)i, i=0,1,, and the premium rate equals two, i.e., κ=2. We find the ultimate time survival probability-generating function Ξ(s) and calculate φ(0) and φ(1) when the net profit condition is satisfied, i.e., EX<2.

    We start with an observation on the net profit condition:

    EX=1pp<213<p<1.

    Then, according to Theorem 3.1 and the description (i) in Section 3,

    GX(s)=p1(1p)s=s2α:=s=p4p3p22(1p)(1,0)

    when 1/3<p<1, and by Corollary 3.2 with κ=2 and x0=p>0,

    Ξ(s)=(3p1)(p4p3p2)p(3p24p3p2)1(1p)s(1s)s2+(s31)p,13<p<1.

    For κ=2, u=0 and 1/3<p<1, the recurrence equation (1.3) or Theorem 3.3 yields

    φ(0)=x1φ(1)+x0φ(2)=(1p)pΞ(0)+pΞ(0)=2EX1α=3p2+4p3p22p,φ(1)=Ξ(0)=2EXx0αα1=3p4p3p22p2.

    One may check that, for p=101/300,

    φ(0)=90597297202=0.0197691,φ(1)=454501509059710201=0.0295066,

    and that coincides with the approximate values of φ(0) and φ(1) in [2,page 12] obtained via recurrent sequences.

    Example 6.3. Suppose that X attains the natural values only, i.e., x0=0, x1>0, κ=2 and the net profit condition is satisfied EX<2. We provide the ultimate time survival probability φ(u) formulas for all uN0.

    Let us recall that

    ˜GX(s)=i=0xi+1si,|s|1.

    The recurrence equation (1.3) and Corollary 3.2 for x0=0 and x1>0 imply

    φ(0)=x1φ(1)=2EX,φ(1)=Ξ(0)=2EXx1,φ(u)=2EX(u1)!du1dsu1(1˜GX(s)s)|s=0=1x1(φ(u1)u1i=1xui+1φ(i)),u2,

    which echoes and widens the statement of Theorem 3 in [2], providing another method of φ(u),u2 calculation.

    Example 6.4. Suppose that the random claim amount X is distributed geometrically with the parameter p=101/300, i.e., P(X=i)=p(1p)i, i=0,1,, and the premium rate equals three, i.e., κ=3. We set up the ultimate time survival probability-generating function Ξ(s) and calculate or provide formulas for φ(u),uN0.

    First, we observe that the net profit condition is satisfied, i.e., EX=199/101<3. We now follow the statement of Theorem 3.1 and the surrounding comments beneath it. Then, for p=101/300, the equation

    GX(s)=p1(1p)s=s3

    has two complex conjugate solutions α1:=0.368094+0.522097i and α2:=0.3680940.522097i inside the unit circle |s|<1. Then, by Theorem 3.2,

    Ξ(s)=2i=0πi2ij=0si+jFX(j)s3GX(s),

    where (π0,π1,π2)=(0.582072,0.0818989,0.0658497) is the unique solution of

    (x0+FX(1)α1+FX(2)α21x0α1+FX(1)α21x0α21x0+FX(1)α2+FX(2)α22x0α2+FX(1)α22x0α223x0+2x1+x22x0+x1x0)(π0π1π2)=(003EX),

    with appropriate numerical characteristics of the provided distribution. Theorem 3.3 and the recurrence equation (1.3) imply

    φ(0)=3i=1x3iφ(i)=3EX(1α1)(1α2)=0.480212,φ(1)=π0=(3EX)α1α2x0(1α1)(1α2)=0.582072,φ(2)=π0+π1=(3EX)x0(α1+α2)+x1α1α2x20(α11)(1α2)=0.663971,φ(3)=π0+π1+π2=0.729821,whereπ2=3EXx0(1(α11)(α21)1i=0πiFX(2i)),φ(u)=1x0(φ(u3)u1i=1xuiφ(i))=du1dsu1Ξ(s)(u1)!|s=0,u3.

    The provided values of φ(0),φ(1),φ(2) and φ(3) coincide with the ones given in [2,page 14], where they are obtained approximately from certain recurrent sequences.

    Example 6.5. Suppose x0=0.128, x1=0.576, x2=0.264, x3=0.032, 3i=0xi=1 and κ=3. We set up the ultimate time survival probability-generating function Ξ(s) and calculate φ(u),uN0.

    For the provided distribution EX=1.2<3, the equation

    0.128+0.576s+0.264s2+0.032s2=s3

    has one root s=4/11=:α of multiplicity two. Then, according to Theorem 3.1 and the comments (i)–(iv) beneath it, we create the modified system, replacing the second line with the corresponding derivatives:

    (x0+FX(1)α+FX(2)α2x0α+FX(1)α2x0α2FX(1)+2FX(2)αx0+2FX(1)α2x0α3x0+2x1+x22x0+x1x0)(π0π1π2)=(003EX),

    which implies (π0,π1,π2)=(1,0,0); consequently

    φ(0)=0.968,φ(u)=1,uN,Ξ(s)=11s,|s|<1.

    One may observe that the obtained result is expected, because u+3nni=1Xi>0 for all nN, except when u=0 and Xi attains the value of 3.

    This work shows that, if certain conditions are met, there exist exact closed-form expressions of the ultimate time survival probability

    φ(u)=P(n=1{u+κnni=1Xi>0}),

    where the roots of sκ=GX(s), the distance κEX>0 and the distribution function FX(s) are involved; see Theorem 3.3. Moreover, having the values of the probability mass function

    P(supn1(ni=1Xiκ)+=u),u=0,1,,κ1,

    we can get the exact expression of the survival probability-generating function

    Ξ(s)=i=0φ(i+1)si,|s|<1;

    see Theorem 3.2. As mentioned in Section 1, the expression of survival φ(u) or ruin 1φ(u) probability is heavily dependent on what type of random variables generate the random walk

    {ni=1(Xiκ),nN}. (7.1)

    The random variables Xi in the sequence (7.1) can be discrete/continuous or dependent/independent, and their quantity for each nN can be deterministic/random, etc. Equally, the premium, or just an intercept technically, κN in (7.1) influences the sequence's distribution, too. As demonstrated, the ultimate time survival probabilities φ(u) are solutions of systems of linear equations, which are based on the roots of sκ=GX(s). The recent work in [22] shows that similar systems can be used to find φ(u) when Xi are distributed differently. Thus, it is of interest to study the sequence (7.1), assuming various other mentioned options for Xi and κ; see, for instance, [26,27,28,29]. Also, the broadness of a random walk's occurrence in mathematics and other applied sciences indicates that this work and referenced research should not be applicable to ruin theory only.

    The author appreciates any kind of constructive criticism of the paper and thanks professor Jonas Šiaulys for the detailed reading of the first draft of the manuscript. Sincere thanks to the anonymous reviewers for their work.

    The author declares no conflict of interest.



    [1] S. M. Li, Y. Lu, J. Garrido, A review of discrete-time risk models, Rev. R. Acad. Cienc. Exactas Fis. Nat. Ser. A Mat., 103 (2009), 321–337. https://doi.org/10.1007/BF03191910 doi: 10.1007/BF03191910
    [2] A. Grigutis, J. Šiaulys, Recurrent sequences play for survival probability of discrete time risk model, Symmetry, 12 (2020), 2111. https://doi.org/10.3390/sym12122111 doi: 10.3390/sym12122111
    [3] A. Grigutis, J. Jankauskas, On 2×2 determinants originating from survival probabilities in homogeneous discrete time risk model, Results Math., 77 (2022), 204. https://doi.org/10.1007/s00025-022-01736-y doi: 10.1007/s00025-022-01736-y
    [4] W. Feller, An introduction to probability theory and its applications. Vol. Ⅱ., 2 Ed., New York: John Wiley & Sons, 1971.
    [5] E. S. Andersen, On the collective theory of risk in case of contagion between the claims, Trans. XVth Int. Congr. Actuaries, 2 (1957), 219–229.
    [6] F. Spitzer, Principles of random walk, 2 Ed., New York: Springer, 2001.
    [7] H. U. Gerber, Mathematical fun with the compound binomial process, ASTIN Bull., 18 (1988), 161–168. https://doi.org/10.2143/AST.18.2.2014949 doi: 10.2143/AST.18.2.2014949
    [8] H. U. Gerber, Mathematical fun with ruin theory, Insur. Math. Econ., 7 (1988), 15–23. https://doi.org/10.1016/0167-6687(88)90091-1 doi: 10.1016/0167-6687(88)90091-1
    [9] E. S. W. Shiu, Calculation of the probability of eventual ruin by Beekman's convolution series, Insur. Math. Econ., 7 (1988), 41–47. https://doi.org/10.1016/0167-6687(88)90095-9 doi: 10.1016/0167-6687(88)90095-9
    [10] E. S. W. Shiu, Ruin probability by operational calculus, Insur. Math. Econ., 8 (1989), 243–249. https://doi.org/10.1016/0167-6687(89)90060-7. doi: 10.1016/0167-6687(89)90060-7
    [11] F. De Vylder, M. J. Goovaerts, Recursive calculation of finite-time ruin probabilities, Insur. Math. Econ., 7 (1988), 1–7. https://doi.org/10.1016/0167-6687(88)90089-3 doi: 10.1016/0167-6687(88)90089-3
    [12] P. Picard, C. Lefèvre, Probabilité de ruine éventuelle dans un modèle de risque à temps discret, J. Appl. Probab., 40 (2003), 543–556. https://doi.org/10.1239/jap/1059060887 doi: 10.1239/jap/1059060887
    [13] S. M. Li, F. J. Huang, C. Jin, Joint distributions of some ruin related quantities in the compound binomial risk model, Stoch. Models, 29 (2013), 518–539. https://doi.org/10.1080/15326349.2013.847610 doi: 10.1080/15326349.2013.847610
    [14] L. Rincón, D. J. Santana, Ruin probability for finite negative binomial mixture claims via recurrence sequences, Comm. Statist. Theory Methods, 2022. https://doi.org/10.1080/03610926.2022.2087091 doi: 10.1080/03610926.2022.2087091
    [15] Y. Q. Cang, Y. Yang, X. X. Shi, A note on the uniform asymptotic behavior of the finite-time ruin probability in a nonstandard renewal risk model, Lith. Math. J., 60 (2020), 161–172. https://doi.org/10.1007/s10986-020-09473-x doi: 10.1007/s10986-020-09473-x
    [16] C. Lefèvre, M. Simon, Schur-constant and related dependence models, with application to ruin probabilities, Methodol. Comput. Appl. Probab., 23 (2021), 317–339. https://doi.org/10.1007/s11009-019-09744-2 doi: 10.1007/s11009-019-09744-2
    [17] D. G. Kendall, The genealogy of genealogy branching processes before (and after) 1873, Bull. London Math. Soc., 7 (1975), 225–253. https://doi.org/10.1112/blms/7.3.225 doi: 10.1112/blms/7.3.225
    [18] L. Arguin, D. Belius, A. J. Harper, Maxima of a randomized Riemann zeta function, and branching random walks, Ann. Appl. Probab., 27 (2017), 178–215. https://doi.org/10.1214/16-AAP1201 doi: 10.1214/16-AAP1201
    [19] L. Arguin, D. Belius, P. Bourgade, M. Radziwiłł, K. Soundararajan, Maximum of the Riemann zeta function on a short interval of the critical line, Commun. Pure Appl. Math., 72 (2019), 500–535. https://doi.org/10.1002/cpa.21791 doi: 10.1002/cpa.21791
    [20] E. Rawashdeh, A simple method for finding the inverse matrix of Vandermonde matrix, Mat. Vesn., 71 (2019), 207–213.
    [21] W. Rudin, Real and complex analysis, 3 Eds., New York: McGraw-Hill, 1987.
    [22] A. Grigutis, J. Jaunkauskas, J. Šiaulys, Multi seasonal discrete time risk model revisited, 2022, In press. https://doi.org/10.48550/arXiv.2207.03196
    [23] R. A. Horn, C. R. Johnson, Topics in matrix analysis, Cambridge University Press, 1991. https://doi.org/10.1017/CBO9780511840371
    [24] J. Kiefer, J. Wolfowitz, On the characteristics of the general queueing process, with applications to random walk, Ann. Math. Statist., 27 (1956), 147–161. https://doi.org/10.1214/aoms/1177728354 doi: 10.1214/aoms/1177728354
    [25] Mathematica (Version 9.0), Champaign, Illinois: Wolfram Research, Inc., 2012. Available from: https://www.wolfram.com/mathematica.
    [26] O. Navickienė, J. Sprindys, J. Šiaulys, Ruin probability for the bi-seasonal discrete time risk model with dependent claims, Modern Stoch. Theory Appl., 6 (2018), 133–144. https://doi.org/10.15559/18-VMSTA118 doi: 10.15559/18-VMSTA118
    [27] Y. Miao, K. P. Sendova, B. L. Jones, On a risk model with dual seasonalities, N. Am. Actuar. J., 2022. https://doi.org/10.1080/10920277.2022.2068611 doi: 10.1080/10920277.2022.2068611
    [28] A. Alencenovič, A. Grigutis, Bi-seasonal discrete time risk model with income rate two, Commun. Statist. Theory Methods, 2022. https://doi.org/10.1080/03610926.2022.2026962 doi: 10.1080/03610926.2022.2026962
    [29] A. Grigutis, A. Nakliuda, Note on the bi-risk discrete time risk model with income rate two, Modern Stoch. Theory Appl., 9 (2022), 401–412. https://doi.org/10.15559/22-VMSTA209 doi: 10.15559/22-VMSTA209
  • This article has been cited by:

    1. David J. Santana, Luis Rincón, Ruin probabilities as functions of the roots of a polynomial, 2023, 2351-6046, 1, 10.15559/23-VMSTA226
    2. Ernesto Cruz, Luis Rincón, David J. Santana, Ruin Probabilities as Recurrence Sequences in a Discrete-Time Risk Process, 2024, 26, 1387-5841, 10.1007/s11009-024-10102-0
    3. Andrius Grigutis, Artur Nakliuda, The limit law of maximum of discrete partial-sums distribution, 2024, 0363-1672, 10.1007/s10986-024-09650-2
    4. Simonas Gervė, Andrius Grigutis, Distribution of shifted discrete random walk generated by distinct random variables and applications in ruin theory, 2024, 2351-6046, 323, 10.15559/24-VMSTA249
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1743) PDF downloads(137) Cited by(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog