Next Article in Journal
Secretary Problem with Possible Errors in Observation
Next Article in Special Issue
Distributed Fusion Estimation with Sensor Gain Degradation and Markovian Delays
Previous Article in Journal
Solving Integro-Differential Boundary Value Problems Using Sinc-Derivative Collocation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotically Normal Estimators of the Gerber-Shiu Function in Classical Insurance Risk Model

School of Insurance, Shandong University of Finance and Economics, Jinan 250014, China
*
Author to whom correspondence should be addressed.
Submission received: 8 July 2020 / Revised: 1 August 2020 / Accepted: 10 August 2020 / Published: 23 September 2020
(This article belongs to the Special Issue Stochastic Statistics and Modeling)

Abstract

:
Nonparametric estimation of the Gerber-Shiu function is a popular topic in insurance risk theory. Zhang and Su (2018) proposed a novel method for estimating the Gerber-Shiu function in classical insurance risk model by Laguerre series expansion based on the claim number and claim sizes of sample. However, whether the estimators are asymptotically normal or not is unknown. In this paper, we give the details to verify the asymptotic normality of these estimators and present some simulation examples to support our result.

1. Introduction

In this paper, we consider the classical risk model
U t = u + c t i = 1 N t X i , t 0 ,
where u 0 denotes the initial surplus level and c > 0 denotes the premium rate. The claim number { N t } t 0 follows a Poisson process with intensity λ > 0 . X 1 , X 2 , denote the claim sizes which are independent and identically distributed random variables with density f and mean μ . Meanwhile, suppose that { N t } and { X i } are mutually independent. Let τ = inf t 0 : U t < 0 denote ruin time and we set τ = when U t 0 for all t 0 .
Gerber and Shiu [1] considered the expected discounted penalty function about the ruin time τ , the surplus before ruin time U τ and the deficit at ruin time | U τ | which is defined
ϕ ( u ) = E [ e δ τ w ( U τ ) , | U τ | 1 ( τ < ) | U 0 = u ] , u 0 ,
where w : R + × R + R + is a measurable penalty function of U τ and | U τ | , 1 { · } is the indicator function, δ is the interest force. In this paper, we verify asymptotically normal estimators of the Gerber-Shiu function when δ = 0 . Many researchers have made remarkable contributions to this model and its generalizations have been made in various risk model. Asmussen and Albercher [2] considered the classical risk model; Gerber and Landry [3], Tsai [4,5] considered the compound Poisson risk model perturbed by diffusion; Yu [6] and Yu et al. [7] investigated the absolute ruin risk model. Zhao and Yin [8] studied the Gerber-Shiu function in Lévy risk model; Xie and Zou [9] and Zhu et al. [10] studied the compound Poisson with delayed claims. The Gerber-Shiu function has also been studied by Yin and Wang [11], Shen et al. [12], Yin and Yuen [13], Cai et al. [14], Deng et al. [15], Dong et al. [16], Wang and Zhang [17], Yu et al. [18], Peng et al. [19] among others.
In all of these mentioned papers, they assumed that the parameter λ and the claim size density f are known. However, the probability characteristics of the surplus process are usually unknown. Recently, statistical estimation of risk measures based on the observed data information of claim number and claim sizes has become a popular topic, see e.g., Shimizu [20] applied regularized Laplace transform to estimate the Gerber-Shiu function in Lévy risk model; Zhang and Yang [21,22] estimated the ruin probability based on high-frequency observation and low-frequency observation, respectively; Zhang [23] construted an estimator of the Gerber-Shiu by Fourier-Sinc series expansion and Zhang [24] estimated the finite time ruin probability by double Fourier transform. Later, Zhang and Su [25] proposed a new method for estimating the Gerber-Shiu function by Laguerre basis and derived the convergence rate of the estimate. In this paper, we show that the above estimator for the Gerber-Shiu function is asymptotically normal. Asymptotic normality results allow us to establish confidence intervals. For more study on the statistical estimation of risk model, the interested readers are referred to the work of Shimizu and Zhang [26] Yu et al. [27], Peng and Wang [28], Yu et al. [29] and Zhang et al. [30].
The remainder of this paper is organized as follows. We present some preliminaries on Laguerre expansion of Gerber-Shiu function and the estimator of Gerber-Shiu function in Section 2. We give the details to verify the asymptotic normality of these estimators in Section 3. Finally, we present some numerical examples to verify the proposed estimator is asymptotically normal in Section 4.

2. Preliminaries

First, we introduce the following notation:
  • R + = [ 0 , ) , N = { 0 , 1 , 2 , } , N 0 = { 1 , 2 , } .
  • L 2 ( R + ) = f : R + R | R + | f ( x ) | 2 d x < .
  • L f ( s ) denotes the Laplace transform of claim size density f,
    L f ( s ) = R + e s x f ( x ) d x .
  • Φ X = E [ e x p ( i t X ) ] denotes the characteristic function of random vector X R k , where t is a random vector in R k , k N 0 .
  • For f , g L 2 ( R + ) , let f , g = R + f ( x ) g ( x ) d x be the scalar product and f 2 = R + f 2 ( x ) d x be L 2 -norm.
  • For positive function f 1 , f 2 , let f 1 f 2 be f 1 ( x ) C f 2 ( x ) , where C is a positive constant.
  • ⊤ means the transpose of matrix.
  • 0 is a zero vector.
  • Let P be convergence in probability and D be convergence in distribution.
  • f = O ( g ) means that | f | C | g | for some constant C.
Meanwhile, we need the following conditions, which have also been considered in Zhang and Su [25]:
  • Condition 1 The premium rate c > λ μ ;
  • Condition 2 Suppose that R + R + ( 1 + x ) w ( x , y ) f ( x + y ) d y d x < ;
  • Condition 3 For some k 1 , k 2 R + , suppose that the penalty function w ( x 1 , x 2 ) 1 + x 1 k 1 1 + x 2 k 2 .

2.1. Laguerre Expansion of Gerber-Shiu Function

The Laguerre functions { ψ k } k N are defined by
ψ k ( x ) = 2 L k ( 2 x ) e x , x R + , k N .
where L k ( x ) = j = 0 k ( 1 ) j k j x j j ! , x 0 . It follows from Abramowitz and Stegun [31] that
  • { ψ k } k N form a complete orthogonal basis over L 2 ( R + ) . Then, when f L 2 ( R + ) , it can be expanded by the Laguerre function: f ( x ) = k N f , ψ k ψ k ( x ) ;
  • { ψ k } k N are uniformly bounded, i.e. | ψ k ( x ) | 2 , x R + , k N ;
  • k , j N , 0 x ψ k ( x y ) ψ j ( y ) d y = 1 2 [ ψ k + j ( x ) ψ k + j + 1 ( x ) ] , ψ k , ψ k = 1 and ψ j , ψ k = 0 for k j .
According to Gerber and Shiu [1], the Gerber-Shiu function satisfies
ϕ ( u ) = 0 u ϕ ( u x ) g ( x ) d x + h ( u ) , u 0 ,
where g ( x ) = λ c x f ( y ) d y , h ( u ) = λ c u x w ( x , y x ) f ( y ) d y d x . Since ϕ , g , h L 2 ( R + ) by the condition 1, 2 of Zhang and Su [25], ϕ , g , h can be expressed by Laguerre basis as
ϕ ( x ) = k P k ψ k ( x ) , g ( x ) = k Q k ψ k ( x ) , h ( x ) = k R k ψ k ( x ) , x R + , k N ,
where P k = ϕ , ψ k , Q k = g , ψ k , R k = h , ψ k . It follows from Zhang and Su [25], we obtain the approximation of the Gerber-Shiu function
ϕ ( u ) = k = 1 P k ψ k ( x ) k = 0 K P k ψ k ( x ) : = ϕ K ( u ) = p K ψ K ( u ) ,
where K N 0 denotes truncation parameter, ψ K ( u ) = ( ψ 0 ( u ) , ψ 1 ( u ) , , ψ K ( u ) ) and p K = ( P 0 , P 1 , , P K ) . p K can be expressed as p K = A K 1 r K , where r K = ( R 0 , R 1 , , R K ) and A K = ( a i j ) 1 i , j K + 1 is a lower triangular invertible Toeplitz matrix, whose components are given by
a i j = 1 1 2 Q 0 i = j ; 1 2 ( Q i j 1 Q i j ) i > j ; 0 i < j .

2.2. Coefficient Q k and R k

By changing the order of integrals, we have
Q k = 0 g ( x ) ψ k ( x ) d x = λ c 0 x f X ( y ) d y ψ k ( x ) d x = λ c 0 0 y ψ k ( x ) d x f X ( y ) d y = λ c E [ H k Q ( X ) ] , k N
and
R k = 0 h ( u ) ψ k ( u ) d u = λ c 0 u x w ( x , y x ) f X ( y ) d y d x ψ k ( u ) d u = λ c 0 0 y u y w ( x , y x ) d x ψ k ( u ) d u f X ( y ) d y = λ c E [ H k R ( X ) ] , k N
where H k Q ( z ) = 0 z ψ k ( x ) d x , H k R ( z ) = 0 z u z w ( x , z x ) d x ψ k ( u ) d u . Meanwhile, by Condition 3 and { ψ k } k N are uniformly bounded, we have
sup k | H k Q ( z ) | 2 0 z d x z , z R +
and
sup k | H k R ( z ) | 2 0 z u z w ( x , z x ) d x d u z k 1 + k 2 + 2 , z R + .

2.3. Statistical Inference

For insurer, the parameter λ and the claim size density f are usually unknown. But, they can be obtained by the following data information,
N T , X 1 , X 2 , , X N T ,
where N T is the claim number over [ 0 , T ] and { X j } 1 j N T are individual claim sizes. We can estimate Q k and R k by
Q ^ k = λ ^ c 1 N t i = 1 N t H k Q ( X i ) , R ^ k = λ ^ c 1 N t i = 1 N t H k R ( X i ) ,
where λ ^ = N T T is the estimator of λ . Furthermore, we can estimate the Gerber-Shiu function by
ϕ ^ K ( u ) : = p ^ K ψ K ( u ) ,
where p ^ K = ( P ^ 0 , P ^ 1 , , P ^ K ) are the estimators of p K . p ^ K holds that p ^ K = A ^ K 1 r ^ K , where r ^ K = ( R ^ 0 , R ^ 1 , , R ^ K ) and A ^ K = ( a ^ i j ) 1 i , j K + 1 are given by
a ^ i j = 1 1 2 Q ^ 0 i = j ; 1 2 ( Q ^ i j 1 Q ^ i j ) i > j ; 0 i < j .

3. Asymptotically Normality

In this part, we show that the estimator is asymptotically normal. For this purpose, we introduce some lemmas for the asymptotic normality of Laguerre coefficients.
Lemma 1.
Suppose E [ X 2 ( 2 + k 1 + k 2 ) ] < and sup k | H k ( x ) | x k 1 + k 2 + 2 , let
H = λ ^ c 1 N T i = 1 N T H 0 ( X i ) , , λ ^ c 1 N T i = 1 N T H d ( X i ) and u = λ c E [ H 0 ( X ) ] , , λ c E [ H d ( X ) ] .
Then
N T H u D N d + 1 0 , λ 2 c 2 Σ , T ,
where Σ = ( σ i j 2 ) 0 i , j d , σ i j 2 : = E ( H i ( X ) E [ H i ( X ) ] ) ( H j ( X ) E [ H j ( X ) ] ) < .
Proof of Lemma 1.
Since N T , which is independent of X j , follows Poisson distribution with intensity λ T , we have λ c E [ H i ( X ) ] = λ c 1 N T j = 1 N T E [ H i ( X j ) ] and
H i u i = λ ^ c λ c 1 N T j = 1 N T H i ( X j ) + λ c 1 N T j = 1 N T H i ( X j ) E [ H i ( X j ) ] , 0 i d .
Furthermore, since λ ^ c P λ c when T and E [ sup k | H k ( X ) | ] < , we have
λ ^ c λ c 1 N T j = 1 N T H i ( X j ) P 0 , 0 i d .
For convenience, we set
Z 1 : = 1 N T j = 1 N T ( H 0 ( X j ) E [ H 0 ( X j ) ] ) , , 1 N T j = 1 N T ( H d ( X j ) E [ H d ( X j ) ] ) = 1 N T j = 1 N T Z 2 ; Z 2 : = H 0 ( X j ) E [ H 0 ( X j ) ] , H 1 ( X ) E [ H 1 ( X j ) ] , , H d ( X j ) E [ H d ( X j ) ] ; Z 3 : = H 0 ( X ) E [ H 0 ( X ) ] , H 1 ( X ) E [ H 1 ( X ) ] , , H d ( X ) E [ H d ( X ) ] .
We can obtain E Z 3 = 0 and for every random vector t R d + 1 ,
E e i N T t Z 1 = E E e x p i 1 N T j = 1 N T t Z 2 | N T = E Φ Z 3 t N T N T .
It follows from (VI 2.13) in Cinlar [32] that N T almost everywhere when T . Furthermore, by dominated convergence theorem, we obtain
lim T E e i N T t Z 1 = E lim T Φ Z 3 t N T N T = e 1 2 t Σ t .
According to Theorem 2.13 in van der Vaart [33], we derive that N T H u D N d + 1 0 , λ 2 c 2 Σ , T . This completes the proof. □
Next, we have a weak consistency for Laguerre coefficients Q ^ k , R ^ k .
Lemma 2.
If E X 2 ( k 1 + k 2 + 2 ) < , then ( Q ^ k , R ^ k ) P ( Q k , R k ) for all k N when T .
Proof of Lemma 2.
For ε > 0 ,
P ( | Q ^ k Q k | > ε ) E [ Q ^ k Q k ] 2 ε 2 = 1 c 2 T 2 ε 2 V a r j = 1 N T 0 X i ψ k ( x ) d x = λ c 2 T ε 2 E 0 X ψ k ( x ) d x 2 2 λ c 2 T ε 2 E X 2 1 ε 2 T 0 , T ,
and
P ( | R ^ k R k | > ε ) E [ R ^ k R k ] 2 ε 2 = 1 c 2 T 2 ε 2 V a r j = 1 N T 0 X i u X i w ( x , X i x ) d x ψ k ( u ) d u = λ c 2 T ε 2 E 0 X u X w ( x , X x ) d x ψ k ( u ) d u 2 2 λ c 2 T ε 2 E X 2 ( k 1 + k 2 + 2 ) 1 ε 2 T 0 , T .
It implies that Q ^ k P Q k and R ^ k P R k . By Theorem 2.7 (vi) in van der Vaart [33], we have that ( Q ^ k , R ^ k ) P ( Q k , R k ) , T .  □
Then, we prove that the Laguerre coefficients are asymptotically normal.
Lemma 3.
Suppose that E X 2 ( k 1 + k 2 + 2 ) < . Then, for all K N we have
N T q ^ K q K r ^ K r K D N 2 ( K + 1 ) 0 , Σ Q , R , T ,
where Σ Q , R = σ Q Q σ Q R σ Q R σ R R , σ X Y = ( σ i , j X Y ) 0 i , j K , σ i , j X Y = 0 ( H i X ( z ) X i ) ( H i Y ( z ) Y j ) f ( z ) d z < , X = Q or R, Y = Q or R.
Proof of Lemma 3.
According to Lemma 1, N T q ^ K q K r ^ K r K converges in distribution to a normal variable with mean 0 . Then, we study its convariance matrix. For convenience, we set Y = ( Y 0 , Y 1 , , Y 2 K + 1 ) and Y ^ = ( Y ^ 0 , Y ^ 1 , , Y ^ 2 K + 1 ) with
Y i = Q i 0 i K ; R i ( K + 1 ) K + 1 i 2 K + 1 , Y ^ i = Q ^ i , 0 i K ; R ^ i ( K + 1 ) , K + 1 i 2 K + 1 .
where
N T Y ^ Y D N d + 1 ( 0 , Σ Q , R ) , T .
The components of Σ Q , R = ( σ i j ) 0 i , j d are
σ i j = E λ c H i Q ( X ) Q i λ c H j Q ( X ) Q j , i f 0 i K , 0 j K ; E λ c H i Q ( X ) Q i λ c H j ( K + 1 ) R ( X ) R j ( K + 1 ) , i f 0 i K , K + 1 j 2 K + 1 ; E λ H i ( K + 1 ) R ( X ) c R i ( K + 1 ) λ c H j Q ( X ) Q j , i f K + 1 i 2 K + 1 , 0 j K ; E λ H i ( K + 1 ) R ( X ) c R i ( K + 1 ) λ H j ( K + 1 ) R ( X ) c R j ( K + 1 ) , i f K + 1 i 2 K + 1 , K + 1 j 2 K + 1 .
Following from Lemma 1, we have N T q ^ K q K r ^ K r K D N 2 ( K + 1 ) 0 , Σ Q , R , T . This completes the proof. □
Finally, we derive the asymptotic normality of Laguerre estimators.
Theorem 1.
Suppose that E X 2 ( k 1 + k 2 + 2 ) < , ϕ W ( R + , r , B ) and K = O ( T 1 r 2 ) for r > 2 , then
N T ϕ ^ K ( x ) ϕ ( x ) D N 0 , Σ ( x ) , T ,
where Σ ( x ) = Ψ K ( x ) A K 1 P K * Σ Q , R ( Ψ K ( x ) A K 1 P K * ) , Ψ K ( x ) = ( ψ 0 ( x ) , ψ 1 ( x ) , , ψ K ( x ) ) , P K * = ( P ¯ K * , I K + 1 ) and W ( R + , r , B ) denotes Sobolev-Laguerre space (see Bongioanni and Torrea [34]). P ¯ K * is given by P ¯ K * = ( P ¯ K * ) i j , where
( P ¯ K * ) i j = 1 2 P 0 , i f i = j , 1 2 ( P i j P i j 1 ) , i f i > j , 0 , i f i < j . w h e r e i , j = 0 , 1 , 2 , , K .
Proof of Theorem 1.
Since
N T ϕ ^ K ( x ) ϕ ( x ) = N T ϕ ^ K ( x ) ϕ K ( x ) + N T ϕ K ( x ) ϕ ( x ) ,
according to Shimizu and Zhang [26], we have
N T ( p ^ K p K ) = A K 1 ( P ^ K * I K + 1 ) N T q ^ K q K r ^ K r K ,
where
P ^ K * = ( P ^ K * ) i j D P ¯ K * , w h e r e ( P ^ K * ) i j = 1 2 P ^ 0 , i f i = j , 1 2 ( P ^ i j P ^ i j 1 ) , i f i > j , 0 , i f i < j . w h e r e i , j = 0 , 1 , 2 , , K .
Then
N T ( p ^ K p K ) D N ( 0 , A K 1 P K * Σ Q , R ( A K 1 P K * ) ) ,
and
N T ( ϕ ^ K ( x ) ϕ K ( x ) ) = Ψ K ( x ) N T ( p ^ K p K ) D N 0 , Ψ K ( x ) A K 1 P K * Σ Q , R ( Ψ K ( x ) A K 1 P K * ) .
For N T ϕ K ( x ) ϕ ( x ) , by Cauchy-Schwarz inequality,
sup x R + N T ϕ K ( x ) ϕ ( x ) = sup x R + N T k K + 1 ϕ , ψ k ψ k ( x ) 2 N T k K + 1 k r 2 ϕ , ψ k k r 2 2 N T k K + 1 k r ϕ , ψ k 2 1 2 k K + 1 k r 1 2 2 N T B K x r d x 1 2 = 2 B λ ^ r 1 1 2 T 1 2 K r 1 2 P 0 , T , r > 2 .
It means that N T ϕ K ( x ) ϕ ( x ) P 0 . This completes the proof. □

4. Simulation

In this part, we provide numerical examples to verify the asymptotic normality of Laguerre estimators. We set the premium rate c = 1.5 , the Poisson density λ = 1 and we consider three claim size densities:
  • Exponential: f ( x ) = e x , x R + ;
  • Erlang(2): f ( x ) = 4 x e 2 x , x R + ;
  • Combination-of-exponentials: f ( x ) = 3 e 1.5 x 3 e 3 x , x R + .
Meanwhile, we estimate three special Gerber-Shiu functions:
  • Ruin probability (RP);
  • Expected claim size causing ruin (ECS);
  • Expected deficit at ruin (ED).
By Asmussen and Albrecher [2] and Dickson [35], the explicit formulas for those function can be obtained by Laplace transform method. Note that under the above claim size densities and special Gerber-Shiu function assumptions, Conditions 1–3 hold true. We set T = 720 , 1440 and cut-off parameter K = 5 T 1 10 , where · means the integer part.
To check the asymptotic normality of the estimators, we present quantile-quantile plot (QQ-plot). QQ-plot displays each data point by ( + ) marks, if the data points display a linear trend, then the distribution of the data is normal. To conclude this discussion, we show QQ-plot for different claim size densities and Gerber-Shiu functions with T = 720 , 1440 . In Figure 1 and Figure 2, we plot QQ-plot of Ruin probability for exponential claim size density with different u and T. The plot display an approximately straight line mean that the estimators are asymptotic normal. Then, we display the QQ-plot of ECS for Erlang(2) density in Figure 3 and Figure 4. The data point produces a straight line means that the estimators are asymptotically normal. In the end, we present QQ-plot of ED for combination-of-exponentials density in Figure 5 and Figure 6 with different u and T. The results manifest the Laguerre estimators are asymptotically normal. In addition, it can be observed from Figure 1, Figure 3 and Figure 5 with T = 720 that the tails seem not converge to the normal distribution. Nevertheless, it is observed from Figure 2, Figure 4 and Figure 6 with T = 1440 that the results manifest asymptotic normality. In conclusion, we prove the asymptotic normality of our estimator as the value of T becomes large.

5. Conclusions

Recently, Zhang and Su [25] proposed a novel method for estimating the Gerber-Shiu function by Laguerre series expansion. Based on the observed information, we obtain the estimator of the Gerber-Shiu function. In this paper, we analyze the asymptotic normality of the Laguerre estimator mathematically and present some simulation experiments to support our result.

Author Contributions

Data curation, W.S.; formal analysis, W.S. and W.Y.; methodology, W.Y.; writing—original draft, W.S. All authors have read and agreed to the published version of the manuscript.

Funding

The research of Wenguang Yu was supported by the National Natural Science Foundation of China (No. 11301303), the National Social Science Foundation of China (No. 15BJY007), the Taishan Scholars Program of Shandong Province (No. tsqn20161041), the Humanities and Social Sciences Project of the Ministry Education of China (No. 19YJA910002), the Natural Science Foundation of Shandong Province (No. ZR2018MG002), the Fostering Project of Dominant Discipline and Talent Team of Shandong Province Higher Education Institutions (No. 1716009), Shandong Provincial Social Science Project Planning Research Project(No. 19CQXJ08), the Risk Management and Insurance Research Team of Shandong University of Finance and Economics, Excellent Talents Project of Shandong University of Finance and Economics, the Collaborative Innovation Center Project of the Transformation of New and old Kinetic Energy and Government Financial Allocation.

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. Asmussen, S.; Albrecher, H. Ruin Probabilities, 2nd ed.; World Scientific: Singapore, 2010. [Google Scholar]
  3. Gerber, H.U.; Landry, B. On the discounted penalty at ruin in a jump-diffusion and the perpetual put option. Insur. Math. Econ. 1998, 22, 263–276. [Google Scholar] [CrossRef]
  4. Tsai, C.C.L. On the discounted distribution functions of the surplus process perturbed by diffusion. Insur. Math. Econ. 2001, 28, 401–419. [Google Scholar] [CrossRef]
  5. Tsai, C.C.L. On the expectations of the present values of the time of ruin perturbed by diffusion. Insur. Math. Econ. 2003, 32, 413–429. [Google Scholar] [CrossRef]
  6. 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]
  7. 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] [Green Version]
  8. Zhao, X.H.; Yin, C.C. The Gerber-Shiu expected discounted penalty function for Lévy insurance risk processes. Acta Math. Appl. Sin. 2010, 26, 575–586. [Google Scholar] [CrossRef]
  9. Xie, J.H.; Zou, W. On the expected discounted penalty function for the compound Poisson risk model with delayed claims. J. Comput. Appl. Math. 2011, 235, 2392–2404. [Google Scholar] [CrossRef]
  10. Zhu, H.M.; Huang, Y.; Yang, X.Q.; Zhou, J.M. On the expected discounted penalty function for the classical risk model with potentially delayed claims and random incomes. J. Appl. Math. 2014, 2014, 717269. [Google Scholar] [CrossRef] [Green Version]
  11. 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]
  12. 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. 2013, 29, 705–716. [Google Scholar] [CrossRef] [Green Version]
  13. 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] [Green Version]
  14. Cai, C.; Chen, N.; You, H. Nonparametric estimation for a spectrally negative Lévy risk process based on low-frequency observation. J. Comput. Appl. Math. 2018, 328, 432–442. [Google Scholar] [CrossRef]
  15. 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]
  16. 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]
  17. Wang, W.Y.; Zhang, Z.M. Computing the Gerber-Shiu function by frame duality projection. Scand. Actuar. J. 2019, 2019, 291–307. [Google Scholar] [CrossRef]
  18. Yu, W.G.; Guo, P.; Wang, Q.; Guan, G.F.; Yang, Q.; Huang, Y.J.; Yu, X.L.; Jin, B.Y.; Cui, C.R. On a periodic capital injection and barrier dividend strategy in the compound Poisson risk model. Mathematics 2020, 8, 511. [Google Scholar] [CrossRef] [Green Version]
  19. Peng, X.H.; Su, W.; Zhang, Z.M. On a perturbed compound Poisson risk model under a periodic threshold-type dividend strategy. J. Ind. Manag. Optim. 2020, 16, 1967–1986. [Google Scholar] [CrossRef]
  20. Shimizu, Y. Nonparametric estimation of the Gerber-Shiu function for the Winer-Poisson risk model. Scand. Actuar. J. 2012, 1, 56–69. [Google Scholar] [CrossRef]
  21. 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]
  22. 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]
  23. Zhang, Z.M. Estimating the Gerber-Shiu function by Fourier-Sinc series expansion. Scand. Actuar. J. 2017, 10, 898–919. [Google Scholar] [CrossRef]
  24. Zhang, Z.M. Nonparametric estimation of the finite time ruin probability in the classical risk model. Scand. Actuar. J. 2017, 5, 452–469. [Google Scholar] [CrossRef]
  25. Zhang, Z.M.; Su, W. A new efficient method for estimating the Gerber-Shiu function in the classical risk model. Scand. Actuar. J. 2018, 5, 426–449. [Google Scholar] [CrossRef]
  26. Shimizu, Y.; Zhang, Z.M. Asymptotically normal estimators of the ruin probability for Lévy insurance surplus from discrete samples. Risks 2019, 7, 37. [Google Scholar] [CrossRef] [Green Version]
  27. Yu, W.G.; Wang, F.; Huang, Y.J.; Liu, H.D. Social optimal mean field control problem for population growth model. Asian J. Control 2019, 1–8. [Google Scholar] [CrossRef]
  28. 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]
  29. Yu, W.G.; Yong, Y.D.; Guan, G.F.; Huang, Y.J.; Su, W.; Cui, C.R. Valuing guaranteed minimum death benefits by cosine series expansion. Mathematics 2019, 7, 835. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, Z.M.; Yong, Y.D.; Yu, W.G. Valuing equity-linked death benefits in general exponential Lévy models. J. Comput. Appl. Math. 2020. [Google Scholar] [CrossRef]
  31. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; National Bureau of Standards Applied Mathematics Series; US Government Printing Office: Washington, DC, USA, 1964.
  32. Cinlar, E. Probability and Stochastics; Springer: New York, NY, USA, 2015. [Google Scholar]
  33. Van Der Vaart, A.W. Asymptotic Statistics; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
  34. Bongioanni, B.; Torrea, J.L. What is a Sobolev space for the Laguerre function system? Studia Math. 2009, 192, 147–172. [Google Scholar] [CrossRef] [Green Version]
  35. Dickson, D.C.M. Insurance Risk and Ruin; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
Figure 1. QQ-plot of RP for Exponential claim size density with T = 720 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Figure 1. QQ-plot of RP for Exponential claim size density with T = 720 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Mathematics 08 01638 g001
Figure 2. QQ-plot of RP for Exponential claim size density with T = 1440 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Figure 2. QQ-plot of RP for Exponential claim size density with T = 1440 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Mathematics 08 01638 g002
Figure 3. QQ-plot of ECS for Erlang(2) density with T = 720 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Figure 3. QQ-plot of ECS for Erlang(2) density with T = 720 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Mathematics 08 01638 g003
Figure 4. QQ-plot of ECS for Erlang(2) density with T = 1440 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Figure 4. QQ-plot of ECS for Erlang(2) density with T = 1440 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Mathematics 08 01638 g004
Figure 5. QQ-plot of ED for Combination-of-exponentials claim size density with T = 720 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Figure 5. QQ-plot of ED for Combination-of-exponentials claim size density with T = 720 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Mathematics 08 01638 g005
Figure 6. QQ-plot of ED for Combination-of-exponentials claim size density with T = 1440 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Figure 6. QQ-plot of ED for Combination-of-exponentials claim size density with T = 1440 : (a) u = 1 ; (b) u = 3 ; (c) u = 5 .
Mathematics 08 01638 g006

Share and Cite

MDPI and ACS Style

Su, W.; Yu, W. Asymptotically Normal Estimators of the Gerber-Shiu Function in Classical Insurance Risk Model. Mathematics 2020, 8, 1638. https://0-doi-org.brum.beds.ac.uk/10.3390/math8101638

AMA Style

Su W, Yu W. Asymptotically Normal Estimators of the Gerber-Shiu Function in Classical Insurance Risk Model. Mathematics. 2020; 8(10):1638. https://0-doi-org.brum.beds.ac.uk/10.3390/math8101638

Chicago/Turabian Style

Su, Wen, and Wenguang Yu. 2020. "Asymptotically Normal Estimators of the Gerber-Shiu Function in Classical Insurance Risk Model" Mathematics 8, no. 10: 1638. https://0-doi-org.brum.beds.ac.uk/10.3390/math8101638

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