Probabilistic Meshless Methods for Bayesian Inverse Problems. Jon Cockayne July 8, 2016
|
|
- Walter Phillips
- 5 years ago
- Views:
Transcription
1 Probabilistic Meshless Methods for Bayesian Inverse Problems Jon Cockayne July 8,
2 Co-Authors Chris Oates Tim Sullivan Mark Girolami 2
3 What is PN? Many problems in mathematics have no analytical solution, and must be solved numerically. In Probabilistic Numerics we phrase such problems as inference problems and construct a probabilistic description of this error. 3
4 What is PN? Many problems in mathematics have no analytical solution, and must be solved numerically. In Probabilistic Numerics we phrase such problems as inference problems and construct a probabilistic description of this error. This is not a new idea 1! Lots of recent development on Integration, Optimization, ODEs, PDEs... see 2 [Kadane, 1985, Diaconis, 1988, O Hagan, 1992, Skilling, 1991] 3
5 What does this mean for PDEs? A prototypical linear PDE. Given g, κ, b find u (κ(x) u(x)) = g(x) u(x) = b(x) in D on D For general D, κ(x) this cannot be solved analytically. 4
6 What does this mean for PDEs? A prototypical linear PDE. Given g, κ, b find u (κ(x) u(x)) = g(x) u(x) = b(x) in D on D For general D, κ(x) this cannot be solved analytically. The majority of PDE solvers produce an approximation like: û(x) = N w i φ i (x) i=1 We want to quantify the error from finite N probabilistically. 4
7 What does this mean for PDEs? Inverse Problem: Given partial information of g, b, u find κ (κ(x) u(x)) = g(x) u(x) = b(x) in D on D 5
8 What does this mean for PDEs? Inverse Problem: Given partial information of g, b, u find κ (κ(x) u(x)) = g(x) u(x) = b(x) in D on D Bayesian Inverse Problem: κ Π κ Data u(x i ) = y i κ y Π y κ We want to account for an inaccurate forward solver in the inverse problem. 5
9 Why do this? Using an inaccurate forward solver in an inverse problem can produce biased and overconfident posteriors. 6
10 Why do this? Using an inaccurate forward solver in an inverse problem can produce biased and overconfident posteriors Probabilistic Standard ma = 4 ma = 8 ma = 16 ma = 32 ma = θ θ Comparison of inverse problem posteriors produced using the Probabilistic Meshless Method (PMM) vs. symmetric collocation. 6
11 Forward Problem 7
12 Abstract Formulation Au(x) = g(x) in D Forward inference procedure: u Π u Data Au(x i ) = g(x i ) u g Π g u 8
13 Posterior for the forward problem Use a Gaussian Process prior u Π u = GP(0, k). Assuming linearity, the posterior Π g u is available in closed-form 2. 2 [Cockayne et al., 2016, Särkkä, 2011, Cialenco et al., 2012, Owhadi, 2014] 9
14 Posterior for the forward problem Use a Gaussian Process prior u Π u = GP(0, k). Assuming linearity, the posterior Π g u is available in closed-form 2. Π g u GP(m 1, Σ 1 ) m 1 (x) = ĀK(x, X ) [ AĀK(X, X ) ] 1 g Σ 1 (x, x ) = k(x, x ) ĀK(x, X ) [ AĀK(X, X ) ] 1 AK(X, x ) Ā the adjoint of A Observation: The mean function is the same as in symmetric collocation! 2 [Cockayne et al., 2016, Särkkä, 2011, Cialenco et al., 2012, Owhadi, 2014] 9
15 Theoretical Results Theorem (Forward Contraction) For a ball B ɛ (u 0 ) of radius ɛ centered on the true solution u 0 of the PDE, we have ( ) h 1 Π g 2β 2ρ d u [B ɛ (u 0 )] = O ɛ h the fill distance β the smoothness of the prior ρ < β d/2 the order of the PDE d the input dimension 10
16 Toy Example 2 u(x) = g(x) x (0, 1) u(x) = 0 x = 0, 1 To associate with the notation from before... Π u GP(0, k(x, y)) A = d 2 dx 2 Ā = d 2 dy 2 11
17 Forward problem: posterior samples g(x) = sin(2πx) Mean Truth Samples 0.02 u x 12
18 Forward problem: convergence µ u Tr(Σ1) ma ma (a) Mean error from truth (b) Trace of posterior covariance Figure 2: Convergence 13
19 Inverse Problem 14
20 Recap (κ(x) u(x)) = g(x) u(x) = b(x) in D on D Now we need to incorporate the forward posterior measure Π g u into the posterior measure for the inverse problem, κ 15
21 Incorporation of Forward Measure Assuming the data in the inverse problem is: y i = u(x i ) + ξ i ξ N(0, Γ) i = 1,..., n implies the standard likelihood: But we don t know u p(y κ, u) N(y; u, Γ) 16
22 Incorporation of Forward Measure Assuming the data in the inverse problem is: y i = u(x i ) + ξ i ξ N(0, Γ) i = 1,..., n implies the standard likelihood: p(y κ, u) N(y; u, Γ) But we don t know u Marginalise the forward posterior Π g u to obtain a PN likelihood: p PN (y κ) p(y κ, u)dπ g u N(y; m 1, Γ + Σ 1 ) 16
23 Inverse Contraction Denote by Π y κ the posterior for κ from likelihood p, and by Π y κ,pn the posterior for κ from likelihood p PN. Theorem (Inverse Contraction) Assume Π y κ δ(κ 0 ) as n. Then Π y κ,pn δ(κ 0) provided h = o(n 1/(β ρ d/2) ) 17
24 Back to the Toy Example (κ u(x)) = sin(2πx) x (0, 1) u(x) = 0 x = 0, 1 Infer κ R + ; data generated for κ = 1 at x = 0.25, Corrupted with independent Gaussian noise ξ N(0, ) 18
25 θ θ Posteriors for κ Probabilistic Standard ma = 4 ma = 8 ma = 16 ma = 32 ma = θ θ (a) Posterior Distributions for different numbers of design points Probabilistic Standard n n (b) Convergence of posterior distributions with number of design points. 19
26 Nonlinear Example: Steady-State Allen Cahn 20
27 Allen Cahn A prototypical nonlinear model. δ 2 u(x) + δ 1 (u(x) 3 u(x)) = 0 x (0, 1) 2 u(x) = 1 x 1 {0, 1} ; 0 < x 2 < 1 u(x) = 1 x 2 {0, 1} ; 0 < x 1 < 1 Goal: infer δ from 16 equally spaced observations of u(x) in the interior of the domain. 1.0 Negative Stable 1.0 Unstable 1.0 Positive Stable
28 Allen Cahn: Forward Solutions Nonlinear PDE - non-gp posterior sampling schemes required, see [Cockayne et al., 2016]. 1.0 Initial Design Optimal design
29 Allen Cahn: Inverse Problem PMM l = 5 PMM l = 10 PMM l = 20 PMM l = 40 PMM l = FEM 5x5 FEM 10x10 FEM 25x25 FEM 50x θ (a) PMM θ (b) FEA Comparison of posteriors for δ with different solver resolutions, when using the PMM forward solver with PN likelihood, vs. FEA forward solver with Gaussian likelihood. 23
30 Conclusions 24
31 We have shown... How to build probability measures for the forward solution of PDEs. How to use this to make rhobust inferences in PDE inverse problems, even with inaccurate forward solvers. 25
32 We have shown... How to build probability measures for the forward solution of PDEs. How to use this to make rhobust inferences in PDE inverse problems, even with inaccurate forward solvers. Coming soon... Evolutionary problems ( / t). More profound nonlinearities. Non-Gaussian priors. 25
33 I. Cialenco, G. E. Fasshauer, and Q. Ye. Approximation of stochastic partial differential equations by a kernel-based collocation method. Int. J. Comput. Math., 89(18): , J. Cockayne, C. J. Oates, T. Sullivan, and M. Girolami. Probabilistic Meshless Methods for Partial Differential Equations and Bayesian Inverse Problems. arxiv preprint arxiv: , P. Diaconis. Bayesian numerical analysis. Statistical decision theory and related topics IV, 1: , J. B. Kadane. Parallel and Sequential Computation: A Statistician s view. Journal of Complexity, 1: , A. O Hagan. Some Bayesian numerical analysis. Bayesian Statistics, 4: , H. Owhadi. Bayesian numerical homogenization. arxiv preprint arxiv: , S. Särkkä. Linear operators and stochastic partial differential equations in Gaussian process regression. In Artificial Neural Networks and Machine Learning ICANN 2011: 21 st International Conference on Artificial Neural Networks, Espoo, Finland, June 14-17, 2011, Proceedings, Part II, pages Springer, J. Skilling. Bayesian solution of ordinary differential equations. Maximum Entropy and Bayesian Methods, 50: 23 37,
Machine Learning for Quantitative Finance
Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing
More informationKernel Conditional Quantile Estimation via Reduction Revisited
Kernel Conditional Quantile Estimation via Reduction Revisited Novi Quadrianto Novi.Quad@gmail.com The Australian National University, Australia NICTA, Statistical Machine Learning Program, Australia Joint
More informationImplementing an Agent-Based General Equilibrium Model
Implementing an Agent-Based General Equilibrium Model 1 2 3 Pure Exchange General Equilibrium We shall take N dividend processes δ n (t) as exogenous with a distribution which is known to all agents There
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationPDE Project Course 1. Adaptive finite element methods
PDE Project Course 1. Adaptive finite element methods Anders Logg logg@math.chalmers.se Department of Computational Mathematics PDE Project Course 03/04 p. 1 Lecture plan Introduction to FEM FEM for Poisson
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationarxiv: v2 [q-fin.pr] 23 Nov 2017
VALUATION OF EQUITY WARRANTS FOR UNCERTAIN FINANCIAL MARKET FOAD SHOKROLLAHI arxiv:17118356v2 [q-finpr] 23 Nov 217 Department of Mathematics and Statistics, University of Vaasa, PO Box 7, FIN-6511 Vaasa,
More informationChapter 7: Estimation Sections
1 / 31 : Estimation Sections 7.1 Statistical Inference Bayesian Methods: 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods: 7.5 Maximum Likelihood
More informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More informationA THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES
Proceedings of ALGORITMY 01 pp. 95 104 A THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES BEÁTA STEHLÍKOVÁ AND ZUZANA ZÍKOVÁ Abstract. A convergence model of interest rates explains the evolution of the
More informationAccelerated Stochastic Gradient Descent Praneeth Netrapalli MSR India
Accelerated Stochastic Gradient Descent Praneeth Netrapalli MSR India Presented at OSL workshop, Les Houches, France. Joint work with Prateek Jain, Sham M. Kakade, Rahul Kidambi and Aaron Sidford Linear
More informationIntroduction to Sequential Monte Carlo Methods
Introduction to Sequential Monte Carlo Methods Arnaud Doucet NCSU, October 2008 Arnaud Doucet () Introduction to SMC NCSU, October 2008 1 / 36 Preliminary Remarks Sequential Monte Carlo (SMC) are a set
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationApplication of MCMC Algorithm in Interest Rate Modeling
Application of MCMC Algorithm in Interest Rate Modeling Xiaoxia Feng and Dejun Xie Abstract Interest rate modeling is a challenging but important problem in financial econometrics. This work is concerned
More informationPart 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL)
Part 3: Trust-region methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective
More informationLecture 2: The Neoclassical Growth Model
Lecture 2: The Neoclassical Growth Model Florian Scheuer 1 Plan Introduce production technology, storage multiple goods 2 The Neoclassical Model Three goods: Final output Capital Labor One household, with
More informationBuilding Infinite Processes from Regular Conditional Probability Distributions
Chapter 3 Building Infinite Processes from Regular Conditional Probability Distributions Section 3.1 introduces the notion of a probability kernel, which is a useful way of systematizing and extending
More informationPICOF, Palaiseau, April 2-4, 2012
The Sobolev gradient regularization strategy for optical tomography coupled with a finite element formulation of the radiative transfer equation Fabien Dubot, Olivier Balima, Yann Favennec, Daniel Rousse
More informationPhys. Lett. A, 372/17, (2008),
Phys. Lett. A, 372/17, (2008), 3064-3070. 1 Wave scattering by many small particles embedded in a medium. A. G. Ramm (Mathematics Department, Kansas State University, Manhattan, KS66506, USA and TU Darmstadt,
More informationThe data-driven COS method
The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica Reading group, March 13, 2017 Reading group, March
More informationThe data-driven COS method
The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica CMMSE 2017, July 6, 2017 Álvaro Leitao (CWI & TUDelft)
More informationCalibration Lecture 4: LSV and Model Uncertainty
Calibration Lecture 4: LSV and Model Uncertainty March 2017 Recap: Heston model Recall the Heston stochastic volatility model ds t = rs t dt + Y t S t dw 1 t, dy t = κ(θ Y t ) dt + ξ Y t dw 2 t, where
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationSimulation Wrap-up, Statistics COS 323
Simulation Wrap-up, Statistics COS 323 Today Simulation Re-cap Statistics Variance and confidence intervals for simulations Simulation wrap-up FYI: No class or office hours Thursday Simulation wrap-up
More informationSTART HERE: Instructions. 1 Exponential Family [Zhou, Manzil]
START HERE: Instructions Thanks a lot to John A.W.B. Constanzo and Shi Zong for providing and allowing to use the latex source files for quick preparation of the HW solution. The homework was due at 9:00am
More informationHigh Dimensional Bayesian Optimisation and Bandits via Additive Models
1/20 High Dimensional Bayesian Optimisation and Bandits via Additive Models Kirthevasan Kandasamy, Jeff Schneider, Barnabás Póczos ICML 15 July 8 2015 2/20 Bandits & Optimisation Maximum Likelihood inference
More informationResearch Statement. Dapeng Zhan
Research Statement Dapeng Zhan The Schramm-Loewner evolution (SLE), first introduced by Oded Schramm ([12]), is a oneparameter (κ (0, )) family of random non-self-crossing curves, which has received a
More informationPricing Barrier Options under Local Volatility
Abstract Pricing Barrier Options under Local Volatility Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 16 November 2002 We study pricing under the local volatility. Our research is mainly
More informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationMulti-period Portfolio Choice and Bayesian Dynamic Models
Multi-period Portfolio Choice and Bayesian Dynamic Models Petter Kolm and Gordon Ritter Courant Institute, NYU Paper appeared in Risk Magazine, Feb. 25 (2015) issue Working paper version: papers.ssrn.com/sol3/papers.cfm?abstract_id=2472768
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationSparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models
Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models Norbert Hilber Seminar of Applied Mathematics ETH Zürich Workshop on Financial Modeling with Jump Processes p. 1/18 Outline
More informationFinite Additivity in Dubins-Savage Gambling and Stochastic Games. Bill Sudderth University of Minnesota
Finite Additivity in Dubins-Savage Gambling and Stochastic Games Bill Sudderth University of Minnesota This talk is based on joint work with Lester Dubins, David Heath, Ashok Maitra, and Roger Purves.
More informationDerivation Of The Capital Asset Pricing Model Part I - A Single Source Of Uncertainty
Derivation Of The Capital Asset Pricing Model Part I - A Single Source Of Uncertainty Gary Schurman MB, CFA August, 2012 The Capital Asset Pricing Model CAPM is used to estimate the required rate of return
More informationFinite Difference Approximation of Hedging Quantities in the Heston model
Finite Difference Approximation of Hedging Quantities in the Heston model Karel in t Hout Department of Mathematics and Computer cience, University of Antwerp, Middelheimlaan, 22 Antwerp, Belgium Abstract.
More informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
More informationFINITE DIFFERENCE METHODS
FINITE DIFFERENCE METHODS School of Mathematics 2013 OUTLINE Review 1 REVIEW Last time Today s Lecture OUTLINE Review 1 REVIEW Last time Today s Lecture 2 DISCRETISING THE PROBLEM Finite-difference approximations
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationStability in geometric & functional inequalities
Stability in geometric & functional inequalities A. Figalli The University of Texas at Austin www.ma.utexas.edu/users/figalli/ Alessio Figalli (UT Austin) Stability in geom. & funct. ineq. Krakow, July
More informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationSensitivity Analysis on Long-term Cash flows
Sensitivity Analysis on Long-term Cash flows Hyungbin Park Worcester Polytechnic Institute 19 March 2016 Eastern Conference on Mathematical Finance Worcester Polytechnic Institute, Worceseter, MA 1 / 49
More informationImproving Stock Price Prediction with SVM by Simple Transformation: The Sample of Stock Exchange of Thailand (SET)
Thai Journal of Mathematics Volume 14 (2016) Number 3 : 553 563 http://thaijmath.in.cmu.ac.th ISSN 1686-0209 Improving Stock Price Prediction with SVM by Simple Transformation: The Sample of Stock Exchange
More informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More informationLecture 7: Bayesian approach to MAB - Gittins index
Advanced Topics in Machine Learning and Algorithmic Game Theory Lecture 7: Bayesian approach to MAB - Gittins index Lecturer: Yishay Mansour Scribe: Mariano Schain 7.1 Introduction In the Bayesian approach
More informationOn Risk Measures, Market Making, and Exponential Families
On Risk Measures, Market Making, and Exponential Families JACOB D. ABERNETHY University of Michigan and RAFAEL M. FRONGILLO Harvard University and SINDHU KUTTY University of Michigan In this note we elaborate
More informationQuarterly Storage Model of U.S. Cotton Market: Estimation of the Basis under Rational Expectations. Oleksiy Tokovenko 1 Lewell F.
Quarterly Storage Model of U.S. Cotton Market: Estimation of the Basis under Rational Expectations Oleksiy Tokovenko 1 Lewell F. Gunter Selected Paper prepared for presentation at the American Agricultural
More informationProject 1: Double Pendulum
Final Projects Introduction to Numerical Analysis II http://www.math.ucsb.edu/ atzberg/winter2009numericalanalysis/index.html Professor: Paul J. Atzberger Due: Friday, March 20th Turn in to TA s Mailbox:
More informationNonlinear Black-Scholes Equation. Through Radial Basis Functions
Journal of Applied Mathematics & Bioinformatics, vol.4, no.3, 2014, 75-86 ISSN: 1792-6602 (print), 1792-6939 (online) Scienpress Ltd, 2014 Nonlinear Black-Scholes Equation Through Radial Basis Functions
More informationSharing the Burden: Monetary and Fiscal Responses to a World Liquidity Trap David Cook and Michael B. Devereux
Sharing the Burden: Monetary and Fiscal Responses to a World Liquidity Trap David Cook and Michael B. Devereux Online Appendix: Non-cooperative Loss Function Section 7 of the text reports the results for
More informationRegret-based Selection
Regret-based Selection David Puelz (UT Austin) Carlos M. Carvalho (UT Austin) P. Richard Hahn (Chicago Booth) May 27, 2017 Two problems 1. Asset pricing: What are the fundamental dimensions (risk factors)
More informationAmerican Option Pricing: A Simulated Approach
Utah State University DigitalCommons@USU All Graduate Plan B and other Reports Graduate Studies 5-2013 American Option Pricing: A Simulated Approach Garrett G. Smith Utah State University Follow this and
More informationarxiv: v1 [math.st] 18 Sep 2018
Gram Charlier and Edgeworth expansion for sample variance arxiv:809.06668v [math.st] 8 Sep 08 Eric Benhamou,* A.I. SQUARE CONNECT, 35 Boulevard d Inkermann 900 Neuilly sur Seine, France and LAMSADE, Universit
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationPh.D. Preliminary Examination MICROECONOMIC THEORY Applied Economics Graduate Program June 2017
Ph.D. Preliminary Examination MICROECONOMIC THEORY Applied Economics Graduate Program June 2017 The time limit for this exam is four hours. The exam has four sections. Each section includes two questions.
More informationOn the pricing equations in local / stochastic volatility models
On the pricing equations in local / stochastic volatility models Hao Xing Fields Institute/Boston University joint work with Erhan Bayraktar, University of Michigan Kostas Kardaras, Boston University Probability
More informationZürich Spring School on Lévy Processes. Poster abstracts
Zürich Spring School on Lévy Processes Poster abstracts 31 March 2015 Akhlaque Ahmad Option Pricing Using Fourier Transforms: An Integrated Approach In this paper, we model stochastic volatility using
More informationECON 815. A Basic New Keynesian Model II
ECON 815 A Basic New Keynesian Model II Winter 2015 Queen s University ECON 815 1 Unemployment vs. Inflation 12 10 Unemployment 8 6 4 2 0 1 1.5 2 2.5 3 3.5 4 4.5 5 Core Inflation 14 12 10 Unemployment
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More informationFinancial Mathematics and Supercomputing
GPU acceleration in early-exercise option valuation Álvaro Leitao and Cornelis W. Oosterlee Financial Mathematics and Supercomputing A Coruña - September 26, 2018 Á. Leitao & Kees Oosterlee SGBM on GPU
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationAnalysis of the Bitcoin Exchange Using Particle MCMC Methods
Analysis of the Bitcoin Exchange Using Particle MCMC Methods by Michael Johnson M.Sc., University of British Columbia, 2013 B.Sc., University of Winnipeg, 2011 Project Submitted in Partial Fulfillment
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationModeling of Price. Ximing Wu Texas A&M University
Modeling of Price Ximing Wu Texas A&M University As revenue is given by price times yield, farmers income risk comes from risk in yield and output price. Their net profit also depends on input price, but
More informationSYLLABUS AND SAMPLE QUESTIONS FOR MSQE (Program Code: MQEK and MQED) Syllabus for PEA (Mathematics), 2013
SYLLABUS AND SAMPLE QUESTIONS FOR MSQE (Program Code: MQEK and MQED) 2013 Syllabus for PEA (Mathematics), 2013 Algebra: Binomial Theorem, AP, GP, HP, Exponential, Logarithmic Series, Sequence, Permutations
More informationAveraged bond prices for Fong-Vasicek and the generalized Vasicek interest rates models
MATHEMATICAL OPTIMIZATION Mathematical Methods In Economics And Industry 007 June 3 7, 007, Herl any, Slovak Republic Averaged bond prices for Fong-Vasicek and the generalized Vasicek interest rates models
More informationThe Game-Theoretic Framework for Probability
11th IPMU International Conference The Game-Theoretic Framework for Probability Glenn Shafer July 5, 2006 Part I. A new mathematical foundation for probability theory. Game theory replaces measure theory.
More informationIdentifying Long-Run Risks: A Bayesian Mixed-Frequency Approach
Identifying : A Bayesian Mixed-Frequency Approach Frank Schorfheide University of Pennsylvania CEPR and NBER Dongho Song University of Pennsylvania Amir Yaron University of Pennsylvania NBER February 12,
More informationEstimation Appendix to Dynamics of Fiscal Financing in the United States
Estimation Appendix to Dynamics of Fiscal Financing in the United States Eric M. Leeper, Michael Plante, and Nora Traum July 9, 9. Indiana University. This appendix includes tables and graphs of additional
More informationWKB Method for Swaption Smile
WKB Method for Swaption Smile Andrew Lesniewski BNP Paribas New York February 7 2002 Abstract We study a three-parameter stochastic volatility model originally proposed by P. Hagan for the forward swap
More informationCalibration of Interest Rates
WDS'12 Proceedings of Contributed Papers, Part I, 25 30, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Calibration of Interest Rates J. Černý Charles University, Faculty of Mathematics and Physics, Prague,
More informationA Stochastic Levenberg-Marquardt Method Using Random Models with Application to Data Assimilation
A Stochastic Levenberg-Marquardt Method Using Random Models with Application to Data Assimilation E Bergou Y Diouane V Kungurtsev C W Royer July 5, 08 Abstract Globally convergent variants of the Gauss-Newton
More information13.3 A Stochastic Production Planning Model
13.3. A Stochastic Production Planning Model 347 From (13.9), we can formally write (dx t ) = f (dt) + G (dz t ) + fgdz t dt, (13.3) dx t dt = f(dt) + Gdz t dt. (13.33) The exact meaning of these expressions
More informationThe Neoclassical Growth Model
The Neoclassical Growth Model 1 Setup Three goods: Final output Capital Labour One household, with preferences β t u (c t ) (Later we will introduce preferences with respect to labour/leisure) Endowment
More informationComprehensive Exam. August 19, 2013
Comprehensive Exam August 19, 2013 You have a total of 180 minutes to complete the exam. If a question seems ambiguous, state why, sharpen it up and answer the sharpened-up question. Good luck! 1 1 Menu
More informationWhat can we do with numerical optimization?
Optimization motivation and background Eddie Wadbro Introduction to PDE Constrained Optimization, 2016 February 15 16, 2016 Eddie Wadbro, Introduction to PDE Constrained Optimization, February 15 16, 2016
More informationTop-down particle filtering for Bayesian decision trees
Top-down particle filtering for Bayesian decision trees Balaji Lakshminarayanan 1, Daniel M. Roy 2 and Yee Whye Teh 3 1. Gatsby Unit, UCL, 2. University of Cambridge and 3. University of Oxford Outline
More informationHeterogeneous Firm, Financial Market Integration and International Risk Sharing
Heterogeneous Firm, Financial Market Integration and International Risk Sharing Ming-Jen Chang, Shikuan Chen and Yen-Chen Wu National DongHwa University Thursday 22 nd November 2018 Department of Economics,
More informationEstimating Macroeconomic Models of Financial Crises: An Endogenous Regime-Switching Approach
Estimating Macroeconomic Models of Financial Crises: An Endogenous Regime-Switching Approach Gianluca Benigno 1 Andrew Foerster 2 Christopher Otrok 3 Alessandro Rebucci 4 1 London School of Economics and
More informationApproximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications
Approximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications Anna Timonina University of Vienna, Abraham Wald PhD Program in Statistics and Operations
More informationNotes on Estimating the Closed Form of the Hybrid New Phillips Curve
Notes on Estimating the Closed Form of the Hybrid New Phillips Curve Jordi Galí, Mark Gertler and J. David López-Salido Preliminary draft, June 2001 Abstract Galí and Gertler (1999) developed a hybrid
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationChapter 8. Markowitz Portfolio Theory. 8.1 Expected Returns and Covariance
Chapter 8 Markowitz Portfolio Theory 8.1 Expected Returns and Covariance The main question in portfolio theory is the following: Given an initial capital V (0), and opportunities (buy or sell) in N securities
More informationA Schauder estimate for stochastic PDEs
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 2016 A Schauder estimate for stochastic PDEs Kai
More informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationMonte Carlo Simulations in the Teaching Process
Monte Carlo Simulations in the Teaching Process Blanka Šedivá Department of Mathematics, Faculty of Applied Sciences University of West Bohemia, Plzeň, Czech Republic CADGME 2018 Conference on Digital
More informationDecision Analysis. Carlos A. Santos Silva June 5 th, 2009
Decision Analysis Carlos A. Santos Silva June 5 th, 2009 What is decision analysis? Often, there is more than one possible solution: Decision depends on the criteria Decision often must be made in uncertain
More informationPublic Information and Effi cient Capital Investments: Implications for the Cost of Capital and Firm Values
Public Information and Effi cient Capital Investments: Implications for the Cost of Capital and Firm Values P O. C Department of Finance Copenhagen Business School, Denmark H F Department of Accounting
More informationSocially-Optimal Design of Crowdsourcing Platforms with Reputation Update Errors
Socially-Optimal Design of Crowdsourcing Platforms with Reputation Update Errors 1 Yuanzhang Xiao, Yu Zhang, and Mihaela van der Schaar Abstract Crowdsourcing systems (e.g. Yahoo! Answers and Amazon Mechanical
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationOptimal switching problems for daily power system balancing
Optimal switching problems for daily power system balancing Dávid Zoltán Szabó University of Manchester davidzoltan.szabo@postgrad.manchester.ac.uk June 13, 2016 ávid Zoltán Szabó (University of Manchester)
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose
More informationThe Correlation Smile Recovery
Fortis Bank Equity & Credit Derivatives Quantitative Research The Correlation Smile Recovery E. Vandenbrande, A. Vandendorpe, Y. Nesterov, P. Van Dooren draft version : March 2, 2009 1 Introduction Pricing
More informationRISK-NEUTRAL VALUATION AND STATE SPACE FRAMEWORK. JEL Codes: C51, C61, C63, and G13
RISK-NEUTRAL VALUATION AND STATE SPACE FRAMEWORK JEL Codes: C51, C61, C63, and G13 Dr. Ramaprasad Bhar School of Banking and Finance The University of New South Wales Sydney 2052, AUSTRALIA Fax. +61 2
More informationOptimal exercise price of American options near expiry
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 2009 Optimal exercise price of American options near expiry W.-T. Chen
More informationIntroduction to Affine Processes. Applications to Mathematical Finance
and Its Applications to Mathematical Finance Department of Mathematical Science, KAIST Workshop for Young Mathematicians in Korea, 2010 Outline Motivation 1 Motivation 2 Preliminary : Stochastic Calculus
More informationRisk minimizing strategies for tracking a stochastic target
Risk minimizing strategies for tracking a stochastic target Andrzej Palczewski Abstract We consider a stochastic control problem of beating a stochastic benchmark. The problem is considered in an incomplete
More informationStochastic Volatility (SV) Models
1 Motivations Stochastic Volatility (SV) Models Jun Yu Some stylised facts about financial asset return distributions: 1. Distribution is leptokurtic 2. Volatility clustering 3. Volatility responds to
More information