Unbiased and efficient Greeks of financial options

The price of a derivative security equals the discounted expected payoff of the security under a suitable measure, and Greeks are price sensitivities with respect to parameters of interest. When closed-form formulas do not exist, Monte Carlo simulation has proved very useful for computing the prices and Greeks of derivative securities. Although finite difference with resimulation is the standard method for estimating Greeks, it is in general biased and suffers from erratic behavior when the payoff function is discontinuous. Direct methods, such as the pathwise method and the likelihood ratio method, are proposed to differentiate the price formulas directly and hence produce unbiased Greeks (Broadie and Glasserman, Manag. Sci. 42:269–285, 1996). The pathwise method differentiates the payoff function, whereas the likelihood ratio method differentiates the densities. When both methods apply, the pathwise method generally enjoys lower variances, but it requires the payoff function to be Lipschitz-continuous. Similarly to the pathwise method, our method differentiates the payoff function but lifts the Lipschitz-continuity requirements on the payoff function. We build a new but simple mathematical formulation so that formulas of Greeks for a broad class of derivative securities can be derived systematically. We then present an importance sampling method to estimate the Greeks. These formulas are the first in the literature. Numerical experiments show that our method gives unbiased Greeks for several popular multi-asset options (also called rainbow options) and a path-dependent option.     

The stochastic collocation Monte Carlo sampler: highly efficient sampling from ‘expensive’ distributions

In this article, we propose an efficient approach for inverting computationally expensive cumulative distribution functions. A collocation method, called the Stochastic Collocation Monte Carlo sampler (SCMC sampler), within a polynomial chaos expansion framework, allows us the generation of any number of Monte Carlo samples based on only a few inversions of the original distribution plus independent samples from a standard normal variable. We will show that with this path-independent collocation approach the exact simulation of the Heston stochastic volatility model, as proposed in Broadie and Kaya [Oper. Res., 2006, 54, 217–231], can be performed efficiently and accurately. We also show how to efficiently generate samples from the squared Bessel process and perform the exact simulation of the SABR model.     

Pricing and hedging guaranteed minimum withdrawal benefits under a general Lévy framework using the COS method

This paper extends the Fourier-cosine (COS) method to the pricing and hedging of variable annuities embedded with guaranteed minimum withdrawal benefit (GMWB) riders. The COS method facilitates efficient computation of prices and hedge ratios of the GMWB riders when the underlying fund dynamics evolve under the influence of the general class of Lévy processes. Formulae are derived to value the contract at each withdrawal date using a backward recursive dynamic programming algorithm. Numerical comparisons are performed with results presented in Bacinello et al. [Scand. Actuar. J., 2014, 1–20], and Luo and Shevchenko [Int. J. Financ. Eng., 2014, 2, 1–24], to confirm the accuracy of the method. The efficiency of the proposed method is assessed by making comparisons with the approach presented in Bacinello et al. [op. cit.]. We find that the COS method presents highly accurate results with notably fast computational times. The valuation framework forms the basis for GMWB hedging. A local risk minimisation approach to hedging intra-withdrawal date risks is developed. A variety of risk measures are considered for minimisation in the general Lévy framework. While the second moment and variance have been considered in existing literature, we show that the Value-at-Risk (VaR) may also be of interest as a risk measure to minimise risk in variable annuities portfolios.     

Tail mutual exclusivity and Tail-VaR lower bounds

In this paper, we extend the concept of mutual exclusivity proposed by [Dhaene, J. & Denuit, M. (1999). The safest dependence structure among risks. Insurance: Mathematics and Economics 25, 11–21] to its tail counterpart and baptize this new dependency structure as tail mutual exclusivity. Probability levels are first specified for each component of the random vector. Under this dependency structure, at most one exceedance over the corresponding Value-at-Risks (VaRs) is possible, the other components being zero in such a case. No condition is imposed when all components stay below the VaRs. Several properties of this new negative dependence concept are derived. We show that this dependence structure gives rise to the smallest value of Tail-VaR (TVaR) of a sum of risks within a given Fréchet space, provided that the probability level of the TVaR is close enough to one.     

Static hedging and pricing American options

This paper utilizes the static hedge portfolio (SHP) approach of Derman et al. [Derman, E., Ergener, D., Kani, I., 1995. Static options replication. Journal of Derivatives 2, 78–95] and Carr et al. [Carr, P., Ellis, K., Gupta, V., 1998. Static hedging of exotic options. Journal of Finance 53, 1165–1190] to price and hedge American options under the Black-Scholes (1973) model and the constant elasticity of variance (CEV) model of Cox [Cox, J., 1975. Notes on option pricing I: Constant elasticity of variance diffusion. Working Paper, Stanford University]. The static hedge portfolio of an American option is formulated by applying the value-matching and smooth-pasting conditions on the early exercise boundary. The results indicate that the numerical efficiency of our static hedge portfolio approach is comparable to some recent advanced numerical methods such as Broadie and Detemple [Broadie, M., Detemple, J., 1996. American option valuation: New bounds, approximations, and a comparison of existing methods. Review of Financial Studies 9, 1211–1250] binomial Black-Scholes method with Richardson extrapolation (BBSR). The accuracy of the SHP method for the calculation of deltas and gammas is especially notable. Moreover, when the stock price changes, the recalculation of the prices and hedge ratios of the American options under the SHP method is quick because there is no need to solve the static hedge portfolio again. Finally, our static hedging approach also provides an intuitive derivation of the early exercise boundary near expiration.     

Singular Fourier–Padé series expansion of European option prices

We apply a new numerical method, the singular Fourier–Padé (SFP) method invented by Driscoll and Fornberg [Numer. Algorithms, 2001, 26, 77–92; The Gibbs Phenomenon in Various Representations and Applications, 2011], to price European-type options in Lévy and affine processes. The motivation behind this application is to reduce the inefficiency of current Fourier techniques when they are used to approximate piecewise continuous (non-smooth) probability density functions. When techniques such as fast Fourier transforms and Fourier series are applied to price and hedge options with non-smooth probability density functions, they cause the Gibbs phenomenon; accordingly, the techniques converge slowly for density functions with jumps in value or derivatives. This seriously adversely affects the efficiency and accuracy of these techniques. In this paper, we derive pricing formulae and their option Greeks using the SFP method to resolve the Gibbs phenomenon and restore the global spectral convergence rate. Moreover, we show that our method requires a small number of terms to yield fast error convergence, and it is able to accurately price any European-type option deep in/out of the money and with very long/short maturities. Furthermore, we conduct an error-bound analysis of the SFP method in option pricing. This new method performs favourably in numerical experiments compared with existing techniques.     

Skewness premium with Lévy processes

We study the skewness premium (SK) introduced by Bates [J. Finance, 1991, 46(3), 1009–1044] in a general context using Lévy processes. Under a symmetry condition, Fajardo and Mordecki [Quant. Finance, 2006, 6(3), 219–227] obtained that SK is given by Bates' x% rule. In this paper, we study SK in the absence of that symmetry condition. More exactly, we derive sufficient conditions for the excess of SK to be positive or negative, in terms of the characteristic triplet of the Lévy process under a risk-neutral measure.     

A novel Monte Carlo approach to hybrid local volatility models

We present in a Monte Carlo simulation framework, a novel approach for the evaluation of hybrid local volatility [Risk, 1994, 7, 18–20], [Int. J. Theor. Appl. Finance, 1998, 1, 61–110] models. In particular, we consider the stochastic local volatility model—see e.g. Lipton et al. [Quant. Finance, 2014, 14, 1899–1922], Piterbarg [Risk, 2007, April, 84–89], Tataru and Fisher [Quantitative Development Group, Bloomberg Version 1, 2010], Lipton [Risk, 2002, 15, 61–66]—and the local volatility model incorporating stochastic interest rates—see e.g. Atlan [ArXiV preprint math/0604316, 2006], Piterbarg [Risk, 2006, 19, 66–71], Deelstra and Rayée [Appl. Math. Finance, 2012, 1–23], Ren et al. [Risk, 2007, 20, 138–143]. For both model classes a particular (conditional) expectation needs to be evaluated which cannot be extracted from the market and is expensive to compute. We establish accurate and ‘cheap to evaluate’ approximations for the expectations by means of the stochastic collocation method [SIAM J. Numer. Anal., 2007, 45, 1005–1034], [SIAM J. Sci. Comput., 2005, 27, 1118–1139], [Math. Models Methods Appl. Sci., 2012, 22, 1–33], [SIAM J. Numer. Anal., 2008, 46, 2309–2345], [J. Biomech. Eng., 2011, 133, 031001], which was recently applied in the financial context [Available at SSRN 2529691, 2014], [J. Comput. Finance, 2016, 20, 1–19], combined with standard regression techniques. Monte Carlo pricing experiments confirm that our method is highly accurate and fast.     

A Note on Gaussian Estimation of the CKLS and CIR Models with Feedback Effects for Japan

In this note we extend the Gaussian estimation of two factor CKLS and CIR models recently considered in Nowman, K. B. (2001, Gaussian estimation and forecasting of multi-factor term structure models with an application to Japan and the United Kingdom, Asia Pacif. Financ. Markets 8, 23–34) to include feedback effects in the conditional mean as was originally formulated in general continuous time models by Bergstrom, A. R. (1966, Non-recursive models as discrete approximations to systems of stochastic differential equations, Econometrica 34, 173–182) with constant volatility. We use the exact discrete model of Bergstrom, A. R. (1966, Non-recursive models as discrete approximations to systems of stochastic differential equations, Econometrica 34, 173–182) to estimate the parameters which was first used by Brennan, M. J. and Schwartz, E. S. (1979, A continuous time approach to the pricing of bonds, J. Bank. Financ. 3, 133–155) to estimate their two factor interest model but incorporating the assumption of Nowman, K. B. (1997, Gaussian estimation of single-factor continuous time models of the term structure of interest rates, J. Financ. 52, 1695–1706; 2001, Gaussian estimation and forecasting of multi-factor term structure models with an application to Japan and the United Kingdom, Asia Pacif. Financ. Markets 8, 23–34). An application to monthly Japanese Euro currency rates indicates some evidence of feedback from the 1-year rate to the 1-month rate in both the CKLS and CIR models. We also find a low level-volatility effect supporting Nowman, K. B. (2001, Gaussian estimation and forecasting of multi-factor term structure models with an application to Japan and the United Kingdom, Asia Pacif. Financ. Markets 8, 23–34).     

Learning minimum variance discrete hedging directly from the market

Option hedging is a critical risk management problem in finance. In the Black–Scholes model, it has been recognized that computing a hedging position from the sensitivity of the calibrated model option value function is inadequate in minimizing variance of the option hedge risk, as it fails to capture the model parameter dependence on the underlying price (see e.g. Coleman et al., J. Risk, 2001, 5(6), 63–89; Hull and White, J. Bank. Finance, 2017, 82, 180–190). In this paper, we demonstrate that this issue can exist generally when determining hedging position from the sensitivity of the option function, either calibrated from a parametric model from current option prices or estimated nonparametricaly from historical option prices. Consequently, the sensitivity of the estimated model option function typically does not minimize variance of the hedge risk, even instantaneously. We propose a data-driven approach to directly learn a hedging function from the market data by minimizing variance of the local hedge risk. Using the S&P 500 index daily option data for more than a decade ending in August 2015, we show that the proposed method outperforms the parametric minimum variance hedging method proposed in Hull and White [J. Bank. Finance, 2017, 82, 180–190], as well as minimum variance hedging corrective techniques based on stochastic volatility or local volatility models. Furthermore, we show that the proposed approach achieves significant gain over the implied BS delta hedging for weekly and monthly hedging.     

Real-time market microstructure analysis: online transaction cost analysis

Motivated by the practical challenge in monitoring the performance of a large number of algorithmic trading orders, this paper provides a methodology that leads to automatic discovery of causes that lie behind poor trading performance. It also gives theoretical foundations to a generic framework for real-time trading analysis. The common acronym for investigating the causes of bad and good performance of trading is transaction cost analysis Rosenthal [Performance Metrics for Algorithmic Traders, 2009]). Automated algorithms take care of most of the traded flows on electronic markets (more than 70% in the US, 45% in Europe and 35% in Japan in 2012). Academic literature provides different ways to formalize these algorithms and show how optimal they can be from a mean-variance (like in Almgren and Chriss [J. Risk, 2000, 3(2), 5–39]), a stochastic control (e.g. Guéant et al. [Math. Financ. Econ., 2013, 7(4), 477–507]), an impulse control (see Bouchard et al. [SIAM J. Financ. Math., 2011, 2(1), 404–438]) or a statistical learning (as used in Laruelle et al. [Math. Financ. Econ., 2013, 7(3), 359–403]) viewpoint. This paper is agnostic about the way the algorithm has been built and provides a theoretical formalism to identify in real-time the market conditions that influenced its efficiency or inefficiency. For a given set of characteristics describing the market context, selected by a practitioner, we first show how a set of additional derived explanatory factors, called anomaly detectors, can be created for each market order (following for instance Cristianini and Shawe-Taylor [An Introduction to Support Vector Machines and Other Kernel-based Learning Methods, 2000]). We then will present an online methodology to quantify how this extended set of factors, at any given time, predicts (i.e. have influence, in the sense of predictive power or information defined in Basseville and Nikiforov [Detection of Abrupt Changes: Theory and Application, 1993], Shannon [Bell Syst. Tech. J., 1948, 27, 379–423] and Alkoot and Kittler [Pattern Recogn. Lett., 1999, 20(11), 1361–1369]) which of the orders are underperforming while calculating the predictive power of this explanatory factor set. Armed with this information, which we call influence analysis, we intend to empower the order monitoring user to take appropriate action on any affected orders by re-calibrating the trading algorithms working the order through new parameters, pausing their execution or taking over more direct trading control. Also we intend that use of this method can be taken advantage of to automatically adjust their trading action in the post trade analysis of algorithms.     

Sensitivity analysis of long-term cash flows

This paper conducts a sensitivity analysis of long-term cash flows. The price of the cash flow at time zero is given by the pricing operator of a Markov diffusion acting on the cash flow function. We study the extent to which the price of the cash flow is affected by small perturbations of the underlying Markov diffusion. The main tool is the Hansen–Scheinkman decomposition, which is a method to express the cash flow in terms of eigenvalues and eigenfunctions of the pricing operator. By incorporating techniques developed by Fournié et al. (Finance Stoch. 3:391–412, 1999), the sensitivities of long-term cash flows can be represented via simple expressions in terms of eigenvalues and eigenfunctions.     

Weighing asset pricing factors: a least squares model averaging approach

Empirical evidence has demonstrated that certain factors in asset pricing models are more important than others for explaining specific portfolio returns. We propose a technique that evaluates the factors included in popular linear asset pricing models. Our method has the advantage of simultaneously ranking the relative importance of those pricing factors through comparing their model weights. As an empirical verification, we apply our method to portfolios formed following Fama and French [A five-factor asset pricing model. J. Financ. Econ., 2015, 116, 1–22] and demonstrate that models accommodated to our factor rankings do improve their explanatory power in both in-sample and out-of-sample analyses.     

High-frequency trading in a limit order book

Abstract

We examine two performance measures advocated for asymmetric return distributions: the Sortino ratio—originally introduced by Sortino and Price (Sortino F and Price L 1994 J. Investing 59–65)—and a measure based on power utility introduced in Leland (Leland H 1999 Financial Analysts J. 27–36). In particular, we investigate the role of the maximum principle in this context, and assess the conditions under which the measures satisfy it. Our results add further motivation for the use of a modified Sortino ratio, by placing it on a sound theoretical foundation. In this light, we discuss its relative merits compared with alternative approaches.     

Marginal consistent dependence modelling using weak subordination for Brownian motions

We present an approach for modelling dependencies in exponential Lévy market models with arbitrary margins originated from time changed Brownian motions. Using weak subordination of Buchmann et al. [Bernoulli, 2017], we face a new layer of dependencies, superior to traditional approaches based on pathwise subordination, since weakly subordinated processes are not required to have independent components considering multivariate stochastic time changes. We apply a subordinator being able to incorporate any joint or idiosyncratic information arrivals. We emphasize multivariate variance gamma and normal inverse Gaussian processes and state explicit formulae for the Lévy characteristics. Using maximum likelihood, we estimate multivariate variance gamma models on various market data and show that these models are highly preferable to traditional approaches. Consistent values of basket-options under given marginal pricing models are achieved using the Esscher transform, generating a non-flat implied correlation surface.     

Optimal portfolio positioning within generalized Johnson distributions

Many empirical studies have shown that financial asset returns do not always exhibit Gaussian distributions, for example hedge fund returns. The introduction of the family of Johnson distributions allows a better fit to empirical financial data. Additionally, this class can be extended to a quite general family of distributions by considering all possible regular transformations of the standard Gaussian distribution. In this framework, we consider the portfolio optimal positioning problem, which has been first addressed by Brennan and Solanki [J. Financial Quant. Anal., 1981, 16, 279–300], Leland [J. Finance, 1980, 35, 581–594] and further developed by Carr and Madan [Quant. Finance, 2001, 1, 9–37] and Prigent [Generalized option based portfolio insurance. Working Paper, THEMA, University of Cergy-Pontoise, 2006]. As a by-product, we introduce the notion of Johnson stochastic processes. We determine and analyse the optimal portfolio for log return having Johnson distributions. The solution is characterized for arbitrary utility functions and illustrated in particular for a CRRA utility. Our findings show how the profiles of financial structured products must be selected when taking account of non Gaussian log-returns.     

The performance of enhanced-return index funds: evidence from bootstrap analysis

We apply the bootstrap technique proposed by Kosowski et al. [J. Finance, 2006, 61, 2551–2595] in conjunction with Carhart's [J. Finance, 1997, 52, 57–82] unconditional and Ferson and Schadt's [J. Finance, 1996, 51, 425–461] conditional four-factor models of performance to examine whether the performances of enhanced-return index funds over the 1996 to 2007 period are based on luck or superior ‘enhancing’ skills. The advantages of using the bootstrap to rank fund performance are many. It eliminates the need to specify the exact shape of the distribution from which returns are drawn and does not require estimating correlations between portfolio returns. It also eliminates the need to explicitly control for potential ‘data snooping’ biases that arise from an ex-post sort. Our results show evidence of enhanced-return index funds with positive and significant alphas after controlling for luck and sampling variability. The results are robust to both stock-only and derivative-enhanced index funds, although the spread of cross-sectional alphas for derivative-enhanced funds is slightly more pronounced. The study also examines various sub-periods within the sample horizon.     

Le calcul de la prime en réassurance d'excédent de sinistres

Abstract

Deux des rapports traitant la réassurance en excédent de sinistres présentés au XI^ème Congrès International d'Actuaires à Paris, par M. I var Hesselberg et par M. P. Dubois, ont abouti à des conclusions si divergentes au sujet de l'applicabilité de la méthode qu'il peut y avoir intérêt à les examiner de plus près.     

Market Value of Liabilities Mortality Risk

Abstract

Market values of the invested assets are frequently published. For most insurance liabilities, there are no published market values and, therefore, these have to be constructed. This construction can be based on a best estimate and a price for the risks in the liabilities. This paper presents a model explaining how the best estimate and the price of mortality risk can be constructed. Several methods to describe the risks are already known. The purpose of this paper is to describe a method to determine the mortality risk in a practical way.     

Fraud Detection Using a Multinomial Logit Model With Missing Information

Recently, Artís, Ayuso, and Guillén (2002, Journal of Risk and Insurance 69: 325–340; henceforth AAG) estimate a logit model using claims data. Some of the claims are categorized as “honest” and other claims are known to be fraudulent. Using the approach of Hausman, Abrevaya, and Scott‐Morton (1998 Journal of Econometrics 87: 239‐269), AAG estimate a modified logit model allowing for the possibility that some claims classified as honest might actually be fraudulent. Applying this model to data on Spanish automobile insurance claims, AGG find that 5 percent of the fraudulent claims go undetected. The purpose of this article is to estimate the model of AAG using a logit model with missing information. A constrained version of this model is used to reexamine the Spanish insurance claim data. The results indicate how to identify misclassified claims. We also show how misclassified claims can be identified using the AAG approach. We show that both approaches can be used to probabilistically identify misclassified claims.