In this paper, we define weighted failure rate and their means from the stand point of an application. We begin by emphasizing that the formation of n independent component series system having weighted failure rates with sum of weight functions being unity is same as a mixture of n distributions. We derive some parametric and non-parametric characterization results. We discuss on the form invariance property of baseline failure rate for a specific choice of weight function. Some bounds on means of aging functions are obtained. Here, we establish that weighted increasing failure rate average (IFRA) class is not closed under formation of coherent systems unlike the IFRA class. An interesting application of the present work is credited to the fact that the quantile version of means of failure rate is obtained as a special case of weighted means of failure rate.

This paper investigates the precise large deviations of the net loss process in a two-dimensional risk model with consistently varying tails and dependence structures, and gives some asymptotic formulas which hold uniformly for all x varying in t-intervals. The study is among the initial efforts to analyze potential risk via large deviation results for the net loss process of the two-dimensional risk model, and can provide a novel insight to assess the operation risk in a long run by fully considering the premium income factors of the insurance company.

In this paper, a new multivariate counting process model (called Multivariate Poisson Generalized Gamma Process) is developed and its main properties are studied. Some basic stochastic properties of the number of events in the new multivariate counting process are initially derived. It is shown that this new multivariate counting process model includes the multivariate generalized Pólya process as a special case. The dependence structure of the multivariate counting process model is discussed. Some results on multivariate stochastic comparisons are also obtained.

Aggregate implements an efficient fast Fourier transform (FFT)-based algorithm to approximate compound probability distributions. Leveraging FFT-based methods offers advantages over recursion and simulation-based approaches, providing speed and accuracy to otherwise time-consuming calculations. Combining user-friendly features and an expressive domain-specific language called DecL, Aggregate enables practitioners and nonprogrammers to work with complex distributions effortlessly. The software verifies the accuracy of its FFT-based numerical approximations by comparing their first three moments to those calculated analytically from the specified frequency and severity. This moment-based validation, combined with carefully chosen default parameters, allows users without in-depth knowledge of the underlying algorithm to be confident in the results. Aggregate supports a wide range of frequency and severity distributions, policy limits and deductibles, and reinsurance structures and has applications in pricing, reserving, risk management, teaching, and research. It is written in Python.

Generalized additive models (GAMs) are a leading model class for interpretable machine learning. GAMs were originally defined with smooth shape functions of the predictor variables and trained using smoothing splines. Recently, tree-based GAMs where shape functions are gradient-boosted ensembles of bagged trees were proposed, leaving the door open for the estimation of a broader class of shape functions (e.g. Explainable Boosting Machine (EBM)). In this paper, we introduce a competing three-step GAM learning approach where we combine (i) the knowledge of the way to split the covariates space brought by an additive tree model (ATM), (ii) an ensemble of predictive linear scores derived from generalized linear models (GLMs) using a binning strategy based on the ATM, and (iii) a final GLM to have a prediction model that ensures auto-calibration. Numerical experiments illustrate the competitive performances of our approach on several datasets compared to GAM with splines, EBM, or GLM with binarsity penalization. A case study in trade credit insurance is also provided.

We investigate bottom-up risk aggregation applied by insurance companies facing reserve risk from multiple lines of business. Since risk capitals should be calculated in different time horizons and calendar years, depending on the regulatory or reporting regime (Solvency II vs IFRS 17), we study correlations of ultimate losses and correlations of one-year losses in future calendar years in lines of business. We consider a multivariate version of a Hertig’s lognormal model and we derive analytical formulas for the ultimate correlation and the one-year correlations in future calendar years. Our main conclusion is that the correlation coefficients that should be used in a bottom-up aggregation formula depend on the time horizon and the future calendar year where the risk emerges. We investigate analytically and numerically properties of the ultimate and the one-year correlations, their possible values observed in practice, and the impact of misspecified correlations on the diversified risk capital.