Next Article in Journal
Consideration of Risk Factors in Corporate Property Portfolio Management
Next Article in Special Issue
Natural Disasters and Economic Growth: A Semiparametric Smooth Coefficient Model Approach
Previous Article in Journal
Effect of Fisheries Subsidies Negotiations on Fish Production and Interest Rate
Previous Article in Special Issue
The Environmental Kuznets Curve: A Semiparametric Approach with Cross-Sectional Dependence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonparametric Estimation of the Ruin Probability in the Classical Compound Poisson Risk Model

1
School of Mathematical Sciences, Qufu Normal University, Qufu 273165, China
2
College of Mathematics and Data Science, Minjiang University, Fuzhou 350108, China
3
Department of Mathematics and Statistics, School of Data Science, University of North Carolina, Charlotte, NC 28223, USA
4
School of Statistics, Qufu Normal University, Qufu 273165, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Risk Financial Manag. 2020, 13(12), 298; https://doi.org/10.3390/jrfm13120298
Submission received: 18 October 2020 / Revised: 16 November 2020 / Accepted: 23 November 2020 / Published: 29 November 2020
(This article belongs to the Special Issue Nonparametric Econometric Methods and Application II)

Abstract

:
In this paper we study estimating ruin probability which is an important problem in insurance. Our work is developed upon the existing nonparametric estimation method for the ruin probability in the classical risk model, which employs the Fourier transform but requires smoothing on the density of the sizes of claims. We propose a nonparametric estimation approach which does not involve smoothing and thus is free of the bandwidth choice. Compared with the Fourier-transformation-based estimators, our estimators have simpler forms and thus are easier to calculate. We establish asymptotic distributions of our estimators, which allows us to consistently estimate the asymptotic variances of our estimators with the plug-in principle and enables interval estimates of the ruin probability.

1. Introduction

It is well-known that estimating the ruin probability is an important problem in economics and insurance. There are many works on this topic in the literature, and they can be divided into two groups, parametric and nonparametric. See (Baumgartner and Gatto 2010); (Croux and Veraverbeke 1990), (Feng et al. 2020); (Frees 1986); (Mnatsakanov et al. 2008); (Pitts 1994); (Politis 2003); (You et al. 2020) and (Zhang et al. 2014), among others. In this paper we propose a nonparametric approach to estimating the ruin probability. Unlike the existing nonparametric methods, our procedure does not involve smoothing. To illustrate our idea, we begin with the classical risk model
U t = u + c t i = 1 N t X i , t 0 ,
where u 0 is the initial surplus, c > 0 is the constant premium rate, and N = { N t } t 0 is a Poisson process with intensity λ > 0 , and X 1 , X 2 , X 3 , are independent and identically distributed random variables with density function f and distribution function F supported on ( 0 , ) . The corresponding process N = { N t , t 0 } is called the claim number process, and X = { X i , i = 1 , 2 , 3 } are the sizes of claims. Further, N and { X i } i = 1 , 2 , are independent. Suppose the mean and variance of claim sizes are finite, i.e., μ = 0 x F ( d x ) < , σ 2 = 0 x 2 F ( d x ) μ 2 < .
Let τ = inf { t > 0 : U t < 0 } , which is the first time when U t < 0 . Then the ruin probability is
Ψ ( u ) = P τ < | U 0 = u .
In order to prevent ruin happening with probability one, we assume that the risk process has a relative safety loading (Gerber 1979) and (Bower et al. 1997), that is,
θ = c / ( λ μ ) 1 > 0 .
As illustrated in (Asmussen and Albrecher 2010) and (Landriault and Willmot 2008), there are many methods for estimating the ruin probability. For example, the integro-differential equation technique, renewal theory, Laplace transform, martingale theory, and so on. Since explicit formulas for ruin probabilities are usually not available, some authors roughly estimate some bounds of ruin probabilities (Cai and Wu 1997) and (Dickson 1994). However, these methods assume knowing the inter-claim distribution and the claim size distribution.
Naturally, some authors seek to nonparametrically estimate the ruin probability which does not require knowing the distribution of the claim size. Examples include but not limited to the sample reuse method in (Frees 1986), which was extended by (Croux and Veraverbeke 1990) by virtue of the Pollaczeck–Khinchin formula and U-statistics, the Laplace transform method in (Mnatsakanov et al. 2008), and the inverse Fourier method in (Zhang et al. 2014). These approaches assumed that the distribution of the individual claim size F is unknown, but the Poisson intensity λ is known. Other authors contributed their works with both F and λ estimated from the data, for example, (Pitts 1994); (Politis 2003); (Baumgartner and Gatto 2010) and (You et al. 2020).
In this paper, we make our contribution to statistical inference for Ψ ( u ) . Our approach is based on the Fourier inversion in (Zhang et al. 2014) and (You et al. 2020). Using the Fourier inversion and kernel density, (Zhang et al. 2014) obtained an explicit expression of Ψ ( u ) , with which they proposed a nonparametric estimator of the ruin probability based on the plug-in principle (see Section 2 for details).
The problem of estimating the probability of ruin when the Poisson intensity λ are unknown was not studied in (Zhang et al. 2014). Based on (Zhang et al. 2014) and (Politis 2003), (You et al. 2020) established asymptotically normal distributions of the nonparametric estimators of Ψ ( u ) with known λ and with unknown λ . Both (Zhang et al. 2014) and (You et al. 2020) employed the kernel density estimation of f ( x ) , which requires choice of the smoothing parameter. Specially, given a sample of claims X = { X i , i = 1 , 2 , , n } , (You et al. 2020) and (Zhang et al. 2014) estimated the claim size density f by the kernel density estimator
f ^ n ( x ) = 1 n h n j = 1 n K ( x X j h n ) ,
where K is a kernel function and h n is the bandwidth. Instead of estimating f ( x ) , in this paper we estimate the distribution function F by the empirical distribution function
F n ( x ) = n 1 j = 1 n I { X j x } ,
where I { · } is the indicator function. Therefore, our estimator does not involve nonparametric smoothing and avoids choice of the bandwidth h n . Given the expression of
Ψ ( u ) = λ μ c ( 1 λ μ c ) 1 2 π 1 e i s u i s λ 0 e i s u f ( u ) d u 1 i c s 1 λ 0 e i s u f ( u ) d u 1 i c s d s
from (Zhang et al. 2014), we will construct an estimator of Ψ ( u ) by the plug-in device, which replaces λ and μ by their empirical estimators and f ( u ) d u by d F n ( u ) . Similar to (You et al. 2020), we will establish asymptotically normal distributions of the nonparametric estimators of Ψ ( u ) with known λ and with unknown λ . Since our estimators do not involve the kernel smoothing, they are easier to calculate than those in (Zhang et al. 2014) and (You et al. 2020). Furthermore, since there is no bandwidth choice for our procedure, our procedure is expected to have more stable performance. The asymptotic distributions of our estimators allow us to consistently estimate the asymptotic variances of our estimators with the plug-in principle and thus enable interval estimates of the ruin probability.

2. Main Results

The ruin probability Ψ ( u ) in (3) can be rewritten as
Ψ ( u ) = λ μ c ( 1 λ μ c ) φ ( u ) ,
where
φ ( u ) = 1 2 π 1 e i s u i s λ ϕ f ( s ) 1 λ ϕ f ( s ) d s ,
with ϕ f ( s ) = ϕ f ( s ) 1 i c s and ϕ f ( s ) = 0 e i s u f ( u ) d u being the inverse Fourier transform of density f. In the following we consider estimating the ruin probability when λ is known or unknown.

2.1. Estimation with Known Intensity λ

Suppose that a sample { X 1 , X 2 , , X n } of the claim size are observed. When λ is known, a natural estimation method for Ψ ( u ) is to use the plug-in device, based on (4) and (5). This requires us to estimate the inverse Fourier transform ϕ f ( · ) of density f ( · ) .
Note that ϕ f ( s ) = 0 e i s u d F ( u ) . It follows that an empirical type estimator of is given by ϕ e m p ( s ) = 0 e i s x d F n ( s ) = 1 n j = 1 n e i s X j . By (5), a plug-in estimator of φ ( u ) is
φ ˜ n ( u ) = 1 2 π 1 e i s u i s λ ϕ n ( s ) 1 λ ϕ n ( s ) d s ,
where ϕ n ( s ) = ϕ e m p ( s ) 1 i c s . It is easy to see that φ ˜ n ( u ) must be real, because the integrand is conjugate symmetric. Since the integral in (6) is possibly infinite, following (Zhang et al. 2014) we estimate Ψ ( u ) by
Ψ ^ n ( u ) = λ μ ^ n c ( 1 λ μ ^ n c ) φ ^ n ( u ) , u 0 ,
where μ ^ n = 1 n j = 1 n X j , φ ^ n ( u ) = max { min ( M , φ ˜ n ( u ) ) , M } , and M is a large constant. Since the above estimate does not involve the kernel density estimate f ^ n ( x ) , it is simpler and easy to calculate than that of (Zhang et al. 2014). In particular, we do not need to select the bandwidth h n , and thus it is expected that our estimator is more stable than Zhang et al.’s and You et al.’s.
In the following we illustrate the above estimation method is feasible through its asymptotic distribution. To this end, we propose the following assumptions:
(A1)
The density function f ( x ) satisfies that
(1) function f ( x ) is right-continuous at zero and continuously differentiable in ( 0 , ) , and f ( x ) exists almost everywhere. Further, 0 | f ( x ) | d x < , 0 | f ( x ) | d x < , and 0 ( f ( x ) ) 2 d x < ;
(2) E ( X 4 ) < , and for some 0 < α < 1 2 , x 1 + α F ¯ ( x α ) 0 as x , where F ¯ ( x ) = 1 F ( x ) .
(A2)
There exists a constant ρ > 0 such that 1 λ μ c 1 ρ .
The above assumptions are wild and were used in (Zhang et al. 2014) and (You et al. 2020). The following result demonstrates that Ψ ^ n ( u ) is n -consistent and asymptotically normal.
Theorem 1.
Suppose that Assumptions (A1)–(A2) hold, then
n { Ψ ^ n ( u ) Ψ ( u ) } D N ( 0 , S 2 ) ,
where S 2 = [ λ c { 1 + φ ( u ) } ] 2 σ 2 + ( λ μ 1 ) 2 D 2 + C 0 , with C 0 = 1 π λ 2 c 2 { 1 + φ ( u ) } ( λ c μ 1 ) ( 1 e i s u ) { i ϕ f ( s ) μ ϕ f ( s ) } s 2 { 1 λ ϕ f ( s ) } 2 d s and D 2 = λ 2 4 π 2 c 2 E [ 1 e i s u s 2 ϕ f ( s ) e i s X j ( 1 λ ϕ f ( s ) ) 2 d s ] 2 .
Remark 1.
Compared with Theorem 1 in (You et al. 2020), we find that the difference of asymptotic variances between our estimator and their estimator is
n 1 ( λ μ 1 ) 2 ( D 2 D 2 ) ,
where D 2 = λ 2 c 2 { j = 1 3 Q j ( u ) + 2 j = 4 6 Q j ( u ) } (for details, see page 11 of (You et al. 2020)). Hence, the asymptotic variance of our estimator has a much simpler form, which allows us a plug-in estimate, but this is contrary to that of (You et al. 2020) which is hard to estimate directly. It fact, they proposed ro estimate the variance with resampling methods such as the bootstrap.

2.2. Estimation with Unknown Intensity λ

Theorem 1 is useful if λ is known. However, the true value of λ is not available in practice. Suppose there are κ n contiguous observational time intervals with same length of d. For j = 1 , , κ n , let { N j d } be the total number of claims from the starting time to the jth observational time interval with length d > 0 . Then an unbiased estimator of Poisson intensity parameter λ is
λ ^ κ n = N κ n d / ( κ n d ) ,
where N 0 d = 0 . By (7), a natural estimator of Ψ ( u ) (see Remark 1 of (Zhang et al. 2014)) is
Ψ ^ n ( u ) = λ ^ κ n μ ^ n c ( 1 λ ^ κ n μ ^ n c ) φ ^ n ( u ) , u 0 ,
where φ ^ n ( u ) = ( M φ ˜ n ( u ) ) ( M ) for a large constant M , and
φ ˜ n ( u ) = 1 2 π 1 e i s u i s λ ^ κ n ϕ n ( s ) 1 λ ^ κ n ϕ n ( s ) d s .
Additionally, we need the following assumption for establishing the asymptotic distribution of Ψ ^ n ( u ) .
(A3)
For integer m 2 , E | X | m m ! H m 2 σ 2 / 2 , where H is a positive constant.
Assumption (A3) is the classical Cramér condition for all moments of X. It was used in (You et al. 2020). The following theorem indicates Ψ ^ n ( u ) converges to the ruin probability at the best parametric rate of n .
Theorem 2.
Suppose Assumptions (A1)–(A3) hold. If 0 λ lim n n / κ n < , then
n { Ψ ^ n ( u ) Ψ ( u ) } D N ( 0 , S 2 ) ,
where S 2 = S 2 + μ c ( 1 + φ ( u ) ) ( 1 λ μ c ) a ( f , λ ) 2 λ λ d and
a ( f , λ ) = 1 2 π 1 e i s u i s ϕ f ( s ) { 1 λ ϕ f ( s ) } 2 d s .
Remark 2.
Both Ψ ^ n ( u ) and Ψ ^ n ( u ) are n -consistent, but the latter has larger variance, which is understandable because estimating λ brings an additional variation.
Given an initial surplus u, Theorem 2 allows us to construct an 100 ( 1 α ) % confidence interval for Ψ ( u ) :
Ψ ^ ( u ) ± n 1 / 2 z 1 α / 2 S ^ ,
where z 1 α / 2 is the 100 ( 1 α / 2 ) th percentile of the standard normal distribution, and S ^ is a consistent estimator of S , for example, from the plug-in device.

2.3. Computation on Ruin Probability

To obtain the estimate Ψ ^ n ( u ) of ruin probability in (10), one needs to calculate φ ˜ n ( u ) in (11). Let G n ( s ) = λ ^ κ n ϕ n ( s ) i s { 1 λ ^ κ n ϕ n ( s ) } . Then
φ ˜ n ( u ) = 1 2 π G n ( s ) d s 1 2 π e i s u G n ( s ) d s .
The above equation is well defined, since the integrals are finite. Let H n ( u ) = 1 2 π e i s u G n ( s ) d s be the Fourier transform of G n ( s ) . Then φ ˜ n ( u ) = H n ( 0 ) H n ( u ) . Therefore, our major task is to calculate H n ( u ) , which is easy to implement in standard software.

2.4. Proofs of Theorems

To facilitate the proofs, we first introduce some technical lemmas. Then we give the proofs of theorems.
Lemma 1.
If Assumption (A1) holds. Then, for any ε > 0 , P ( sup s | ϕ f ( s ) ϕ e m p ( s ) i s | > ε ) = o ( n 1 ) .
Proof of Lemma 1.
By the inequality | e i s x 1 | | s x | , we have
ϕ f ϕ e m p ( s ) i s = 0 ( e i s x 1 ) d F ( x ) 0 ( e i s x 1 ) d F n ( x ) i s = 0 0 x e i s y d y d F ( x ) 0 0 x e i s y d y d F n ( x ) i s = 0 ( 1 F ( x ) ) e i s x d x 0 ( 1 F n ( x ) ) e i s x d x i s .
Hence, for ε > 0 , by Lemma 1 of (Zhang et al. 2014), we obtain that
P sup s ϕ f ϕ e m p ( s ) i s > ε P 0 | F ( x ) F n ( x ) | d x > ε = o ( n 1 ) .
Lemma 2.
Let A n ( s ) = λ c ϕ f ( s ) ϕ e m p ( s ) i s . If Assumptions (A1)–(A2) hold, then E | A n ( s ) | 2 d s = O ( n 1 ) and E | A n ( s ) | 4 d s = O ( n 2 ) .
Proof of Lemma 2.
Note that ϕ e m p ( s ) = 1 n j = 1 n e i s X j , it follows that
ϕ f ( s ) ϕ e m p ( s ) i s = 1 n j = 1 n ϕ f ( s ) e i s X j i s .
Taking expectation, we obtain that E { ϕ e m p ( s ) ϕ f ( s ) / ( i s ) } = 0 and
V a r { ϕ e m p ( s ) ϕ f ( s ) / ( i s ) } = n 1 V a r { ϕ f ( s ) e i s X j i s } = ϕ f 2 ( s ) ϕ f ( 2 s ) n s 2 .
By Fubini’s theorem, we have
E | A n ( s ) | 2 d s = E λ c ϕ f ( s ) ϕ e m p ( s ) i s 2 d s = ( λ c ) 2 1 n ϕ f 2 ( s ) ϕ f ( 2 s ) s 2 d s = O ( n 1 ) .
Again, using Fubini’s theorem, we establish that
E | A n ( s ) | 4 d s = E λ c ϕ f ( s ) ϕ e m p ( s ) i s 4 d s = V a r ( λ c ϕ f ( s ) ϕ e m p ( s ) i s ) 2 + ( V a r ( λ c ϕ f ( s ) ϕ e m p ( s ) i s ) ) 2 d s = ( λ c ) 4 1 n 2 V a r ( ϕ f ( s ) e i s X j i s ) 2 + ( V a r ( ϕ f ( s ) e i s X j i s ) ) 2 d s .
By independence of the sample points, we have
E | A n ( s ) | 4 d s = ( λ c ) 4 1 n 2 E ϕ f ( s ) e i s X j i s 4 d s = ( λ c ) 4 1 n 2 6 ϕ f 2 ( s ) ϕ f ( 2 s ) 4 ϕ f ( s ) ϕ f ( 3 s ) + ϕ f ( 4 s ) 3 ϕ f 4 ( s ) s 4 d s = O ( n 2 ) .
Lemma 3.
Suppose that Assumptions (A1)–(A2) hold. Then
(i) 
(Consistency) μ ^ n μ P 0 , φ ^ n ( u ) φ ( u ) P 0 , and λ ^ κ n λ P 0 .
(ii) 
(Asymptotic normality) n ( μ ^ n μ ) D N ( 0 , σ 2 ) , n ( λ ^ κ n λ ) D N ( 0 , λ d ) , and n ( φ ^ n ( u ) φ ( u ) ) D N ( 0 , D 2 ) .
Proof of Lemma 3.
By the property of the Poisson process N and the size of claim X, it is easy to obtain μ ^ n μ P 0 , λ ^ κ n λ P 0 , n ( μ ^ n μ ) D N ( 0 , σ 2 ) , and n ( λ ^ κ n λ ) D N ( 0 , λ d ) . In the following we show that φ ^ n ( u ) φ ( u ) P 0 .
First, combining (5) and (6) leads to
φ ˜ n ( u ) φ ( u ) = 1 2 π 1 e i s u i s A n ( s ) ( 1 λ ϕ f ( s ) ) 2 d s + 1 2 π 1 e i s u i s A n 2 ( s ) ( 1 λ ϕ f ( s ) + A n ( s ) ) ( 1 λ ϕ f ( s ) ) 2 d s I n ( u ) + II n ( u ) .
Let ϵ 1 and ϵ 2 be two small positive constants such that φ ( u ) + ϵ 1 + ϵ 2 < M . Define V n = { | I n ( u ) | < ϵ 1 , | II n ( u ) | < ϵ 2 } and its complementary set V ¯ n . Then
φ ^ n φ ( u ) = { φ ^ n φ ( u ) } I V n + { φ ^ n φ ( u ) } I V ¯ n = { φ ˜ n φ ( u ) } I V n + { φ ^ n φ ( u ) } I V ¯ n = i = 1 4 H n , i ( u ) ,
where H n , 1 ( u ) = I n ( u ) , H n , 2 ( u ) = II n ( u ) I V n , H n , 3 = I n ( u ) I V ¯ n and H n , 4 = ( φ ^ n ( u ) φ ( u ) ) I V ¯ n . Using the inequality | ϕ f 1 | | i μ s | , Hölder’s inequality and Assumption (A2), we obtain that
| I n ( u ) | 1 2 π | 1 e i s u i s | | A n ( s ) | | 1 λ ϕ f ( s ) | 2 d s 1 2 ρ 2 π | 1 e i s u i s | | A n ( s ) | d s 1 2 ρ 2 π ( | 1 e i s u i s | 4 3 d s ) 3 4 ( | A n ( s ) | 4 d s ) 1 4 C ( | A n ( s ) | 4 d s ) 1 4 .
By Markov’s inequality and Lemma 2, we have
P ( | I n ( u ) | ϵ 1 ) E | I n ( u ) | 4 ϵ 1 4 C · E | A n ( s ) | 4 d s = O ( 1 n 2 ) .
For 0 < δ < ρ , let B n , δ = { sup s | A n ( s ) | δ } , and let B ¯ n , δ be its complementary set. On the set B n , δ , we have
| II n ( u ) | 1 2 π | 1 e i s u i s | A n 2 ( s ) ( 1 λ ϕ f ( s ) + A n ( s ) ) ( 1 λ ϕ f ( s ) ) 2 d s C 1 e i s u i s | A n ( s ) | 2 d s .
Then, the Cauchy-Schwarz inequality, | II n ( u ) | C | A n ( s ) | 2 d s and | II n ( u ) | C ( | A n ( s ) | 4 d s ) 1 2 . Therefore,
E | H n , 2 | = E [ | H n , 2 | I ( V n B n , δ ) ] + E [ | H n , 2 | I ( V n B ¯ n , δ ) ] C · E | A n ( s ) | 2 d s + ε 2 P ( B ¯ n , δ ) = O ( n 1 ) .
In order to get the bounds of H n , 3 and H n , 4 , we need to study probability P ( V ¯ n ) . By Markov’s inequality, we have
P ( { | II n ( u ) | ε 2 } B n , δ ) E [ | II n ( u ) | 2 I ( B n , δ ) ] ε 2 2 C · E | A n ( s ) | 4 d s = O ( 1 n 2 ) ,
which, combined with Lemma 1, leads to
P ( | II n ( u ) | ε 2 ) = P ( { | II n ( u ) | ε 2 } B n , δ ) + P ( { | II n ( u ) | ε 2 } B ¯ n , δ ) P ( { | II n ( u ) | ε 2 } B n , δ ) + P ( B ¯ n , δ ) = O ( 1 n 2 ) + o ( 1 n ) .
This, combined with (15) and (18), gives us that
P ( V ¯ n ) P ( | I n ( u ) | ε 1 ) + P ( | II n ( u ) | ε 2 ) P ( | I n ( u ) | ε 1 ) + P ( { | II n ( u ) | ε 2 } B n , δ ) + P ( B ¯ n , δ ) = o ( 1 n ) + O ( 1 n 2 ) .
Hence, applying Hölder’s inequality and Lemma 2, we obtain that
E | H n , 3 ( u ) | E [ | I n ( u ) | I V ¯ n | ] = ( E | I n ( u ) | 4 ) 1 4 P ( V ¯ n ) 3 4 = o ( 1 n ) + O ( 1 n 2 ) .
By (19), we get
E | H n , 4 ( u ) | ( M + φ ( u ) ) P ( V ¯ n ) = o ( 1 n ) + O ( 1 n 2 ) .
It follows from (14), (15), (17), (20) and (21) and Markov’s inequality that φ ^ n ( u ) φ ( u ) P 0 .
Next, we show that n ( φ ^ n ( u ) φ ( u ) ) D N ( 0 , D 2 ) . By (14), n { φ ^ n ( u ) φ ( u ) } = n i = 1 4 H n , i . It follows from Markov’s inequality that for any ε > 0 ,
P ( | n i = 2 4 H n , i | > ε ) n i = 2 4 E | H n , i | ε ,
which together with (15), (17), (20) and (21) implies that n i = 2 4 H n , i P 0 . Therefore,
n { φ ^ n ( u ) φ ( u ) } = n H n , 1 + o P ( 1 ) = n 1 2 π 1 e i s u i s A n ( s ) ( 1 λ ϕ f ( s ) ) 2 d s + o P ( 1 ) = λ 2 c n π j = 1 n 1 e i s u s 2 ϕ f ( s ) e i s X j ( 1 λ ϕ f ( s ) ) 2 d s + o P ( 1 ) .
Applying Fubini’s theorem, we establish that
E [ n H n , 1 ] = E 1 2 n π λ c j = 1 n 1 e i s u s 2 ϕ f ( s ) e i s X j ( 1 λ ϕ f ( s ) ) 2 d s = 1 2 n π λ c j = 1 n 1 e i s u s 2 E [ ϕ f ( s ) e i s X j ] ( 1 λ ϕ f ( s ) ) 2 d s = 0
and
V a r ( n H n , 1 ) = ( 1 2 n π λ c ) 2 j = 1 n V a r 1 e i s u s 2 ϕ f ( s ) e i s X j ( 1 λ ϕ f ( s ) ) 2 d s = λ 2 4 π 2 c 2 V a r 1 e i s u s 2 ϕ f ( s ) e i s X j ( 1 λ ϕ f ( s ) ) 2 d s = λ 2 4 π 2 c 2 E 1 e i s u s ϕ f ( s ) e i s X j s ( 1 λ ϕ f ( s ) ) 2 d s 2 λ 2 4 π 2 c 2 E 1 e i s u s 2 d s ϕ f ( s ) e i s X j s ( 1 λ ϕ f ( s ) ) 2 2 d s ,
where the last inequality is from the H o ¨ lder inequality. Therefore,
V a r ( n H n , 1 ) λ 2 4 π 2 c 2 1 e i s u s 2 d s 1 s 2 ( 1 λ ϕ f ( s ) ) 4 E [ ϕ f ( s ) e i s X j ] 2 d s = λ 2 4 π 2 c 2 ( 1 e i s u ) 2 s 2 d s ϕ f ( 2 s ) ϕ f 2 ( s ) s 2 ( 1 λ ϕ f ( s ) ) 4 d s .
Using the inequality | ϕ f ( s ) 1 | | i μ s | and Assumption (A2), we obtain that
( 1 λ ϕ f ( s ) ) 4 ( 1 | λ ϕ f ( s ) | ) 4 = ( 1 | λ c ϕ f ( s ) 1 i s | ) 4 ( 1 μ λ c ) 4 ρ 4 .
Then
V a r ( n H n , 1 ) λ 2 4 ρ 4 π 2 c 2 ( 1 e i s u ) 2 s 2 d s ϕ f ( 2 s ) ϕ f 2 ( s ) s 2 d s < .
Therefore, by the CLT and Slutsky’s theorem, we obtain that
n ( φ ^ n ( u ) φ ( u ) ) N ( 0 , D 2 ) ,
where D 2 = λ 2 4 π 2 c 2 E [ 1 e i s u s 2 ϕ f ( s ) e i s X j ( 1 λ ϕ f ( s ) ) 2 d s ] 2 . □
Lemma 4.
Let ξ n 1 = ( λ ^ κ n λ ) a ( f , λ ) , ξ n 2 = n 1 j = 1 n Z j and ξ n = ξ n 1 + ξ n 2 , where
Z j = λ 2 π c 1 e i s u s 2 [ 1 λ ϕ f ( s ) ] 2 [ e i s X j ϕ f ( s ) ] d s .
Suppose that assumptions (A1)–(A3) hold. Then P { | ξ n | > C + o p ( n 1 / 2 ) } = o ( n 1 / 2 ) , where C is any given constant.
Proof of Lemma 4.
The proof is similar to Lemma A.2 in (You et al. 2020). □
Proof of Theorem 1.
Our argument for the proof is to first establish the Bahardur representation of our estimator, based on its explicit formula in (7). Then the result of theorem follows from the Bahardur’s representation. By (4) and (7), we have
n Ψ ^ n ( u ) Ψ ( u ) = n λ c { 1 + φ ( u ) } ( μ ^ n μ ) + ( λ c μ ^ n 1 ) { φ ^ n ( u ) φ ( u ) } = λ c { 1 + φ ( u ) } n ( μ ^ n μ ) + ( λ c μ 1 ) n { φ ^ n ( u ) φ ( u ) } + λ c n ( μ ^ n μ ) { φ ^ n ( u ) φ ( u ) } .
By Lemma 3 and Slutsky’s theorem, the last term of (23) converges to zero in probability as n . Thus,
n { Ψ ^ n ( u ) Ψ ( u ) } = λ c { 1 + φ ( u ) } n ( μ ^ n μ ) + ( λ c μ 1 ) n { φ ^ n ( u ) φ ( u ) } + o p ( 1 ) .
This, combined with (22), leads to the following Bahardur representation:
n { Ψ ^ n ( u ) Ψ ( u ) } = λ c { 1 + φ ( u ) } n ( μ ^ n μ ) + ( λ c μ 1 ) n H n , 1 + o p ( 1 ) = 1 n j = 1 n λ c [ ( 1 + φ ( u ) ) ( X j μ ) + ( λ c μ 1 ) 1 2 π ( 1 e i s u ) ( e i s X j ϕ f ( s ) ) s 2 ( 1 λ ϕ f ( s ) ) 2 d s ] + o p ( 1 ) = 1 n j = 1 n B j + o p ( 1 ) ,
where B j = λ c ( 1 + φ ( u ) ) ( X j μ ) + ( λ c μ 1 ) 1 2 π ( 1 e i s u ) ( e i s X j ϕ f ( s ) ) s 2 ( 1 λ ϕ f ( s ) ) 2 d s . Note that E [ B 1 ] = E [ 1 n j = 1 n B j ] = ( λ c μ 1 ) E [ H n , 1 ] = 0 and
Var ( B 1 ) = λ 2 c 2 Var { 1 + φ ( u ) } ( X j μ ) + λ 2 c 2 Var ( λ c μ 1 ) 1 2 π ( 1 e i s u ) ( e i s X j ϕ f ( s ) ) s 2 ( 1 λ ϕ f ( s ) ) 2 d s + 2 λ 2 c 2 C o v { 1 + φ ( u ) } ( X j μ ) , ( λ c μ 1 ) 1 2 π ( 1 e i s u ) ( e i s X j ϕ f ( s ) ) s 2 ( 1 λ ϕ f ( s ) ) 2 d s = ( λ c ( 1 + φ ( u ) ) ) 2 σ 2 + n ( λ c μ 1 ) 2 Var 1 2 π 1 e i s u i s A n ( s ) ( 1 λ ϕ f ( s ) ) 2 d s + 1 π λ 2 c 2 ( 1 + φ ( u ) ) ( λ c μ 1 ) E ( X j μ ) ( 1 e i s u ) ( e i s X j ϕ f ( s ) ) s 2 ( 1 λ ϕ f ( s ) ) 2 d s .
Then, by simple algebra, we arrive at
Var ( B 1 ) = λ c { 1 + φ ( u ) } 2 σ 2 + ( λ c μ 1 ) 2 Var ( n H n , 1 ) + 1 π λ 2 c 2 { 1 + φ ( u ) } ( λ c μ 1 ) ( 1 e i s u ) ( i ϕ f ( s ) μ ϕ f ( s ) ) s 2 ( 1 λ ϕ f ( s ) ) 2 d s . = λ c { 1 + φ ( u ) } 2 σ 2 + ( λ c μ 1 ) 2 D 2 + 1 π λ 2 c 2 { 1 + φ ( u ) } ( λ c μ 1 ) ( 1 e i s u ) ( i ϕ f ( s ) μ ϕ f ( s ) ) s 2 ( 1 λ ϕ f ( s ) ) 2 d s = S 2 .
It follow from the central limit theorem that 1 n j = 1 n B j N ( 0 , S 2 ) . Applying Slutsky’s theorem and (25), we obtain that n { Ψ ^ n ( u ) Ψ ( u ) } D N ( 0 , S 2 ) , as n .
Proof of Theorem 2.
Our idea of the proof is to apply Lemma 4, and to follow the argument for the proof of Theorem 2.2 in (You et al. 2020). By (4) and (10), we have
Ψ ^ n ( u ) Ψ ( u ) = λ ^ κ n c ( μ ^ n μ ) + μ c ( λ ^ κ n λ ) 1 + φ ( u ) ( 1 λ ^ κ n μ ^ n c ) { φ ^ n ( u ) φ ( u ) } ,
where φ ^ n ( u ) = ( M φ ˜ n ( u ) ) ( M ) . To establish the asymptotic distribution of Ψ ^ n ( u ) Ψ ( u ) , we need to study φ ˜ n ( u ) . By (6) and (11),we have
φ ˜ n ( u ) φ ˜ n ( u ) = ( λ ^ κ n λ ) 1 2 π 1 e i s u i s ϕ n ( s ) ( 1 λ ϕ n ( s ) ) ( 1 λ ^ κ n ϕ n ( s ) ) d s .
Applying Lemmas 1 and 3, we obtain that sup s | ϕ n ( s ) ϕ f ( s ) | = o p ( 1 ) and
1 2 π 1 e i s u i s ϕ n ( s ) ( 1 λ ϕ n ( s ) ) ( 1 λ ^ κ n ϕ n ( s ) ) d s = a ( f , λ ) + o p ( 1 ) ,
where a ( f , λ ) = 1 2 π 1 e i s u i s ϕ f ( s ) { 1 λ ϕ f ( s ) } 2 d s . Therefore,
φ ˜ n ( u ) φ ˜ n ( u ) = ( λ ^ κ n λ ) a ( f , λ ) + o p ( n 1 / 2 ) .
Combining (13), (16) and (26) and Lemma 2 leads to
φ ˜ n ( u ) φ ( u ) = ( λ ^ κ n λ ) a ( f , λ ) I n ( u ) + r n ,
where r n = o p ( n 1 / 2 ) . By (15), we see that I n ( u ) = O p ( n 2 ) . Hence, φ ˜ n ( u ) φ ( u ) = o p ( 1 ) . Let ξ n = ( λ ^ κ n λ ) a ( f , λ ) I n ( u ) . Then
ξ n = ( λ ^ κ n λ ) a ( f , λ ) λ 2 π c n 1 j = 1 n 1 e i s u s 2 [ 1 λ ϕ f ( s ) ] 2 [ e i s X j ϕ f ( s ) ] d s ξ n 1 + ξ n 2 ,
φ ˜ n ( u ) = ξ n + φ ( u ) + r n , and P ( | φ ˜ n ( u ) | > M ) P { | ξ n | > ( M | φ ( u ) | | r n | } . Then, by Lemma 4 and the proof of Theorem 2.2 in (You et al. 2020), it is easy to obtain n ( Ψ ^ n ( u ) Ψ ( u ) ) D N ( 0 , S 2 ) .

Author Contributions

Conceptualization, L.C.; methodology, J.J.; validation, L.C. and J.J.; formal analysis, J.J., Y.G. and H.Y.; investigation, L.C.; resources, L.C.; writing—original draft preparation, J.J.; writing—review and editing, L.C., Y.G., H.Y. and J.J.; visualization, L.C.; supervision, J.J.; project administration, J.J.; funding acquisition, H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Asmussen, Søren, and Hansjörg Albrecher. 2010. Ruin Probabilities, 2nd ed. Singapore: World Scientific. [Google Scholar]
  2. Baumgartner, Benjamin, and Riccardo Gatto. 2010. A bootstrap test for the probability of ruin in the compound poisson risk process. ASTIN Bulletin 40: 241–55. [Google Scholar] [CrossRef] [Green Version]
  3. Bowers, Newton, Hans Gerber, James Hickmann, Donald Jones, and Cecil Nesbitt. 1997. Actuarial Mathematics, 2nd ed. Schaumburg: The Society of Actuaries. [Google Scholar]
  4. Cai, Jun, and Yanhong Wu. 1997. Some improvements on the lundberg bound for the ruin probability. Statistics and Probability Letters 33: 395–403. [Google Scholar] [CrossRef]
  5. Croux, Kristof, and Noel Veraverbeke. 1990. Nonparametric estimators for the probability of ruin. Insurance: Mathematics and Economics 9: 127–30. [Google Scholar] [CrossRef]
  6. Dickson, David. 1994. An upper bound for the probability of ultimate ruin. Scandinavian Actuarial Journal 2: 131–8. [Google Scholar] [CrossRef]
  7. Feng, Minyu, Liang-Jian Deng, Feng Chen, Matjaž Perc, and Jürgen Kurths. 2020. The accumulative law and its probability model: An extension of the Pareto distribution and the log-normal distribution. Proceedings of the Royal Society A, 47620200019. [Google Scholar] [CrossRef] [PubMed]
  8. Frees, Edward W. 1986. Nonparametric estimation of the probability of ruin. ASTIN Bulletin 16: 81–90. [Google Scholar] [CrossRef] [Green Version]
  9. Gerber, Hans U. 1979. An Introduction to Mathematical Risk Theory, no.8. In Huebner Foundation Monograph Series. Philadelphia: S. S. Huebner Foundation for Insurance Education. [Google Scholar]
  10. Landriault, David, and Gordon Willmot. 2008. On the gerber-shiu discounted penalty function in the sparre andersen model with an arbitrary interclaim time distribution. Insurance: Mathematics and Economics 42: 600–8. [Google Scholar] [CrossRef]
  11. Mnatsakanov, Robert, L. L. Ruymgaart, and Frits H. Ruymgaart. 2008. Nonparametric estimation of ruin probabilities given a random sample of claims. Mathematical Methods of Statistics 17: 35–43. [Google Scholar] [CrossRef]
  12. Politis, Konstadinos. 2003. Semiparametric estimation for non-ruin probabilities. Scandinavian Actuarial Journal 1: 75–96. [Google Scholar] [CrossRef]
  13. Pitts, Susan M. 1994. Nonparametric estimation of compound distributions with applications in insurance. Annals of the Institute of Statistical Mathematics 46: 537–55. [Google Scholar]
  14. You, Honglong, Junyi Guo, and Jiancheng Jiang. 2020. Interval estimation of the ruin probability in the classical compound Poisson risk model. Computational Statistics and Data Analysis 144: 106890. [Google Scholar] [CrossRef]
  15. Zhang, Zhimin, Hailiang Yang, and Hu Yang. 2014. On a nonparametric estimator for ruin probability in the classical risk model. Scandinavian Actuarial Journal 4: 309–38. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, Y.; Chen, L.; Jiang, J.; You, H. Nonparametric Estimation of the Ruin Probability in the Classical Compound Poisson Risk Model. J. Risk Financial Manag. 2020, 13, 298. https://doi.org/10.3390/jrfm13120298

AMA Style

Gao Y, Chen L, Jiang J, You H. Nonparametric Estimation of the Ruin Probability in the Classical Compound Poisson Risk Model. Journal of Risk and Financial Management. 2020; 13(12):298. https://doi.org/10.3390/jrfm13120298

Chicago/Turabian Style

Gao, Yuan, Lingju Chen, Jiancheng Jiang, and Honglong You. 2020. "Nonparametric Estimation of the Ruin Probability in the Classical Compound Poisson Risk Model" Journal of Risk and Financial Management 13, no. 12: 298. https://doi.org/10.3390/jrfm13120298

APA Style

Gao, Y., Chen, L., Jiang, J., & You, H. (2020). Nonparametric Estimation of the Ruin Probability in the Classical Compound Poisson Risk Model. Journal of Risk and Financial Management, 13(12), 298. https://doi.org/10.3390/jrfm13120298

Article Metrics

Back to TopTop