Source author record

Pavel V. Shevchenko

Pavel V. Shevchenko appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

Cyber Risk Frequency, Severity and Insurance Viability

In this study an exploration of insurance risk transfer is undertaken for the cyber insurance industry in the United States of America, based on the leading industry dataset of cyber events provided by Advisen. We seek to address two core unresolved questions. First, what factors are the most significant covariates that may explain the frequency and severity of cyber loss events and are they heterogeneous over cyber risk categories? Second, is cyber risk insurable in regards to the required premiums, risk pool sizes and how would this decision vary with the insured companies industry sector and size? We address these questions through a combination of regression models based on the class of Generalised Additive Models for Location Shape and Scale (GAMLSS) and a class of ordinal regressions. These models will then form the basis for our analysis of frequency and severity of cyber risk loss processes. We investigate the viability of insurance for cyber risk using a utility modelling framework with premium calculated by classical certainty equivalence analysis utilising the developed regression models. Our results provide several new key insights into the nature of insurability of cyber risk and rigorously address the two insurance questions posed in a real data driven case study analysis.

preprint2022arXiv

The Nature of Losses from Cyber-Related Events: Risk Categories and Business Sectors

In this study we examine the nature of losses from cyber related events across different risk categories and business sectors. Using a leading industry dataset of cyber events, we evaluate the relationship between the frequency and severity of individual cyber-related events and the number of affected records. We find that the frequency of reported cyber related events has substantially increased between 2008 and 2016. Furthermore, the frequency and severity of losses depend on the business sector and type of cyber threat: the most significant cyber loss event categories, by number of events, were related to data breaches and the unauthorized disclosure of data, while cyber extortion, phishing, spoofing and other social engineering practices showed substantial growth rates. Interestingly, we do not find a distinct pattern between the frequency of events, the loss severity, and the number of affected records as often alluded to in the literature. We also analyse the severity distribution of cyber related events across all risk categories and business sectors. This analysis reveals that cyber risks are heavy-tailed, i.e., cyber risk events have a higher probability to produce extreme losses than events whose severity follows an exponential distribution. Furthermore, we find that the frequency and severity of cyber related losses exhibits a very dynamic and time varying nature.

preprint2019arXiv

Optimal life-cycle consumption and investment decisions under age-dependent risk preferences

In this article we solve the problem of maximizing the expected utility of future consumption and terminal wealth to determine the optimal pension or life-cycle fund strategy for a cohort of pension fund investors. The setup is strongly related to a DC pension plan where additionally (individual) consumption is taken into account. The consumption rate is subject to a time-varying minimum level and terminal wealth is subject to a terminal floor. Moreover, the preference between consumption and terminal wealth as well as the intertemporal coefficient of risk aversion are time-varying and therefore depend on the age of the considered pension cohort. The optimal consumption and investment policies are calculated in the case of a Black-Scholes financial market framework and hyperbolic absolute risk aversion (HARA) utility functions. We generalize Ye (2008) (2008 American Control Conference, 356-362) by adding an age-dependent coefficient of risk aversion and extend Steffensen (2011) (Journal of Economic Dynamics and Control, 35(5), 659-667), Hentschel (2016) (Doctoral dissertation, Ulm University) and Aase (2017) (Stochastics, 89(1), 115-141) by considering consumption in combination with terminal wealth and allowing for consumption and terminal wealth floors via an application of HARA utility functions. A case study on fitting several models to realistic, time-dependent life-cycle consumption and relative investment profiles shows that only our extended model with time-varying preference parameters provides sufficient flexibility for an adequate fit. This is of particular interest to life-cycle products for (private) pension investments or pension insurance in general.

preprint2016arXiv

A unified approach to mortality modelling using state-space framework: characterisation, identification, estimation and forecasting

This paper explores and develops alternative statistical representations and estimation approaches for dynamic mortality models. The framework we adopt is to reinterpret popular mortality models such as the Lee-Carter class of models in a general state-space modelling methodology, which allows modelling, estimation and forecasting of mortality under a unified framework. Furthermore, we propose an alternative class of model identification constraints which is more suited to statistical inference in filtering and parameter estimation settings based on maximization of the marginalized likelihood or in Bayesian inference. We then develop a novel class of Bayesian state-space models which incorporate apriori beliefs about the mortality model characteristics as well as for more flexible and appropriate assumptions relating to heteroscedasticity that present in observed mortality data. We show that multiple period and cohort effect can be cast under a state-space structure. To study long term mortality dynamics, we introduce stochastic volatility to the period effect. The estimation of the resulting stochastic volatility model of mortality is performed using a recent class of Monte Carlo procedure specifically designed for state and parameter estimation in Bayesian state-space models, known as the class of particle Markov chain Monte Carlo methods. We illustrate the framework we have developed using Danish male mortality data, and show that incorporating heteroscedasticity and stochastic volatility markedly improves model fit despite an increase of model complexity. Forecasting properties of the enhanced models are examined with long term and short term calibration periods on the reconstruction of life tables.

preprint2016arXiv

Crunching Mortality and Life Insurance Portfolios with extended CreditRisk+

Using an extended version of the credit risk model CreditRisk+, we develop a flexible framework with numerous applications amongst which we find stochastic mortality modelling, forecasting of death causes as well as profit and loss modelling of life insurance and annuity portfolios which can be used in (partial) internal models under Solvency II. Yet, there exists a fast and numerically stable algorithm to derive loss distributions exactly, even for large portfolios. We provide various estimation procedures based on publicly available data. Compared to the Lee-Carter model, we have a more flexible framework, get tighter bounds and can directly extract several sources of uncertainty. Straight-forward model validation techniques are available.

preprint2016arXiv

Optimal Consumption, Investment and Housing with Means-tested Public Pension in Retirement

In this paper, we develop an expected utility model for the retirement behavior in the decumulation phase of Australian retirees with sequential family status subject to consumption, housing, investment, bequest and government provided means-tested Age Pension. We account for mortality risk and risky investment assets, and introduce a health proxy to capture the decreasing level of consumption for older retirees. Then we find optimal housing at retirement, and optimal consumption and optimal risky asset allocation depending on age and wealth. The model is solved numerically as a stochastic control problem, and is calibrated using the maximum likelihood method on empirical data of consumption and housing from the Australian Bureau of Statistics 2009-2010 Survey. The model fits the characteristics of the data well to explain the behavior of Australian retirees. The key findings are the following: First, the optimal policy is highly sensitive to means-tested Age Pension early in retirement but this sensitivity fades with age. Secondly, the allocation to risky assets shows a complex relationship with the means-tested Age Pension that disappears once minimum withdrawal rules are enforced. As a general rule, when wealth decreases the proportion allocated to risky assets increases, due to the Age Pension working as a buffer against investment losses. Finally, couples can be more aggressive with risky allocations due to their longer life expectancy compared with singles.

preprint2016arXiv

Should the advanced measurement approach be replaced with the standardized measurement approach for operational risk?

Recently, Basel Committee for Banking Supervision proposed to replace all approaches, including Advanced Measurement Approach (AMA), for operational risk capital with a simple formula referred to as the Standardised Measurement Approach (SMA). This paper discusses and studies the weaknesses and pitfalls of SMA such as instability, risk insensitivity, super-additivity and the implicit relationship between SMA capital model and systemic risk in the banking sector. We also discuss the issues with closely related operational risk Capital-at-Risk (OpCar) Basel Committee proposed model which is the precursor to the SMA. In conclusion, we advocate to maintain the AMA internal model framework and suggest as an alternative a number of standardization recommendations that could be considered to unify internal modelling of operational risk. The findings and views presented in this paper have been discussed with and supported by many OpRisk practitioners and academics in Australia, Europe, UK and USA, and recently at OpRisk Europe 2016 conference in London.

preprint2016arXiv

The 2015-2017 policy changes to the means-tests of Australian Age Pension: implication to decisions in retirement

The Australian Government uses the means-test as a way of managing the pension budget. Changes in Age Pension policy impose difficulties in retirement modelling due to policy risk, but any major changes tend to be `grandfathered' meaning that current retirees are exempt from the new changes. In 2015, two important changes were made in regards to allocated pension accounts -- the income means-test is now based on deemed income rather than account withdrawals, and the income-test deduction no longer applies. We examine the implications of the new changes in regards to optimal decisions for consumption, investment, and housing. We account for regulatory minimum withdrawal rules that are imposed by regulations on allocated pension accounts, as well as the 2017 asset-test rebalancing. The new policy changes are modelled in a utility maximizing lifecycle model and solved as an optimal stochastic control problem. We find that the new rules decrease the benefits from planning the consumption in relation to the means-test, while the housing allocation increases slightly in order to receive additional Age Pension. The difference in optimal drawdown between the old and new policy are only noticeable early in retirement until regulatory minimum withdrawal rates are enforced. However, the amount of extra Age Pension received for many households is now significantly different due to the new deeming income rules, which benefit slightly wealthier households who previously would receive no Age Pension due to the income-test and minimum withdrawals.

preprint2015arXiv

A State-Space Estimation of the Lee-Carter Mortality Model and Implications for Annuity Pricing

In this article we investigate a state-space representation of the Lee-Carter model which is a benchmark stochastic mortality model for forecasting age-specific death rates. Existing relevant literature focuses mainly on mortality forecasting or pricing of longevity derivatives, while the full implications and methods of using the state-space representation of the Lee-Carter model in pricing retirement income products is yet to be examined. The main contribution of this article is twofold. First, we provide a rigorous and detailed derivation of the posterior distributions of the parameters and the latent process of the Lee-Carter model via Gibbs sampling. Our assumption for priors is slightly more general than the current literature in this area. Moreover, we suggest a new form of identification constraint not yet utilised in the actuarial literature that proves to be a more convenient approach for estimating the model under the state-space framework. Second, by exploiting the posterior distribution of the latent process and parameters, we examine the pricing range of annuities, taking into account the stochastic nature of the dynamics of the mortality rates. In this way we aim to capture the impact of longevity risk on the pricing of annuities. The outcome of our study demonstrates that an annuity price can be more than 4% under-valued when different assumptions are made on determining the survival curve constructed from the distribution of the forecasted death rates. Given that a typical annuity portfolio consists of a large number of policies with maturities which span decades, we conclude that the impact of longevity risk on the accurate pricing of annuities is a significant issue to be further researched. In addition, we find that mis-pricing is increasingly more pronounced for older ages as well as for annuity policies having a longer maturity.

preprint2015arXiv

Forecasting Leading Death Causes in Australia using Extended CreditRisk$+$

Recently we developed a new framework in Hirz et al (2015) to model stochastic mortality using extended CreditRisk$^+$ methodology which is very different from traditional time series methods used for mortality modelling previously. In this framework, deaths are driven by common latent stochastic risk factors which may be interpreted as death causes like neoplasms, circulatory diseases or idiosyncratic components. These common factors introduce dependence between policyholders in annuity portfolios or between death events in population. This framework can be used to construct life tables based on mortality rate forecast. Moreover this framework allows stress testing and, therefore, offers insight into how certain health scenarios influence annuity payments of an insurer. Such scenarios may include improvement in health treatments or better medication. In this paper, using publicly available data for Australia, we estimate the model using Markov chain Monte Carlo method to identify leading death causes across all age groups including long term forecast for 2031 and 2051. On top of general reduced mortality, the proportion of deaths for certain certain causes has changed massively over the period 1987 to 2011. Our model forecasts suggest that if these trends persist, then the future gives a whole new picture of mortality for people aged above 40 years. Neoplasms will become the overall number-one death cause. Moreover, deaths due to mental and behavioural disorders are very likely to surge whilst deaths due to circulatory diseases will tend to decrease. This potential increase in deaths due to mental and behavioural disorders for older ages will have a massive impact on social systems as, typically, such patients need long-term geriatric care.

preprint2015arXiv

Sequential Monte Carlo Samplers for capital allocation under copula-dependent risk models

In this paper we assume a multivariate risk model has been developed for a portfolio and its capital derived as a homogeneous risk measure. The Euler (or gradient) principle, then, states that the capital to be allocated to each component of the portfolio has to be calculated as an expectation conditional to a rare event, which can be challenging to evaluate in practice. We exploit the copula-dependence within the portfolio risks to design a Sequential Monte Carlo Samplers based estimate to the marginal conditional expectations involved in the problem, showing its efficiency through a series of computational examples.

preprint2015arXiv

Valuation of Barrier Options using Sequential Monte Carlo

Sequential Monte Carlo (SMC) methods have successfully been used in many applications in engineering, statistics and physics. However, these are seldom used in financial option pricing literature and practice. This paper presents SMC method for pricing barrier options with continuous and discrete monitoring of the barrier condition. Under the SMC method, simulated asset values rejected due to barrier condition are re-sampled from asset samples that do not breach the barrier condition improving the efficiency of the option price estimator; while under the standard Monte Carlo many simulated asset paths can be rejected by the barrier condition making it harder to estimate option price accurately. We compare SMC with the standard Monte Carlo method and demonstrate that the extra effort to implement SMC when compared with the standard Monte Carlo is very little while improvement in price estimate can be significant. Both methods result in unbiased estimators for the price converging to the true value as $1/\sqrt{M}$, where $M$ is the number of simulations (asset paths). However, the variance of SMC estimator is smaller and does not grow with the number of time steps when compared to the standard Monte Carlo. In this paper we demonstrate that SMC can successfully be used for pricing barrier options. SMC can also be used for pricing other exotic options and also for cases with many underlying assets and additional stochastic factors such as stochastic volatility; we provide general formulas and references.

preprint2015arXiv

Valuation of Variable Annuities with Guaranteed Minimum Withdrawal and Death Benefits via Stochastic Control Optimization

In this paper we present a numerical valuation of variable annuities with combined Guaranteed Minimum Withdrawal Benefit (GMWB) and Guaranteed Minimum Death Benefit (GMDB) under optimal policyholder behaviour solved as an optimal stochastic control problem. This product simultaneously deals with financial risk, mortality risk and human behaviour. We assume that market is complete in financial risk and mortality risk is completely diversified by selling enough policies and thus the annuity price can be expressed as appropriate expectation. The computing engine employed to solve the optimal stochastic control problem is based on a robust and efficient Gauss-Hermite quadrature method with cubic spline. We present results for three different types of death benefit and show that, under the optimal policyholder behaviour, adding the premium for the death benefit on top of the GMWB can be problematic for contracts with long maturities if the continuous fee structure is kept, which is ordinarily assumed for a GMWB contract. In fact for some long maturities it can be shown that the fee cannot be charged as any proportion of the account value -- there is no solution to match the initial premium with the fair annuity price. On the other hand, the extra fee due to adding the death benefit can be charged upfront or in periodic instalment of fixed amount, and it is cheaper than buying a separate life insurance.

preprint2014arXiv

Fast and Simple Method for Pricing Exotic Options using Gauss-Hermite Quadrature on a Cubic Spline Interpolation

There is a vast literature on numerical valuation of exotic options using Monte Carlo, binomial and trinomial trees, and finite difference methods. When transition density of the underlying asset or its moments are known in closed form, it can be convenient and more efficient to utilize direct integration methods to calculate the required option price expectations in a backward time-stepping algorithm. This paper presents a simple, robust and efficient algorithm that can be applied for pricing many exotic options by computing the expectations using Gauss-Hermite integration quadrature applied on a cubic spline interpolation. The algorithm is fully explicit but does not suffer the inherent instability of the explicit finite difference counterpart. A `free' bonus of the algorithm is that it already contains the function for fast and accurate interpolation of multiple solutions required by many discretely monitored path dependent options. For illustrations, we present examples of pricing a series of American options with either Bermudan or continuous exercise features, and a series of exotic path-dependent options of target accumulation redemption note (TARN). Results of the new method are compared with Monte Carlo and finite difference methods, including some of the most advanced or best known finite difference algorithms in the literature. The comparison shows that, despite its simplicity, the new method can rival with some of the best finite difference algorithms in accuracy and at the same time it is significantly faster. Virtually the same algorithm can be applied to price other path-dependent financial contracts such as Asian options and variable annuities.

preprint2014arXiv

Holder-extendible European option: corrections and extensions

Financial contracts with options that allow the holder to extend the contract maturity by paying an additional fixed amount found many applications in finance. Closed-form solutions for the price of these options have appeared in the literature for the case when the contract underlying asset follows a geometric Brownian motion with the constant interest rate, volatility, and non-negative "dividend" yield. In this paper, the option price is derived for the case of the underlying asset that follows a geometric Brownian motion with the time-dependent drift and volatility which is important to use the solutions in real life applications. The formulas are derived for the drift that may include non-negative or negative "dividend" yield. The latter case results in a new solution type that has not been studied in the literature. Several typographical errors in the formula for the holder-extendible put, typically repeated in textbooks and software, are corrected.

preprint2014arXiv

Markov chain Monte Carlo estimation of default and recovery: dependent via the latent systematic factor

It is a well known fact that recovery rates tend to go down when the number of defaults goes up in economic downturns. We demonstrate how the loss given default model with the default and recovery dependent via the latent systematic risk factor can be estimated using Bayesian inference methodology and Markov chain Monte Carlo method. This approach is very convenient for joint estimation of all model parameters and latent systematic factors. Moreover, all relevant uncertainties are easily quantified. Typically available data are annual averages of defaults and recoveries and thus the datasets are small and parameter uncertainty is significant. In this case Bayesian approach is superior to the maximum likelihood method that relies on a large sample limit Gaussian approximation for the parameter uncertainty. As an example, we consider a homogeneous portfolio with one latent factor. However, the approach can be easily extended to deal with non-homogenous portfolios and several latent factors.

preprint2013arXiv

Loss Distribution Approach for Operational Risk Capital Modelling under Basel II: Combining Different Data Sources for Risk Estimation

The management of operational risk in the banking industry has undergone significant changes over the last decade due to substantial changes in operational risk environment. Globalization, deregulation, the use of complex financial products and changes in information technology have resulted in exposure to new risks very different from market and credit risks. In response, Basel Committee for banking Supervision has developed a regulatory framework, referred to as Basel II, that introduced operational risk category and corresponding capital requirements. Over the past five years, major banks in most parts of the world have received accreditation under the Basel II Advanced Measurement Approach (AMA) by adopting the loss distribution approach (LDA) despite there being a number of unresolved methodological challenges in its implementation. Different approaches and methods are still under hot debate. In this paper, we review methods proposed in the literature for combining different data sources (internal data, external data and scenario analysis) which is one of the regulatory requirement for AMA.

preprint2013arXiv

Optimal insurance purchase strategies via optimal multiple stopping times

In this paper we study a class of insurance products where the policy holder has the option to insure $k$ of its annual Operational Risk losses in a horizon of $T$ years. This involves a choice of $k$ out of $T$ years in which to apply the insurance policy coverage by making claims against losses in the given year. The insurance product structure presented can accommodate any kind of annual mitigation, but we present three basic generic insurance policy structures that can be combined to create more complex types of coverage. Following the Loss Distributional Approach (LDA) with Poisson distributed annual loss frequencies and Inverse-Gaussian loss severities we are able to characterize in closed form analytical expressions for the multiple optimal decision strategy that minimizes the expected Operational Risk loss over the next $T$ years. For the cases where the combination of insurance policies and LDA model does not lead to closed form expressions for the multiple optimal decision rules, we also develop a principled class of closed form approximations to the optimal decision rule. These approximations are developed based on a class of orthogonal Askey polynomial series basis expansion representations of the annual loss compound process distribution and functions of this annual loss.

preprint2013arXiv

Understanding Operational Risk Capital Approximations: First and Second Orders

We set the context for capital approximation within the framework of the Basel II / III regulatory capital accords. This is particularly topical as the Basel III accord is shortly due to take effect. In this regard, we provide a summary of the role of capital adequacy in the new accord, highlighting along the way the significant loss events that have been attributed to the Operational Risk class that was introduced in the Basel II and III accords. Then we provide a semi-tutorial discussion on the modelling aspects of capital estimation under a Loss Distributional Approach (LDA). Our emphasis is to focus on the important loss processes with regard to those that contribute most to capital, the so called high consequence, low frequency loss processes. This leads us to provide a tutorial overview of heavy tailed loss process modelling in OpRisk under Basel III, with discussion on the implications of such tail assumptions for the severity model in an LDA structure. This provides practitioners with a clear understanding of the features that they may wish to consider when developing OpRisk severity models in practice. From this discussion on heavy tailed severity models, we then develop an understanding of the impact such models have on the right tail asymptotics of the compound loss process and we provide detailed presentation of what are known as first and second order tail approximations for the resulting heavy tailed loss process. From this we develop a tutorial on three key families of risk measures and their equivalent second order asymptotic approximations: Value-at-Risk (Basel III industry standard); Expected Shortfall (ES) and the Spectral Risk Measure. These then form the capital approximations.

preprint2011arXiv

Bayesian Model Choice of Grouped t-copula

One of the most popular copulas for modeling dependence structures is t-copula. Recently the grouped t-copula was generalized to allow each group to have one member only, so that a priori grouping is not required and the dependence modeling is more flexible. This paper describes a Markov chain Monte Carlo (MCMC) method under the Bayesian inference framework for estimating and choosing t-copula models. Using historical data of foreign exchange (FX) rates as a case study, we found that Bayesian model choice criteria overwhelmingly favor the generalized t-copula. In addition, all the criteria also agree on the second most likely model and these inferences are all consistent with classical likelihood ratio tests. Finally, we demonstrate the impact of model choice on the conditional Value-at-Risk for portfolios of six major FX rates.

preprint2011arXiv

Calibration and filtering for multi factor commodity models with seasonality: incorporating panel data from futures contracts

We examine a general multi-factor model for commodity spot prices and futures valuation. We extend the multi-factor long-short model in Schwartz and Smith (2000) and Yan (2002) in two important aspects: firstly we allow for both the long and short term dynamic factors to be mean reverting incorporating stochastic volatility factors and secondly we develop an additive structural seasonality model. Then a Milstein discretized non-linear stochastic volatility state space representation for the model is developed which allows for futures and options contracts in the observation equation. We then develop numerical methodology based on an advanced Sequential Monte Carlo algorithm utilising Particle Markov chain Monte Carlo to perform calibration of the model jointly with the filtering of the latent processes for the long-short dynamics and volatility factors. In this regard we explore and develop a novel methodology based on an adaptive Rao-Blackwellised version of the Particle Markov chain Monte Carlo methodology. In doing this we deal accurately with the non-linearities in the state-space model which are therefore introduced into the filtering framework. We perform analysis on synthetic and real data for oil commodities.

preprint2011arXiv

Dependent default and recovery: MCMC study of downturn LGD credit risk model

There is empirical evidence that recovery rates tend to go down just when the number of defaults goes up in economic downturns. This has to be taken into account in estimation of the capital against credit risk required by Basel II to cover losses during the adverse economic downturns; the so-called "downturn LGD" requirement. This paper presents estimation of the LGD credit risk model with default and recovery dependent via the latent systematic risk factor using Bayesian inference approach and Markov chain Monte Carlo method. This approach allows joint estimation of all model parameters and latent systematic factor, and all relevant uncertainties. Results using Moody's annual default and recovery rates for corporate bonds for the period 1982-2010 show that the impact of parameter uncertainty on economic capital can be very significant and should be assessed by practitioners.

preprint2010arXiv

A Short Tale of Long Tail Integration

Integration of the form $\int_a^\infty {f(x)w(x)dx} $, where $w(x)$ is either $\sin (ω{\kern 1pt} x)$ or $\cos (ω{\kern 1pt} x)$, is widely encountered in many engineering and scientific applications, such as those involving Fourier or Laplace transforms. Often such integrals are approximated by a numerical integration over a finite domain $(a,\,b)$, leaving a truncation error equal to the tail integration $\int_b^\infty {f(x)w(x)dx} $ in addition to the discretization error. This paper describes a very simple, perhaps the simplest, end-point correction to approximate the tail integration, which significantly reduces the truncation error and thus increases the overall accuracy of the numerical integration, with virtually no extra computational effort. Higher order correction terms and error estimates for the end-point correction formula are also derived. The effectiveness of this one-point correction formula is demonstrated through several examples.

preprint2010arXiv

Calculation of aggregate loss distributions

Estimation of the operational risk capital under the Loss Distribution Approach requires evaluation of aggregate (compound) loss distributions which is one of the classic problems in risk theory. Closed-form solutions are not available for the distributions typically used in operational risk. However with modern computer processing power, these distributions can be calculated virtually exactly using numerical methods. This paper reviews numerical algorithms that can be successfully used to calculate the aggregate loss distributions. In particular Monte Carlo, Panjer recursion and Fourier transformation methods are presented and compared. Also, several closed-form approximations based on moment matching and asymptotic result for heavy-tailed distributions are reviewed.

preprint2010arXiv

Chain ladder method: Bayesian bootstrap versus classical bootstrap

The intention of this paper is to estimate a Bayesian distribution-free chain ladder (DFCL) model using approximate Bayesian computation (ABC) methodology. We demonstrate how to estimate quantities of interest in claims reserving and compare the estimates to those obtained from classical and credibility approaches. In this context, a novel numerical procedure utilising Markov chain Monte Carlo (MCMC), ABC and a Bayesian bootstrap procedure was developed in a truly distribution-free setting. The ABC methodology arises because we work in a distribution-free setting in which we make no parametric assumptions, meaning we can not evaluate the likelihood point-wise or in this case simulate directly from the likelihood model. The use of a bootstrap procedure allows us to generate samples from the intractable likelihood without the requirement of distributional assumptions, this is crucial to the ABC framework. The developed methodology is used to obtain the empirical distribution of the DFCL model parameters and the predictive distribution of the outstanding loss liabilities conditional on the observed claims. We then estimate predictive Bayesian capital estimates, the Value at Risk (VaR) and the mean square error of prediction (MSEP). The latter is compared with the classical bootstrap and credibility methods.

preprint2010arXiv

Computing Tails of Compound Distributions Using Direct Numerical Integration

An efficient adaptive direct numerical integration (DNI) algorithm is developed for computing high quantiles and conditional Value at Risk (CVaR) of compound distributions using characteristic functions. A key innovation of the numerical scheme is an effective tail integration approximation that reduces the truncation errors significantly with little extra effort. High precision results of the 0.999 quantile and CVaR were obtained for compound losses with heavy tails and a very wide range of loss frequencies using the DNI, Fast Fourier Transform (FFT) and Monte Carlo (MC) methods. These results, particularly relevant to operational risk modelling, can serve as benchmarks for comparing different numerical methods. We found that the adaptive DNI can achieve high accuracy with relatively coarse grids. It is much faster than MC and competitive with FFT in computing high quantiles and CVaR of compound distributions in the case of moderate to high frequencies and heavy tails.

preprint2010arXiv

Impact of Insurance for Operational Risk: Is it worthwhile to insure or be insured for severe losses?

Under the Basel II standards, the Operational Risk (OpRisk) advanced measurement approach allows a provision for reduction of capital as a result of insurance mitigation of up to 20%. This paper studies the behaviour of different insurance policies in the context of capital reduction for a range of possible extreme loss models and insurance policy scenarios in a multi-period, multiple risk settings. A Loss Distributional Approach (LDA) for modelling of the annual loss process, involving homogeneous compound Poisson processes for the annual losses, with heavy tailed severity models comprised of alpha-stable severities is considered. There has been little analysis of such models to date and it is believed, insurance models will play more of a role in OpRisk mitigation and capital reduction in future. The first question of interest is when would it be equitable for a bank or financial institution to purchase insurance for heavy tailed OpRisk losses under different insurance policy scenarios? The second question then pertains to Solvency II and addresses what the insurers capital would be for such operational risk scenarios under different policy offerings. In addition we consider the insurers perspective with respect to fair premium as a percentage above the expected annual claim for each insurance policy. The intention being to address questions related to VaR reduction under Basel II, SCR under Solvency II and fair insurance premiums in OpRisk for different extreme loss scenarios. In the process we provide closed form solutions for the distribution of loss process and claims process in an LDA structure as well as closed form analytic solutions for the Expected Shortfall, SCR and MCR under Basel II and Solvency II. We also provide closed form analytic solutions for the annual loss distribution of multiple risks including insurance mitigation.

preprint2010arXiv

The t copula with Multiple Parameters of Degrees of Freedom: Bivariate Characteristics and Application to Risk Management

The t copula is often used in risk management as it allows for modelling tail dependence between risks and it is simple to simulate and calibrate. However, the use of a standard t copula is often criticized due to its restriction of having a single parameter for the degrees of freedom (dof) that may limit its capability to model the tail dependence structure in a multivariate case. To overcome this problem, grouped t copula was proposed recently, where risks are grouped a priori in such a way that each group has a standard t copula with its specific dof parameter. In this paper we propose the use of a grouped t copula, where each group consists of one risk factor only, so that a priori grouping is not required. The copula characteristics in the bivariate case are studied. We explain simulation and calibration procedures, including a simulation study on finite sample properties of the maximum likelihood estimators and Kendall's tau approximation. This new copula can be significantly different from the standard t copula in terms of risk measures such as tail dependence, value at risk and expected shortfall. Keywords: grouped t copula, tail dependence, risk management.

preprint2009arXiv

Implementing Loss Distribution Approach for Operational Risk

To quantify the operational risk capital charge under the current regulatory framework for banking supervision, referred to as Basel II, many banks adopt the Loss Distribution Approach. There are many modeling issues that should be resolved to use the approach in practice. In this paper we review the quantitative methods suggested in literature for implementation of the approach. In particular, the use of the Bayesian inference method that allows to take expert judgement and parameter uncertainty into account, modeling dependence and inclusion of insurance are discussed.