Next Article in Journal
A Study on Balanced Scorecard and Its Impact on Sustainable Development of Renewable Energy Organizations; A Mediating Role of Political and Regulatory Institutions
Previous Article in Journal
Merton Investment Problems in Finance and Insurance for the Hawkes-Based Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monte Carlo Simulation of the Moments of a Copula-Dependent Risk Process with Weibull Interwaiting Time

by
Sharifah Farah Syed Yusoff Alhabshi
,
Zamira Hasanah Zamzuri
and
Siti Norafidah Mohd Ramli
*
Department of Mathematical Sciences, Universiti Kebangsaan Malaysia, Bandar Baru Bangi 43600, Selangor, Malaysia
*
Author to whom correspondence should be addressed.
Risks 2021, 9(6), 109; https://doi.org/10.3390/risks9060109
Submission received: 15 March 2021 / Revised: 30 April 2021 / Accepted: 5 May 2021 / Published: 3 June 2021

Abstract

:
The widely used Poisson count process in insurance claims modeling is no longer valid if the claims occurrences exhibit dispersion. In this paper, we consider the aggregate discounted claims of an insurance risk portfolio under Weibull counting process to allow for dispersed datasets. A copula is used to define the dependence structure between the interwaiting time and its subsequent claims amount. We use a Monte Carlo simulation to compute the higher-order moments of the risk portfolio, the premiums and the value-at-risk based on the New Zealand catastrophe historical data. The simulation outcomes under the negative dependence parameter θ , shows the highest value of moments when claims experience exhibit overdispersion. Conversely, the underdispersed scenario yields the highest value of moments when θ is positive. These results lead to higher premiums being charged and more capital requirements to be set aside to cope with unfavorable events borne by the insurers.

1. Introduction

The primary objective of an insurance company is to ensure that it can pay its promised obligations and remain solvent in the business. This can be achieved by managing surplus processes effectively and charging the appropriate premium amount, which will then guarantee an adequate reserve and capital requirement to cope with any unfavorable events. The ability to estimate accurately the aggregate discounted claims would help insurers manage their future liabilities and estimate appropriate premiums to be charged for an insurance contract under a particular line of business. The aggregate discounted claims, or the present value of the total claims amount paid out by an insurance company, are based on the number of claims occurrences up to time t, claims arrival time and claims amount, as well as the discount factor. Past studies related to the first two moments of the aggregate discounted claims were seen in Léveillé and Garrido (2001a) using renewal theory argument while, per Jang (2004), obtaining the Laplace transformation of the distribution using a different martingale approach. Léveillé and Garrido (2001b) applied the renewal theory arguments and conditioning on the first claim arrival to derive the m-th recursive aggregate moment, assuming independence of claim arrival time from its severity.
The independence assumption between claims size and inter-claim arrival time in a classical risk model, as seen in Jang (2004), Waters (1983) and Yang and Zhang (2001) may no longer be appropriate for insurance risk portfolio modeling. With the increased frequency of catastrophic events (MunichRE 2018), the independence assumption will underestimate the risks faced by insurers, especially the values related to pricing, reserving and ruin measures of a risk portfolio. Past literature has relaxed the independence assumption used in the classical risk theory to include dependency elements between claims amounts and claim occurrences that follows a Poisson process (Barges et al. 2011; Mohd Ramli and Jang 2014; Sun et al. 2020; Woo and Cheung 2013). Barges et al. (2011) and Woo and Cheung (2013) were among the first to derive the m-th moment of the aggregate discounted claims using the Laplace transformation and relax the independence assumptions through a copula. Mohd Ramli and Jang (2014) extended the works by explicitly finding the Neumann series expressions and numerically solving the first two moments, where the dependence structure is captured by Farlie–Gumbel–Morgenstern (FGM), Gumbel, and Gaussian copulae. Mohd Ramli et al. (2018), on the other hand solved the recursive moments which take the form of a Volterra integral equation (VIE) using the Laplace transform, assuming an exponential claims size under the FGM, Frank and the HRT copulae. Recent study by Sun et al. (2020) derived the explicit expressions for the first three moments using a Laplace transformation by adopting the Spearman copula to cover a wide range of positive dependency, and developed the convex approximation of the copula. The studies mentioned above illustrated the use of the Poisson counting process to represent event occurrence.
Owing to the practical limitation in the Poisson process which requires an equidispersed set of data, researchers has also adopted the Poisson-gamma or also known as negative binomial count process that only allows for an overdispersed dataset (Lora and Singer 2011). The Weibull count process, on the other hand, is a good alternative to allow for both underdispersed and overdispersed datasets, as commonly seen in applied data analysis with heterogeneous populations (Kreer et al. 2015; McShane et al. 2008; Winkelmann 1995). McShane et al. (2008) derived a closed form of the renewal process based on Weibull inter-arrival times using a Taylor series expansion. Additionally, as the constant hazard rate assumption is not valid in real practice, the duration dependence or time-varying function under a renewal process replaces the constant hazard function. The non-constant hazard rates that vary according to the duration of the interwaiting time (IWT) would be useful in modeling the unexpected random shocks that lead to the breakdown of many engineering systems, as seen in the work of Liu (2019). Compared with a simple Poisson count model with independent marginals, the Weibull count model with a copula is found to be more useful for predicting the number of goals scored in a football match (Boshnakov et al. 2017).
The Weibull process has also been widely used in insurance and financial studies. Kreer et al. (2015) and Hasumi et al. (2009) mentioned that the probability distribution of the inter-occurrence time of the earthquake events and the automotive claim sizes distribution can be described by the Weibull distribution. It has also been applied in ruin theory see (Albrecher et al. 2011) to establish explicit formulas for ruin probabilities under a Sparre–Andersen risk process. Past researchers worked with explicit solutions of the moments derived from the Laplace transform (Barges et al. 2011; Li and Lu 2018; Mohd Ramli et al. 2018) or Neumann series (Mohd Ramli and Jang 2014), which may require heavy computational capacity resulting from lengthy expressions. The variance expression derived from the Neumann series (Mohd Ramli and Jang 2014) has up to sextuple numerical integration even with a simple Poisson counting process. Thus, Mohd Ramli et al. (2019) work with numerical inversion of the Laplace transform via Gaver–Stehfest algorithm to compute the first two moments of Weibull risk process. However, the slow convergence of the numerical solutions under copulae other than FGM may be improved with better computational techniques. This study seeks to overcome the mathematical complexity under a Weibull risk process by using the Monte Carlo simulation to allow for faster computation and more flexibility to compute the higher-order moments through a wider choice of copula.
We organize the remainder of this article as follows: In Section 2, we introduce the aggregate discounted claims model under Weibull counting process together with copulae to represent the dependency between the IWT and the subsequent claim sizes. In Section 3, we first show the comparison of the simulated moments under Monte Carlo simulation for a special case when λ = 1 , with results obtained in previous studies to ensure consistency. We then fit the distributions and estimate the associated parameters to the New Zealand catastrophe occurrences and losses amount data to illustrate the dispersion in the dataset as allowed by a Weibull process. We also find the best copula to accommodate the dependency in the dataset. We then perform a scenario analysis to examine the characteristics of the resulting risk portfolio through the moments, the VaR and the premium computation. Section 4 concludes the article.

2. Materials and Methods

2.1. Aggregate Risk Model

We define the aggregate discounted claims as in Léveillé and Garrido (2001b), whereby Z ̲ = Z t t 0 with:
Z t = i = 1 N t e δ T i X i if N t > 0 0 if N t = 0 ,
where X i , i = 1 , 2 , is a continuous, non-negative, independent and identically distributed (i.i.d.) random variable that represents claims sizes occurring at random times T i , i = 1 , 2 , , N t . The instantaneous rate of net interest δ is assumed to be deterministic.

2.1.1. Weibull Counting Process

In this study, we let the counting process N ̲ = N t t 0 follow a basic Weibull counting process as in Mohd Ramli et al. (2019), whereby the corresponding continuous random variable of the inter-claim arrival time, W i , is defined as:
W i = T i for i = 1 T i T i 1 for i = 2 , 3 , . . . ,
which follows a Weibull distribution. The time elapsed between two successive claims arrival is called the interwaiting time (henceforth IWT) of the counting process N t . Each pair of the joint variables forms a sequence of i.i.d. random vectors X i , W i i N . The mean and variance of a Weibull count model with shape parameter λ and scale parameter β , derived by McShane et al. (2008) are given by:
E N t = n = 1 j = n n 1 j + n t λ / β j τ j n Γ λ j + 1 , Var N t = E N 2 t E 2 N t
= n = 2 j = n n 2 1 j + n t λ / β j τ j n Γ λ j + 1 n = 1 j = n n 1 j + n t λ / β j τ j n Γ λ j + 1 2 ,
where
τ j n = Γ λ j + 1 Γ j + 1 , for n = 0 , j = 0 , 1 , 2 , . . . m = n 1 j 1 τ j n 1 Γ λ j λ m + 1 Γ j m + 1 , for n = 1 , 2 , 3 , . . . , j = n , n + 1 , n + 2 , . . . ,
in which the parameter n denotes the number of events that occur in the interval [ 0 , t ] .
Although the memoryless property of the exponentially distributed IWT allows for mathematical convenience, the Poisson counting process is only adequate if the data satisfy the restrictive assumption of equidispersion in which the variance of the data is equal to their mean. Hence, the Weibull count model is a better alternative for its ability to accommodate both overdispersed and underdispersed datasets. It also nests commonly used count models, including the Poisson and the negative binomial distributions. Additionally, the Weibull count model also allows for non-constant hazard rates that vary according to the duration of the IWT (McShane et al. 2008). The hazard function characterizes the IWT distribution and relates them to the type of dispersion observed in the corresponding count data. Past studies (Jose and Abraham 2011; McShane et al. 2008; Winkelmann 1995) have verified that the underlying IWT displays negative (positive) duration dependence under a decreasing (increasing) hazard for shape parameter of 0 < λ < 1 ( λ > 1 ) which causes overdispersion (underdispersion). The lack of duration dependence for λ = 1 leads to the Poisson distribution with a constant hazard function, as shown in Figure 1.
We also note that there are two types of dependency applied when modeling the IWT under the Weibull counting process. First dependency lying between the duration of the IWT and the frequency of the claims arrival represented by non-constant hazard rates. As mentioned in McShane et al. (2008) and Winkelmann (1995), the negative (positive) duration dependence indicates a higher (lower) probability of claims occurring immediately after the occurrence of previous claims and decreases (increases) steadily as the IWT increases. This implies that the overdispersed (underdispersed) cases produce higher claims frequencies under short (long) IWT and vice versa.

2.1.2. Copula

The second dependency is between the IWT and its subsequent claims amount, which is captured by a copula. A copula provides a more convenient way to model any joint distribution of two or more random variables. We can model the marginal distributions of each variable by itself, and the copula can link these into a joint distribution. As we allow dependency between the claims size X i , and the inter-arrival time W i , we have a dependent Sparre–Andersen risk process with a dependence structure defined by a copula or other joint probability functions. Please note that when the marginals are negatively (positively) correlated, a large claim amount X i will occur following a short (long) inter-arrival claim time W i , and vice versa. In this study, the dependence structure between the two marginals of X and W are described by a few copulae including the Clayton, Frank and FGM copulae.
Both the Clayton and Frank copulae belong to the same Archimedean family, in which the cumulative distribution function (CDF) of each copula is given by:
C θ C l a y t o n F X x , F W w = m a x F X x θ + F W w θ 1 , 0 1 θ ,
where θ 1 , ,
C θ F r a n k F X x , F W w = 1 θ l n 1 + ( e θ F X x 1 ) ( e θ F W w 1 ) ( e θ 1 ) ,
where θ R . The CDF of an FGM copula, which is from a different class of copula family is given by:
C θ F G M F X x , F W w = F X x F W w + θ 1 F X x 1 F W s ,
where θ 1 , 1 represents the dependence parameter (Klugman and Parsa 1999; Ly et al. 2019). The copulae return to an independence copula when θ = 0 . The FGM copula is a popular choice in extreme value analysis and risk management (Mao and Yang 2015) due to its simplicity and analytical tractability, despite its ability to capture only moderate dependency. The moderate dependence structure for FGM copula can be seen in Figure 2k–o, in which θ 1 , 1 is restricted to Kendall’s tau, τ 2 / 9 , 2 / 9 . In contrast, the Clayton copula (Figure 2a,e) and Frank copula (Figure 2f,j) allow for a wider range of dependency with τ 1 , 1 , as shown by the almost linear scatter plots of strong negative and positive dependency. However, the Frank copula could neither capture lower nor upper tail dependence, whereas the Clayton copula allows for lower tail dependence, or dependency on small values (Ly et al. 2019). The lower tail dependency can be seen clearly in its ’ice cream cone’ shapes as in Figure 2b,d.

2.2. Recursive Moment Expressions

The general form of the m-th moment of aggregate discounted claims Z t , as given here,
μ z m T = E Z m T = 0 T f w s e m δ s E X m | W = s d s + 0 T f w s e m δ s μ z m T s d s + j = 1 m 1 m j 0 T f w s e m δ s E X j | W = s μ z m j T s d s ,
has a recursive nature (Barges et al. 2011; Mohd Ramli and Jang 2014). Working with explicit solutions of Equation (8) require heavy computational capacity due to lengthy expressions as seen in Mohd Ramli and Jang (2014) and Mohd Ramli et al. (2018). We therefore use the Monte Carlo simulation to compute the higher order of moments in a much shorter period with a common computational capacity. Hence, the mathematical properties of the insurance risk portfolio can be examined in a more efficient and viable manner despite the complicated expressions, which would be welcomed by the industrial community (Driels and Shin 2004).

2.3. Monte Carlo Simulation

Monte Carlo simulation is a popular and powerful quantitative tool often used in risk analysis such as in evaluation of insurers’ capital requirements (Casarano et al. 2017), estimation of claim size distribution (Bar-Lev and Ridder 2019), mortality projection (Zamzuri and Hui 2020) and bootstrapping approach in Structural Equation Modeling (Razak et al. 2019). We apply the Monte Carlo simulations to incorporate and understand the impact of risk and uncertainty in our model by computing the probability of different outcomes in a risk process (Zamzuri and Hui 2020).
In this study, we simulate the aggregate discounted claims process using the estimated parameters and distributions from a real insurance dataset. In simulating the copula function, we adopt the same approach as Kelly (2007), whereby the first stage involves modeling the unidimensional marginal distributions while the second stage involves modeling the dependence structure. The steps taken in the simulation process are the following:
  • Generate pairs of dependent random variates X i , W i from multivariate distributions constructed from the chosen copula. The multivariate distributions are based on the best-fit distribution obtained from the insurance dataset.
  • Compute the random time T i , from the accumulated W i as in Equation (2).
  • Compute the aggregate discounted claims Z ( T ) , by assuming deterministic δ .
  • Stop running the iterations once the T i is above a pre-determined term of a policy contract.
  • Repeat the process from step 1 to 4 for n simulations.
  • Determine the moments, premium and VaR from the simulated risk process Z ( T ) .

3. Results and Discussion

In this section, we organize the results and discussion by first showing the consistency of results returned by the Monte Carlo simulation method with results obtained in previous studies (Section 3.1). Next, we fit the aggregate discounted claims model to the New Zealand catastrophe data in Section 3.2. Using the estimated parameters of the claims size distribution, the IWT distribution and the dependency between these marginals, we then examine the risk portfolio of the New Zealand catastrophe dataset. This includes computing the value of moments, premium and VaR in Section 3.3. Finally, we perform scenario analysis in Section 3.4 by varying the parameters of the Weibull process and the copulae to examine the impact on the risk portfolio.

3.1. Results Verification

We first consider the special case of a Weibull IWT with the shape parameter λ = 1 to ensure consistency of the Monte Carlo simulations with the numerical values obtained through Laplace transform as in Barges et al. (2011), Mohd Ramli and Jang (2014) and Mohd Ramli et al. (2018). We simulate pairs of the joint variables of random vectors ( X i , W i ) ( i N ) with an FGM copula defining the dependency between the variables. The values were computed using the same parameters in which X Exp ( α = 0.01 ) , W∼ Exp ( β = 1 ) , δ = 0.04 and T = 5 for different values of θ . The average number of iterations for each simulation is 5, which is consistent with the average of the exponentially distributed IWT. The process is then repeated for n = 100 , 000 simulations.
The mean and variance of the aggregate discounted claims obtained via Monte Carlo simulation are presented in Table 1 and we compare the results with the explicit solution under the Laplace transform as in Mohd Ramli et al. (2018). The simulation results show only slight deviation from the values computed using the Laplace transform, of less than 1% differences for the mean and 2 % for the variance.

3.2. Fitting Distribution and Parameter Estimation of Insurance Datasets

In this section, we fit the aggregate discounted claim models on the New Zealand catastrophe occurrences and losses amount data from 1968 to 2014 which was retrieved from the CASdatasets package in R Statistical Software (Dutang and Charpentier 2020). We estimate the parameters of the fitted distributions of the claims size and the IWT, as well as the dependency between the marginals, using the maximum likelihood estimation method.

3.2.1. The Claim Sizes Distribution

The Cullen and Frey (Cullen et al. 1999) graph in Figure 3 illustrates the skewness-kurtosis plot of the empirical distribution of the New Zealand catastrophe loss claims. The values are computed on 1000 bootstrap samples to allow for the uncertainty of the estimated values of kurtosis and skewness from data. This graph exhibits some common distributions on the plot to help with the choice of distributions to fit to data. The normal, uniform, and logistic distributions are not good fits for the data since the skewness-kurtosis points of the distribution are far from the empirical value. The descriptive statistic in Table 2 shows a positive skewness and a heavy-tailed kurtosis. Thus, the fit of three common right-skewed distributions could be considered, which are the Weibull, Gamma and Lognormal distributions, together with other heavy-tailed distributions such as Pareto and Burr. The results in Table 3 and Figure 4 show that the lognormal distribution, with a meanlog of 1.35969 and an sdlog of 1.52856, provides the best fit for the claims size distribution (values in million USD). The fitted lognormal distribution has the lowest negative log likelihood, Akaike information criterion (AIC) and Bayesian information criterion (BIC) values, with significant p-values under the Kolmogorov–Smirnov (KS) and Anderson-Darling (AD) test.

3.2.2. The Interwaiting Time Distribution

The descriptive statistics of the IWT of the New Zealand catastrophe claims arrival in Table 4 imply that the claims arrive on average every 138.29 days, with a maximum IWT of 7.39 years. The IWT exhibits positive skewness with high kurtosis. Table 5 and Figure 5 show that the Weibull distribution is the best fit for the IWT when compared with exponential distribution. The Weibull shape parameter value of 0.700152 < 1 indicates an overdispersed claims occurrences, implying that the claims arrive at times that vary more than expected.

3.2.3. The Dependency between the Claim Sizes and the IWT

Table 6 exhibits that the Clayton copula is the best fit for the dependence structure between the IWT and its subsequent claim amount embedded in the New Zealand catastrophe historical data. Since Clayton copula allows for lower tail dependence, this indicates that the IWT and the subsequent claims amount show more dependency on small values. The positive dependency between the marginals as shown by parameter θ = 0.32 or the equivalent Kendall’s τ of 0.14, implies that the large (small) claim amount would occur following the long (short) IWT.

3.3. Risk Characteristic of the Risk Process with Overdispersed Claim Arrival

We will examine the effect of overdispersion in the claims arrival of the New Zealand catastrophe historical data on the moments, premium and VaR in comparison to the equidispersed counting process through the Monte Carlo simulation. We are setting δ = 0.025 , X i lognormal (1.3597,1.5286), with W i following the Weibull and exponential distributions from the estimated parameters as in Table 5. The dependency of X i and W i is captured through a Clayton copula with θ = 0.32 .
The simulation results as illustrated in Table 7 show that both risk processes of Z ( t ) are a right-skewed distribution based on the positive skewness. The high value of kurtosis from the results indicates a heavy-tailed distribution. However, the exponential IWT underestimates the value of the mean and variance, which results in a lower value of premium and VaR compared with the overdispersed Weibull IWT as seen in Table 8.

3.4. Scenario Analysis on Risk Characteristic of the Risk Process under Various Dispersion Effect on Claims Arrival

In this section, we perform a scenario analysis on the moments, premium and VaR by varying the shape λ and the scale β parameters under the Weibull IWT distribution. This is done to allow for various dispersion effect of claims occurrences, across different dependence parameter θ . The parameters of the Weibull count process are chosen so that they have an equal mean of claim frequencies but with different variances, obtained through Equation (4). Table 9 illustrates the chosen parameter in which the equidispersed case has an equal mean and variance of 4.243, whereas the variance under the overdispersed (underdispersed) case is higher (lower) than the mean. This is then reflected in the probability distribution of the claims frequencies as seen in Figure 6. There is an increased chance of more claims occurring under the overdispersed case, while the claim frequencies are more concentrated around the mean for the underdispersed case in a fixed time interval.
Please note that the simulation of the aggregate discounted claims under Weibull process adopt two types of dependency. The first type is the duration dependence between the duration of the IWT, W i , and the frequency of claims arrivals, which causes dispersed dataset. The second type is the dependency between the W i and its subsequent claims amount X i , which is captured by a dependence parameter θ of a copula.
Table 10 shows that the highest mean and variance are obtained when the X i and W i exhibit strong negative correlation as seen when θ = 1 . This is intuitively reasonable; short IWT (or frequent claim occurrence) followed by large claims amount X i , will produce higher value of moments. However, as θ changes from negative to positive dependency, the mean and variance are decreasing under all dispersion scenarios.
In comparison to an equidispersed process, the case of overdispersed claims occurrences produces the highest first two moments when θ is negative. Conversely, the underdispersed case shows the highest values of moments when θ is positive. These outcomes could be related to both aforementioned types of dependencies under Weibull process. Higher claim frequencies under short (long) IWT will result in many large claims when θ is negative (positive).
Additionally, the positive skewness and the high value of kurtosis imply that the risk process exhibits a right-skewed with a heavy-tailed risk distribution. However, as θ changes to a strong positive dependency, the tail of the risk process under the overdispersed (underdispersed) case is less (more) skewed with less (more) extreme values on the tail distribution as shown in the value for skewness and kurtosis. This is then reflected in Table 11, whereby the risk process with heavier-tailed distribution leads to a greater difference between the value of 95% and the 99.5% VaR. The 99.5% VaR has more than tripled in value than the 95% VaR under the most heavy-tailed risk process shown by the underdispersed case with a strong positive θ .

Premium Computation and VaR of the Risk Portfolio

The solvency capital requirement (SCR), which is the additional amount of funds an insurer is required to hold, should correspond to the 99.5% VaR as outlined under the Solvency II regime (Christiansen and Niemeyer 2014). In doing so, the premium should be charged appropriately, taking into account risks related to the claims experience. The computation of loaded premium, Π according to the mean and standard deviation premium principle are given by the following:
Π T = E Z T + κ E Z T
Π T = E Z T + κ Var Z T .
with κ [ 0 , 1 ] as the loading factor.
The results presented in Table 11 and Figure 7 show that an insurer needs to charge higher premiums and set aside more capital requirement if the claims experience exhibits overdispersion (underdispersion) with negative (positive) θ , in comparison to equidispersed claims experience. These results are consistent with a Frank copula which is seen in Table 12 and Figure 8. However, the claims modeling based on FGM copula as represented in Table 13 and Figure 9, do not produce the same pattern because of that copula’s weak dependence structure as compared with the other two Archimedean copulae with greater dependency.
We plotted again the simulated premium and VaR across different Kendall’s τ to standardize the value on x-axis for comparison. The VaR and premium for each dispersion case exhibit the same value, respectively when τ = 0 as seen in Figure 10 and Figure 11. As τ changes from −1 to 1, the overdispersed (underdispersed) case illustrates the widest (narrowest) spread of values among all. On average, Frank (Clayton) copula illustrates higher value of premiums and VaR than the Clayton (Frank) copula when τ is negative (positive), while the values under FGM copula are illustrated between the two copulae.
In comparison to Mohd Ramli et al. (2019) which only illustrated the computations of the first two moments of a Weibull risk process with λ = 2 under an FGM copula, we used the Monte Carlo simulation technique and obtained the higher-order moments of the Weibull risk process with wider choice of copulae. With the additional information on the asymmetry and the tail distribution of the aggregate discounted claims, insurance providers can allow for the dispersion in the claims occurrences and dependency between variables when estimating the premium and capital requirement for an insurance risk portfolio.

4. Conclusions

The ability of an insurance firm to manage surplus effectively is important to ensure capital adequacy, as required by the Solvency II standard. To allow for dispersed datasets, we propose the Weibull counting process to represent the claims arrival process of an insurance risk portfolio and an arbitrary continuous claims size. Under this model, two types of dependency were accommodated which were captured by the hazard function and the copula. We then fit the model to the historical data of New Zealand catastrophe events occurrences and losses, which was best represented by an overdispersed Weibull process with claim amounts following a lognormal distribution and a Clayton copula with a weak dependency. The moments estimation via Monte Carlo simulations shows the highest mean and variance of the aggregate discounted claims were obtained under an overdispersed claim arrivals with a strong negative correlation. The positive skewness and the high value of kurtosis of the aggregate discounted claims with Weibull IWT also implied an asymmetric and a heavy-tailed risk distribution. The scenario analysis conducted indicated that insurers will need to charge higher premiums and set aside more capital requirements as the claims experience exhibits overdispersed (underdispersed) IWT with negative (positive) dependency captured by copulae, in comparison to claims experience with equidispersed IWT. Other than computing the respective ruin probability under the Weibull risk process and applying commonly used copulae in finance and insurance, we may also examine the applicability of the model on a random annuity plan with random payments being paid at random times.

Author Contributions

Conceptualization, S.N.M.R.; methodology, S.N.M.R. and Z.H.Z.; software, Z.H.Z.; validation, S.N.M.R. and Z.H.Z.; formal analysis, S.F.S.Y.A.; investigation, S.F.S.Y.A.; resources, S.N.M.R.; data curation, Z.H.Z. and S.F.S.Y.A.; writing–original draft preparation, S.F.S.Y.A.; writing–review and editing, S.N.M.R. and Z.H.Z.; visualization, S.F.S.Y.A.; supervision, S.N.M.R. and Z.H.Z.; project administration, S.N.M.R. and Z.H.Z.; funding acquisition, S.N.M.R. and Z.H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universiti Kebangsaan Malaysia Grant Number GUP-2020-033.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: http://cas.uqam.ca/, http://dutangc.perso.math.cnrs.fr/RRepository/ (accessed on 15 March 2021).

Conflicts of Interest

The authors declare no conflict of interest. The funder had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

ADAnderson-Darling
AICAkaike information criterion
BICBayesian information criterion
CDFCumulative distribution function
FGMFarlie-Gumbel-Mogenstern
GOFGoodness-of-fit
IWTInterwaiting time
KSKolmogorov–Smirnov
MDPIMultidisciplinary Digital Publishing Institute
SCRSolvency capital requirement
Std DevStandard deviation
VaRValue-at-Risk

References

  1. Albrecher, Hansjörg, Corina Constantinescu, and Stéphane Loisel. 2011. Explicit ruin formulas for models with dependence among risks. Insurance: Mathematics and Economics 48: 265–70. [Google Scholar] [CrossRef] [Green Version]
  2. Bar-Lev, Shaul, and Ad Ridder. 2019. Monte Carlo methods for insurance risk computation. International Journal of Statistics and Probability 8: 54–74. [Google Scholar] [CrossRef]
  3. Barges, Mathieu, Hélene Cossette, Stéphane Loisel, and Etienne Marceau. 2011. On the moments of aggregate discounted claims with dependence introduced by a FGM copula. ASTIN Bulletin: The Journal of the IAA 41: 215–38. [Google Scholar]
  4. Boshnakov, Georgi, Tarak Kharrat, and Ian G. McHale. 2017. A bivariate Weibull count model for forecasting association football scores. International Journal of Forecasting 33: 458–66. [Google Scholar] [CrossRef] [Green Version]
  5. Casarano, Giuseppe, Gilberto Castellani, Luca Passalacqua, Francesca Perla, and Paolo Zanetti. 2017. Relevant applications of Monte Carlo simulation in Solvency II. Soft Computing 21: 1181–92. [Google Scholar] [CrossRef]
  6. Christiansen, Marcus C., and Andreas Niemeyer. 2014. Fundamental definition of the solvency capital requirement in Solvency II. Astin Bulletin 44: 501–33. [Google Scholar] [CrossRef] [Green Version]
  7. Cullen, Alison C., H. Christopher Frey, and Christopher H. Frey. 1999. Probabilistic Techniques in Exposure Assessment: A Handbook for Dealing with Variability and Uncertainty in Models and Inputs. Berlin: Springer Science & Business Media. [Google Scholar]
  8. Driels, Morris R., and Young S. Shin. 2004. Determining the Number of Iterations for Monte Carlo Simulations of Weapon Effectiveness. Technical Report. Monterey: Naval Postgraduate School Monterey. [Google Scholar]
  9. Dutang, Christophe, and Arthur Charpentier. 2020. CASdatasets: Insurance Datasets. R Package Version 1.0-11. Vienna: R Core Team. [Google Scholar]
  10. Hasumi, Tomohiro, Takuma Akimoto, and Yoji Aizawa. 2009. The Weibull–log Weibull distribution for inter-occurrence times of earthquakes. Physica A: Statistical Mechanics and Its Applications 388: 491–98. [Google Scholar] [CrossRef] [Green Version]
  11. Jang, Jiwook. 2004. Martingale approach for moments of discounted aggregate claims. Journal of Risk and Insurance 71: 201–11. [Google Scholar] [CrossRef]
  12. Jose, Kanichukattu K., and Bindu Abraham. 2011. A count model based on Mittag-Leffler interarrival times. Statistica 71: 501–14. [Google Scholar]
  13. Kelly, Dana L. 2007. Using copulas to model dependence in simulation risk assessment. Paper presented at the ASME 2007 International Mechanical Engineering Congress and Exposition, Seattle, WA, USA, November 11–15; vol. 43084, pp. 81–89. [Google Scholar]
  14. Klugman, Stuart A., and Rahul Parsa. 1999. Fitting bivariate loss distributions with copulas. Insurance: Mathematics and Economics 24: 139–48. [Google Scholar] [CrossRef]
  15. Kreer, Markus, Ayşe Kızılersü, Anthony W. Thomas, and Alfredo D Egídio dos Reis. 2015. Goodness-of-fit tests and applications for left-truncated Weibull distributions to non-life insurance. European Actuarial Journal 5: 139–63. [Google Scholar] [CrossRef]
  16. Léveillé, Ghislain, and Jose Garrido. 2001a. Moments of compound renewal sums with discounted claims. Insurance: Mathematics and Economics 28: 217–31. [Google Scholar]
  17. Léveillé, Ghislain, and José Garrido. 2001b. Recursive moments of compound renewal sums with discounted claims. Scandinavian Actuarial Journal 2001: 98–110. [Google Scholar] [CrossRef]
  18. Li, Shuanming, and Yi Lu. 2018. On the moments and the distribution of aggregate discounted claims in a Markovian environment. Risks 6: 59. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, Hanlin. 2019. Reliability and maintenance modeling for competing risk processes with Weibull inter-arrival shocks. Applied Mathematical Modelling 71: 194–207. [Google Scholar] [CrossRef]
  20. Lora, Mayra Ivanoff, and Julio M. Singer. 2011. Beta-binomial/gamma-Poisson regression models for repeated counts with random parameters. Brazilian Journal of Probability and Statistics 25: 218–35. [Google Scholar] [CrossRef]
  21. Ly, Sel, Kim-Hung Pho, Sal Ly, and Wing-Keung Wong. 2019. Determining distribution for the product of random variables by using copulas. Risks 7: 23. [Google Scholar] [CrossRef] [Green Version]
  22. Mao, Tiantian, and Fan Yang. 2015. Risk concentration based on expectiles for extreme risks under FGM copula. Insurance: Mathematics and Economics 64: 429–39. [Google Scholar] [CrossRef]
  23. McShane, Blake, Moshe Adrian, Eric T. Bradlow, and Peter S. Fader. 2008. Count models based on Weibull interarrival times. Journal of Business & Economic Statistics 26: 369–78. [Google Scholar]
  24. Mohd Ramli, Siti Norafidah, and Jiwook Jang. 2014. Neumann series on the recursive moments of copula-dependent aggregate discounted claims. Risks 2: 195–210. [Google Scholar] [CrossRef] [Green Version]
  25. Mohd Ramli, Siti Norafidah, Nur Atikah Mohamed Rozali, Sharifah Farah Syed Yusoff Alhabshi, and Ishak Hashim. 2018. Laplace transform on the recursive moments of copula-dependent aggregate discounted claims. In AIP Conference Proceedings. Melville: AIP Publishing LLC, vol. 1974, p. 020110. [Google Scholar]
  26. Mohd Ramli, Siti Norafidah, Nur Atikah Mohamed Rozali, Sharifah Farah Syed Yusoff Alhabshi, and Ishak Hashim. 2019. Laplace transform on the recursive moments of aggregate discounted claims with Weibull interwaiting time. In AIP Conference Proceedings. Melville: AIP Publishing LLC, vol. 2184, p. 050017. [Google Scholar]
  27. MunichRE. 2018. NatCatSERVICE. Available online: http://natcatservice.munichre.com (accessed on 14 February 2018).
  28. Razak, Nor Iza Anuar, Zamira Hasanah Zamzuri, and Nur Riza Mohd Suradi. 2019. The Implementation Of Double Bootstrap Method In Structural Equation Modeling. Asm Science Journal 12: 8–14. [Google Scholar]
  29. Sun, Weiwei, Xiang Hu, and Lianzeng Zhang. 2020. Moments of discounted aggregate claims with dependence based on Spearman copula. Journal of Computational and Applied Mathematics 377: 112889. [Google Scholar] [CrossRef]
  30. Waters, Howard R. 1983. Probability of ruin for a risk process with claims cost inflation. Scandinavian Actuarial Journal 1983: 148–64. [Google Scholar] [CrossRef]
  31. Winkelmann, Rainer. 1995. Duration dependence and dispersion in count-data models. Journal of Business & Economic Statistics 13: 467–74. [Google Scholar]
  32. Woo, Jae-Kyung, and Eric C. K. Cheung. 2013. A note on discounted compound renewal sums under dependency. Insurance: Mathematics and Economics 52: 170–79. [Google Scholar] [CrossRef]
  33. Yang, Hailiang, and Lihong Zhang. 2001. On the distribution of surplus immediately after ruin under interest force. Insurance: Mathematics and Economics 29: 247–55. [Google Scholar] [CrossRef]
  34. Zamzuri, Zamira Hasanah, and Gwee Jia Hui. 2020. Comparing and forecasting using stochastic mortality models: A Monte Carlo simulation. Sains Malaysiana 49: 2013–22. [Google Scholar] [CrossRef]
Figure 1. Weibull hazard function with different shape parameters to indicate different dispersions.
Figure 1. Weibull hazard function with different shape parameters to indicate different dispersions.
Risks 09 00109 g001
Figure 2. The simulated scatter plots of Clayton, Frank and FGM copulae under different dependence parameters. (a) Clayton with θ = 1 ; (b) Clayton with θ = 0.7 ; (c) Clayton with θ = 0 ; (d) Clayton with θ = 2 ; (e) Clayton with θ = 100 ; (f) Frank with θ = 100 ; (g) Frank with θ = 6 ; (h) Frank with θ = 0 ; (i) Frank with θ = 6 ; (j) Frank with θ = 100 ; (k) FGM with θ = 1 ; (l) FGM with θ = 0.5 ; (m) FGM with θ = 0 ; (n) FGM with θ = 0.5 ; (o) FGM with θ = 1 .
Figure 2. The simulated scatter plots of Clayton, Frank and FGM copulae under different dependence parameters. (a) Clayton with θ = 1 ; (b) Clayton with θ = 0.7 ; (c) Clayton with θ = 0 ; (d) Clayton with θ = 2 ; (e) Clayton with θ = 100 ; (f) Frank with θ = 100 ; (g) Frank with θ = 6 ; (h) Frank with θ = 0 ; (i) Frank with θ = 6 ; (j) Frank with θ = 100 ; (k) FGM with θ = 1 ; (l) FGM with θ = 0.5 ; (m) FGM with θ = 0 ; (n) FGM with θ = 0.5 ; (o) FGM with θ = 1 .
Risks 09 00109 g002
Figure 3. Skewness-kurtosis plot for the New Zealand catastrophe claims sizes.
Figure 3. Skewness-kurtosis plot for the New Zealand catastrophe claims sizes.
Risks 09 00109 g003
Figure 4. The theoretical and empirical plots for various distributions fitted to the New Zealand catastrophe claim sizes.
Figure 4. The theoretical and empirical plots for various distributions fitted to the New Zealand catastrophe claim sizes.
Risks 09 00109 g004
Figure 5. The theoretical and empirical plots for Weibull and exponential distributions fitted to the New Zealand catastrophe IWT.
Figure 5. The theoretical and empirical plots for Weibull and exponential distributions fitted to the New Zealand catastrophe IWT.
Risks 09 00109 g005
Figure 6. Probability distribution of claims frequencies under three dispersion cases.
Figure 6. Probability distribution of claims frequencies under three dispersion cases.
Risks 09 00109 g006
Figure 7. Premium amount and VaR under a Clayton copula across different values for θ .
Figure 7. Premium amount and VaR under a Clayton copula across different values for θ .
Risks 09 00109 g007
Figure 8. Premium amount and VaR under a Frank copula across different values for θ .
Figure 8. Premium amount and VaR under a Frank copula across different values for θ .
Risks 09 00109 g008
Figure 9. Premium amount and VaR under a FGM copula across different values for θ .
Figure 9. Premium amount and VaR under a FGM copula across different values for θ .
Risks 09 00109 g009
Figure 10. Premium amount under all copulae across different values for τ .
Figure 10. Premium amount under all copulae across different values for τ .
Risks 09 00109 g010aRisks 09 00109 g010b
Figure 11. VaR under all copulae across different values for τ .
Figure 11. VaR under all copulae across different values for τ .
Risks 09 00109 g011
Table 1. Verification of the Monte Carlo Simulation.
Table 1. Verification of the Monte Carlo Simulation.
Moments θ Monte CarloLaplace TransformRelative Deviation
−0.999479.23477.660.330%
−0.9474.83475.230.084%
−0.5469.16465.430.803%
Mean0452.79453.170.084%
0.5440.59440.920.074%
0.9433.43431.120.537%
0.999430.23428.690.359%
−0.999106,554.00106,351.840.190%
−0.9105,554.20103,929.501.563%
−0.594,099.8994,253.780.163%
Variance080,630.8582,420.232.171%
0.569,601.1170,874.441.797%
0.961,217.1461,845.861.017%
0.99959.182.1359,638.740.766%
Table 2. Descriptive statistics for New Zealand catastrophe claims sizes.
Table 2. Descriptive statistics for New Zealand catastrophe claims sizes.
MinMaxMedianMeanEstimated Std DevEstimated SkewnessEstimated Kurtosis
0.011123.610.6725217.145092.99779214.29496
Table 3. Goodness-of-fit (GOF) measures for fitted distributions on the claim sizes.
Table 3. Goodness-of-fit (GOF) measures for fitted distributions on the claim sizes.
GOF CriterionGOF Test
Distribution−2 Log LikelihoodAICBICKS Test
(p-Value)
AD Test
(p-Value)
Exponential828.447830.447833.2600.00000.0000
Lognormal787.928791.928797.5520.99800.9775
Gamma804.195808.195813.8200.09120.0389
Weibull796.037800.037805.6620.42910.1797
Pareto788.541792.541798.1650.90960.8201
Burr788.420794.421802.8570.93390.8521
Table 4. Descriptive statistics for the IWT (in years) of claims arrival.
Table 4. Descriptive statistics for the IWT (in years) of claims arrival.
MinMaxMedianMeanEstimated
Std Dev
Estimated
Skewness
Estimated
Kurtosis
0.00277.39220.14780.37860.81256.081949.9694
Table 5. Goodness-of-fit (GOF) measures for fitted distributions on the IWT.
Table 5. Goodness-of-fit (GOF) measures for fitted distributions on the IWT.
GOF CriterionGOF Test
DistributionEstimated
Parameters
AICBICKS Test
(p-Value)
AD Test
(p-Value)
Weibull λ = 0.700152
β = 0.282022
−27.4824−21.85800.53830.2877
Exponential α = 2.6411389.0824111.89460.00010.0000
Table 6. Goodness-of-fit criterion for fitted copula on the dependency between the claim sizes and the IWT.
Table 6. Goodness-of-fit criterion for fitted copula on the dependency between the claim sizes and the IWT.
CopulaEstimated Parameter, θ Log LikelihoodAICBIC
Independence0000
Gaussian0.192.02−2.040.77
T0.192.12−0.255.38
Clayton0.323.57−5.13−2.32
Gumbel1.090.920.172.98
Frank1.11.91−1.811
Joe1.060.171.654.47
Table 7. Simulated moments of the New Zealand catastrophe risk process for T = 5.
Table 7. Simulated moments of the New Zealand catastrophe risk process for T = 5.
IWT DistributionMeanVarianceSkewnessKurtosis
Weibull169.34321313.5406.296145.444
Exponential154.03418239.7907.242150.931
Table 8. Simulated VaR and premium of the New Zealand catastrophic risk process for T = 5.
Table 8. Simulated VaR and premium of the New Zealand catastrophic risk process for T = 5.
IWTVaR ($m)Premium Amount ($m)
Distribution95%99.50%Mean PrincipleStd Dev Principle
Weibull406.757846.901186.278183.942
Exponential363.706809.141169.438167.540
Table 9. Marginal mean and variance of overdispersed, equidispersed and underdispersed Weibull count models, Weibull C when T = 5.
Table 9. Marginal mean and variance of overdispersed, equidispersed and underdispersed Weibull count models, Weibull C when T = 5.
OverdispersedEquidispersedUnderdispersed
Weibull C
( λ = 0.5 , β = 0.945 )
Weibull C
( λ = 1 , β = 1.17835 )
Weibull C
( λ = 2 , β = 1.5 )
Mean4.2434.2434.243
Variance9.9034.2431.375
Table 10. Simulated moments under various dispersion of claims occurrences across different θ under a Clayton copula when T = 5.
Table 10. Simulated moments under various dispersion of claims occurrences across different θ under a Clayton copula when T = 5.
Overdispersed Weibull C ( λ = 0.5 , β = 0.945 )
θ MeanVarianceSkewnessKurtosis
−166.35910202.7509.423276.883
−0.558.9337595.1397.010128.587
053.7647168.82410.130294.856
0.3251.2667538.55622.3921599.203
1026.7221205.30316.188714.862
5019.757109.5212.620105.380
10019.45490.367−0.336−0.839
Equidispersed Weibull C ( λ = 1 , β = 1.17835 )
θ MeanVarianceSkewnessKurtosis
−160.6638240.1787.779162.108
−0.554.6576749.6118.895190.364
050.7515835.95010.186315.298
0.3248.3895557.51213.696576.241
1034.7682378.91813.372.840
5028.978855.98321.9681236.261
10027.540533.32225.7541589.694
Underdispersed Weibull C ( λ = 2 , β = 1.5 )
θ MeanVarianceSkewnessKurtosis
−147.4146147.8938.180148.807
−0.542.0144487.10011.205347.476
039.5873991.3589.331190.339
0.3238.6724477.05213.927431.441
1032.0962714.67414.606497.202
5030.0362288.30422.5671324.124
10029.6232289.05725.2801591.957
Table 11. Simulated VaR and premium ( κ = 0.1 ) under various dispersion of claims occurrences across different θ under a Clayton copula when T = 5.
Table 11. Simulated VaR and premium ( κ = 0.1 ) under various dispersion of claims occurrences across different θ under a Clayton copula when T = 5.
Overdispersed Weibull C ( λ = 0.5 , β = 0.945 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−1225.744573.57772.99476.459
−0.5197.729498.65864.82767.649
0180.155473.74759.14062.231
0.32167.568454.26456.39259.948
1063.731199.87629.39430.194
5034.14741.33521.73320.804
10032.85537.79321.39920.405
Equidispersed Weibull C ( λ = 1 , β = 1.17835 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−1202.152514.11866.72969.741
−0.5176.053466.06760.12262.872
0161.369443.09455.82758.391
0.32151.449424.32653.22855.844
1088.749288.37938.24539.646
5058.319170.82331.87631.904
10052.508120.99130.29329.849
Equidispersed Weibull C ( λ = 2 , β = 1.5 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−1163.186463.17652.15555.254
−0.5138.329380.70846.21548.712
0128.152374.84143.54645.905
0.32121.294364.80242.54045.363
1091.612298.82435.30637.307
5080.044269.07833.03934.819
10077.195256.81232.58534.407
Table 12. Simulated VaR and premium ( κ = 0.1 ) under various dispersion of claims occurrences across different θ under a Frank copula when T = 5.
Table 12. Simulated VaR and premium ( κ = 0.1 ) under various dispersion of claims occurrences across different θ under a Frank copula when T = 5.
Overdispersed Weibull C ( λ = 0.5 , β = 0.945 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−100225.816566.75173.36376.605
−50227.023554.46673.35676.737
−10222.041544.25372.45175.522
0179.770475.84359.57762.559
1074.381223.73631.72232.858
5034.97945.03821.94021.018
10033.04538.07221.42720.435
Equidispersed Weibull C ( λ = 1 , β = 1.17835 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−100203.215536.92567.42071.469
−50204.435532.85467.10070.079
−10200.478504.21466.21569.807
0161.197429.60855.87658.620
1094.393300.31939.51241.085
5058.900177.42231.99032.082
10052.942124.69630.38329.798
Underdispersed Weibull C ( λ = 2 , β = 1.5 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−100162.916459.14152.38755.664
−50161.985445.69152.35555.655
−10153.385432.80850.56053.662
0128.944368.50144.23247.049
1092.540297.59035.40537.388
5081.145275.84933.25435.059
10077.572252.96732.44533.838
Table 13. Simulated VaR and premium ( κ = 0.1 ) under various dispersion of claims occurrences across different θ under a FGM copula when T = 5.
Table 13. Simulated VaR and premium ( κ = 0.1 ) under various dispersion of claims occurrences across different θ under a FGM copula when T = 5.
Overdispersed Weibull C ( λ = 0.5 , β = 0.945 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−0.999204.390515.81466.56469.693
−0.5190.822509.91462.89166.048
0178.699480.09259.10662.356
0.5165.488448.83755.45158.282
0.999154.182426.93551.78954.361
Equidispersed Weibull C ( λ = 1 , β = 1.17835 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−0.999178.922492.66460.98864.017
−0.5172.109458.16358.60861.381
0161.823436.83856.02858.640
0.5151.357423.80953.24756.144
0.999140.298411.16850.55653.115
Underdispersed Weibull C ( λ = 2 , β = 1.5 )
θ VaR ($m)Premium Amount ($m)
95%99.50%Mean PrincipleStd Dev Principle
−0.999140.496403.16946.94049.749
−0.5133.912392.82645.23147.897
0126.499372.11543.53846.174
0.5120.662367.65142.26144.904
0.999115.166349.60240.98243.725
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Syed Yusoff Alhabshi, S.F.; Zamzuri, Z.H.; Mohd Ramli, S.N. Monte Carlo Simulation of the Moments of a Copula-Dependent Risk Process with Weibull Interwaiting Time. Risks 2021, 9, 109. https://doi.org/10.3390/risks9060109

AMA Style

Syed Yusoff Alhabshi SF, Zamzuri ZH, Mohd Ramli SN. Monte Carlo Simulation of the Moments of a Copula-Dependent Risk Process with Weibull Interwaiting Time. Risks. 2021; 9(6):109. https://doi.org/10.3390/risks9060109

Chicago/Turabian Style

Syed Yusoff Alhabshi, Sharifah Farah, Zamira Hasanah Zamzuri, and Siti Norafidah Mohd Ramli. 2021. "Monte Carlo Simulation of the Moments of a Copula-Dependent Risk Process with Weibull Interwaiting Time" Risks 9, no. 6: 109. https://doi.org/10.3390/risks9060109

APA Style

Syed Yusoff Alhabshi, S. F., Zamzuri, Z. H., & Mohd Ramli, S. N. (2021). Monte Carlo Simulation of the Moments of a Copula-Dependent Risk Process with Weibull Interwaiting Time. Risks, 9(6), 109. https://doi.org/10.3390/risks9060109

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