Error analysis of finite difference and Markov chain approximations for option pricing

Mijatovi? and Pistorius proposed an efficient Markov chain approximation method for pricing European and barrier options in general one‐dimensional Markovian models. However, sharp convergence rates of this method for realistic financial payoffs, which are nonsmooth, are rarely available. In this paper, we solve this problem for general one‐dimensional diffusion models, which play a fundamental role in financial applications. For such models, the Markov chain approximation method is equivalent to the method of lines using the central difference. Our analysis is based on the spectral representation of the exact solution and the approximate solution. By establishing the convergence rate for the eigenvalues and the eigenfunctions, we obtain sharp convergence rates for the transition density and the price of options with nonsmooth payoffs. In particular, we show that for call‐/put‐type payoffs, convergence is second order, while for digital‐type payoffs, convergence is generally only first order. Furthermore, we provide theoretical justification for two well‐known smoothing techniques that can restore second‐order convergence for digital‐type payoffs and explain oscillations observed in the convergence for options with nonsmooth payoffs. As an extension, we also establish sharp convergence rates for European options for a rich class of Markovian jump models constructed from diffusions via subordination. The theoretical estimates are confirmed using numerical examples.     

Spatial statistical modelling of insurance risk: a spatial epidemiological approach to car insurance

Spatial models, such as the Besag, York and Mollie (BYM) model, have long been used in epidemiology and disease mapping. A common research question in these subjects is modelling the number of disease events per region; here the BYM models provides a holistic framework for both covariates and dependencies between regions. We use these tools to assess the relative insurance risk associated with the policyholders geographical location. A Bayesian modelling approach is presented and an elastic net is used to reduce the large number of possible geographic covariates. The final inference is performed using Integrated Nested Laplace Approximation. The model is applied to car insurance data from If P&C Insurance together with spatially referenced covariate data of high resolution, provided by Insightone. The entire analysis is performed using freely available R-packages. Including spatial dependence when modelling the number of claims significantly improves on the result obtained using ordinary generalised linear models. However, the support for adding a spatial component to the model for claims cost is weaker.     

Conditional risk measures in a bipartite market structure

In this paper, we study the effect of network structure between agents and objects on measures for systemic risk. We model the influence of sharing large exogeneous losses to the financial or (re)insurance market by a bipartite graph. Using Pareto-tailed losses and multivariate regular variation, we obtain asymptotic results for conditional risk measures based on the Value-at-Risk and the Conditional Tail Expectation. These results allow us to assess the influence of an individual institution on the systemic or market risk and vice versa through a collection of conditional risk measures. For large markets, Poisson approximations of the relevant constants are provided. Differences of the conditional risk measures for an underlying homogeneous and inhomogeneous random graph are illustrated by simulations.     

Saddlepoint approximations to the mean and variance of the extended hypergeometric distribution

Conditional inference on 2 x 2 tables with fixed margins and unequal probabilities is based on the extended hypergeometric distribution. If the support of the distribution is large, exact calculation of the conditional mean and variance of the table entry may be computationally demanding. This paper proposes a single‐saddlepoint approximation to the mean and variance. While the approximation achieves acceptable accuracy for ordinary practical purposes, an alternative saddlepoint approximation is provided that gives much closer to exact results. It improves the accuracy of current approximations to up to more than four powers of ten.     

A saddlepoint approximation to the distribution of the sum of independent non‐identically uniform random variables

Calculating the probability of the corresponding significance point is important for finite sample sizes. However, it is difficult to evaluate this probability when the sample sizes are moderate to large. Under these circumstances, consideration of a more accurate approximation for the distribution function is extremely important. Herein, we performed a saddlepoint approximation in the upper tails for the distribution of the sum of independent non‐identically uniform random variables under finite sample sizes. Saddlepoint approximation results were compared with those for a normal approximation. Additionally, the order of errors of the saddlepoint approximation was derived. © 2014 The Authors. Statistica Neerlandica © 2014 VVS.     

FIRE SALES FORENSICS: MEASURING ENDOGENOUS RISK

We propose a tractable framework for quantifying the impact of loss‐triggered fire sales on portfolio risk, in a multi‐asset setting. We derive analytical expressions for the impact of fire sales on the realized volatility and correlations of asset returns in a fire sales scenario and show that our results provide a quantitative explanation for the spikes in volatility and correlations observed during such deleveraging episodes. These results are then used to develop an econometric framework for the forensic analysis of fire sales episodes, using observations of market prices. We give conditions for the identifiability of model parameters from time series of asset prices, propose a statistical test for the presence of fire sales, and an estimator for the magnitude of fire sales in each asset class. Pathwise consistency and large sample properties of the estimator are studied in the high‐frequency asymptotic regime. We illustrate our methodology by applying it to the forensic analysis of two recent deleveraging episodes: the Quant Crash of August 2007 and the Great Deleveraging following the default of Lehman Brothers in Fall 2008.     

Barrier present value maximization for a diffusion model of insurance surplus

In this paper, we study a barrier present value (BPV) maximization problem for an insurance entity whose surplus process follows an arithmetic Brownian motion. The BPV is defined as the expected discounted value of a payment made at the time when the surplus process reaches a high barrier level. The insurance entity buys proportional reinsurance and invests in a Black–Scholes market to maximize the BPV. We show that the maximal BPV function is a classical solution to the corresponding Hamilton–Jacobi–Bellman equation and is three times continuously differentiable using a novel operator. Its associated optimal reinsurance-investment control policy is determined by verification techniques.     

Valuation of forward start options under affine jump-diffusion models

Under the general affine jump-diffusion framework of Duffie et al. [Econometrica, 2000, 68, 1343–1376], this paper proposes an alternative pricing methodology for European-style forward start options that does not require any parallel optimization routine to ensure square integrability. Therefore, the proposed methodology is shown to possess a better accuracy–efficiency trade-off than the usual and more general approach initiated by Hong [Forward Smile and Derivative Pricing. Working paper, UBS, 2004] that is based on the knowledge of the forward characteristic function. Explicit pricing solutions are also offered under the nested jump-diffusion setting proposed by Bakshi et al. [J. Finance, 1997, 52, 2003–2049], which accommodates stochastic volatility and stochastic interest rates, and different integration schemes are numerically tested.     

聚合信用风险模型的改进和研究

近年来的研究发现,违约损失率的分布呈现一种双峰特征。文章对传统聚合信用风险模型进行改进,用具有双峰特征的双beta分布来刻画违约损失率的变化,给出全部贷款组合信用风险概率生成函数的解析式,利用数值模拟的结果验证了模型的有效性。