This paper introduces gemact, a Python package for actuarial modeling based on the collective risk model. The library supports applications to risk costing and risk transfer, loss aggregation, and loss reserving. We add new probability distributions to those available in scipy, including the (a, b, 0) and (a, b, 1) discrete distributions, copulas of the Archimedean family, the Gaussian, the Student t and the Fundamental copulas. We provide an implementation of the AEP algorithm for calculating the cumulative distribution function of the sum of dependent, nonnegative random variables, given their dependency structure specified with a copula. The theoretical framework is introduced at the beginning of each section to give the reader with a sufficient understanding of the underlying actuarial models.

For some time now, Solvency II requires that insurance companies calculate minimum capital requirements to face the risk of insolvency, either in accordance with the Standard Formula or using a full or partial Internal Model. An Internal Model must be based on a market-consistent valuation of assets and liabilities at a 1-year time span, where a real-world probabilistic structure is used for the first year of projection. In this paper, we describe the major risks of a non-life insurance company, i.e. the non-life underwriting risk and market risk, and their interactions, focusing on the non-life premium risk, equity risk, and interest rate risk. This analysis is made using some well-known stochastic models in the financial-actuarial literature and practical insurance business, i.e. the Collective Risk Model for non-life premium risk, the Geometric Brownian Motion for equity risk, and a real-world version of the G2++ Model for interest rate risk, where parameters are calibrated on current and real market data. Finally, we illustrate a case study on a single-line and a multi-line insurance company in order to see how the risk drivers behave in both a stand-alone and an aggregate framework.

The distribution of the sum of dependent risks is a crucial aspect in actuarial sciences, risk management and in many branches of applied probability. In this paper, we obtain analytic expressions for the probability density function (pdf) and the cumulative distribution function (cdf) of aggregated risks, modelled according to a mixture of exponential distributions. We first review the properties of the multivariate mixture of exponential distributions, to then obtain the analytical formulation for the pdf and the cdf for the aggregated distribution. We study in detail some specific families with Pareto (Sarabia et al., 2016), gamma, Weibull and inverse Gaussian mixture of exponentials (Whitmore and Lee, 1991) claims. We also discuss briefly the computation of risk measures, formulas for the ruin probability (Albrecher et al., 2011) and the collective risk model. An extension of the basic model based on mixtures of gamma distributions is proposed, which is one of the suggested directions for future research.

We bring the recently developed framework of dependence uncertainty into collective risk models, one of the most classic models in actuarial science. We study the worst-case values of the Value-at-Risk (VaR) and the Expected Shortfall (ES) of the aggregate loss in collective risk models, under two settings of dependence uncertainty: (i) the counting random variable (claim frequency) and the individual losses (claim sizes) are independent, and the dependence of the individual losses is unknown; (ii) the dependence of the counting random variable and the individual losses is unknown. Analytical results for the worst-case values of ES are obtained. For the loss from a large portfolio of insurance policies, an asymptotic equivalence of VaR and ES is established. Our results can be used to provide approximations for VaR and ES in collective risk models with unknown dependence. Approximation errors are obtained in both cases.

This paper reviews Panjer's recursion formula for evaluation of compound distributions and illustrates how it can be applied to life and general insurance problems.

We propose practical solutions for the determination of optimal retentions in a stop-loss reinsurance. We develop two new optimization criteria for deriving the optimal retentions by, respectively, minimizing the value-at-risk (VaR) and the conditional tail expectation (CTE) of the total risks of an insurer. We establish necessary and sufficient conditions for the existence of the optimal retentions for two risk models: individual risk model and collective risk model. The resulting optimal solution of our optimization criterion has several important characteristics: (i) the optimal retention has a very simple analytic form; (ii) the optimal retention depends only on the assumed loss distribution and the reinsurer’s safety loading factor; (iii) the CTE criterion is more applicable than the VaR criterion in the sense that the optimal condition for the former is less restrictive than the latter; (iv) if optimal solutions exist, then both VaR- and CTE-based optimization criteria yield the same optimal retentions. In terms of applications, we extend the results to the individual risk models with dependent risks and use multivariate phase type distribution, multivariate Pareto distribution and multivariate Bernoulli distribution to illustrate the effect of dependence on optimal retentions. We also use the compound Poisson distribution and the compound negative binomial distribution to illustrate the optimal retentions in a collective risk model.

Recursions are derived for a class of compound distributions having a claim frequency distribution of the well known (a,b)-type. The probability mass function on which the recursions are usually based is replaced by the distribution function in order to obtain increasing iterates. A monotone transformation is suggested to avoid an underflow in the initial stages of the iteration. The faster increase of the transformed iterates is diminished by use of a scaling function. Further, an adaptive weighting depending on the initial value and the increase of the iterates is derived. It enables us to manage an arbitrary large portfolio. Some numerical results are displayed demonstrating the efficiency of the different methods. The computation of the stop-loss premiums using these methods are indicated. Finally, related iteration schemes based on the cumulative distribution function are outlined.