Journal of Applied Mathematics
Volume 2019 (2019), Article ID 2191509, 9 pages
https://doi.org/10.1155/2019/2191509

Research Article

An Approximation of Minimum Initial Capital of Investment Discrete Time Surplus Process with Weibull Distribution in a Reinsurance Company

Soontorn Boonta¹ and Somchit Boonthiem²

¹Department of General Science, Faculty of Science and Engineering, Kasetsart University, Chalermphrakiat Sakon Nakhon Province Campus, Sakon Nakhon 47000, Thailand
²Mathematics and Statistics Program, Faculty of Science, Sakon Nakhon Rajabhat University, Sakon Nakhon 47000, Thailand

Correspondence should be addressed to Soontorn Boonta; soontorn.bo@ku.th

Received 31 January 2019; Accepted 23 May 2019; Published 11 June 2019

Academic Editor: Ying Hu

Copyright © 2019 Soontorn Boonta and Somchit Boonthiem. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Catastrophe is a loss that has a low probability of occurring but can lead to high-cost claims. This paper uses the data of fire accidents from a reinsurance company in Thailand for an experiment. Our study is in two parts. First, we approximate the parameters of a Weibull distribution. We compare the parameter estimation using a direct search method with other frequently used methods, such as the least squares method, the maximum likelihood estimation, and the method of moments. The results show that the direct search method approximates the parameters more precisely than other frequently used methods (to four-digit accuracy). Second, we approximate the minimum initial capital (MIC) a reinsurance company has to hold under a given ruin probability (insolvency probability) by using parameters from the first part. Finally, we show MIC with varying the premium rate.

1. Introduction

The risks of an insurance company can be assessed based on disasters of varying severity. The insurance company evaluates its risks in order to maintain consistency. Insolvency cannot occur if the company knows how to manage the risk process. For instance, if an insurance company does not have sufficient initial capital to pay some claims, then the company can share some of the risks by transferring them to reinsurance. Parameter estimation is an important method to construct a risk model in an insurance business. The well-known methods are the least squares method, the maximum likelihood estimation, and the method of moments. This research is interested in the estimation of the parameters of a Weibull distribution which are represented by fire accident data. Many authors have studied the different aspects of Weibull parameters. Bergman [1] and Sullivan and Lauzon [2] proposed four probability estimators which were frequently applied in the least squares method. Based on a Monte Carlo simulation, Khalili and Kromp [3] and Trustrum and Jayatilaka [4] compared an estimation of Weibull parameters by using the least squares method, the maximum likelihood estimation, and the method of moments. Boonta et al. [5] proposed a direct search technique to estimate the parameters of Weibull distribution. They compared the Chi-squared value of the direct search technique to the least squares method, the maximum likelihood estimation, and the method of moments. The results showed that the direct search technique gave a more precise estimation than the least squares method, the maximum likelihood estimation, and the method of moments.

In this paper, we start by introducing the surplus of nonlife insurance. The surplus can be described as

Lundberg [6] was the first actuary who considered the surplus process of nonlife insurance under three assumptions in his model.

(1) Claims happening at times satisfying are called claim arrivals or claim times.

(2) The -th claim arriving at time causes the claim size or claim severity . The sequence of claim sizes constitutes an independent and identically distributed (i.i.d.) sequence of nonnegative random variables.

(3) The claim severity process and the claim arrival process are mutually independent.

Next, we define the claim number process

Thus, is the number of claims in .

Next, we denote as the i.i.d. process of the claims, as the initial capital, and as the premium rate for one unit of time. There are many research studies of ruin probability in terms of initial capital (see in Pavlova and Willmot [7], Dickson [8] Li [9, 10], and Rongming and Haifeng [11]). Chan and Zhang [12] have studied a discrete time surplus process such as claim time (claims which occur every day). They proposed recursive and explicit formulae of the ruin probabilities which are in the form

Sattayatham et al. [13] generalized the results of Chan and Zhang for claim times which do not occur every day. Their model is of the form

This can be rewritten as

where is the inter arrival time process, assuming i.i.d. such that Since the formula of ruin probability is difficult to find explicitly, they proposed the ruin probability in the recursive form:

where . Chengguo Weng et al. [14] studied a model of ruin probability by adding some investments as

where is the constant interest rate for a period of time . They considered as a sequence of dependent individuals that have a regular variation distribution and zero index of upper tail dependence. They also established some asymptotic results for both finite ruin probability and ultimate ruin probability.

Reinsurance and investment are a normal activity of insurance companies because reinsurance can reduce the risk (ruin probability) arising from claims, and the investment can make more profit for the company. The process can be controlled by reinsurance, i.e., by choosing the retention level (or risk exposure) of a reinsurance for one period. The (measurable) function specifies the part of the claim paid by the insurer. Then depends on the retention level (fixed in the risk sharing contract) at the beginning of the respective period where . Hence is the part paid by the reinsurer. It is natural to assume that is increasing in .

In the case of an excess of loss reinsurance we have

with retention level .

In case of a proportional reinsurance we have

with retention level 1.

Therefore, the retention level stands for the control action “no reinsurance” which explains the property “”. The smallest retention level may be chosen in such a way that the condition is satisfied. Then may be calculated according to the expected value principle with safety loading of the reinsurer:

Recently, Luesamai and Chongcharoen [15] expanded the risk model by adding proportional reinsurance and investment. Insurers can invest in the bond and stock markets, and they assume that the interest rates of the bonds have a finite number of possible values and follow a time-homogenous Markov chain. Moreover, they assume that the controlling reinsurance and stock investment values in each time period are constant values. For every time period unit , the risk model is formulated as

where is computed by a proportional reinsurance function, and the sequences , and are mutually independent. The interest rate is assumed to follow a time-homogeneous Markov chain . is the gross return and is assumed to be a sequence of i.i.d. nonnegative random variables. is the amount of stock investment. The results of the study led to the proposal of two upper bounds of ruin probability under a discrete time risk model for reinsurance by generalizing the classical model in terms of two controlling factors: proportional reinsurance and investment.

In this research, we study a risk model of reinsurance by adding investment (buying bonds or fixed accounts). We present two parts consisting of an approximation of the parameters of Weibull distribution and the calculation of the minimum initial capital of investment discrete time surplus process with Weibull distribution. The first part is the estimation of the Weibull parameters using the direct search technique, the least squares method, the maximum likelihood estimation, and the method of moments. The selected method gives the minimum KS statistic value (to four-digit accuracy). The second part is a simulation to calculate the ruin probability of the surplus process under the condition that a reinsurance company can invest in risk-free assets (bonds or fixed accounts). The surplus process is of the form

where is the i.i.d. process of the claims, is the initial capital, and is the daily interest rate which can calculated by per annum. is the premium rate for one unit of time which can be computed by where is a safety loading and is the inter arrival time process, and assuming i.i.d. such that . We approximate the minimum initial capital (MIC) an insurance company has to hold under a given ruin probability (insolvency probability) by using parameters from the first part.

2. Materials and Methods

2.1. Weibull Distribution

The classical Weibull distribution is useful for reliability engineering. Moreover, it can be extended to the various families of probability distributions which deal with the estimation of model parameters by maximum likelihood and it can also be used to illustrate the potentiality of the extended family with two applications to real data [16]. Furthermore, a nonclassical Weibull distribution can be used to estimate the statistical characteristics in a cellular automaton such that a cell’s yield stress is assumed to be a Weibull distribution [17].

Normally, claims that occur infrequently but have high costs will be called catastrophe losses. For example, a fire accident is a type catastrophe loss. Furthermore, a Weibull distribution that shape parameter being less than one and scale parameter being greater than zero is also an example of catastrophe loss. The probability density function of three parameters of a Weibull distribution is of the form

and the cumulative distribution function is of the form

for all , where is a positive shape parameter, is a positive scale parameter, and is a positive location parameter, respectively.

In our work, the costs of all claims are greater than twenty million Baht. We establish and . Thus we revise (13) and (14) as (15) and (16)

and

2.2. Weibull Parameter Estimation

2.2.1. Least Squares Method

Bergman, Sullivan, and Lauzon proposed the probability estimator for the th ranked and is the sample size as shown in the following equations:

We take the natural logarithm to (16),

Then we set and . Therefore, (18) is the linear equation which is of the form .

2.2.2. Maximum Likelihood Estimation

Let be the samples from a Weibull distribution. A log-likelihood function is defined by

By partial derivative the log-likelihood function , , and setting to zero, we have

Thus

2.2.3. Method of Moments

The k th moment of the Weibull distribution, , is defined by

where is a gamma function which is given by

Setting , we obtain

Since ,

We calculate a coefficient of variation (CV) of the Weibull distribution from the formula

If we apply the bisection method to (27), we get the parameter . Thereafter, we substitute in (25) which obtains the parameter .

2.2.4. Kolmogorov-Simirnov Test

The Kolmogorov-Simirnov (KS) test is a distance test. The KS test works well with small samples. Let be the order of statistics of a random sample . The empirical distribution function is of the form

Let be a hypothesized distribution function, so the KS statistic is defined by .

Throughout this research, we use the data of the cost of all claim sizes that are greater than twenty million Baht from the fire insurance of Thai Reinsurance Public. The amount is shown in Table 1 which is greater than twenty million Baht, i.e., . Boonta et al. [5] estimated the shape parameters and scale parameters from a variety of estimation methods by means of a minimum Chi-squared test as shown in Table 2 and using the data according to Table 1.

Table 1: Claims size (millions Baht).

15.50
6.70
49.90
6.40
12.40
102.70
44.90

56.50
138.90
107.30
13.20
13.10
37.70
5.70

40.00
1.80
40.20
84.30
47.30
112.20
9.20

28.50
0.90
43.10
3.60
45.80
70.00
2.30

35.30
7.20
64.60
13.00
2.40
1.40
2.10

7.50
31.50
4.20
37.20
0.70
24.40
14.20

20.10
0.40
33.20
9.30
10.80

Table 2: Shape and scale parameters for a variety of estimation methods.

Method
Type

1
Least squares method 1
5.1
0.8405
28.7721

2
Least squares method 2
5.2
0.8310
28.8602

3
Least squares method 3
5.3
0.7984
29.1888

4
Least squares method 4
5.4
0.8580
28.6168

5
Maximum likelihood estimation
-
0.8633
28.8668

6
Method of moments
-
0.9286
30.0055

We perform the KS test according to the six methods as shown in Table 2. The hypothesis : the set data in Table 1 assumed the Weibull distribution with parameters and . The result shows that the set of data cannot reject the Weibull distribution at 95% confidence. For example, the hypothesis : the set of data assumed the Weibull distribution with and . The KS critical value at 95% confidence is whereas the KS statistic value is 0.0709. Thus the set of data accepts the Weibull distribution at 95% confidence. Figure 1 shows a plot of the cumulative empirical distribution function with the cumulative Weibull distribution function .

Figure 1: The empirical CDF (blue color) and the Weilbull CDF (red color).

3. Results

3.1. Part I: The Direct Search Technique

The direct search technique is a numerical optimization method. The principle of the direct search technique is iterative and random shifting from the beginning solution to a better solution. The direct search technique process can be described as follows.

Step 1. Begin from the simple initial parameters and , i.e., , ( is mean of 47 observe data) and compute the KS statistic value for and

Step 2. Calculate value to and to by using the formulas:

where is randomly chosen from .

Step 3. Calculate the KS statistic value by using and .

Step 4. Compare the KS value of Step 1 and Step 3.

(4.1)If the KS statistic value of Step 1 is less than Step 3, then we do Step 2 again.(4.2)If (4.1) is false we set , and then we perform Step 2 again.

Step 5. Iterate until the process is complete (four-decimal place accuracy).

If value for the KS test is approximated to four-decimal place accuracy, then the direct search technique is completed. We perform the simulation 1000 times. The results are shown in Table 3.

Table 3: Parameters , and the KS statistic value by the direct search technique (Initial parameters ).

Times

Value for the KS test

1
0.5049
31.0142
0.1262

10
0.5489
30.6153
0.1088

100
0.7652
29.5447
0.0544

1000
0.7652
29.5451
0.0544

Table 4 shows that the value for the KS test of the direct search technique equals 0.0544 which is less than that of the other methods. Therefore, the direct search technique with given simple initial parameters and affect the preciseness of the parameter estimation (four-digit accuracy) more than the least squares method, the maximum likelihood estimation, and the method of moments. Consequently, the set data in Table 1 accepts the Weibull distribution at 95% confidence with parameters and . These parameters will be used to approximate the minimum initial capital of the surplus process (12) which will be explained in Part II.

Table 4: The value for the KS test from a variety of estimation methods.

Method
Type

Value for the KS test

1
Least squares method 1
0.8405
28.7721
0.0649

2
Least squares method 2
0.8310
28.8602
0.0636

3
Least squares method 3
0.7984
29.1888
0.0591

4
Least squares method 4
0.8580
28.6168
0.0674

5
Maximum likelihood estimation
0.8633
28.8668
0.0709

6
Method of moments
0.9286
30.0055
0.0964

7
Direct search technique
0.7652
29.5450
0.0544

3.2. Part II: The Simulation of the Ruin Probability and the Minimum Initial Capital

3.2.1. The Simulation of the Ruin Probability

We consider the surplus process (12) under the condition of the claim time for with when days. A simulation method is used to compute the ruin probability , where . The flowchart is shown in Figure 2.

Figure 2: The flow chart of a simulation of the ruin probability.

In Figure 2, the premium rate for one unit of time, , can be calculated by the expected value principle, i.e.,

where is the sequence of claims which are assumed to be i.i.d. and is the sequence of the Poisson random variable which is also assumed to be i.i.d. We estimate that the parameter is equal to 7.766. Then we get the premium rate as the following equation:

Thereafter we perform the simulation of the ruin probability by setting the initial capital million Baht. The simulations are computed 10,000 iterations (see Figure 2). The results are shown in Figure 3.

Figure 3: The relation between the initial capital and the ruin probability in the case of .

3.2.2. Minimum Initial Capital

We consider the relationship between the ruin probability and the initial capital in Figure 3 which is the quadratic function,

Next, we use the quadratic regression method in (32); we obtain the parameters , and as shown in Table 5.

Table 5: The parameters , , and of the quadratic function by varying safety loading .

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

Finally, we consider that the ruin probability is not greater than , i.e.,

For example, we can apply (33) to calculate the MIC in the case of , and . Consequently, we obtain

and we get . Therefore, the MIC of a reinsurance company has to hold for ensuring that the ruin probability is not greater than , i.e., the MIC is equal to 1642.3413 million Baht. In the other cases of and , we use a similar method. Then, we obtain the MIC shown in the Table 6.

Table 6: The MIC of a reinsurance company has to hold for ensuring that the ruin probability is not greater than .

Premium rate
MIC million Baht
MIC million Baht

0.1
4.9027
1642.3413
1622.0112

0.2
5.3484
1627.1924
1606.6316

0.3
5.7941
1615.2543
1594.3263

0.4
6.2398
1602.4812
1581.1251

0.5
6.6855
1589.8425
1567.9948

0.6
7.1312
1582.7270
1560.1270

0.7
7.5769
1574.8999
1551.4706

0.8
8.0226
1563.0873
1538.9719

0.9
8.4683
1555.8489
1530.7069

1.0
8.9140
1546.0777
1520.0044

Table 6 shows that if increases, it affects the premiums which also increase and then the MIC decreases under given regulation.

4. Conclusion and Suggestion

This work presents the use of the direct search technique to estimate the parameters of a Weibull distribution. The direct search technique is compared with the least squares method, the maximum likelihood estimation, and the method of moments. The results show that the direct search technique is more accurate than the other methods (accuracy of four decimal points). Fire accident data is used for special cases exceeding 20 million Baht. The results show that if the safety loading is increased, the MIC decreases under the given regulation; i.e., reinsurance will enable the financial management plan to prevent insolvency.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This paper is supported by the Faculty of Science and Engineering, Kasetsart University, Chalermphrakiat Sakon Nakhon Province Campus, and Mathematics and Statistics Program, Faculty of Science and Technology, Sakon Nakhon Rajabhat University, Sakon Nakhon, 47000, Thailand.

References

B. Bergman, “On the estimation of the Weibull modulus,” Journal of Materials Science Letters, vol. 3, no. 8, pp. 689–692, 1984.

J. D. Sullivan and P. H. Lauzon, “Experimental probability estimators for Weibull plots,” Journal of Materials Science Letters, vol. 5, no. 12, pp. 1245–1247, 1986.

A. Khalili and K. Kromp, “Statistical properties of Weibull estimators,” Journal of Materials Science, vol. 26, no. 24, pp. 6741–6752, 1991.

K. Trustrum and A. D. S. Jayatilaka, “On estimating the Weibull modulus for a brittle material,” Journal of Materials Science, vol. 14, no. 5, pp. 1080–1084, 1979.

S. Boonta, A. Sattayatham, and P. Sattayatham, “Estimation of Weibull parameters using a randomized neighborhood search for the severity of fire accidents,” Journal of Mathematics and Statistics, vol. 9, no. 1, pp. 12–17, 2013.

F. Lundberg, Approximerad Framställning av Sannolikhetsfunktionen, II. Återförsäkering av Kollektivrisker , Almqvist & Wiksell, Uppsala, 1903 (Chinese).

K. P. Pavlova and G. E. Willmot, “The discrete stationary renewal risk model and the gerber – shui discounted penalty function,” Insurance: Mathematics & Economics, vol. 35, no. 2, pp. 267–277, 2004.

D. C. M. Dickson, Insurance Risk and Ruin, Cambridge University Press, New York, NY, USA, 2005.

S. Li, “Distributions of the surplus before ruin, the deficit at ruin and the claim causing ruin in a class of direct time renewal risk models,” Scandinavian Actuarial Journal, vol. 4, pp. 271–284, 2005.

S. Li, “On a class of direct time renewal risk models,” Scandinavian Actuarial Journal, vol. 4, pp. 241–260, 2005.

R. Wang and H. Liu, “On the ruin probability under a class of risk process,” The Journal of the International Actuarial Association, vol. 32, no. 1, pp. 81–90, 2002.

W.-S. Chan and L. Zhang, “Direct derivation of finite-time ruin probabilities in the discrete risk model with exponential or geometric claims,” North American Actuarial Journal, vol. 10, no. 4, pp. 269–279, 2006.

P. Sattayatham, K. Sangaroon, and W. Klongdee, “Ruin probability-based initial capital of the discrete-time surplus process,” Variance: Advancing the Science of Risk, vol. 7, pp. 74–81, 2013.

C. Weng, Y. Zhang, and K. S. Tan, “Ruin probabilities in a discrete time risk model with dependent risks of heavy tail,” Scandinavian Actuarial Journal, no. 3, pp. 205–218, 2009.

A. Luesamai and S. Chongcharoen, “Upper bound of ruin probability for an insurance discrete-time risk model with proportional reinsurance and investment,” Pakistan Journal of Statistics and Operation Research, vol. 14, no. 3, pp. 595–614, 2018.

M. Bourguignon, B. S. Rodrigo, and M. C. Gauss, “The Weibull – G family of probability distributions,” Journal of Data Science, vol. 12, pp. 53–68, 2014.

K. E. Aifantis, A. Konstantinidis, and S. Forest, “Modeling strain localization bands in metal foams,” Journal of Computational and Theoretical Nanoscience, vol. 7, no. 2, pp. 1–7, 2010.