# Computational Finance (q-fin.CP)

• The QLBS model is a discrete-time option hedging and pricing model that is based on Dynamic Programming (DP) and Reinforcement Learning (RL). It combines the famous Q-Learning method for RL with the Black-Scholes (-Merton) model's idea of reducing the problem of option pricing and hedging to the problem of optimal rebalancing of a dynamic replicating portfolio for the option, which is made of a stock and cash. Here we expand on several NuQLear (Numerical Q-Learning) topics with the QLBS model. First, we investigate the performance of Fitted Q Iteration for a RL (data-driven) solution to the model, and benchmark it versus a DP (model-based) solution, as well as versus the BSM model. Second, we develop an Inverse Reinforcement Learning (IRL) setting for the model, where we only observe prices and actions (re-hedges) taken by a trader, but not rewards. Third, we outline how the QLBS model can be used for pricing portfolios of options, rather than a single option in isolation, thus providing its own, data-driven and model independent solution to the (in)famous volatility smile problem of the Black-Scholes model.
• The present paper aims to demonstrate the usage of Convolutional Neural Networks as a generative model for stochastic processes, enabling researchers from a wide range of fields (such as quantitative finance and physics) to develop a general tool for forecasts and simulations without the need to identify/assume a specific system structure or estimate its parameters.
• Deep learning is an effective approach to solving image recognition problems. People draw intuitive conclusions from trading charts; this study uses the characteristics of deep learning to train computers in imitating this kind of intuition in the context of trading charts. The three steps involved are as follows: 1. Before training, we pre-process the input data from quantitative data to images. 2. We use a convolutional neural network (CNN), a type of deep learning, to train our trading model. 3. We evaluate the model's performance in terms of the accuracy of classification. A trading model is obtained with this approach to help devise trading strategies. The main application is designed to help clients automatically obtain personalized trading strategies.
• As machine learning ascends the peak of computer science zeitgeist, the usage and experimentation with sentiment analysis using various forms of textual data seems pervasive. The effect is especially pronounced in formulating securities trading strategies, due to a plethora of reasons including the relative ease of implementation and the abundance of academic research suggesting automated sentiment analysis can be productively used in trading strategies. The source data for such analyzers ranges a broad spectrum like social media feeds, micro-blogs, real-time news feeds, ex-post financial data etc. The abstract technique underlying these analyzers involve supervised learning of sentiment classification where the classifier is trained on annotated source corpus, and accuracy is measured by testing how well the classifiers generalizes on unseen test data from the corpus. Post training, and validation of fitted models, the classifiers are used to execute trading strategies, and the corresponding returns are compared with appropriate benchmark returns (for e.g., the S&P500 returns). In this paper, we introduce $\underline{a\ novel\ technique\ of\ using\ price\ volatilities\ to\ empirically\ determine\ the\ sentiment\ in\ news\ data}$, instead of the traditional reverse approach. We also perform meta sentiment analysis by evaluating the efficacy of existing sentiment classifiers and the precise definition of sentiment from securities trading context. We scrutinize the efficacy of using human-annotated sentiment classification and the tacit assumptions that introduces subjective bias in existing financial news sentiment classifiers.
• Time-series calibrations often suggest that the GARCH diffusion model could also be a suitable candidate for option (risk-neutral) calibration. But unlike the popular Heston model, it lacks a fast, semi-analytic solution for the pricing of vanilla options, perhaps the main reason why it is not used in this way. In this paper we show how an efficient finite difference-based PDE solver can effectively replace analytical solutions, enabling accurate option calibrations in less than a minute. The proposed pricing engine is shown to be robust under a wide range of model parameters and combines smoothly with black-box optimizers. We use this approach to produce a first PDE calibration of the GARCH diffusion model to SPX options and present some benchmark results for future reference.
• We perform a large-scale simulation of an Ising-based financial market model that includes 300 asset time series. The financial system simulated by the model shows a fat-tailed return distribution and volatility clustering and exhibits unstable periods indicated by the volatility index measured as the average of absolute-returns. Moreover, we determine that the cumulative risk fraction, which measures the system risk, changes at high volatility periods. We also calculate the inverse participation ratio (IPR) and its higher-power version, IPR6, from the absolute-return cross-correlation matrix. Finally, we show that the IPR and IPR6 also change at high volatility periods.
• The authors aim to develop numerical schemes of the two representative quadratic hedging strategies: locally risk minimizing and mean-variance hedging strategies, for models whose asset price process is given by the exponential of a normal inverse Gaussian process, using the results of Arai et al. \citeAIS, and Arai and Imai. Here normal inverse Gaussian process is a framework of Lévy processes frequently appeared in financial literature. In addition, some numerical results are also introduced.
• In this paper, I discuss a method to tackle the issues arising from the small data-sets available to data-scientists when building price predictive algorithms that use monthly/quarterly macro-financial indicators. I approach this by training separate classifiers on the equivalent dataset from a range of countries. Using these classifiers, a three level meta learning algorithm (MLA) is developed. I develop a transform, ASG, to create a country agnostic proxy for the macro-financial indicators. Using these proposed methods, I investigate the degree to which a predictive algorithm for the US 5Y bond price, predominantly using macro-financial indicators, can outperform an identical algorithm which only uses statistics deriving from previous price.
• We introduce the simulation tool SABCEMM (Simulate Agent-Based Computational Economic Market Models) for agent-based computational economic market (ABCEM) models. Our simulation tool is implemented in C++ and we can easily run ABCEM models with up to several million agents. Thanks to the object-oriented software design, this tool enables the user to design and compare multiple ABCEM models from a unified perspective. Thus, one can easily change the market mechanism or agent types. This makes it possible to quantitatively compare ABCEM models e.g. regarding the ability of each model to reproduce stylized facts. We present a qualitative study of three known ABCEM models and several variants of those. Furthermore, we discuss finite-size effects and time discretizations of ABCEM models. Finally, we show the great impact of different random number generators on the run time of ABCEM models and even on the qualitative output of the model. The code can be downloaded from GitHub https://github.com/SABCEMM/SABCEMM, such that all results can be reproduced by the reader
• This paper considers the problem of computing Bayesian estimates of system parameters and functions of them on the basis of observed system performance data. This is a previously studied issue where stochastic simulation approaches have been examined using the popular Metropolis--Hastings (MH) algorithm. This prior study has identified a recognised difficulty of tuning the proposal distribution so that the MH method provides realisations with sufficient mixing to deliver efficient convergence. This paper proposes and empirically examines a method of tuning the proposal using ideas borrowed from the numerical optimisation literature around efficient computation of Hessians so that gradient and curvature information of the target posterior can be incorporated in the proposal.
• We derive a semi-analytic formula for the transition probability of three-dimensional Brownian motion in the positive octant with absorption at the boundaries. Separation of variables in spherical coordinates leads to an eigenvalue problem for the resulting boundary value problem in the two angular components. The main theoretical result is a solution to the original problem expressed as an expansion into special functions and an eigenvalue which has to be chosen to allow a matching of the boundary condition. We discuss and test several computational methods to solve a finite-dimensional approximation to this nonlinear eigenvalue problem. Finally, we apply our results to the computation of default probabilities and credit valuation adjustments in a structural credit model with mutual liabilities.
• In this paper, we analyse piecewise deterministic Markov processes, as introduced in Davis (1984). Many models in insurance mathematics can be formulated in terms of the general concept of piecewise deterministic Markov processes. In this context, one is interested in computing certain quantities of interest such as the probability of ruin of an insurance company, or the insurance company's value, defined as the expected discounted future dividend payments until the time of ruin. Instead of explicitly solving the integro-(partial) differential equation related to the quantity of interest considered (an approach which can only be used in few special cases), we adapt the problem in a manner that allows us to apply deterministic numerical integration algorithms such as quasi-Monte Carlo rules; this is in contrast to applying random integration algorithms such as Monte Carlo. To this end, we reformulate a general cost functional as a fixed point of a particular integral operator, which allows for iterative approximation of the functional. Furthermore, we introduce a smoothing technique which is applied to the integrands involved, in order to use error bounds for deterministic cubature rules. On the analytical side, we prove a convergence result for our PDMP approximation, which is of independent interest as it justifies phase-type approximations on the process level. We illustrate the smoothing technique for a risk-theoretic example, and provide a comparative study of deterministic and Monte Carlo integration.
• In financial markets, low prices are generally associated with high volatilities and vice-versa, this well known stylized fact usually being referred to as leverage effect. We propose a local volatility model, given by a stochastic differential equation with piecewise constant coefficients, which accounts of leverage and mean-reversion effects in the dynamics of the prices. This model exhibits a regime switch in the dynamics accordingly to a certain threshold. It can be seen as a continuous time version of the Self-Exciting Threshold Autoregressive (SETAR) model. We propose an estimation procedure for the volatility and drift coefficients as well as for the threshold level. Tests are performed on the daily prices of 21 assets. They show empirical evidence for leverage and mean-reversion effects, consistent with the results in the literature.
• This paper develops a novel analytically tractable Neumann series of Bessel functions representation for pricing (and hedging) European-style double barrier knock-out options, which can be applied to the whole class of one-dimensional time-homogeneous diffusions even for the cases where the corresponding transition density is not known. The proposed numerical method is shown to be efficient and simple to implement. To illustrate the flexibility and computational power of the algorithm, we develop an extended jump to default model that is able to capture several empirical regularities commonly observed in the literature.
• When the underlying asset displays oscillations, spikes or heavy-tailed distributions, the lognormal diffusion process (for which Black and Scholes developed their momentous option pricing formula) is inadequate: in order to overcome these real world difficulties many models have been developed. Merton proposed a jump-diffusion model, where the dynamics of the price of the underlying are subject to variations due to a Brownian process and also to possible jumps, driven by a compound Poisson process. Merton's model admits a series solution for the European option price, and there have been a lot of attempts to obtain a discretisation of the Merton model with tree methods in order to price American or more complex options, e. g. Amin, the $O(n^3)$ procedure by Hilliard and Schwartz and the $O(n^{2.5})$ procedure by Dai et al. Here, starting from the implementation of the seven-nodes procedure by Hilliard and Schwartz, we prove theoretically that it is possible to reduce the complexity to $O(n \ln n)$ in the European case and $O(n^2 \ln n)$ in the American put case. These theoretical results can be obtained through suitable truncation of the lattice structure and the proofs provide closed formulas for the truncation limitations.
• We extend the viscosity solution characterization proved in [5] for call/put American option prices to the case of a general payoff function in a multi-dimensional setting: the price satisfies a semilinear re-action/diffusion type equation. Based on this, we propose two new numerical schemes inspired by the branching processes based algorithm of [8]. Our numerical experiments show that approximating the discontinu-ous driver of the associated reaction/diffusion PDE by local polynomials is not efficient, while a simple randomization procedure provides very good results.