Next Article in Journal
A Fully Resolved Computational Fluid Dynamics Study of the Boundary Layer Flow of an Aqueous Nanoliquid Comprising Gyrotactic Microorganisms over a Stretching Sheet: The Validity of Conventional Similarity Models
Previous Article in Journal
Improved Mixed-Integer Linear Programming Model for Short-Term Scheduling of the Pressing Process in Multi-Layer Printed Circuit Board Manufacturing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Speed of Convergence of the Threshold Estimator of Ruin Probability under the Tempered α-Stable Lévy Subordinator

1
School of Mathematical Sciences, Qufu Normal University, Jining 273165, China
2
School of Statistics, Qufu Normal University, Jining 273165, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(21), 2654; https://doi.org/10.3390/math9212654
Submission received: 14 September 2021 / Revised: 13 October 2021 / Accepted: 18 October 2021 / Published: 20 October 2021

Abstract

:
In this paper, a nonparametric estimator of ruin probability is introduced in a spectrally negative Lévy process where the jump component is a tempered α -stable subordinator. Given a discrete record of high-frequency data, a threshold technique is proposed to estimate the mean of the jump size and use the Fourier transform and the Pollaczek–Khinchin formula to construct the estimator of ruin probability. The convergence rate of the integrated squared error for the estimator is studied.

1. Introduction

In actuarial science, the statistical inference of ruin probability has received much attention from scholars. Many works have been contributed to parametric and nonparametric estimation of ruin probability. See, for example, Croux and Veraverbeke [1], Frees [2], Mnatsakanov et al. [3], Pitts [4], Politis [5], and Veraverbeke [6]. In recent years, many nice results have been obtained by actuarial scholars, such as Huang et al. [7], Li et al. [8], You et al. [9], Zhang and Yang [10], Zhang and Yang [11], Zhang [12], and Zhang and Yang [13]. As an extension of ruin probability, the Gerber–Shiu function has been introduced and studied for its statistical properties. Interested readers can refer to Su and Yu [14,15], Yang et al. [16], Zhang and Su [17], Su et al. [18], Zhang and Su [19], Zhang [20], Shimizu [21], and Shimizu and Zhang [22], among others.
In Asmussen and Albrecher [23], an analytic (or probabilistic) approach was suggested, and it needs much more detailed information about the risk model, such as the claim size distribution. However, in practical situations, it is not easy to obtain the specific distribution information. Instead, one observes the surplus process at some discrete time points. Then, a statistical methodology can be directly used to estimate the claim size distribution with the observed data. In Zhang and Yang [10], a nonparametric estimator of ruin probability was proposed, based on the Pollaczek–Khinchin formula and the Fourier transform in a pure-jump Lévy risk model. This estimation approach was extended by Zhang and Yang [11] to a spectrally negative Lévy risk model. Subsequently, Shimizu and Zhang [22] estimated the Gerber–Shiu function for an insurance surplus process driven by a Lévy subordinator. In Zhang and Yang [10] and Comte and Genon-Catalot [24], they considered high-frequency sampling with n discrete time observations of step width h n > 0 and derived asymptotics under the framework that h n 0 and n h n .
In the present work, our interest is to estimate ruin probability for a spectrally negative Lévy risk model under the above framework. Assume that the surplus of the risk model can be observed at a sequence of discrete time points { t k n = k h n , k = 0 , 1 , 2 , 3 , } with h n 0 being the length of the sampling interval. Without observing the jump and diffusion parts of the risk model, it is challenging to estimate ruin probability, since it depends on both parts in a spectrally negative Lévy risk model. In Mancini [25,26,27,28] and Shimizu [29,30], they developed a threshold technique for identifying the times when jump sizes exceed a suitably defined threshold. Using the threshold technique and the Fourier transform, an estimator of ruin probability is constructed, and the convergence rate of its integrated squared error is obtained.
The remainder of this paper is organized as follows. In Section 2, the risk model, as well as some assumptions for the asymptotic theory are introduced. In Section 3, some estimators are suggested, based on the Fourier transform and the threshold technique. In Section 4, the convergence rate of our estimators is established. In Section 5, we conclude this paper.

2. Preliminaries

2.1. Risk Model and Some Assumptions

A spectrally negative Lévy process is specified by:
Y t = c t + σ W t J t , t 0 ,
where c > 0 is a parameter, σ > 0 represents the perturbation coefficient, W t is a standard Brownian motion, and J t is a subordinator. Suppose that W t and J t are independent of each other. Then, the characteristic exponent of Y t is given by:
ψ Y ( s ) = 1 t ln ( E [ e i s Y t ] ) = i c s 1 2 σ 2 s 2 0 ( 1 e i s x ) ν ( d x ) ,
where ν is the Lévy measure on ( 0 , ) . By Sato [31], it can be rewritten as:
J t = L t + M t * ,
where L t is the sum of jumps over [ 0 , t ] with the jump size larger than one, and M t * is the sum of jumps over [ 0 , t ] with the jump size less than one. Specifically, L t = 0 t x > 1 x μ ( d s , d x ) = k = 1 N t γ k , where μ is the Poisson random measure of J t such that E [ μ ( d s , d x ) ] = ν ( d x ) d s , N t is a Poisson process, and γ 1 , γ 2 , γ 3 , are i.i.d. random variables, that is L t is a compound Poisson process representing the jumps of J t with the jump size larger than one. Process M t * admits decomposition, M t * = b t + M t , where b = x 1 x ν ( d x ) and M t = 0 t x 1 x μ ˜ ( d s , d x ) is a martingale with μ ˜ ( d s , d x ) = μ ( d s , d x ) ν ( d x ) d s being the compensated measure of μ ( d s , d x ) . It is known that M t is a square integrable martingale with infinite activity of the jump such that E [ M t ] = 0 and V a r [ M t ] = t x 1 x 2 ν ( d x ) < . Suppose that γ k , N t , and M t are independent of each other.
Let u > 0 be the initial surplus of an insurance company. Then, the surplus at time t can be modeled by:
U t = u + Y t = u + c t + σ W t J t , t 0 ,
where c is the rate of the premium, σ represents the perturbation coefficient, J t denotes the claim payments and other expenses in insurance businesses, and W t is a perturbation.

2.2. Ruin Probability and Its Fourier Transform

The infinite-time horizon ruin probability Φ ( u ) is defined as:
Φ ( u ) = P inf 0 t < U t 0 | U 0 = u .
By Equation (1) in Zhang and Yang [11], Φ ( u ) admits the following representation:
Φ ( u ) = 1 ( 1 ρ ) i = 0 ρ i ( G ( i + 1 ) * * H ( i ) * ) ( u ) , u > 0 ,
where ρ = μ 1 c , μ 1 = 0 x ν ( d x ) , H ( x ) = 1 μ 1 0 x ν ( y , ) d y , and G is determined by the Fourier transform 0 e i s x d G ( x ) = c { c σ 2 2 i s } 1 . Setting Φ ( u ) 0 for u < 0 , Zhang and Yang [11] obtained the Fourier transform of Φ ( u ) :
F Φ ( s ) = 0 e i s u Φ ( u ) d u = σ 2 2 i s + 1 i s 0 ( e i s x 1 ) ν ( d x ) μ 1 i c s + σ 2 2 s 2 0 ( e i s x 1 ) ν ( d x ) = i c s ( 1 ρ ) + ψ Y ( s ) i s ψ Y ( s ) ,
Once an estimator of F Φ ( s ) is available, Φ ( u ) can be estimated by the inverse Fourier transform.

3. Estimation of Ruin Probability

Suppose that a discrete sample Y n = { Y t i n | t i n = i h n ; i = 0 , 1 , 2 , n } can be observed. Let Z i = Y t i n Y t i 1 n , h n = t i n t i 1 n > 0 , lim n h n = 0 , and lim n n h n = . Our interest is to estimate Φ ( u ) by Z 1 , Z 2 , , Z n when Lévy measure ν and perturbation coefficient σ are unknown.
If one can estimate ρ and ψ Y ( s ) in (6), then F Φ ( s ) can be estimated with the plug-in device. Inspired by Zhang and Yang [10,11] and You and Yin [32], we define the estimator of ψ Y ( s ) :
ψ ^ Y ( s ) = 1 h n ( 1 n k = 1 n e i s Z k 1 ) .
To estimate ρ = μ 1 c , we need to estimate μ 1 , the mean of J 1 . Zhang and Yang [11] proposed to estimate μ 1 , by:
μ ^ 1 * = 1 n h n k = 1 n ( c h n Z k ) .
Note that c h n Z k = ( J t k n J t k 1 n ) σ ( W t k n W t k 1 n ) . Ideally, we hope that the estimator of μ 1 is 1 n h n k = 1 n ( J t k n J t k 1 n ) , but we cannot observe a discrete sample J n = { J t i n | t i n = i h n ; i = 0 , 1 , 2 , n } . To this end, we introduce a threshold technique. Motivated by Shimizu [29,30] and Mancini [26,27], we introduce the filter:
D k n : = { ω Ω : ( c h n Z k ) > r n } ,
where r n > 0 is a suitable threshold parameter dependent on n such that lim h n 0 r n = 0 . Let C k n : = { ω Ω : ( c h n Z k ) r n } be the complement of D k n . By (9), if c h n Z k > r n , we can detect the existence of a jump in an interval ( t k 1 n , t k n ] , and then, we take c h n Z k as an approximation to J t k n J t k 1 n . This leads to a natural estimate of μ 1 :
μ ^ 1 = k = 1 n ( c h n Z k ) I D k n n h n .
Then, ρ is estimated by:
ρ ^ = k = 1 n ( c h n Z k ) I D k n c n h n .
Combining (6), (7), and (11) leads to our estimator:
F ^ Φ ( s ) = i c s ( 1 ρ ^ ) + ψ ^ Y ( s ) i s ψ ^ Y ( s ) .
Note that the above estimate has no definition at s = 0 . When s 0 , ψ ^ Y ( s ) 0 , and thus, F ^ Φ ( s ) may behave erratically. Applying the inverse Fourier transform and removing a small neighborhood of s = 0 , we propose to estimate Φ ( u ) by:
Φ ^ ( u ) = 1 2 π m n M n e i u s F ^ Φ ( s ) d s + 1 2 π M n m n e i u s F ^ Φ ( s ) d s , u > 0 ,
where m n and M n are positive threshold numbers such that m n 0 and M n as n .

4. Asymptotic Properties of Estimators

In this section, the asymptotic properties of ρ ^ and Φ ^ ( u ) are studied. For the ease of exposure, we first introduce some notations. For integer k = 1 , 2 , μ k : = 0 x k ν ( d x ) . For any two positive sequences { x n } n = 1 and { y n } n = 1 , x n y n means that x n C y n for some constant C and large index n. For any function f ( x ) with support ( 0 , ) , define f 2 = 0 f 2 ( x ) d x . Let μ L = 1 x ν ( d x ) , μ γ = E [ γ k ] , σ γ 2 = V a r ( γ k ) , and σ M 2 = 0 1 x 2 ν ( d x ) . Next, we make the following assumptions for our theoretical results:
Assumption 1.
The safety loading condition holds, i.e., c μ 1 < .
Assumption 2.
J t is the tempered α-stable subordinator.
Assumption 3.
h n = n κ 1 , m n = O ( n κ 2 ) , α [ 0 , 1 ) , and θ ( 0 , 1 2 ) , where κ 1 , κ 2 > 0 , κ 1 + 3 κ 2 < 1 , 1 κ 1 ( 2 2 θ ) + κ 2 < 0 , and κ 2 2 θ κ 1 ( 1 α ) < 0 .
Assumption 1 guarantees that the ruin is not a certain event. Assumption 2 means that ν has a density of the form ν ( x ) = z e λ x x 1 + α I x > 0 , where z > 0 , λ > 0 , and α [ 0 , 1 ) . Assumption 2 implies that μ 1 < , μ 2 < ,
0 r n x k ν ( d x ) ( r n ) k α , k = 1 , 2 , 3 , 4 ,
and r n 1 x ν ( d x ) ( m + r n 1 α ) ,
where m indicates a generic constant.
To establish the convergence rate of Φ ^ ( u ) , we need to calibrate the estimation errors of ρ ^ . The following Theorem 1 gives the rate of convergence of ρ ^ .
Theorem 1.
Let r n = h n θ with θ ( 0 , 1 / 2 ) . Then, under Assumptions 1 and 2,
ρ ^ ρ = 1 n h n Q h n + O P ( h n θ ( 1 α ) + n h n 2 2 θ ) ,
where Q h n D N 0 , 1 c 2 ( μ L μ γ + μ L σ γ 2 μ γ + σ M 2 ) .
Proof. 
Let Δ k J = ( J t k n J t k 1 n ) , Δ k M = ( M t k n M t k 1 n ) , Δ k L = ( L t k n L t k 1 n ) , Δ k W = σ ( W t k n W t k 1 n ) , and Δ k N = ( N t k n N t k 1 n ) . By (9) and (11), we have:
ρ ^ ρ = k = 1 n ( c h n Z k ) I D k n c n h n 1 c μ 1 = 1 c k = 1 n ( Δ k J Δ k W ) I { Δ k J Δ k W > h n θ } n h n μ 1 = 1 c k = 1 n ( b h n + Δ k L Δ k W + Δ k M ) I { b h n + Δ k L Δ k W + Δ k M > h n θ } n h n μ 1 = j = 1 4 I j ,
where:
I 1 = 1 c k = 1 n ( Δ k L Δ k W ) I { Δ k L Δ k W > h n θ } n h n 1 x ν ( d x ) ,
I 2 = 1 c k = 1 n ( Δ k L Δ k W ) n h n ( I { b h n + Δ k L Δ k W + Δ k M > h n θ } I { Δ k L Δ k W > h n θ } ) ,
I 3 = 1 c k = 1 n ( b h n + Δ k M ) n h n 0 1 x ν ( d x ) ,
I 4 = 1 c k = 1 n ( b h n + Δ k M ) n h n I { b h n + Δ k L Δ k W + Δ k M h n θ } .
In the following, we study each of I 1 to I 4 .
(i) Note that:
I 1 = 1 c k = 1 n ( Δ k L Δ k W ) I { Δ k L Δ k W > h n θ , Δ k N = 1 } n h n μ L + 1 c k = 1 n ( Δ k L Δ k W ) I { Δ k L Δ k W > h n θ , Δ k N = 0 } n h n + 1 c k = 1 n ( Δ k L Δ k W ) I { Δ k L Δ k W > h n θ , Δ k N 2 } n h n I 11 + I 12 + I 13 .
By Mancini [25], { Δ k L Δ k W > h n θ } is equal to { Δ k N = 1 } almost surely for small h n . Thus, for small h n ,
I 11 = 1 c k = 1 n ( Δ k L Δ k W ) I { Δ k L Δ k W > h n θ , Δ k N = 1 } n h n μ L = 1 c k = 1 n ( γ k Δ k W ) I { Δ k N = 1 } n h n μ L .
Rewrite:
I 11 = 1 n h n k = 1 n { ( γ k Δ k W ) I { Δ k N = 1 } h n μ L } c n h n Z h n n h n .
Then, due to the independence of γ k , N t , and W t , by the central limit theorem, we have Z h n D N ( 0 , 1 c 2 ( μ L μ γ + μ L σ γ 2 μ γ ) ) .
By (A.2)–(A.4) in You and Yin [32], we have:
I 12 = 1 c k = 1 n ( Δ k L Δ k W ) I { Δ k L Δ k W > h n θ , Δ k N = 0 } n h n = o P ( 1 n h n )
and:
I 13 = 1 c k = 1 n ( Δ k L Δ k W ) I { Δ k L Δ k W > h n θ , Δ k N 2 } n h n = o P ( 1 n h n ) .
Therefore, I 1 = Z h n n h n + o P ( 1 n h n ) .
(ii) Next, we show that I 2 = o P ( h n 1 α θ log ( 1 h n ) ) . Since b h n + Δ k M 0 , we have { Δ k L Δ k W > h n θ } { b h n + Δ k M + Δ k L Δ k W > h n θ } . Thus,
I 2 = 1 c k = 1 n ( Δ k L Δ k W ) n h n ( I { b h n + Δ k L Δ k W + Δ k M > h n θ } I { Δ k L Δ k W > h n θ } ) = 1 c k = 1 n ( Δ k L Δ k W ) n h n I { b h n + Δ k L Δ k W + Δ k M > h n θ , Δ k L Δ k W h n θ } .
For small h n , { Δ k L Δ k W h n θ } = { Δ k N = 0 } almost surely. As a result, we have:
I 2 = 1 c k = 1 n ( Δ k W ) n h n I { Δ k M > h n θ , Δ k L Δ k W h n θ } ( 1 + o p ( 1 ) ) .
By Mancini [27], we obtain that:
1 c [ k = 1 n Δ k W n h n I { Δ k M > h n θ } ] = o P ( h n 1 α θ log ( 1 h n ) ) .
(iii) Applying the central limit theorem, we obtain:
I 3 = 1 c [ k = 1 n ( b h n + Δ k M ) n h n 0 1 x ν ( d x ) ] T h n n h n ,
where T h n D N ( 0 , σ M 2 c 2 ) .
(iv) Now, let us consider the last term of (17).
I 4 = 1 c k = 1 n ( b h n + Δ k M ) n h n I { b h n + Δ k L Δ k W + Δ k M h n θ } = 1 c k = 1 n ( b h n + Δ k M ) n h n I { b h n + Δ k L Δ k W + Δ k M h n θ , b h n + Δ k M > 2 h n θ } + 1 c k = 1 n ( b h n + Δ k M ) n h n I { b h n + Δ k L Δ k W + Δ k M h n θ , b h n + Δ k M 2 h n θ } B 1 + B 2 .
By (A.29)–(A.32) in You and Yin [32], we have
B 1 = 1 c k = 1 n ( b h n + Δ k M ) n h n I { b h n + Δ k L Δ k W + Δ k M h n θ , b h n + Δ k M > 2 h n θ } O p ( n h n 1 + ( 1 2 θ ) ) .
Next, we show that the second term of (23) is of order O ( r n 1 α ) . In fact,
B 2 = 1 c k = 1 n ( b h n + Δ k M ) n h n I { b h n + Δ k L Δ k W + Δ k M h n θ , b h n + Δ k M 2 h n θ } 1 c k = 1 n ( b h n + Δ k M ) n h n I { b h n + Δ k M 2 h n θ } = 1 c k = 1 n ( t k 1 n t k n 0 2 h n θ x μ ( d x , d t ) ) n h n I { b h n + Δ k M 2 h n θ } 1 c k = 1 n ( t k 1 n t k n 0 2 h n θ x μ ( d x , d t ) ) n h n .
Using Assumption 2 and the law of large numbers, we establish that:
E [ 1 c k = 1 n ( t k 1 n t k n 0 2 h n θ x μ ( d x , d t ) ) n h n ] = 1 c 0 2 h n θ x ν ( d x ) h n θ ( 1 α )
and:
B 2 = O ( h n θ ( 1 α ) ) .
Thus, I 4 = O p ( n h n 1 + ( 1 2 θ ) + h n θ ( 1 α ) ) .
Finally, combining (i)–(iv) leads to ρ ^ ρ = 1 n h n ( Z h n + T h n ) + O p ( n h n 1 + ( 1 2 θ ) + h n θ ( 1 α ) ) . Note that Z h n and T h n are independent. It follows that the result of the theorem holds. □
The convergence rate of Φ ^ ( u ) depends on the choice of h n , m n , and M n . The following theorem establishes the convergence rate of the integrated squared error of Φ ^ ( u ) .
Theorem 2.
Under Assumptions 1–3, as n ,
Φ ^ ( u ) Φ ( u ) 2 = O P max { n κ 1 + 3 κ 2 1 , n κ 2 2 κ 1 θ ( 1 α ) , n 1 κ 1 ( 2 2 θ ) + κ 2 , 1 M n } .
Proof. 
By (13), we have:
Φ ^ ( u ) Φ ( u ) 2 = 0 1 2 π m n M n e i u s F ^ Φ ( s ) d s + 1 2 π M n m n e i u s F ^ Φ ( s ) d s Φ ( u ) 2 d u .
Using Parseval’s identity, we obtain that:
Φ ^ ( u ) Φ ( u ) 2 1 2 π m n < | s | < M n c ( 1 ρ ^ ) ψ ^ Y ( s ) c ( 1 ρ ) ψ Y ( s ) 2 d s + 1 2 π | s | M n F Φ ( s ) 2 d s + 1 2 π | s | m n F Φ ( s ) 2 d s = ⨿ 1 + ⨿ 2 + ⨿ 3 .
Note that:
⨿ 1 = 1 2 π m n < | s | < M n ( c μ 1 ) [ ψ ^ Y ( s ) ψ Y ( s ) ] + c ψ Y ( s ) [ ρ ^ ρ ] ψ ^ Y ( s ) ψ Y ( s ) 2 d s m n < | s | < M n ( c μ 1 ) [ ψ ^ Y ( s ) ψ Y ( s ) ] ψ ^ Y ( s ) ψ Y ( s ) 2 d s + m n < | s | < M n c [ ρ ^ ρ ] ψ ^ Y ( s ) 2 d s .
Using the fact that | ψ Y ( s ) | > s ( c μ 1 ) and Theorem 4.1 in You and Cai [33], we establish that:
m n < | s | < M n ( c μ 1 ) [ ψ ^ Y ( s ) ψ Y ( s ) ] ψ ^ Y ( s ) ψ Y ( s ) 2 d s = O P ( 1 n h n m n 3 )
and:
m n < | s | < M n c [ ρ ^ ρ ] ψ ^ Y ( s ) 2 d s = O P ( 1 n h n 1 m n ) + O P ( ( n h n 2 2 θ + h n 2 θ ( 1 α ) ) 1 m n ) .
By Lemma 1 in Zhang and Yang [11], we have ⨿ 2 = O ( 1 M n ) and ⨿ 3 = O ( m n ) . Then, combining (30), (31), and Assumption 3 leads to:
Φ ^ ( u ) Φ ( u ) 2 = O P max { n κ 1 + 3 κ 2 1 , n κ 2 2 κ 1 θ ( 1 α ) , n 1 κ 1 ( 2 2 θ ) + κ 2 , 1 M n } .

5. Conclusions

In this paper, the threshold and Fourier transform (inversion) techniques were employed to construct a new estimator of ruin probability for the spectrally negative Lévy process. The convergence rate of the integrated squared error (ISE) of the estimator was obtained when the jump component was the tempered α -stable subordinator. This shows that the ISE of the estimated ruin probability function is well controlled. Further work includes, but is not limited to deriving the asymptotic distribution of the proposed estimator and making statistical inference for ruin probability, under the framework that the risk model is a spectrally negative Lévy process with dividend strategy and investment. Furthermore, statistical inference for the Gerber–Shiu function and the dividend function are worthy of study.

Author Contributions

Formal analysis, H.Y.; Methodology, H.Y.; Writing, Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of Shandong Province (No. ZR2020QA029).

Informed Consent Statement

Not applicable.

Acknowledgments

The authors express many thanks to the anonymous referees for their helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Croux, K.; Veraverbeke, N. Non-parametric estimators for the probability of ruin. Insur. Math. Econ. 1990, 9, 127–130. [Google Scholar] [CrossRef]
  2. Frees, E.W. Nonparametric estimation of the probability of ruin. Astin Bull. 1986, 16, 81–90. [Google Scholar] [CrossRef] [Green Version]
  3. Mnatsakanov, R.; Ruymgaart, L.L.; Ruymgaart, F.H. Nonparametric estimation of ruin probabilities given a random sample of claims. Math. Methods Stat. 2008, 17, 35–43. [Google Scholar] [CrossRef]
  4. Pitts, S.M. Nonparametric estimation of compound distributions with applications in insurance. Ann. Inst. Stat. Math. 1994, 46, 537–555. [Google Scholar]
  5. Politis, K. Semiparametric estimation for non-ruin probabilities. Scand. Actuar. J. 2003, 2003, 75–96. [Google Scholar] [CrossRef]
  6. Veraverbeke, N. Asymptotic estimates for the probability of ruin in a Poisson model with diffusion. Insur. Math. Econ. 1993, 13, 57–62. [Google Scholar] [CrossRef]
  7. Huang, Y.; Li, J.; Liu, H.; Yu, W. Estimating ruin probability in an insurance risk model with stochastic premium income based on the CFS method. Mathematics 2021, 9, 982. [Google Scholar] [CrossRef]
  8. Li, J.; Yu, W.; Liu, C. Nonparametric estimation of ruin probability by complex Fourier series expansion in the compound Poisson mode. Commun. Stat.-Theory Methods 2020. [Google Scholar] [CrossRef]
  9. You, H.; Guo, J.; Jiang, J. Interval estimation of the ruin probability in the classical compound Poisson risk model. Comput. Stat. Data Anal. 2020, 144, 106890. [Google Scholar] [CrossRef]
  10. Zhang, Z.; Yang, H. 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]
  11. Zhang, Z.; Yang, H. 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]
  12. Zhang, Z. Nonparametric estimation of the finite time ruin probability in the classical risk model. Scand. Actuar. J. 2017, 5, 452–469. [Google Scholar] [CrossRef]
  13. Zhang, Z.; Yang, H.; Yang, H. On a nonparametric estimator for ruin probability in the classical risk model. Scand. Actuar. J. 2014, 4, 309–338. [Google Scholar] [CrossRef] [Green Version]
  14. Su, W.; Yu, W. Asymptotically normal estimators of the Gerber–Shiu function in classical insurance risk model. Mathematics 2020, 8, 1638. [Google Scholar] [CrossRef]
  15. Su, W.; Yu, W. Estimating the Gerber–Shiu Function in Lévy Insurance Risk Model by Fourier–Cosine Series Expansion. Mathematics 2021, 9, 1402. [Google Scholar] [CrossRef]
  16. Yang, Y.; Su, W.; Zhang, Z. Estimating the discounted density of the deficit at ruin by Fourier cosine series expansion. Stat. Probab. Lett. 2019, 146, 147–155. [Google Scholar] [CrossRef]
  17. Zhang, Z.; 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]
  18. Su, W.; Yong, Y.; Zhang, Z. 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]
  19. Zhang, Z.; 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]
  20. Zhang, Z. Estimating the Gerber–Shiu function by Fourier–Sinc series expansion. Scand. Actuar. J. 2017, 10, 898–919. [Google Scholar] [CrossRef]
  21. Shimizu, Y. Non-parametric estimation of the Gerber–Shiu function for the Wiener–Poisson risk model. Scand. Actuar. J. 2012, 2012, 56–69. [Google Scholar] [CrossRef]
  22. Shimizu, Y.; Zhang, Z. Estimating Gerber–Shiu functions from discretely observed Lévy driven surplus. Insur. Math. Econ. 2017, 74, 84–98. [Google Scholar] [CrossRef]
  23. Asmussen, S.; Albrecher, H. Ruin Probabilities, 2nd ed.; World Scientific: Singapore, 2010. [Google Scholar]
  24. Comte, F.; Genon-Catalot, V. Nonparamtetric estimation for pure jump Lévy processes based on high frequency data. Stoch. Process. Their Appl. 2009, 119, 4088–4123. [Google Scholar] [CrossRef] [Green Version]
  25. Mancini, C. Estimation of the characteristics of the jump of a general Poisson-diffusion model. Scand. Actuar. J. 2004, 1, 42–52. [Google Scholar] [CrossRef]
  26. Mancini, C. Non-parametric threshold estimation for models with stochastic diffusion coefficient and jumps. Scand. J. Stat. 2009, 36, 270–296. [Google Scholar] [CrossRef]
  27. Mancini, C. Non-parametric tests for pathwise properties of semimartingales. Bernoulli 2011, 17, 781–813. [Google Scholar]
  28. Mancini, C. The speed of convergence of the threshold estimator of integrated variance. Stoch. Process. Their Appl. 2011, 121, 845–855. [Google Scholar] [CrossRef] [Green Version]
  29. Shimizu, Y. A new aspect of a risk process and its statistical inference. Insur. Math. Econ. 2009, 44, 70–77. [Google Scholar] [CrossRef]
  30. Shimizu, Y. Functional estimation for Lévy measures of semimartingales with Poissonian jumps. J. Multivar. Anal. 2009, 100, 1073–1092. [Google Scholar] [CrossRef] [Green Version]
  31. Sato, K. Lévy Processes and Infinitely Divisible Distributions; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  32. You, H.; Yin, C. Threshold estimation for a spectrally negative Lévy process. Math. Probl. Eng. 2020, 1, 1–12. [Google Scholar] [CrossRef]
  33. You, H.; Cai, C. Nonparametric estimation for a spectrally negative Lévy process based on high frequency data. J. Comput. Appl. Math. 2019, 345, 196–205. [Google Scholar] [CrossRef]
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.; You, H. The Speed of Convergence of the Threshold Estimator of Ruin Probability under the Tempered α-Stable Lévy Subordinator. Mathematics 2021, 9, 2654. https://doi.org/10.3390/math9212654

AMA Style

Gao Y, You H. The Speed of Convergence of the Threshold Estimator of Ruin Probability under the Tempered α-Stable Lévy Subordinator. Mathematics. 2021; 9(21):2654. https://doi.org/10.3390/math9212654

Chicago/Turabian Style

Gao, Yuan, and Honglong You. 2021. "The Speed of Convergence of the Threshold Estimator of Ruin Probability under the Tempered α-Stable Lévy Subordinator" Mathematics 9, no. 21: 2654. https://doi.org/10.3390/math9212654

APA Style

Gao, Y., & You, H. (2021). The Speed of Convergence of the Threshold Estimator of Ruin Probability under the Tempered α-Stable Lévy Subordinator. Mathematics, 9(21), 2654. https://doi.org/10.3390/math9212654

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