Estimating the Gerber-Shiu Function by Fourier Cosine Series Expansion in the Wiener-Poisson Risk Model ()
1. Introduction
In this paper, we consider that the financial surplus of an insurance company evolves as a compound Poisson risk model perturbed by a Wiener process
(1)
where
is the initial capital and
is the constant premium per time. The aggregate claims
follows a compound Poisson process, where the number of claims
is an homogeneous Poisson process with intensity
, and the individual claim sizes
is a sequence of positive i.i.d. random variables generated by a generic variable X with density
and mean
. Finally,
is a standard Brownian motion with
, and
is the diffusion volatility parameter. We suppose that
,
and
are mutually independent. In this paper, we shall assume throughout the safety loading condition
, so that ruin is an uncertain event.
The ruin time is defined by
with the convention
if
for all
. Problems related to ruin are a hot topic in risk theory literature. Starting with Cramér and Lundberg’s fundamental collective insurance risk model which was established at the turn of the century [1] [2] , actuarial scholars are still investigating into specific examples of it or tweaking some of its aspects to make it more useful. [3] proposed three equivalent equations for ruin probability in a Cramér Lundberg model, and the solutions were means of inverse Laplace transform. To extend the ruin problem, [4] was the first to propose the joint distribution of the time to ruin, the surplus before ruin, and the deficit at ruin, named the Gerber-Shiu function.
(2)
where
is a measurable penalty function of the surplus prior to ruin and the deficit at ruin,
represents the force of interest and
is the indicator function of the event A. Set
. The Gerber-Shiu function is a powerful tool for studying ruin related problems. To make it simple, when
,
becomes the ruin probability when
or the Laplace transform of the time to ruin when
. For
,
becomes the expected discounted deficit at ruin when ruin is caused by a claim. Finally, for
,
becomes the expected discounted claim size causing ruin. Over the last decades, Gerber-Shiu has drawn the attention of numerous actuarial scholars in various risk models. [5] [6] [7] [8] [9]
The Wiener-Poisson risk model was initially proposed in the literature of actuarial science by [10] to extend the classical compound Poisson risk model, where the diffusion perturbation is used to describe the uncertainty of the premium income and aggregate claims. Since then, many scholars have made substantial contributions to this model. For instance, [11] showed that the probability of ruin satisfied certain renewal equations; [12] studied the ruin problems; [13] derived and solved boundary conditions for the n-th moment of the discounted dividend payments; [14] investigated some explicit solutions of the expected discounted penalty function; [15] considered a model with two-sided jumps. For further investigation, we refer the interested readers to [16] [17] and [18] . The explicit formula of the Gerber-Shiu function can be determined when it is assumed that the Poisson intensity, claim size density and diffusion volatility parameter are known. However, insurance companies do not have these probability distributions, instead they have the data information on the surplus flow levels, individual claim sizes and claim numbers. So the big challenge is to figure out how to estimate the Gerber-Shiu function based on the data information.
Recently, non-parametric estimation of risk measures has emerged as a hot topic in risk theory. For the classical model, [19] proposed a semi parametric estimator of the ruin probability, [20] constructed a non parametric estimator for ruin probability by Fourier inversion and kernel density estimation method, [21] estimated the finite time ruin by double Fourier transform. For the Wiener-Poisson risk model, [22] presented the estimator of the Gerber-Shiu function based on a regularized inversion of the Laplace transform, [23] presented an alternative method, the Fourier sinc series expansion to ameliorate the regularized inversion of Laplace transform, [24] proposed a new method for the estimation of the Gerber-Shiu function, Laguerre series expansion method which is not based on the Fourier transform and the Laplace transform.
In this paper, we implement the Fourier cosine series expansion to estimate the Gerber-Shiu function in the risk model (1). The Fourier cosine series expansion has some applications in pricing financial products and it is also named the COS method in the literature. In the field of risk theory, [25] approximated the Gerber-Shiu function under the Levy subordinator model, [26] studied the density of the time to ruin in the classical compound Poisson risk model and [27] estimated the expected discounted penalty function in the Levy risk model based on the Fourier cosine series expansion.
The remainder of this paper is organized as follows. In Section 2, we introduce the Fourier cosine series expansion. In Section 3, we show how to construct our estimator based on the COS method. The convergence rate of the estimator is derived in Section 4. Finally, some numerical simulations are presented in Section 5 to show that the estimator performs well when the sample size is finite and Section 6 is the conclusion.
2. Fourier Cosine Series Expansion
This section presents some known results of the COS method. For an integrable function g with a finite support
, we have the following Fourier Cosine series expansion,
(3)
where
indicates that the first term in the summation is weighted by one-half, and the cosine coefficients are given by
(4)
Let define
as the Fourier transform of g. Assume that
(5)
Then
(6)
Finally, the approximation of our function is given by
(7)
The approximation error for the COS method is provided in [28] by the following Lemma.
Lemma 1. For the real-value integrable function g supported on
, suppose that
,
and
. Then, for some positive constants
and
, we have
(8)
Remark. Lemma 8 shows that COS method has two types of approximation errors depending on the parameter a which is in fact an integration domain truncation parameter in the COS method and the parameter K which is used in the approximation of the of the COS coefficients by the Fourier transform
. The first error term
means that larger K yields better approximation, but large a may slowdown the convergence rate. The second error term
means that larger a can result in more accurate approximation of the COS coefficient. However, this error term is increasing with respect to K since larger K means more COS coefficients have to be approximated.
3. Estimation Procedure
In this section, we shall first approximate the expected discounted penalty function by the COS method. Afterward, we shall replace the estimate values in the approximate function to have our estimator.
3.1. The COS Approximation Method
We consider how to use the COS method to approximate the expected discounted penalty function
. The Laplace transform of the Gerber-Shiu function is given as follow, see [29]
(9)
where
and
.
From the relation between Laplace transform and Fourier transform, it follows that
(10)
By Formula 7, we use the COS method approximation as follow
(11)
where the COS coefficients are given by:
Remark. The approximation error can also be obtained Lemma 8. Furthermore, suppose that
,
,
,
and
,
and let C be a positive constant, then from Lemma 8 and Remark 2, we can obtain
(12)
where
(13)
3.2. The COS Estimation Method
Since we have approximated the expected discounted penalty function, it can now be estimated based of the dataset of the surplus flow level, a random sample on individual claim sizes and claim numbers.
Let us assume that the premium rate c is constant, but the Poisson parameter
and the claim size density function
are both unknown, but as in [23] , we assume that the surplus process can be observed over a long time interval
. Let
be a sampling interval. Without loss of generality, we assume that
is an integer and let
.
Suppose that the insurer can get the following data-set.
• Data-set of surplus level:
Where
is the observed surplus level at time
• Data-set of claim numbers and claim sizes:
where
is the total claim number up to time
.
We shall propose an estimator for the expected discounted penalty function
. Obviously, we need to estimate the following quantities
and
. As in [23] , we can estimate
and
by
It is known that
and
.
Finally, we propose an estimator for the expected discounted penalty function
. It is easily seen that
see [23] . From which we obtain the estimators for
and
Now the Fourier transforms
can be estimated as follows,
(14)
Then
is estimated as follows,
(15)
where the COS coefficients are given by:
4. Consistency Properties
In this section, we derive the consistency properties of our estimators when the observation interval
is very large. First, we know that
,
Lemma 2. Suppose that
,
and
and
,
, then we have
and
Proof. The convergence of
is well known, see [30] .
For the Gerber-Shiu function, we have that
(16)
where
It is shown in [31] that
Hence,
(17)
Then by formula 16 and 17, we can conclude that
5. Numerical Simulation
In this section, we present some numerical results to show the effectiveness of our method. All computations are done in MATLAB on a EliteBook, with Intel(R) Core(TM) i5-6300U CPU@2.40GHz 2.50GHz and a RAM of 8GB. Throughout this section, we set
,
,
and consider claim size density functions, the exponential EXP(1):
and the Erlang (2, 2):
. Then the true value of the Gerber-Shiu function can be found in [23] . Now for these two claim size density functions, closed form of Fourier transforms exist so that the COS method approximation of
can be computed. When using Fourier Cosine Series expansion method to approximate
, we apply the cumulant method given in [32] to determine the parameter a. We take
,
to provide the benchmark.
First of all, we test the effectiveness of the Fourier Cosine method for approximating
. In Table 1 and Table 2, we present some average relative errors and average absolute errors for
respectively for exponential and Erlang distributions which are calculated by the following
Here, we take
, since when
,
is very small. We consider the truncation parameter
,
. It can be observed that in each column of Table 1 and Table 2, both average relative errors and the average absolute errors are decreasing w.r.t. q, which implies that large truncation parameter can reduce the bias under the model setting. In Table 1 and Table 2, we find that the average relative errors are larger than the average absolute errors, which is due to that
is always smaller than 1, and in particular, it is close to zero for large initial surplus.
Next we test the performance of our estimator
. We fixe
. As for the observation interval
, We shall take
for
. For the COS parameter a, we also apply the cumulant method given in [32] . We take
where
. We repeat 300 simulations, and compute empirical average relative errors and average absolute errors for
. Which are defined by
where
denote the jth simulation values of
. In Table 3 and Table 4, we present the empirical estimation errors for
respectively for exponential and Erlang distributions. As expected, both the empirical average errors and empirical absolute errors are decreasing w.r.t. p which is due to that, as p increases (or equivalently T increases), more sample can be used to estimate
. Again, we observe that the average relative errors are larger than the average absolute errors.
Finally, we plot 300 consecutive estimators (green curves) on the same picture together with the true curve (red curve) to illustrate variability bands and show the stability of the procedures. In Figure 1, we observe that the beams of estimators are much closed to the true curves. In particular, it follows that the variances for estimating
are very small for a large observation interval. We find that
is a decreasing function of the initial surplus u, which means that the ruin is more likely to happen when u is small. At the same time, we can also observe that as T increases, the estimator tends to be stable and converges to
. We plot the corresponding curves. Again, we can observe that the estimator becomes better as T becomes larger.
Table 1. Approximation Errors for Gerber-Shiu function
based on EXP (1).
Table 2. Approximation Errors for the Gerber-Shiu function
based on Erlang (2, 2).
Table 3. Estimation Errors Errors for the Gerber-Shiu function
based on EXP (1).
Table 4. Estimation Errors for the Gerber-Shiu function
based on Erlang (2, 2).
6. Conclusion
In this paper, we have proposed a new estimator of the expected discounted penalty function in the perturbed compound Poisson risk model. Our estimator is based on the COS method. The COS coefficient is easily derived based on the Fourier transform of the Gerber-Shiu function which is useful when approximating the function by the COS method. Suppose we have the surplus flow level, the claim number, and the claim sizes over a long-term interval. Then, we construct our estimator based on the COS approximation we have derived, by replacing the different functions with their estimates. We have derived theoretical errors and presented some simulation results to show the effectiveness of our estimator. We have shown that our estimator has an accurate convergence rate. Since the COS method does not utilize the Fast Fourier transform algorithm, our estimator is easier to compute compared to the Fourier sinc method. Since this method and the Fourier sinc method use the Fourier transform of the Gerber-Shiu, the comparison of the two methods could be future research.
Acknowledgements
The authors would like to thank the Pan African University Institute of Basic Sciences Technology and Innovation (PAUSTI) for their financial assistance.