Many pension plans and private retirement products contain annuity factors, converting the funds at some future time into lifelong income. In general model settings like, for example, the Li-Lee mortality model, analytical values for the annuity factors are not available and one has to rely on numerical techniques. Their computation typically requires nested simulations as they depend on the interest rate level and the mortality tables at the time of retirement. We exploit the flexibility and efficiency of feed-forward neural networks (NNs) to value the annuity factors at the time of retirement. In a numerical study, we compare our deep learning approach to (least-squares) Monte-Carlo, which can be represented as a special case of the NN.

Conditional risk measures and their associated risk contribution measures are commonly employed in finance and actuarial science for evaluating systemic risk and quantifying the effects of risk interactions. This paper introduces various types of contribution ratio measures based on the multivariate conditional value-at-risk (MCoVaR), multivariate conditional expected shortfall (MCoES), and multivariate marginal mean excess (MMME) studied in [34] (Ortega-Jiménez, P., Sordo, M., & Suárez-Llorens, A. (2021). Stochastic orders and multivariate measures of risk contagion. Insurance: Mathematics and Economics, vol. 96, 199–207) and [11] (Das, B., & Fasen-Hartmann, V. (2018). Risk contagion under regular variation and asymptotic tail independence. Journal of Multivariate Analysis 165(1), 194–215) to assess the relative effects of a single risk when other risks in a group are in distress. The properties of these contribution risk measures are examined, and sufficient conditions for comparing these measures between two sets of random vectors are established using univariate and multivariate stochastic orders and statistically dependent notions. Numerical examples are presented to validate these conditions. Finally, a real dataset from the cryptocurrency market is used to analyze the spillover effects through our proposed contribution measures.

In this article, we study an optimization problem for a couple including two breadwinners with uncertain life times. Both breadwinners need to choose the optimal strategies for consumption, investment, housing, and life insurance purchasing to maximize the utility. In this article, the prices of housing assets and investment risky assets are assumed to be correlated. These two breadwinners are considered to have dependent mortality rates to include the breaking heart effect. The method of copula functions is used to construct the joint survival functions of two breadwinners. The analytical solutions of optimal strategies can be achieved, and numerical results are demonstrated.

Decentralized consensus protocols have a variety of parameters to be set during their deployment for practical applications in blockchains. The analysis given in most research papers proves the security state of the blockchain, at the same time usually providing a range of acceptable values, thus allowing further tuning of the protocol parameters. In this paper, we investigate Ouroboros Praos, the proof-of-stake consensus protocol deployed in Cardano and other blockchains. In contrast to its predecessor, Praos allows multiple honest slot leaders that lead to fork creation and resolution, consequently decreasing the block rate per time unit. In our analysis of dependence on protocol parameters such as active slot coefficient and p2p network block propagation time, we obtain new theoretical results and explicit formulas for the expectation of the length of the longest chain created during the Praos epoch, the length of the longest unintentional fork created by honest slot leaders, the efficiency of block generation procedure (the ratio of blocks included in the final longest chain vs the total number of created blocks), and other characteristics of the blockchain throughput.

We study these parameters as stochastic characteristics of the block generation process. The model is described in terms of the two-parametric family ξij of independent Bernoulli random variables which generate deformation of the binomial distribution by a positive integer parameter—the delay (deterministic or random). An essential part of our paper is a study of this deformation in terms of denumerable Markov chains and generating functions.

This report explores key considerations in relation to adopting a dynamic discount rate funding approach and the impacts of doing so in a range of areas, including funding volatility, investment strategy and end game objectives. It considers the advantages and disadvantages of this approach from the perspective of a range of stakeholders and the challenges that need overcoming in order to fully implement and support the approach, for example data challenges and the new skills required in the industry. The report includes sample modelling to highlight the practical issues that arise when adopting this approach. It describes a step-by-step approach for assessing the risks to be considered when determining an appropriate level of assets to provide funding for a sample set of pension scheme cash flows, as summarised in the table below.

It also considers how a dynamic discount rate approach fits within the proposed future funding regulations. Finally, the report puts forward recommendations for the IFoA, Scheme Actuaries and TPR.

Consequences of schemes adopting a dynamic discount rate approach could include very different investment strategies with investment in a wider pool of assets, less use of leveraged Liability Driven Investment, fewer schemes targeting buy-out as their end game strategy and an increase in technical work for actuaries in advising on the optimisation of asset and liability cash flows.

We consider the moments and the distribution of hitting times on the lollipop graph which is the graph exhibiting the maximum expected hitting time among all the graphs having the same number of nodes. We obtain recurrence relations for the moments of all order and we use these relations to analyze the asymptotic behavior of the hitting time distribution when the number of nodes tends to infinity.

In the present technological age, where cyber-risk ranks alongside natural and man-made disasters and catastrophes – in terms of global economic loss – businesses and insurers alike are grappling with fundamental risk management issues concerning the quantification of cyber-risk, and the dilemma as to how best to mitigate this risk. To this end, the present research deals with data, analysis, and models with the aim of quantifying and understanding cyber-risk – often described as “holy grail” territory in the realm of cyber-insurance and IT security. Nonparametric severity models associated with cyber-related loss data – identified from several competing sources – and accompanying parametric large-loss components, are determined, and examined. Ultimately, in the context of analogous cyber-coverage, cyber-risk is quantified through various types and levels of risk adjustment for (pure-risk) increased limit factors, based on applications of actuarially founded aggregate loss models in the presence of various forms of correlation. By doing so, insight is gained into the nature and distribution of volatile severity risk, correlated aggregate loss, and associated pure-risk limit factors.

The cumulative residual extropy has been proposed recently as an alternative measure of extropy to the cumulative distribution function of a random variable. In this paper, the concept of cumulative residual extropy has been extended to cumulative residual extropy inaccuracy (CREI) and dynamic cumulative residual extropy inaccuracy (DCREI). Some lower and upper bounds for these measures are provided. A characterization problem for the DCREI measure under the proportional hazard rate model is studied. Nonparametric estimators for CREI and DCREI measures based on kernel and empirical methods are suggested. Also, a simulation study is presented to evaluate the performance of the suggested measures. Simulation results show that the kernel-based estimator performs better than the empirical-based estimator. Finally, applications of the DCREI measure for model selection are provided using two real data sets.

We propose Rényi information generating function (RIGF) and discuss its properties. A connection between the RIGF and the diversity index is proposed for discrete-type random variables. The relation between the RIGF and Shannon entropy of order q > 0 is established and several bounds are obtained. The RIGF of escort distribution is derived. Furthermore, we introduce the Rényi divergence information generating function (RDIGF) and discuss its effect under monotone transformations. We present nonparametric and parametric estimators of the RIGF. A simulation study is carried out and a real data relating to the failure times of electronic components is analyzed. A comparison study between the nonparametric and parametric estimators is made in terms of the standard deviation, absolute bias, and mean square error. We have observed superior performance for the newly proposed estimators. Some applications of the proposed RIGF and RDIGF are provided. For three coherent systems, we calculate the values of the RIGF and other well-established uncertainty measures, and similar behavior of the RIGF is observed. Further, a study regarding the usefulness of the RDIGF and RIGF as model selection criteria is conducted. Finally, three chaotic maps are considered and then used to establish a validation of the proposed information generating function.

This work studies the reliability function of K-out-of-N systems with a general repair time distribution and a single repair facility. It introduces a new repair mechanism using an effort function, described by a nonlinear ordinary differential equation. Three theoretical results are obtained: regularity properties preventing simultaneous failures and repairs, derivation of a Kolmogorov forward system for micro-state and macro-state probabilities, and comparison of reliability functions of two K-out-of-N systems. An additional hypothesis on the model’s parameters allows us to obtain an ordering relation between the reliability functions. A numerical example demonstrates the model’s practical application and confirms the theoretical results.

In this paper, we investigate the number of customers that overlap or coincide with a virtual customer in an Erlang-A queue. Our analysis starts with the fluid and diffusion limit differential equations to obtain the mean and variance of the queue length. We then develop precise approximations for waiting times using fluid limits and the polygamma function. Building on this, we introduce a novel approximation scheme to calculate the mean and variance of the number of overlapping customers. This method facilitates the assessment of transient overlap risks in complex service systems, offering a useful tool for service providers to mitigate significant overlaps during pandemic seasons.