Next Article in Journal
Improving the Computational Efficiency of a Variant of Steffensen’s Method for Nonlinear Equations
Next Article in Special Issue
On the Rate of Convergence for a Characteristic of Multidimensional Birth-Death Process
Previous Article in Journal
Role of Media and Effects of Infodemics and Escapes in the Spatial Spread of Epidemics: A Stochastic Multi-Region Model with Optimal Control Approach
Previous Article in Special Issue
Monte Carlo Algorithms for the Parabolic Cauchy Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating the Expected Discounted Penalty Function in a Compound Poisson Insurance Risk Model with Mixed Premium Income

1
College of Mathematics and Statistics, Chongqing University, Chongqing 401331, China
2
School of Insurance, Shandong University of Finance and Economics, Jinan 250014, China
3
School of Science, Shandong Jiaotong University, Jinan 250357, China
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(3), 305; https://doi.org/10.3390/math7030305
Submission received: 24 February 2019 / Revised: 15 March 2019 / Accepted: 18 March 2019 / Published: 26 March 2019
(This article belongs to the Special Issue Stochastic Processes: Theory and Applications)

Abstract

:
In this paper, we consider an insurance risk model with mixed premium income, in which both constant premium income and stochastic premium income are considered. We assume that the stochastic premium income process follows a compound Poisson process and the premium sizes are exponentially distributed. A new method for estimating the expected discounted penalty function by Fourier-cosine series expansion is proposed. We show that the estimation is easily computed, and it has a fast convergence rate. Some numerical examples are also provided to show the good properties of the estimation when the sample size is finite.

1. Introduction

In this paper, we consider an insurance risk model with mixed premium income defined by
U ( t ) = u + c t + i = 1 M ( t ) Y i j = 1 N ( t ) X j , t 0 ,
where u 0 is the initial surplus, c 0 is the constant premium rate, and U ( t ) denotes the surplus level of an insurance company at time t. The premium number process M ( t ) and the claim number process N ( t ) are homogenous Poisson processes with intensity μ > 0 and λ > 0 , respectively. The individual claim sizes, X 1 , X 2 , , are positive independent and identically distributed (i.i.d) continuous random variables with density function f. The premium sizes, Y 1 , Y 2 , are positive i.i.d. continuous random variables with exponential distribution function g ( y ) = β e β y , y , β > 0 , where β is unknown. Throughout this paper, we assume that { M ( t ) } t 0 , { N ( t ) } t 0 , { X i } i 1 and { Y j } j 1 are mutually independent.
Whenever the surplus process becomes negative, we say that ruin occurs. The ruin time is defined by
τ = inf { t 0 : U ( t ) < 0 } ,
if for all t 0 , U ( t ) > 0 , we denote τ = . To avoid that ruin is a certain event, suppose that the following condition holds throughout this paper.
Net profit condition
c + μ β λ E [ X ] > 0 .
The above condition guarantees that the expectation of the surplus process will always be positive at any time t > 0 .
Let δ 0 be the interest force, defining the expected discounted penalty function by
Φ ( u ) = E [ e δ τ w ( U ( τ ) , | U ( τ ) | ) 1 ( τ < ) | U ( 0 ) = u ] , u 0 ,
where, w : [ 0 , ) × [ 0 , ) [ 0 , ) is a measurable penalty function of the surplus before ruin and the deficit at ruin, and 1 A is an indicator function of the event A. This function was first introduced by Gerber and Shiu [1], and is called Gerber-Shiu function in the literature. It has become an important and standard risk measure in ruin theory since various quantities of interests in ruin theory can be obtained for different values of the discount factor δ and different penalty functions w . For recent research progress on the Gerber-Shiu function, we can refer to work by Lin et al. [2], Yuen et al. [3,4], Zhao and Yin [5], Chi [6], Yin and Wang [7], Chi and Lin [8], Shen et al. [9], Yin and Yuen [10], Zhao and Yao [11], Li et al. [12,13], Dong et al. [14], and Wang and Zhang [15], among others.
For the classical insurance risk model, the premium rate is usually set to a constant. Extensive research has been done on the ruin measures of the model and its extension under the constant premium rate, such as ruin probability, Gerber-Shiu function (see, e.g., Gerber and Shiu [1], Willmot and Dickson [16], Wang et al. [17], Yin and Yuen [18], Dong and Yin [19], Yu [20], Yin et al. [21], Zeng et al. [22], Zhao et al. [23] and Yu et al. [24]). However, in the actual insurance business, the premium income of an insurance company, especially a small one, is sometimes volatile and random. In view of this situation, Boucherie et al. [25] first extended the classical risk model to the risk model with stochastic premium income by replacing the constant premium income with a compound Poisson process. Since then, the risk model with stochastic premium income has been widely studied by many scholars. The non-ruin probability and ruin probability were, respectively, studied by Boikov [26] and Temnov [27]. Bao [28], Bao and Ye [29] and Yang and Zhang [30] studied the Gerber-Shiu function in the classical risk model, delayed renewal risk model and Sparre Andersen risk model by assuming the premium process are Poisson process, respectively. Supposing premium and claim process follow compound Poisson processes, a defective renewal equation satisfied by the Gerber-Shiu function was established by Labbé and Sendova [31]. Yang and Zhang [32] further studied the above model by assuming that there exists a specific dependence structure among the claim sizes, inter-claim times and premium sizes and the individual premium sizes are exponentially distributed. Zhao and Yin [33] considered a renewal risk model where the premiums follow a compound Poisson risk process and the claims follow a generalized Erlang(n) process. The Laplace transform and a defective renewal equation for the Gerber-Shiu function are derived when the premium sizes are exponentially distributed. For more study on the risk model with stochastic premium income, the interested readers are referred to the work by Xu et al. [34], Yu [35,36], Zhou et al. [37], Gao and Wu [38], Zhou et al. [39], Deng et al. [40] and Zeng et al. [41] and the references therein.
The above-mentioned papers assume that some probability characteristics of the surplus process are known, however, which are usually unknown for an insurance company. Actually, some data information on the surplus process, such as surplus levels, claim and premium numbers, and claim and premium sizes, can be obtained. Thus, some semi-parametric and parametric estimations of the ruin probability and Gerber-Shiu function for different risk models are presented in recent literature (e.g., Politis [42], Wang and Yin [43], Yuen and Yin [44], Zhang et al. [45], Wang et al. [46], Zhang [47], Zhang and Yang [48,49] and Shimizu and Zhang [50], Peng and Wang [51,52], Yang et al. [53,54]). Besides, some effective methods, such as Laplace transform, Fourier-Sinc series expansion and Laguerre series expansion, have been applied in estimating the ruin probability and Gerber-Shiu function based on the observed data information. We refer the interested readers to Shimizu [55], Zhang [56], Zhang and Su [57,58] and Su et al. [59], among others.
In recent years, Fang and Oosterlee [60] proposed a novel method for pricing European options by Fourier-cosine series expansion, which is also called the COS method in the literature. It can be easily applied to approximate an integrable function as long as the corresponding Fourier transform has closed-form expression. Now, the COS method has been widely used for pricing options and other financial derivatives (see, e.g., Fang and Oosterlee [61] and Zhang and Oosterlee [62]). Recently, the COS method has also been used in risk theory to compute and estimate some risk measures by some actuarial researchers. For example, Chau et al. [63,64] used the COS method to compute the ruin probability and Gerber-Shiu function in the Lévy risk models. Zhang [65] applied the COS method to compute the density of the time to ruin in the classical risk model. Yang et al. [66] used a two-dimensional COS method to estimate the discounted density function of the deficit at ruin in the classical risk model with stochastic income where the premiums are only described by a compound Poisson process and the premium sizes are exponential distributed. Inspired by Yang et al. [66], in this paper, we extend the risk model of Yang et al. [66] by considering that there is also constant rate premium income process. Then, we use the COS method to estimate the expected discounted penalty function in this risk model with mixed premiums income process. The remainder of this paper is organized as follows. In Section 2, we first briefly introduce the Fourier-cosine series expansion method, and then derive the Fourier transform of the expected discounted penalty function. In Section 3, an estimator of the expected discounted penalty function is proposed by the observed sample of the surplus process. The consistent property is studied in Section 4 under large sample size setting. Finally, in Section 5, we present some simulation results to show that the estimator behaves well under finite sample size setting.

2. Preliminaries on Expected Discounted Penalty Function

2.1. Fourier-Cosine Series Expansion

In this subsection, we present some known results on the Fourier-cosine series expansion method. Let L 1 ( R + ) denote the class of integrable functions on the positive axis, and let F f and L f denote the Fourier transform and Laplace transform of a f L 1 ( R + ) , respectively. For any complex number z, we denote its real part and imaginary part by R e ( z ) and I m ( z ) , respectively.
It is known that an integrable function f defined on [ a 1 , a 2 ] has the following cosine series expansion,
f ( x ) = k = 0 2 a 2 a 1 a 1 a 2 f ( x ) cos k π x a 1 a 2 a 1 d x cos k π x a 1 a 2 a 1 ,
where means the first term of the summation has half weight. For a function f L 1 ( R + ) , we introduce an auxiliary function
f a ( x ) = f ( x ) · 1 ( 0 x a ) , a > 0 ,
then f a has a finite domain [ 0 , a ] and f ( x ) = f a ( x ) when x [ 0 , a ] . By Equation (3), we have
f ( x ) = f a ( x ) = k = 0 2 a 0 a f ( x ) cos k π x a d x cos k π x a , 0 x a .
For large a, we have
2 a 0 a f ( x ) cos k π x a d x = 2 a R e 0 a f ( x ) e i k π a x d x 2 a R e F f k π a ,
thus
f ( x ) k = 0 2 a R e F f k π a cos k π x a , 0 x a .
Furthermore, for a large integer K, the above summation can be truncated as follows:
f ( x ) k = 0 K 1 2 a R e F f k π a cos k π x a , 0 x a .
As for the expected discounted penalty function, it follows from Equation (5) that for 0 u a ,
Φ ( u ) Φ K , a ( u ) : = k = 0 K 1 2 a R e F Φ k π a cos k π x a .
It can be easily seen that the key of approximating the expected discounted penalty function by Fourier-cosine series expansion method is to calculate the Fourier transform F Φ ( s ) for s { k π a : k = 0 , 1 , , K 1 } . Therefore, we derive a specific expression of Fourier transform F Φ ( · ) in the next section.

2.2. The Fourier Transform of Expected Discounted Penalty Function

In this subsection, we derive the Fourier transform of the expected discounted penalty function. For convenience, we introduce the Dickson-Hipp operator T s (see, e.g., Dickson and Hipp [67] and Li and Garrido [68]), which for any integrable function f on ( 0 , ) and any complex number s with R e ( s ) 0 is defined as
T s f ( y ) = y e s ( x y ) f ( x ) d x = 0 e s x f ( x + y ) d x , y 0 .
The Dickson-Hipp operator has been widely used in ruin theory to simplify the expression of ruin related functions. For properties on this operator, we refer the interested readers to Li and Garrido [68].
We call the following equation (with respect to s) the Lundberg’s fundamental equation for the risk model defined in Equation (1),
( β s ) s δ + μ + λ c + μ c · β β s + λ c L f ( s ) = 0 .
By Lemma 2.1 in Zhao and Yin [33], we known that for δ 0 and c > 0 , the above equation has exactly two nonnegative roots, denoted as ρ 1 , ρ 2 in this literature, and one of the roots ρ 1 equals 0 when δ = 0 .
Remark 1.
Let
χ ( s ) = δ + μ + λ c s λ L f ( s ) + μ β s β , s 0 .
It is obvious that χ ( s ) = 0 also has exactly two nonnegative roots ρ 1 , ρ 2 . Under net profit condition, we have
χ ( s ) = c + λ 0 x e s x f ( x ) d x μ β ( s β ) 2 c λ E [ X ] μ β < 0 , 0 s < β ,
which implies that χ ( s ) is decreasing on [ 0 , β ) . In addition, χ ( s ) is continuous on [ 0 , β ) , χ ( 0 ) = δ 0 and χ ( β 0 ) = . Therefore, we conclude that the root ρ 1 of equation χ ( s ) = 0 is located in the interval [ 0 , β ) , and, in particular, ρ 1 = 0 when δ = 0 ; then, ρ 2 is located in the interval ( β , ) .
By Theorem 3.1 in Zhao and Yin [33], when stochastic premium income follows a compound Poisson process, we known that the expected discounted penalty function satisfies the following renewal equation
Φ ( u ) = 0 u Φ ( u x ) H ( x ) d x + K ( u ) ,
where
H ( x ) = λ c β ρ 1 ρ 1 ρ 2 T ρ 1 f ( x ) + β ρ 2 ρ 2 ρ 1 T ρ 2 f ( x ) , K ( u ) = λ c β ρ 1 ρ 1 ρ 2 T ρ 1 ω ( u ) + β ρ 2 ρ 2 ρ 1 T ρ 2 ω ( u ) , ω ( u ) = u w ( u , x u ) f ( x ) d x .
To derive the Fourier transform of Φ ( u ) , we assume the following condition holds true.
Condition 1. The penalty function w satisfies
0 0 ( 1 + x ) w ( x , y ) f ( x + y ) d y d x < .
This condition guarantees Φ L 1 ( R + ) .
Now, we compute the Fourier transform of the expected discounted penalty function. Applying the Fourier transform on both sides of Equation (8) gives
F Φ ( s ) = 0 e i s u Φ ( u ) d u = 0 e i s u 0 u Φ ( u x ) H ( x ) d x d u + 0 e i s u K ( u ) d u = 0 0 e i s ( z + x ) Φ ( z ) d z d x + F K ( s ) = F Φ ( s ) F H ( s ) + F K ( s ) ,
leading to
F Φ ( s ) = F K ( s ) 1 F H ( s ) .
For the Fourier transform F H ( s ) , we have
F H ( s ) = 0 e i s u H ( x ) d x = λ c β ρ 1 ρ 1 ρ 2 0 e i s x T ρ 1 f ( x ) d x + β ρ 2 ρ 2 ρ 1 0 e i s x T ρ 2 f ( x ) d x ,
where
0 e i s u T ρ j f ( x ) d x = 0 e i s x x e ρ j ( y x ) f ( y ) d y d x = 0 x e ( i s + ρ j ) x · e ρ j y f ( y ) d y d x = 1 ρ j + i s 0 e i s y f ( y ) d y 0 e ρ i y f ( y ) d y = 1 ρ j + i s [ F f ( s ) L f ( ρ j ) ] , j = 1 , 2 .
Then, we obtain
F H ( s ) = λ c ( ρ 1 ρ 2 ) ( β ρ 2 ) [ F f ( s ) L f ( ρ 2 ) ] ρ 2 + i s ( β ρ 1 ) [ F f ( s ) L f ( ρ 1 ) ] ρ 1 + i s .
Similarly, for F K ( s ) , we obtain
F K ( s ) = λ c ( ρ 1 ρ 2 ) ( β ρ 2 ) [ F ω ( s ) L ω ( ρ 2 ) ] ρ 2 + i s ( β ρ 1 ) [ F ω ( s ) L ω ( ρ 1 ) ] ρ 1 + i s .

3. Estimation Procedure

In this section, we study how to estimate the expected discounted penalty function by Fourier-cosine series expansion based on the discretely observed information of the surplus process, the aggregate claims and premiums processes. According to Equation (6), we know that the key is to construct an estimation of F Φ ( s ) based on these discrete information.
Assume that we can observe the surplus process over a long time interval [ 0 , T ] . Let Δ > 0 be a fixed inter-observation interval. Furthermore, without loss of generality, we assume T / Δ is an integer denoted as n.
Suppose that the insurer can get the following datasets.
(1)
Dataset of surplus levels:
{ U j Δ : j = 0 , 1 , 2 , , n } ,
where U j Δ is the observed surplus level at time t = j Δ .
(2)
Dataset of claim numbers and claim sizes:
{ N j Δ , X 1 , X 2 , , X N j Δ } , j = 1 , , n ,
where N j Δ is the total claim number up to time t = j Δ .
(3)
Dataset of premium numbers and claim sizes:
{ M j Δ , Y 1 , Y 2 , , Y N j Δ } , j = 1 , , n ,
where M j Δ is the total premium number up to time t = j Δ .
Next, we study how to estimate the Fourier transform F Φ ( s ) based on the above datasets. To estimate F Φ ( s ) by Equations (9)–(11), we should first estimate the following characteristics:
λ , μ , β , F f , F ω , ρ j , L f ( ρ j ) , L ω ( ρ j ) , j = 1 , 2 .
First, we can estimate F f ( s ) by the empirical characteristic function
F f ^ ( s ) = 1 N T j = 1 N T e i s X j .
Similarly, L f ( s ) can be estimated by
L f ^ ( s ) = 1 N T j = 1 N T e s X j .
Next, for the function ω ( u ) , we have
F ω ( s ) = E 0 X e i s u w ( u , X u ) d u , L ω ( s ) = E 0 X e s u w ( u , X u ) d u .
Then, F ω ( s ) and L ω ( s ) can be respectively estimated by
F ω ^ ( s ) = 1 N T j = 1 N T 0 X j e i s u w ( u , X j u ) d u , L ω ^ ( s ) = 1 N T j = 1 N T 0 X j e s u w ( u , X j u ) d u .
According to the property of Poisson distribution, λ and μ can be estimated by
λ ^ = 1 T N T , μ ^ = 1 T M T .
It is easily seen that
λ ^ λ = O p ( T 1 2 ) , μ ^ μ = O p ( T 1 2 ) .
Since the premium size Y follows exponential distribution with parameter β , we have E [ Y ] = 1 β ; then, we can estimate β by
β ^ = 1 1 M T i = 1 M T Y i .
It is also easily seen that β ^ β = O p ( T 1 2 ) . The estimation of ρ j , j = 1 , 2 , denoted as ρ j ^ , j = 1 , 2 , are defined to be the nonnegative roots of the following equation (in s),
δ + μ ^ + λ ^ c s λ ^ L f ^ ( s ) + μ ^ β ^ s β ^ = 0 .
Furthermore, we can, respectively, estimate L f ( ρ j ) , L ω ( ρ j ) , j = 1 , 2 by L f ^ ( ρ j ^ ) , L ω ^ ( ρ j ^ ) , j = 1 , 2 .
Remark 2.
Let
χ ^ ( s ) = δ + μ ^ + λ ^ c s λ ^ L f ^ ( s ) + μ ^ β ^ s β ^ , s 0 .
By the similar arguments of the Lunderg’s fundamental equation in Zhao and Yin [33], we can obtain that equation χ ^ ( s ) = 0 has exactly two nonnegative roots, denoted as ρ ^ 1 , ρ ^ 2 in this literature. It is clear that χ ^ ( 0 ) = δ 0 and χ ^ ( β ^ 0 ) = . Under net profit condition, we have
χ ^ ( s ) = c + λ ^ 1 N T j = 1 N T X j e s X j μ ^ β ^ ( s β ^ ) 2 c + λ ^ 1 N T j = 1 N T X j μ ^ β ^ a . s . c + λ E [ X ] μ β < 0 , 0 s < β ^ .
Therefore, we obtain that the probability that χ ^ ( s ) = 0 has unique root ρ 1 ^ on [ 0 , β ^ ) tends to one as T . We set ρ 1 ^ = 0 when δ = 0 since ρ 1 = 0 in this case. Then, the other root ρ 2 ^ is located in the interval ( β ^ , ) with probability tends to one as T .
Proposition 1.
Suppose that net profit condition holds true, then we have ρ 1 ^ p ρ 1 and ρ 2 ^ p ρ 2 .
Proof of Proposition 1.
By Remark 2, when 0 s < β ^ , χ ^ ( s ) is continuous and the probability that χ ^ ( s ) = 0 has unique nonnegative root tends to one as T ; and when s > β ^ , χ ^ ( s ) is continuous and the probability that χ ^ ( s ) = 0 has unique positive root tends to one as T . By Remark 1, it is easily seen that, for every ε > 0 , χ ( ρ 1 ε ) > 0 > χ ( ρ 1 + ε ) and χ ( ρ 2 ε ) > 0 > χ ( ρ 2 + ε ) . Besides, we find that for any s > 0 , χ ^ ( s ) p χ ( s ) . Thus, it follows from Lemma 5.10 in Van Der Vaart [69] that ρ 1 ^ p ρ 1 and ρ 2 ^ p ρ 2 . ☐
Once we have obtained the estimation of the above characteristics, by Equations (9)–(11), the estimation of Fourier transform F Φ ( s ) , denoted as F Φ ^ ( s ) , can be defined as
F Φ ^ ( s ) = F K ^ ( s ) 1 F H ^ ( s ) ,
where
F K ^ ( s ) = λ ^ c ( ρ 1 ^ ρ 2 ^ ) ( β ^ ρ 2 ^ ) [ F ω ^ ( s ) L ω ^ ( ρ 2 ^ ) ] ρ 2 ^ + i s ( β ^ ρ 1 ^ ) [ F ω ^ ( s ) L ω ^ ( ρ 1 ^ ) ] ρ 1 ^ + i s , F H ^ ( s ) = λ ^ c ( ρ 1 ^ ρ 2 ^ ) ( β ^ ρ 2 ^ ) [ F f ^ ( s ) L f ^ ( ρ 2 ^ ) ] ρ 2 ^ + i s ( β ^ ρ 1 ^ ) [ F f ^ ( s ) L f ^ ( ρ 1 ^ ) ] ρ 1 ^ + i s .
Finally, replacing F Φ ( · ) in Equation (6) by the estimation F Φ ^ ( · ) , the expected discounted penalty function can be estimated by
Φ ^ K , a ( u ) : = k = 0 K 1 2 a R e F Φ ^ k π a cos k π x a , 0 u a .

4. Consistency Properties

In this section, we study the asymptotic properties of the estimation Φ ^ K , a . Let L 2 ( R + ) denote the class of square integrable functions on the positive axis. For any function f L 2 ( R + ) , its L 2 -norm is defined by f = 0 f 2 ( x ) d x 1 2 . Throughout this section, C represents a positive generic constant that may take different values at different steps. In addition, we define
H j ( x ) = 0 x u j w ( u , x u ) d u , j = 0 , 1 , 2 .
It is easy to see that
0 u j ω ( u ) d u = E H j ( X ) .
For reader’s convenience, we introduce some definitions in empirical process theory, which are used to study the asymptotic properties. For any measurable function f, its L r ( P ) -norm is defined by f P , r = | f ( ω ) | r d P ( ω ) 1 r . Given two functions l and u, the bracket [ l , u ] is the set of all functions f with l f u . An ε -bracket in L r ( P ) is a bracket [ l , u ] with u l P , r < ε . For a class G L r ( P ) , the bracketing number N ( ε , G , L r ( P ) ) is the minimum number of ε -brackets needed to cover G . For β > 0 , the bracketing integral is defined by J β , G , L r ( P ) = 0 β N ( ε , G , L r ( P ) ) d ε .
We use the L 2 -norm to study the asymptotic properties of the estimator. The following condition is useful in our discussion, which ensures Φ L 2 ( R + ) .
Condition 2. For the penalty function w , there exist some integers α 1 , α 2 and constant C such that
w ( x , y ) C ( 1 + x ) α 1 ( 1 + y ) α 2 .
Put Φ K , a = Φ ^ K , a = 0 when u > a . By triangle inequality, we have
Φ Φ ^ K , a Φ Φ K , a + Φ K , a Φ ^ K , a ,
where the first term Φ Φ K , a is the bias caused by Fourier cosine series approximation, and the second term Φ K , a Φ ^ K , a is the bias caused by statistical estimation.
For the bias Φ Φ K , a , by similar arguments to those of Zhang [65], we obtain the following result.
Theorem 1.
Suppose that 0 | Φ ( u ) | d u < and for some integer m, Φ ( u ) C u ( m + 1 ) ; then, under net profit condition, Conditions 1 and 2, we have
Φ Φ K , a C K + 1 a 2 m + 1 + a K 1 .
Next, we study the error Φ K , a Φ ^ K , a . For ρ 1 ^ and ρ 2 ^ , we derive the following result.
Theorem 2.
Supposing that net profit condition holds, we have ρ 1 ^ ρ 1 = O p ( T 1 2 ) . Supposing that c > λ E [ X ] , we have ρ 2 ^ ρ 2 = O p ( T 1 2 ) .
Proof of Theorem 2.
By the mean value theorem,
χ ^ ( ρ j ) = χ ^ ( ρ j ^ ) + χ ^ ( ρ j ) ( ρ j ρ j ^ ) = χ ^ ( ρ j ) ( ρ j ρ j ^ ) , j = 1 , 2 ,
where ρ j ( j = 1 , 2 ) is a random number between ρ j ( j = 1 , 2 ) and ρ j ^ ( j = 1 , 2 ) . Since χ ( ρ j ) = 0 , j = 1 , 2 , we obtain
ρ j ^ ρ j = χ ( ρ j ) χ ^ ( ρ j ) χ ^ ( ρ j ) , j = 1 , 2 .
It is easily seen that χ ( ρ j ) χ ^ ( ρ j ) = O p ( T 1 2 ) , j = 1 , 2 for λ ^ λ = O p ( T 1 2 ) , μ ^ μ = O p ( T 1 2 ) , and β ^ β = O p ( T 1 2 ) .
For ρ 1 , under net profit condition, we introduce the following set:
A T , 1 = | χ ^ ( ρ 1 ) | > 1 2 c + μ β λ E [ X ] .
Since c + μ ^ β ^ λ ^ 1 N T j = 1 N T X j p c + μ β λ E [ X ] and | χ ^ ( ρ 1 ) | c + μ ^ β ^ λ ^ 1 N T j = 1 N T X j , we have
P ( A T , 1 ) P c + μ ^ β ^ λ ^ 1 N T j = 1 N T X j 1 2 ( c + μ β λ E [ X ] ) = 1 P c + μ ^ β ^ λ ^ 1 N T j = 1 N T X j < 1 2 ( c + μ β λ E [ X ] ) = 1 P c + μ β λ E [ X ] c μ ^ β ^ + λ ^ 1 N T j = 1 N T X j 1 2 ( c + μ β λ E [ X ] ) 1 , a s T .
For ρ 2 , we have
| χ ^ ( ρ 2 ) | c + μ ^ β ^ ( s β ^ ) 2 λ ^ 1 N T j = 1 N T X j c λ ^ 1 N T j = 1 N T X j .
Under condition c > λ E [ X ] , we introduce the following set:
A T , 2 = | χ ^ ( ρ 2 ) | > 1 2 c λ E [ X ] .
Similarly, we obtain that P ( A T , 2 ) 1 a s T in that c λ ^ 1 N T j = 1 N T X j p c λ E [ X ] .
Furthermore,
P | ρ j ^ ρ j | > C T 1 2 = P | χ ( ρ j ) χ ^ ( ρ j ) | | χ ^ ( ρ j ) | > C T 1 2 P | χ ( ρ j ) χ ^ ( ρ j ) | | χ ^ ( ρ j ) | > C T 1 2 A T , j + P ( A T , j c ) = P | χ ( ρ j ) χ ^ ( ρ j ) | > C T 1 2 + P ( A T , j c ) , j = 1 , 2 .
As a result, since P ( A T , 1 c ) 0 under net profit condition, P ( A T , 2 c ) 0 under c > λ E [ X ] , and χ ( ρ j ) χ ^ ( ρ j ) = O p ( T 1 2 ) , j = 1 , 2 , we derive desired results. ☐
The following two theorems give the uniform convergence rates of F H and F K .
Theorem 3.
Suppose that c > λ E [ X ] , H j ( X ) P , 1 < , j = 0 , 1 , and H j ( X ) P , 2 < , j = 1 , 2 . Then, for large a , K and T, we have
sup s [ 0 , K π / a ] F K ( s ) F K ^ ( s ) = O p log K a / T .
Proof of Theorem 3.
By Equations (11) and (12),
F K ^ ( s ) F K ( s ) = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s F ω ^ ( s ) L ω ^ ( ρ 2 ^ ) λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s F ω ( s ) L ω ( ρ 2 ) + λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 1 ^ ρ 1 ^ + i s F ω ^ ( s ) L ω ^ ( ρ 1 ^ ) λ c ( ρ 1 ρ 2 ) β ρ 1 ρ 1 + i s F ω ( s ) L ω ( ρ 1 ) : = I 1 + I 2 ,
Since I 1 and I 2 have similar formations, we only study I 1 in detail.
I 1 = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s F ω ^ ( s ) L ω ^ ( ρ 2 ^ ) λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s F ω ( s ) L ω ( ρ 2 ) = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T 0 X j ( e i s u e ρ 2 ^ u ) w ( u , X j u ) d u λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s E 0 X ( e i s u e ρ 2 u ) w ( u , X u ) d u = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T 0 X j ( e ρ 2 u e ρ 2 ^ u ) w ( u , X j u ) d u + λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T 0 X j ( e i s u e ρ 2 ^ u ) w ( u , X j u ) d u λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s E 0 X ( e i s u e ρ 2 u ) w ( u , X u ) d u = I I 1 + I I 2 .
where
I I 1 = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T 0 X j ( e ρ 2 u e ρ 2 ^ u ) w ( u , X j u ) d u , I I 2 = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T 0 X j ( e i s u e ρ 2 ^ u ) w ( u , X j u ) d u λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s E 0 X ( e i s u e ρ 2 u ) w ( u , X u ) d u .
For I I 1 , we have
| I I 1 | λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T 0 X j u | ρ 2 ρ 2 ^ | w ( u , X j u ) d u = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s | ρ 2 ρ 2 ^ | 1 N T j = 1 N T H 1 ( X j ) .
Since 1 N T j = 1 N T H 1 ( X j ) p H 1 ( X ) P , 1 < , ρ 2 ρ 2 ^ = O p ( T 1 2 ) , we have
| I I 1 | = O p ( T 1 2 ) .
For I I 2 we derive
I I 2 = λ ^ ( β ^ ρ 2 ^ ) c ( ρ 1 ^ ρ 2 ^ ) ρ 2 + i s ρ 2 ^ + i s 1 N T j = 1 N T 0 X j e i s u e ρ 2 ^ u ρ 2 + i s w ( u , X j u ) d u λ ( β ρ 2 ) c ( ρ 1 ρ 2 ) E 0 X e i s u e ρ 2 u ρ 2 + i s w ( u , X u ) d u = λ ^ ( β ^ ρ 2 ^ ) c ( ρ 1 ^ ρ 2 ^ ) ρ 2 + i s ρ 2 ^ + i s 1 N T j = 1 N T g 1 , s ( X j ) E [ g 1 , s ( X ) ] + ( ρ 2 + i s ) λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s λ c ( ρ 1 ρ 2 ) ( β ρ 2 ) ρ 2 + i s E [ g 1 , s ( X ) ] ,
where
g 1 , s ( x ) = 0 x e i s u e ρ 2 u ρ 2 + i s w ( u , x u ) d u = 0 x 0 u e i s ( u y ) ρ 2 y d y w ( u , x u ) d u , x 0 .
For g 1 , s ( x ) , we have
sup s [ 0 , K π / a ] E g 1 , s ( X ) 0 0 x 0 u e ρ 2 y d y w ( u , x u ) d u f ( x ) d x 0 0 x u w ( u , x u ) d u f ( x ) d x = H 1 ( X ) P , 1 < ,
then
sup s [ 0 , K π / a ] ( ρ 2 + i s ) λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s λ c ( ρ 1 ρ 2 ) ( β ρ 2 ) ρ 2 + i s E [ g 1 , s ( X ) ] ( ρ 2 + i s ) λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s λ c ( ρ 1 ρ 2 ) ( β ρ 2 ) ρ 2 + i s · H 1 ( X ) P , 1 = O p ( T 1 2 )
for λ ^ λ = O p ( T 1 2 ) , μ ^ μ = O p ( T 1 2 ) , ρ 1 ^ ρ 1 = O p ( T 1 2 ) and ρ 2 ^ ρ 2 = O p ( T 1 2 ) .
Let us introduce the following two types of real-valued functions,
G K , R = g : g = R e ( g 1 , s ) , s [ 0 , K π / a ] , G K , I = g : g = I m ( g 1 , s ) , s [ 0 , K π / a ] .
Then, we have
sup s [ 0 , K π / a ] 1 N T j = 1 N T g 1 , s ( X j ) E g 1 , s ( X ) sup g G K , R 1 N T j = 1 N T g 1 , s ( X j ) E g 1 , s ( X ) + sup g G K , I 1 N T j = 1 N T g 1 , s ( X j ) E g 1 , s ( X ) .
We only study the convergence rate of the first term sup g G K , R 1 N T j = 1 N T g 1 , s ( X j ) E g 1 , s ( X ) , since the second term follows similarly.
For any real-valued function g G K , R , we have
| g ( x ) | sup s [ 0 , K π / a ] 0 x 0 u e i s ( u y ) ρ 2 y d y w ( u , x u ) d u 0 x 0 u e ρ 2 y d y w ( u , x u ) d u 0 x u w ( u , x u ) d u = H 1 ( x ) ,
which implies that G K , R is contained in the single bracket [ H 1 , H 1 ] . For two functions g 1 , s 1 , g 1 , s 2 , where s j [ 0 , K π / a ] , j = 1 , 2 , the mean value theorem gives
R e ( g 1 , s 1 ) R e ( g 1 , s 2 ) = 0 x 0 u ( c o s ( s 1 ( u y ) ) c o s ( s 2 ( u y ) ) ) e ρ 2 y d y w ( u , x u ) d u = 0 x 0 u s i n ( s ( u y ) ) · ( u y ) · ( s 1 s 2 ) e ρ 2 y d y w ( u , x u ) d u 0 x 0 u ( u y ) e ρ 2 y d y w ( u , x u ) d u · | s 1 s 2 | H 2 ( x ) | s 1 s 2 | ,
where s is a number between s 1 and s 2 . Under the condition H 2 ( X ) P , 2 < , it follows from Example 19.7 in Van Der Vaart [69] that, for any 0 < ε < K π / a , there exists a constant C such that the bracket number for G K , R satisfies
N ( ε , G K , R , L 2 ( P ) ) C K π ε a H 2 ( x ) P , 2 .
As a result, for every δ > 0 , the bracketing integral
J ( δ , G K , R , L 2 ( P ) ) 0 δ log C K π ε a H 2 ( x ) P , 2 d ε log K a .
Furthermore, by Corollary 19.35 in Van Der Vaart [69], we have
E 1 N T sup g G K , R j = 1 N T g ( X j ) E g ( X ) | N T J ( δ , G K , R , L 2 ( P ) ) 0 δ log C K π ε a H 2 ( x ) P , 2 d ε log K a ,
then
E 1 N T sup g G K , R j = 1 N T g ( X j ) E g ( X ) = E 1 N T E 1 N T sup g G K , R j = 1 N T g ( X j ) E g ( X ) | N T log K a / E [ N T ] log K a / E [ N T ] = log K a / λ T .
Therefore,
sup g G K , R 1 N T j = 1 N T g ( X j ) E g ( X ) = O p log K a / T .
Similarly,
sup g G K , I 1 N T j = 1 N T g ( X j ) E g ( X ) = O p log K a / T .
Combining Equations (18)–(20), we obtain
sup s [ 0 , K π / a ] 1 N T j = 1 N T g 1 , s ( X j ) E g 1 , s ( X ) = O p log K a / T .
As a result, Equations (17) and (21) give
sup s [ 0 , K π / a ] | I I 2 | = O p log K a / T .
By Equations (16) and (22), we have
sup s [ 0 , K π / a ] | I 1 | = O p log K a / T ;
then by the similar arguments of I 1 , we have
sup s [ 0 , K π / a ] | I 2 | = O p log K a / T .
Finally, we can derive the desired result by Equations (23) and (24). ☐
Theorem 4.
Suppose that c > λ E [ X ] , E [ X k ] < , k = 1 , 2 . Then, for large a , K and T, we have
sup s [ 0 , K π / a ] F H ( s ) F H ^ ( s ) = O p log K a / T .
Proof of Theorem 4.
By Equations (11) and (12),
F H ^ ( s ) F H ( s ) = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s F f ^ ( s ) L f ^ ( ρ 2 ^ ) λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s F f ( s ) L f ( ρ 2 ) + λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 1 ^ ρ 1 ^ + i s F f ^ ( s ) L f ^ ( ρ 1 ^ ) λ c ( ρ 1 ρ 2 ) β ρ 1 ρ 1 + i s F f ( s ) L f ( ρ 1 ) : = l 1 + l 2 ,
Since l 1 and l 2 have similar formations, we only study l 1 in detail. For l 1 , we have
l 1 = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T e i s X j e ρ 2 ^ X j λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s E e i s X e ρ 2 X = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T e ρ 2 X j e ρ 2 ^ X j + λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T e i s X j e ρ 2 X j λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s E e i s X e ρ 2 X : = l l 1 + l l 2 ,
where
l l 1 = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T e ρ 2 X j e ρ 2 ^ X j , l l 2 = λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s 1 N T j = 1 N T e i s X j e ρ 2 X j λ c ( ρ 1 ρ 2 ) β ρ 2 ρ 2 + i s E e i s X e ρ 2 X .
For l l 1 , we have
| l l 1 | λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s | ρ 2 ρ 2 ^ | 1 N T j = 1 N T X j .
Since 1 N T j = 1 N T X j p E [ X ] < , ρ 2 ρ 2 ^ = O p ( T 1 2 ) , we have
| l l 1 | = O p ( T 1 2 ) .
For l l 2 , we obtain
l l 2 = λ ^ ( β ^ ρ 2 ^ ) c ( ρ 1 ^ ρ 2 ^ ) ρ 2 + i s ρ 2 ^ + i s 1 N T j = 1 N T g 2 , s ( X j ) E [ g 2 , s ( X ) ] + ( ρ 2 + i s ) λ ^ c ( ρ 1 ^ ρ 2 ^ ) β ^ ρ 2 ^ ρ 2 ^ + i s λ c ( ρ 1 ρ 2 ) ( β ρ 2 ) ρ 2 + i s E [ g 2 , s ( X ) ] ,
where
g 2 , s ( x ) = e i s x , x 0 .
For g 2 , s ( x ) , we have
sup s [ 0 , K π / a ] | E [ g 2 , s ( X ) ] | = sup s [ 0 , K π / a ] 0 e i s x f ( x ) d x 1 < ,
Under condition E [ X k ] < , k = 1 , 2 , by similar arguments of g 1 , s , we conclude that
sup s [ 0 , K π / a ] 1 N T j = 1 N T g 2 , s ( X j ) E g 2 , s ( X ) = O p log K a / T .
In addition, it follows from a similar analysis of I I 2 in the proof of Theorem 3 that
sup s [ 0 , K π / a ] | l l 2 | = O p log K a / T .
Then, Equations (26) and (28) give
sup s [ 0 , K π / a ] | l 1 | sup s [ 0 , K π / a ] | l l 2 | + sup s [ 0 , K π / a ] | l l 2 | = O p log K a / T .
Similarly, we derive
sup s [ 0 , K π / a ] | l 2 | = O p log K a / T .
Finally, we can derive the desired result by Equations (29) and (30). ☐
Based on the above conclusions, we have the following result.
Theorem 5.
Suppose that c > λ E [ X ] , E [ X k ] < , k = 1 , 2 , H j ( X ) P , 1 < , j = 0 , 1 , and H j ( X ) P , 2 < , j = 1 , 2 . Then, for large a , K and T, we have
Φ K , a Φ ^ K , a 2 = O p K a log K a / T .
Proof of Theorem 5.
First, we have
Φ K , a Φ ^ K , a 2 = 0 a Φ K , a ( u ) Φ ^ K , a ( u ) 2 d u k = 0 K 1 4 a 2 Re F Φ k π a F Φ ^ k π a 2 0 a cos k π a u 2 d u = 2 a k = 0 K 1 F Φ k π a F Φ ^ k π a 2 2 K a sup s 0 , K π / a F Φ s F Φ ^ s 2 .
Then, by Equations (9) and (12), we obtain
sup s [ 0 , K π / a ] F Φ ( s ) F Φ ^ ( s ) = sup s [ 0 , K π / a ] F K ( s ) 1 F H ( s ) F K ^ ( s ) 1 F H ^ ( s ) = sup s [ 0 , K π / a ] 1 F H ( s ) F K ( s ) F K ^ ( s ) + F K ( s ) F H ( s ) F H ^ ( s ) 1 F H ( s ) 1 F H ^ ( s ) .
Combining Theorems 3 and 4 gives
sup s [ 0 , K π / a ] | F Φ ( s ) F Φ ^ ( s ) | = O p log K a / T .
Finally, by Equation (31), we derive that
Φ K , a Φ ^ K , a 2 2 K a sup s [ 0 , K π / a ] | F Φ ( s ) F Φ ^ ( s ) | 2 = O p K a log K a / T .
This completes the proof. ☐
Combing Theorems 1 and 5, we finally obtain the following convergence rate:
Φ Φ ^ K , a 2 = O K + 1 a 2 m + 1 + O a K 1 + O p K a log K a / T .
We first find the optimal truncation parameter a = O K ( n + 1 ) to minimize the convergence rate O K + 1 a 2 m + 1 + O a K 1 . Replacing a by a in Equation (33) gives
Φ Φ ^ K , a 2 = O K m m + 1 + O p K m m + 1 log K / T .
In addition, we find the optimal truncation K = O T m + 1 2 m . Thus, we obtain the smallest convergence rate Φ Φ ^ K , a 2 = O p T 2 m 2 ( m + 1 ) 2 .

5. Simulation Studies

In this section, we present some numerical examples to explain the excellent properties of our estimator when the observed sample is finite. We mainly studied the following three kinds of special expected discounted penalty function:
(1) Ruin probability (RP: w 1 , δ = 0 ) ;
(2) Laplace transform of ruin time (LT: w 1 , δ = 0 . 1 ) ;
(3) Expected discounted deficit at ruin (EDD) when ruin is due to a claim ( w ( x , y ) = y , δ = 0 . 1 ) .
Some simulation examples of the above functions are presented for different distributions of claim sizes:
(1) Exponential distribution: f ( x ) = e x , x > 0 ;
(2) Erlang(2) distribution: f ( x ) = 4 x e 2 x , x > 0 ;
(3) Combined exponential distribution: f ( x ) = 3 e 1 . 5 x 3 e 3 x , x > 0 ;
(4) Mixed exponential distribution: f ( x ) = 2 3 e 2 x + 2 3 e x , x > 0 .
We set Δ = 1 , λ = 2 , and μ = 5 , where Δ = 1 can be explained as one week. That is to say, in a long time interval, the insurer will observe the data once a week, the expected claim number is 2 times per week, and the expected premium number is 5 times per week. Furthermore, since there are 52 business weeks every year, we assumed that T = 1 4 × 52 × Δ , T = 1 2 × 52 × Δ , T = 1 × 52 × Δ , T = 5 × 52 × Δ , which means that we observed the surplus process for one quarter, half a year, one year and five years. Then, we used Equation (13) to estimate the above Gerber-Shiu functions, and the corresponding true value obtained by Laplace inversion. In all simulations, we set c = 5 , a = 30 , K = 2 13 , and we carried out the relevant analysis based on 300 simulation experiments. We first introduce several concepts used in this section. The mean value and the mean relative error, which are, respectively, defined by
1 300 j = 1 300 Φ ^ K , a , j ( u ) , 1 300 j = 1 300 Φ ^ K , a , j ( u ) Φ ( u ) 1 ,
and the integrated mean square error (IMSE), which is defined by
1 300 j = 1 300 0 ( Φ ^ K , a , j ( u ) Φ ( u ) ) 2 d u 1 300 j = 1 300 0 30 ( Φ ^ K , a , j ( u ) Φ ( u ) ) 2 d u ,
where Φ ^ K , a , j ( u ) is the estimate of expected discounted penalty function in the jth experiment. For IMSE, we computed the integral on the finite domain [ 0 , 30 ] , as both the true value and the estimator will be very small when u 30 . In reality, both the true value and the estimated value are very close to zero when u > 20 , thus we present all images for u [ 0 , 20 ] to illustrate the performance of our estimators better.
First, we plot the mean curves of the estimated expected discounted penalty functions and compare them with the corresponding true curves. For the above-mentioned four distributions of claim sizes, we show the mean curves and the true curves of RP, LT, EDD due to a claim in Figure 1, Figure 2 and Figure 3, respectively. It is easily observed that, even though we used the quarter book data, the performance of the estimation for each claim distribution was still good. We could hardly distinguish the true curves from the mean curves when we used the five-year book data. Since it is difficult to distinguish the mean curves for different T values when u becomes large, we further show the mean relative error curves of the estimated expected discounted penalty functions for different claim distributions in Figure 4, Figure 5 and Figure 6, respectively. We observed that the mean relative errors became small as T increased, and they were very small when T = 260 . Besides, we found that the mean relative errors were small when u was small, but became very large when u was large. This is because the true values of the expected discounted penalty functions were very small when u was large. All of the above results illustrate the performance of our estimations by images. Finally, we present some values of IMSE for the estimators of RP, LT, and EDD in Table 1 and Table 2 to further show that our estimators perform well. Combining all of the simulation results, we conclude that our estimators could effectively approximate the true values, and they performed well for the large T.

Author Contributions

Data curation, Y.W.; Formal analysis, W.Y., Y.H. and X.Y.; Methodology, W.Y. and H.F.; and Writing—original draft, W.Y.

Funding

This research was partially supported by the National Natural Science Foundation of China (Grant Nos. 11301303 and 71804090), the National Social Science Foundation of China (Grant No. 15BJY007), the Taishan Scholars Program of Shandong Province (Grant No. tsqn20161041), the Humanities and Social Sciences Project of the Ministry Education of China (Grant Nos. 19YJA910002 and 16YJC630070), the Natural Science Foundation of Shandong Province (Grant No. ZR2018MG002), the Fostering Project of Dominant Discipline and Talent Team of Shandong Province Higher Education Institutions (Grant No. 1716009), the Social Science Foundation of Shandong Province (Grant Nos. 18DSHJ04 and 17CQXJ01), the Shandong Jiaotong University “Climbing” Research Innovation Team Program, the Risk Management and Insurance Research Team of Shandong University of Finance and Economics, the 1251 Talent Cultivation Project of Shandong Jiaotong University, and Collaborative Innovation Center Project of the Transformation of New and Old Kinetic Energy and Government Financial Allocation.

Acknowledgments

We thank two anonymous reviewers for their helpful comments on an earlier draft of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gerber, H.U.; Shiu, E.S.W. On the time value of ruin. N. Am. Actuar. J. 1998, 2, 48–78. [Google Scholar] [CrossRef]
  2. Lin, X.S.; Willmot, G.E.; Drekic, S. The classical risk model with a constant dividend barrier: Analysis of the Gerber-Shiu discounted penalty function. Insur. Math. Econ. 2003, 33, 391–408. [Google Scholar]
  3. Yuen, K.C.; Wang, G.J.; Li, W.K. The Gerber-Shiu expected discounted penalty function for risk processes with interest and a constant dividend barrier. Insur. Math. Econ. 2007, 40, 104–112. [Google Scholar] [CrossRef]
  4. Yuen, K.C.; Zhou, M.; Guo, J.Y. On a risk model with debit interest and dividend payments. Stat. Probab. Lett. 2008, 78, 2426–2432. [Google Scholar] [CrossRef]
  5. Zhao, X.H.; Yin, C.C. The Gerber-Shiu expected discounted penalty function for Lévy insurance risk processes. Acta Math. Sin. 2010, 26, 575–586. [Google Scholar] [CrossRef]
  6. Chi, Y.C. Analysis of expected discounted penalty function for a general jump diffusion risk model and applications in finance. Insur. Math. Econ. 2010, 46, 385–396. [Google Scholar] [CrossRef]
  7. Yin, C.C.; Wang, C.W. The perturbed compound Poisson risk process with investment and debit Interest. Methodol. Comput. Appl. 2010, 12, 391–413. [Google Scholar] [CrossRef]
  8. Chi, Y.C.; Lin, X.S. On the threshold dividend strategy for a generalized jump-diffusion risk model. Insur. Math. Econ. 2011, 48, 326–337. [Google Scholar] [CrossRef]
  9. Shen, Y.; Yin, C.C.; Yuen, K.C. Alternative approach to the optimality of the threshold strategy for spectrally negative Lévy processes. Acta Math. Appl. Sin. Engl. 2013, 29, 705–716. [Google Scholar] [CrossRef]
  10. Yin, C.C.; Yuen, K.C. Exact joint laws associated with spectrally negative Lévy processes and applications to insurance risk theory. Front. Math. China 2014, 9, 1453–1471. [Google Scholar] [CrossRef]
  11. Zhao, Y.X.; Yao, D.J. Optimal dividend and capital injection problem with a random time horizon and a ruin penalty in the dual model. Appl. Math. Ser. B 2015, 30, 325–339. [Google Scholar] [CrossRef]
  12. Li, S.M.; Lu, Y.; Sendova, K.P. The expected discounted penalty function: From infinite time to finite time. Scand. Actuar. J. 2019. [Google Scholar] [CrossRef]
  13. Li, Y.Q.; Yin, C.C.; Zhou, X.W. On the last exit times for spectrally negative Lévy processes. J. Appl. Probab. 2017, 54, 474–489. [Google Scholar] [CrossRef] [Green Version]
  14. Dong, H.; Yin, C.C.; Dai, H.S. Spectrally negative Lévy risk model under Erlangized barrier strategy. J. Comput. Appl. Math. 2019, 351, 101–116. [Google Scholar] [CrossRef]
  15. Wang, W.Y.; Zhang, Z.M. Computing the Gerber-Shiu function by frame duality projection. Scand. Actuar. J. 2019. [Google Scholar] [CrossRef]
  16. Willmot, G.E.; Dickson, D.C.M. The Gerber-Shiu discounted penalty function in the stationary renewal risk model. Insur. Math. Econ. 2003, 32, 403–411. [Google Scholar] [CrossRef]
  17. Wang, C.W.; Yin, C.C.; Li, E.Q. On the classical risk model with credit and debit interests under absolute ruin. Stat. Probab. Lett. 2010, 80, 427–436. [Google Scholar] [CrossRef]
  18. Yin, C.C.; Yuen, K.C. Optimality of the threshold dividend strategy for the compound Poisson model. Stat. Probab. Lett. 2011, 81, 1841–1846. [Google Scholar] [CrossRef]
  19. Dong, H.; Yin, C.C. Complete monotonicity of the probability of ruin and DE Finetti’s dividend problem. J. Syst. Sci. Complex. 2012, 25, 178–185. [Google Scholar] [CrossRef]
  20. Yu, W.G. Some results on absolute ruin in the perturbed insurance risk model with investment and debit interests. Econ. Model. 2013, 31, 625–634. [Google Scholar] [CrossRef]
  21. Yin, C.C.; Shen, Y.; Wen, Y.Z. Exit problems for jump processes with applications to dividend problems. J. Comput. Appl. Math. 2013, 245, 30–52. [Google Scholar] [CrossRef]
  22. Zeng, Y.; Li, D.P.; Gu, A.L. Robust equilibrium reinsurance-investment strategy for a mean-variance insurer in a model with jumps. Insur. Math. Econ. 2016, 66, 138–152. [Google Scholar] [CrossRef]
  23. Zhao, Y.X.; Wang, R.M.; Yao, D.J. Optimal dividend and equity issuance in the perturbed dual model under a penalty for ruin. Commun. Stat. Theory Methods 2016, 45, 365–384. [Google Scholar] [CrossRef]
  24. Yu, W.G.; Huang, Y.J.; Cui, C.R. The absolute ruin insurance risk model with a threshold dividend strategy. Symmetry 2018, 10, 377. [Google Scholar] [CrossRef]
  25. Boucherie, R.J.; Boxma, O.J.; Sigman, K. A note on negative customers, GI/G/1 workload, and risk process. Probab. Eng. Inf. Sci. 1997, 11, 305–311. [Google Scholar] [CrossRef]
  26. Boikov, A.V. The Cramér-Lundberg model with stochastic premium process. Theory Probab. Appl. 2003, 47, 489–493. [Google Scholar] [CrossRef]
  27. Temnov, G. Risk process with random income. J. Math. Sci. 2004, 123, 3780–3794. [Google Scholar] [CrossRef]
  28. Bao, Z.H. The expected discounted penalty at ruin in the risk process with random income. Appl. Math. Comput. 2006, 179, 559–566. [Google Scholar] [CrossRef]
  29. Bao, Z.H.; Ye, Z.X. The Gerber-Shiu discounted penalty function in the delayed renewal risk process with random income. Appl. Math. Comput. 2007, 184, 857–863. [Google Scholar] [CrossRef]
  30. Yang, H.; Zhang, Z.M. On a class of renewal risk model with random income. Appl. Stoch. Model. Bus. 2009, 25, 678–695. [Google Scholar] [CrossRef]
  31. Labbé, C.; Sendova, K.P. The expected discounted penalty function under a risk model with sochastic income. Appl. Math. Comput. 2009, 215, 1852–1867. [Google Scholar]
  32. Zhang, Z.M.; Yang, H. On a risk model with stochastic premiums income and dependence between income and loss. J. Comput. Appl. Math. 2010, 234, 44–57. [Google Scholar] [CrossRef] [Green Version]
  33. Zhao, Y.X.; Yin, C.C. The expected discounted penalty function under a renewal risk model with stochastic income. Appl. Math. Comput. 2012, 218, 6144–6154. [Google Scholar] [CrossRef]
  34. Xu, L.; Wang, R.M.; Yao, D.J. On maximizing the expected terminal utility by investment and reinsurance. J. Ind. Manag. Optim. 2008, 4, 801–815. [Google Scholar]
  35. Yu, W.G. On the expected discounted penalty function for a Markov regime switching risk model with stochastic premium income. Discret. Dyn. Nat. Soc. 2013, 2013, 1–9. [Google Scholar] [CrossRef]
  36. Yu, W.G. Randomized dividends in a discrete insurance risk model with stochastic premium income. Math. Probl. Eng. 2013, 2013, 1–9. [Google Scholar] [CrossRef]
  37. Zhou, J.M.; Mo, X.Y.; Ou, H.; Yang, X.Q. Expected present value of total dividends in the compound binomial model with delayed claims and random income. Acta Math. Sci. 2013, 33, 1639–1651. [Google Scholar] [CrossRef]
  38. Gao, J.W.; Wu, L.Y. On the Gerber-Shiu discounted penalty function in a risk model with two types of delayed-claims and random income. J. Comput. Appl. Math. 2014, 269, 42–52. [Google Scholar] [CrossRef]
  39. Zhou, M.; Yuen, K.C.; Yin, C.C. Optimal investment and premium control in a nonlinear diffusion model. Acta Math. Appl. Sin. Engl. 2017, 33, 945–958. [Google Scholar] [CrossRef]
  40. Deng, Y.C.; Liu, J.; Huang, Y.; Li, M.; Zhou, J.M. On a discrete interaction risk model with delayed claims and stochastic incomes under random discount rates. Commun. Stat. Theory Methods 2018, 47, 5867–5883. [Google Scholar] [CrossRef]
  41. Zeng, Y.; Li, D.P.; Chen, Z.; Yang, Z. Ambiguity aversion and optimal derivative-based pension investment with stochastic income and volatility. J. Econ. Dyn. Control 2018, 88, 70–103. [Google Scholar] [CrossRef]
  42. Politis, K. Semiparametric estimation for non-ruin probabilities. Scand. Actuar. J. 2003, 2003, 75–96. [Google Scholar] [CrossRef]
  43. Wang, Y.F.; Yin, C.C. Approximation for the ruin probabilities in a discrete time risk model with dependent risks. Stat. Probab. Lett. 2010, 80, 1335–1342. [Google Scholar] [CrossRef]
  44. Yuen, K.C.; Yin, C.C. Asymptotic results for tail probabilities of sums of dependent and heavy-tailed random variables. Chin. Ann. Math. B 2012, 33, 557–568. [Google Scholar] [CrossRef]
  45. Zhang, Z.M.; Yang, H.L.; Yang, H. On a nonparametric estimator for ruin probability in the classical risk model. Scand. Actuar. J. 2014, 2014, 309–338. [Google Scholar] [CrossRef]
  46. Wang, Y.F.; Yin, C.C.; Zhang, X.S. Uniform estimate for the tail probabilities of randomly weighted sums. Acta Math. Appl. Sin. Engl. 2014, 30, 1063–1072. [Google Scholar] [CrossRef]
  47. Zhang, Z.M. Nonparametric estimation of the finite time ruin probability in the classical risk model. Scand. Actuar. J. 2017, 2017, 452–469. [Google Scholar] [CrossRef]
  48. Zhang, Z.M.; Yang, H.L. Nonparametric estimate of the ruin probability in a pure-jump Lévy risk model. Insur. Math. Econ. 2013, 53, 24–35. [Google Scholar] [CrossRef] [Green Version]
  49. Zhang, Z.M.; Yang, H.L. Nonparametric estimation for the ruin probability in a Lévy risk model under low-frequency observation. Insur. Math. Econ. 2014, 59, 168–177. [Google Scholar] [CrossRef] [Green Version]
  50. Shimizu, Y.; Zhang, Z.M. Estimating Gerber-Shiu functions from discretely observed Lévy driven surplus. Insur. Math. Econ. 2017, 74, 84–98. [Google Scholar] [CrossRef]
  51. Peng, J.Y.; Wang, D.C. Asymptotics for ruin probabilities of a non-standard renewal risk model with dependence structures and exponential lévy process investment returns. J. Ind. Manag. Optim. 2017, 13, 155–185. [Google Scholar] [CrossRef]
  52. Peng, J.Y.; Wang, D.C. Uniform asymptotics for ruin probabilities in a dependent renewal risk model with stochastic return on investments. Stochastics 2018, 90, 432–471. [Google Scholar] [CrossRef]
  53. Yang, Y.; Wang, K.Y.; Liu, J.J.; Zhang, Z.M. Asymptotics for a bidimensional risk model with two geometric Lévy price processes. J. Ind. Manag. Optim. 2019, 15, 481–505. [Google Scholar] [CrossRef]
  54. Yang, Y.; Yuen, K.C.; Liu, J.F. Asymptotics for ruin probabilities in Lévy-driven risk models with heavy-tailed claims. J. Ind. Manag. Optim. 2018, 14, 231–247. [Google Scholar] [CrossRef]
  55. Shimizu, Y. Estimation of the expected discounted penalty function for Lévy insurance risks. Math. Methods Stat. 2011, 20, 125–149. [Google Scholar] [CrossRef]
  56. Zhang, Z.M. Estimating the Gerber-Shiu function by Fourier-Sinc series expansion. Scand. Actuar. J. 2017, 2017, 898–919. [Google Scholar] [CrossRef]
  57. Zhang, Z.M.; Su, W. A new efficient method for estimating the Gerber-Shiu function in the classical risk model. Scand. Actuar. J. 2018, 2018, 426–449. [Google Scholar] [CrossRef]
  58. Zhang, Z.M.; Su, W. Estimating the Gerber-Shiu function in a Lévy risk model by Laguerre series expansion. J. Comput. Appl. Math. 2019, 346, 133–149. [Google Scholar] [CrossRef]
  59. Su, W.; Yong, Y.D.; Zhang, Z.M. Estimating the Gerber-Shiu function in the perturbed compound Poisson model by Laguerre series expansion. J. Math. Anal. Appl. 2019, 469, 705–729. [Google Scholar] [CrossRef]
  60. Fang, F.; Oosterlee, C.W. A novel option pricing method based on Fourier cosine series expansions. SIAM J. Sci. Comput. 2008, 31, 826–848. [Google Scholar] [CrossRef]
  61. Fang, F.; Oosterlee, C.W. Pricing early-exercise and discrete barrier options by Fourier-cosine series expansions. Numer. Math. 2009, 114, 27–62. [Google Scholar] [CrossRef]
  62. Zhang, B.; Oosterlee, C.W. Pricing of early-exercise Asian options under Lévy processes based on Fourier cosine expansions. Appl. Numer. Math. 2014, 78, 14–30. [Google Scholar] [CrossRef] [Green Version]
  63. Chau, K.W.; Yam, S.C.P.; Yang, H.L. Fourier-cosine method for Gerber-Shiu functions. Insur. Math. Econ. 2015, 61, 170–180. [Google Scholar] [CrossRef]
  64. Chau, K.W.; Yam, S.C.P.; Yang, H.L. Fourier-cosine method for ruin probabilities. J. Comput. Appl. Math. 2015, 281, 94–106. [Google Scholar] [CrossRef] [Green Version]
  65. Zhang, Z.M. Approximating the density of the time to ruin via Fourier-cosine series expansion. Astin Bull. 2017, 47, 169–198. [Google Scholar] [CrossRef]
  66. Yang, Y.; Su, W.; Zhang, Z.M. Estimating the discounted density of the deficit at ruin by Fourier cosine series expansion. Stat. Probab. Lett. 2019, 146, 147–155. [Google Scholar] [CrossRef]
  67. Dickson, D.C.M.; Hipp, C. On the time to ruin for Erlang(2) risk processes. Insur. Math. Econ. 2001, 29, 333–344. [Google Scholar] [CrossRef]
  68. Li, S.M.; Garrido, J. On ruin for the Erlang(n) risk process. Insur. Math. Econ. 2004, 34, 391–408. [Google Scholar] [CrossRef]
  69. Van Der Vaart, A.W. Asymptotic Statistics; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
Figure 1. Estimation of the ruin probability for different claim sizes. Mean curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Figure 1. Estimation of the ruin probability for different claim sizes. Mean curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Mathematics 07 00305 g001
Figure 2. Estimation of the Laplace transform of ruin time for different claim sizes. Mean curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Figure 2. Estimation of the Laplace transform of ruin time for different claim sizes. Mean curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Mathematics 07 00305 g002
Figure 3. Estimation of the expected discounted deficit at ruin due to a claim for different claim sizes. Mean curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Figure 3. Estimation of the expected discounted deficit at ruin due to a claim for different claim sizes. Mean curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Mathematics 07 00305 g003
Figure 4. Estimation of the ruin probability for different claim sizes. Mean relative error curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Figure 4. Estimation of the ruin probability for different claim sizes. Mean relative error curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Mathematics 07 00305 g004
Figure 5. Estimation of the Laplace transform of ruin time for different claim sizes. Mean relative error curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Figure 5. Estimation of the Laplace transform of ruin time for different claim sizes. Mean relative error curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Mathematics 07 00305 g005aMathematics 07 00305 g005b
Figure 6. Estimation of the expected discounted deficit at ruin due to a claim for different claim sizes. Mean relative error curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Figure 6. Estimation of the expected discounted deficit at ruin due to a claim for different claim sizes. Mean relative error curves: (a) exponential claim sizes; (b) Erlang claim sizes; (c) combined-exponential claim sizes; and (d) mixed-exponential claim sizes.
Mathematics 07 00305 g006
Table 1. IMSE of Φ ^ K , a .
Table 1. IMSE of Φ ^ K , a .
TExpErlang(2)
RPLTEDDRPLTEDD
130.007170.006440.015860.003510.004140.00412
260.003520.003160.008460.001860.002000.00214
520.001770.001990.005370.001080.000990.00095
2600.000320.000340.000920.000180.000200.00019
Table 2. IMSE of Φ ^ K , a .
Table 2. IMSE of Φ ^ K , a .
TCom-ExpMix-Exp
RPLTEDDRPLTEDD
130.007190.004310.005270.004740.005110.01176
260.002710.001950.002770.002670.002330.00602
520.001600.001140.001240.001170.001070.00257
2600.000460.000210.000240.000250.000180.00053

Share and Cite

MDPI and ACS Style

Wang, Y.; Yu, W.; Huang, Y.; Yu, X.; Fan, H. Estimating the Expected Discounted Penalty Function in a Compound Poisson Insurance Risk Model with Mixed Premium Income. Mathematics 2019, 7, 305. https://doi.org/10.3390/math7030305

AMA Style

Wang Y, Yu W, Huang Y, Yu X, Fan H. Estimating the Expected Discounted Penalty Function in a Compound Poisson Insurance Risk Model with Mixed Premium Income. Mathematics. 2019; 7(3):305. https://doi.org/10.3390/math7030305

Chicago/Turabian Style

Wang, Yunyun, Wenguang Yu, Yujuan Huang, Xinliang Yu, and Hongli Fan. 2019. "Estimating the Expected Discounted Penalty Function in a Compound Poisson Insurance Risk Model with Mixed Premium Income" Mathematics 7, no. 3: 305. https://doi.org/10.3390/math7030305

APA Style

Wang, Y., Yu, W., Huang, Y., Yu, X., & Fan, H. (2019). Estimating the Expected Discounted Penalty Function in a Compound Poisson Insurance Risk Model with Mixed Premium Income. Mathematics, 7(3), 305. https://doi.org/10.3390/math7030305

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