Search Results
The invariant distribution of wealth and employment status in a small open economy with precautionary savings
2019, Bayer, Christian, Rendall, Alan D., Wälde, Klaus
We study optimal savings in continuous time with exogenous transitions between employment and unemployment as the only source of uncertainty in a small open economy. We prove the existence of an optimal consumption path. We exploit that the dynamics of consumption and wealth between jumps can be expressed as a Fuchsian system. We derive conditions under which an invariant joint distribution for the state variables, i.e., wealth and labour market status, exists and is unique. We also provide conditions under which the distribution of these variables converges to the invariant distribution. Our analysis relies on the notion of T-processes and applies results on the stability of Markovian processes from Meyn and Tweedie (1993a, b,c). © 2019 The Author(s)
From rough path estimates to multilevel Monte Carlo
2013, Bayer, Christian, Friz, Peter K., Riedel, Sebastian, Schoenmakers, John G.M.
Discrete approximations to solutions of stochastic differential equations are well-known to converge with strong rate 1=2. Such rates have played a key-role in Giles multilevel Monte Carlo method [Giles, Oper. Res. 2008] which gives a substantial reduction of the computational effort necessary for the evaluation of diffusion functionals. In the present article similar results are established for large classes of rough differential equations driven by Gaussian processes (including fractional Brownian motion with H > 1=4 as special case).
Pricing options under rough volatility with backward SPDEs
2020, Bayer, Christian, Qiu, Jinniao, Yao, Yao
In this paper, we study the option pricing problems for rough volatility models. As the framework is non-Markovian, the value function for a European option is not deterministic; rather, it is random and satisfies a backward stochastic partial differential equation (BSPDE). The existence and uniqueness of weak solutions is proved for general nonlinear BSPDEs with unbounded random leading coefficients whose connections with certain forward-backward stochastic differential equations are derived as well. These BSPDEs are then used to approximate American option prices. A deep learning-based method is also investigated for the numerical approximations to such BSPDEs and associated non-Markovian pricing problems. Finally, the examples of rough Bergomi type are numerically computed for both European and American options.
Numerical smoothing with hierarchical adaptive sparse grids and quasi-Monte Carlo methods for efficient option pricing
2022, Bayer, Christian, Ben Hammouda, Chiheb, Tempone, Raúl F.
When approximating the expectation of a functional of a stochastic process, the efficiency and performance of deterministic quadrature methods, such as sparse grid quadrature and quasi-Monte Carlo (QMC) methods, may critically depend on the regularity of the integrand. To overcome this issue and reveal the available regularity, we consider cases in which analytic smoothing cannot be performed, and introduce a novel numerical smoothing approach by combining a root finding algorithm with one-dimensional integration with respect to a single well-selected variable. We prove that under appropriate conditions, the resulting function of the remaining variables is a highly smooth function, potentially affording the improved efficiency of adaptive sparse grid quadrature (ASGQ) and QMC methods, particularly when combined with hierarchical transformations (i.e., Brownian bridge and Richardson extrapolation on the weak error). This approach facilitates the effective treatment of high dimensionality. Our study is motivated by option pricing problems, and our focus is on dynamics where the discretization of the asset price is necessary. Based on our analysis and numerical experiments, we show the advantages of combining numerical smoothing with the ASGQ and QMC methods over ASGQ and QMC methods without smoothing and the Monte Carlo approach.
A regularity structure for rough volatility
2019, Bayer, Christian, Friz, Peter K., Gassiat, Paul, Martin, Jorg, Stemper, Benjamin
A new paradigm has emerged recently in financial modeling: rough (stochastic) volatility. First observed by Gatheral et al. in high-frequency data, subsequently derived within market microstructure models, rough volatility captures parsimoniously key-stylized facts of the entire implied volatility surface, including extreme skews (as observed earlier by Alòs et al.) that were thought to be outside the scope of stochastic volatility models. On the mathematical side, Markovianity and, partially, semimartingality are lost. In this paper, we show that Hairer's regularity structures, a major extension of rough path theory, which caused a revolution in the field of stochastic partial differential equations, also provide a new and powerful tool to analyze rough volatility models.
Reinforced optimal control
2020, Bayer, Christian, Belomestny, Denis, Hager, Paul, Pigato, Paolo, Schoenmakers, John G. M., Spokoiny, Vladimir
Least squares Monte Carlo methods are a popular numerical approximation method for solving stochastic control problems. Based on dynamic programming, their key feature is the approximation of the conditional expectation of future rewards by linear least squares regression. Hence, the choice of basis functions is crucial for the accuracy of the method. Earlier work by some of us [Belomestny, Schoenmakers, Spokoiny, Zharkynbay, Commun. Math. Sci., 18(1):109?121, 2020] proposes to reinforce the basis functions in the case of optimal stopping problems by already computed value functions for later times, thereby considerably improving the accuracy with limited additional computational cost. We extend the reinforced regression method to a general class of stochastic control problems, while considerably improving the method?s efficiency, as demonstrated by substantial numerical examples as well as theoretical analysis.
Solving linear parabolic rough partial differential equations
2018, Bayer, Christian, Belomestny, Denis, Redmann, Martin, Riedel, Sebastian, Schoenmakers, John
We study linear rough partial differential equations in the setting of [Friz and Hairer, Springer, 2014, Chapter 12]. More precisely, we consider a linear parabolic partial differential equation driven by a deterministic rough path W of Hölder regularity with 1=3 < 1=2. Based on a stochastic representation of the solution of the rough partial differential equation, we propose a regression Monte Carlo algorithm for spatio-temporal approximation of the solution. We provide a full convergence analysis of the proposed approximation method which essentially relies on the new bounds for the higher order derivatives of the solution in space. Finally, a comprehensive simulation study showing the applicability of the proposed algorithm is presented.
Stability of deep neural networks via discrete rough paths
2020, Bayer, Christian, Friz, Peter, Tapia, Nikolas
Using rough path techniques, we provide a priori estimates for the output of Deep Residual Neural Networks. In particular we derive stability bounds in terms of the total p-variation of trained weights for any p ≥ 1.
Hierarchical adaptive sparse grids and quasi Monte Carlo for option pricing under the rough Bergomi model
2019, Bayer, Christian, Hammouda, Chiheb Ben, Tempone, Raúl F.
The rough Bergomi (rBergomi) model, introduced recently in [4], is a promising rough volatility model in quantitative finance. It is a parsimonious model depending on only three parameters, and yet exhibits remarkable fit to empirical implied volatility surfaces. In the absence of analytical European option pricing methods for the model, and due to the non-Markovian nature of the fractional driver, the prevalent option is to use the Monte Carlo (MC) simulation for pricing. Despite recent advances in the MC method in this context, pricing under the rBergomi model is still a timeconsuming task. To overcome this issue, we design a novel, hierarchical approach, based on i) adaptive sparse grids quadrature (ASGQ), and ii) quasi Monte Carlo (QMC). Both techniques are coupled with Brownian bridge construction and Richardson extrapolation. By uncovering the available regularity, our hierarchical methods demonstrate substantial computational gains with respect to the standard MC method, when reaching a sufficiently small relative error tolerance in the price estimates across different parameter constellations, even for very small values of the Hurst parameter. Our work opens a new research direction in this field, i.e., to investigate the performance of methods other than Monte Carlo for pricing and calibrating under the rBergomi model.
Simulation of conditional diffusions via forward-reverse stochastic representations
2013, Bayer, Christian, Schoenmakers, John G.M.
In this paper we derive stochastic representations for the finite dimensional distributions of a multidimensional diffusion on a fixed time interval, conditioned on the terminal state. The conditioning can be with respect to a fixed point or more generally with respect to some subset. The representations rely on a reverse process connected with the given (forward) diffusion as introduced in Milstein et al. [Bernoulli 10(2):281312, 2004] in the context of a forward-reverse transition density estimator. The corresponding Monte Carlo estimators have essentially root-N accuracy, hence they do not suffer from the curse of dimensionality. We provide a detailed convergence analysis and give a numerical example involving the realized variance in a stochastic volatility asset model conditioned on a fixed terminal value of the asset.