On Solving Integral Equations using. Markov Chain Monte Carlo Methods
|
|
- Wilfred Wells
- 5 years ago
- Views:
Transcription
1 On Solving Integral quations using Markov Chain Monte Carlo Methods Arnaud Doucet Department of Statistics and Department of Computer Science, University of British Columbia, Vancouver, BC, Canada mail: Vladislav B. Tadić Department of Mathematics, University of Bristol, Bristol, UK mail: Technical Report CUD-F-INFNG Cambridge University no. 444 First version December, Updated November 4 & January 7
2 ABSTRACT In this report, we propose an original approach to solve Fredholm equations of the second kind. We interpret the standard von Neumann expansion of the solution as an expectation with respect to a probability distribution de ned on an union of subspaces of variable dimension. Based on this representation, it is possible to use trans-dimensional Markov Chain Monte Carlo (MCMC) methods such as Reversible Jump MCMC to approximate the solution numerically. This can be an attractive alternative to standard Sequential Importance Sampling (SIS) methods routinely used in this context. We sketch an application to value function estimation for a Markov decision process. Keywords: Fredholm equation, Trans-dimensional Markov Chain Monte Carlo, Sequential Importance Sampling, Sequential Monte Carlo.
3 I. Fredholm equations and Von Neumann s expansion Fredholm equations of the second kind and related problems appear in many scienti c elds including optimal control [], molecular population genetics [6] and physics [9]. Formally, we are interested in solving the integral equation f (x ) = K (x ; x ) f (x ) dx + g (x ) () where g :! R and K :! R are known and f :! R is unknown. Let us de ne K (x; y), ; K (x; y), K (x; y) and K n (x; y), K (x; z) K n (z; y) dz: If X jk n (x ; x n ) g (x n )j dx n < () n= then the solution of the Fredholm equation () admits the following Von Neumann s series representation; see [], [9] for details: f (x ) = K (x ; x ) f (x ) dx + g (x ) = K (x ; x ) K (x ; x ) f (x ) dx + g (x ) dx + g (x ) = K (x ; x ) K (x ; x ) f (x ) dx dx + K (x ; x ) g (x ) dy + g (x ) and, by iterating X! ny f (x ) = g (x ) + K (x k ; x k ) g (x n ) dx :n (3) n= n k= where we use the notation x i:j, (x i ; : : : ; x j ) for i j: By introducing the notation f (x ) = g (x ) (4) and setting for n ; ny f n (x :n ) = g (x n ) K (x k ; x k ) (5) then it is possible to rewrite (3) as X f (x ) = f (x ) + f n (x :n ) dx :n : n (6) n= We will address two problems in this paper: how to estimate the function f (x ) over the set and how to estimate this function pointwise. k=
4 II. Monte Carlo Methods to Solve Fredholm quations Computing (3) is challenging as it involves an in nite sum of integrals of increasing dimension. Monte Carlo methods are a natural approach to this problem. A. Sequential Importance Sampling. Algorithm The use of Monte Carlo methods to solve problems of this type can be traced back 5 years. The standard approach consists of using Sequential Importance Sampling (SIS) to numerically approximate (3); see for example [], [9]. Consider a Markov chain with initial distribution/density (x) on and a transition kernel M (x; y) which gives the probability or probability density of moving to state y when the current state is x. We select and M such that (x) > over and M (x; y) > if K (x; y) 6=. Moreover, M is chosen to have an absorbing/cemetery state, say fyg =, such that M (x; fyg) = for any x. The algorithm proceeds as follows to approximate the function f: Simulate N independent Markov chain paths y). n o N + until absorption (i.e. X(i) i= k (i) + = Calculate the associated importance weights k(i) Q K W k= M = >: g k (i) ;X(i) k (i) k (i) ;X(i) k (i) A g k (i) if k (i) ; if k (i) = : (7) The empirical measure bf (x ) = N i= W x h R is an unbiased Monte Carlo approximation of the function f (i.e. for any set A, b i f (x A ) dx = R A f (x ) dx ). () If we are interested in estimating the function f (x ) at a given point say x = x, then we n o N simulate paths starting from + = x according to M until absorption/death and i= 3
5 the importance weights become W = >< >: k(i) K k= M k (i) ;X(i) k (i) k (i) ;X(i) k (i) A g k (i) if k (i) ; g(x) if k (i) = : (9) We obtain the following unbiased estimate of f (x) bf (x) = N i= W x; :k (i) : (). Importance Sampling on Path Space To check the unbiasedness of the estimates () and (), we use a slightly non-standard argument which will later prove useful. The rst method to estimate the function f through () can be interpreted as an importance sampling technique using an importance distribution (n; x :n ) de ned on the path space F, U k= fkg k+ where (n; x :n ) = p ;n ;n (x :n ) () with p ;n the probability that the simulated path is of length n + (i.e. X :n n+ and X n+ = fyg) and ;n (x :n ) the probability or probability density of a path conditional upon this path being of length n +. We have p ;n = Pr X :n n+ ; X n+ = fyg = ( ) n ; () and Q (x ) n M (x k ; x k ) k= ;n (x :n ) = ( ) n. (3) Now using (6) and importance sampling, this yields By sampling f (x ) = f (x ) (; x ) (; x ) + = fk (X :k ) (k; X :k ) X n= n f n (x :n) (n; x :n ) (n; x :n ) dx :n (4) n k (i) ; o (i = ; : : : ; N) according to, we can obtain the following approximation bf (x ) = N i= f k (i) k (i) ; (i) X (x ) : (5) 4
6 It is straightforward to check using (4), (5), (7), (), () and (3) that f k (i) thus establishing the unbiasedness of (). k (i) ; = W ; Similarly, the second method to compute f (x) pointwise using () corresponds to an importance sampling method on the space F, U k= fkg k. The importance distribution is given by (; x : ), () = and for n (n; x :n ) = p ;n ;n (x :n ) with p ;n = Pr (X :n n ; X n+ = fg) = ( ) n ; (6) and Q M (x; x ) n M (x k ; x k ) k= ;n (x :n ) = ( ) n. (7) Using the importance sampling identity f (x) = f (x) () () + X n= = fk (x; X ::k ) (k; X :k ) n f n (x; x :n) (n; x :n ) (n; x :n ) dx :n () n o then sampling k (i) ; (i = ; : : : ; N) according to, we obtain the following approximation bf (x) = N i= Using (4), (5), (9), (6) and (7), we have f n x; :k (i) f k (i) k (i) ; :k (i) = W thus establishing the unbiasedness of (). x; :k (i) k (i) ; :k (i) : (9) x; :k (i) 3. Limitations of SIS The estimates (5) and (9) will have a reasonable Monte Carlo variance if the variance of the absolute value of the weights is small. However, this can be di cult to ensure using the standard SIS approach. First, it imposes an arbitrary geometric distribution for the simulated paths length (), 5
7 (6) which might be inappropriate. Second, a product of terms K k ; X(i) k =M k ; X(i) k appears in the expression of the weights if M 6= K ; its variance typically increases approximately exponentially fast with the length of the paths. Third, if we are interested in estimating the function on using (5), the initial distribution appears in the denominator of (7). This distribution has to be selected very carefully to ensure that the variance of the resulting weights will be nite. The performance of SIS algorithms can usually be dramatically improved by introducing a resampling step [], [3]. The basic idea is to monitor the variance of importance weights over time and to discard paths with small weights and multiply paths with high weights; resampled paths are given an equal weight. However, even with an incorporated resampling step, SIS might still be ine cient in this context as we are interested in estimating a function at time. ach time it is used, the resampling step decreases the diversity in the number of paths left from time to the current time index. B. Importance Sampling using Trans-dimensional MCMC In this paper, we propose an alternative approach in which we do not limit ourselves to imputing paths sequentially. The importance sampling identity (4) is valid for any distribution such that R f n n (x :n ) dx :n 6= ) p ;n > and f n (x :n ) 6= ) ;n (x :n ) 6=. Similarly () is valid when R f n n (x; x :n ) dx :n 6= ) p ;n > and f n (x; x :n ) 6= ) ;n (x :n ) 6=. We now show how it is possible to structure e cient importance distributions which can be sampled from using trans-dimensional MCMC methods.. Optimal Importance Distributions We propose selecting importance distributions (n; x :n ) [resp. (n; x :n )] which minimize the variance of the absolute value of the importance weights in (5) [resp. (9)] in order to reduce the Monte Carlo variance of these estimates. Let us rst consider the case (5). We de ne (n; x :n ) on F as follows. The renormalized In many applications K is not a Markov kernel and it is impossible to select M = K: 6
8 version of the absolute value of f n (x :n ) is given by ;n (x :n ) = c ;n jf n (x :n )j () with c ;n = jf n (x :n )j dx :n : n+ Note that if g (x) and K (x; y) for any x; y, then assumption () ensures c n <. However, in the more general case, we need to make the additional assumption c n < for any n. We also consider p ;n = c c ;n () where X c = c ;n : () n= It is assumed here that c < ; this is true if () holds. In this case, where f (x ) = c ; sgn (f (x )) ; (x ) + X c ;n sgn (f n (x :n )) ;n (x :n ) dx :n n n= X = c sgn (f (x )) (; x ) + c Given N random samples approximate (3) by n= >< if u ; sgn (u) = >: if u < : n sgn (f n (x :n )) (n; x :n ) dx :n n k (i) ; o distributed according to, it is possible to bf (x ) = c N i= sgn f k (i) x : This is clearly the optimal importance distribution as the variance of the absolute values of the importance weights is equal to zero. However, it is usually impossible to sample from (n; x :n ) exactly and to compute c in closed-form. We claim that these two problems can be satisfactorily solved in most cases using transdimensional MCMC. To sample from, which is a distribution de ned on a union of subspaces of di erent dimensions, we can use any trans-dimensional MCMC method such as the popular Reversible Jump MCMC (RJMCMC) algorithm [4], [5]. This idea involves building a F -valued 7
9 n ergodic Markov chain k (i) ; which admits :k oi (i) as an invariant distribution. This is a generalization of the standard Metropolis-Hastings algorithm. As i!, one obtains (correlated) samples distributed according to. Moreover, under the standard and realistic assumption that c ; = jg (x)j dx is known or can be estimated numerically; then we can obtain the following estimate of c namely bc = c ; bp ; where bp ; is the proportion of random samples such that k (i) = ; i.e. bp ; = N k (i) : (3) i= Let us now consider the case (9). The importance distribution is de ned on F = U k= fkg k with (n; x :n ) = p ;n ;n (x :n ) (4) where and ;n (x :n ) = c;n n (x; x :n )j ; (5) c ;n = jf n (x; x :n )j dx :n n p ;n = c c ;n; (6) c = X c ;n : (7) n= It is assumed that c < ; this is satis ed if () holds. In this case, f (x) = f (x) + Given N random samples X c ;n sgn (f n (x; x :n )) n (x :n ) dx :n n n= X = f (x) + c n= n sgn (f n (x; x :n )) (n; x :n ) dx :n : n k (i) ; :k (i) o N i= distributed according to, it is possible to approximate (3) by bf (x) = f (x) + c N i= sgn f k (i) x; : :k (i)
10 To sample from, we can use trans-dimensional MCMC. To estimate c, we use the fact that if c ; = jf (x; x )j dx = jg (x ) K (x; x )j dx is known or can be estimated numerically then we can obtain the following estimate of c bc = c ; bp ; where bp ; is the proportion of random samples such that k (i) = ; i.e. bp ; = N k (i) : () i=. A Reversible Jump Markov chain Monte Carlo algorithm For sake of completeness, we describe here a simple RJMCMC algorithm to sample from as de ned by (), () and (). A very similar algorithm could be propose to sample from as de ned by (4), (5), (6) and (7). More elaborate algorithms are discussed in [5]. This algorithm is based on update, birth and death moves. ach move is respectively selected with probability u k (i), b k (i) and d k (i) where u k (i) + b k (i) + d k (i) = at iteration i. We also introduce two proposal distributions on denoted by q u (x; ) and q b () : We denote the uniform distribution on A by U (A). Initialization. Set k () ; X () :k () randomly or deterministically. Iteration i : Sample U U [; ] : If U u k Update move Set k (i) = k, sample J U ; ; : : : ; k (i) and X J q u X J ;. With probability < min : ; set = X :J ; X J k (i) ; :J ; X J X k (i) ; X ) ; X(i q J+:k (i) u q u X ; X(i ) J+:k (i) ; otherwise set J XJ ) ; X(i J ; X J = X : :k 9 = ; (9) 9
11 lse If U u k + b k, Birth move Sample J U ; ; : : : ; k, sample X J q b (). With probability < min : ; k + ; set k (i) = k +, :k = X = X : :k X :J ; X ) J ; X(i k ; X q :k b (XJ ) b k :J ; X J 9 d J:k k + = ; (3) ) ; X(i ; otherwise set k (i) = k, J:k lse Death move Sample J U ; ; : : : ; k. With probability < min : ; k ; set k (i) = k, = = X : :k ) X :J ; X(i q J+:k b X k ; X d :k k X :J ; X(i ) J+:k (i J b k 9 = ; (3) ) ; otherwise set k (i) = k, To compute (9), (3) and (3), one needs to be able to compute ratios of the form (l; x :l ) (k; x :k ) = c l ;l (x :l ) c k ;k (x :k ) = f l (x :l ) f k (x :k ) : This can be performed easily as f l (x :l ) and f k (x :k ) are given by (5). It is easy to check that the invariant distribution of this Markov chain is. stablishing ergodicity has to be performed on a case by case basis. III. Application to Value Function stimation A. Model Our motivating application is related to control. We consider a Markov process fx k g k on with transition kernel P. Let us introduce a reward function r :! R + and a discount factor
12 (; ). When the process is in state x at time k it accumulates a reward k r (x). Thus the expected reward starting from X = x is given by " # X V (x) = X=x k r (X k ) : The expected reward is called the value function in optimal control []. Under standard regularity k= assumptions, it can be established that the value function satis es V (x) = P (x; y) V (y) dy + r (x) ; that is a Fredholm equation of the second kind () where f (x) = V (x), K (x; y) = P (x; y) and g (x) = r (x). B. xample We present here a simple example where all calculations can be performed analytically which emphasizes the limitations of SIS in this context. We denote by N m; the Gaussian distribution of mean m and variance and N x; m; = p exp! (x m) : We set P (x; y) = N y; x; (with jj < ) and r (x) = N x; ; r. In this case, one has X k j (X = x) N m k (x) ; k with m (x) = x; = and for k m k (x) = k x, k =! kx : i= It follows that X f (x) = k N m k (x) ; ; k + r : k= Assume one considers a SIS method to solve this problem. A sensible choice for M is M (x; y) = ( ) P (x; y) + (x ) : If one is interested in estimating the function at a given point x = x, then the importance weights are given by (9); that is >< W = >: k (i) g k (i) ( ) if k (i) ; g(x) if k (i) = :
13 The variance of the importance weights is given by i var hw x; = :k (i) p r This variance (3) will only be nite if X k= k N m k (x) ; ; k + r= f (x) : (3) <. In this case, the optimal importance function ;n can easily be computed in closed-form as p ;n is known and ;n (x :n ) is a Gaussian; the variance of the associated estimate is zero. When estimating the function f (x ), we consider the importance weights (7) given by W = >< >: ( ) g The variance of the importance weights is equal to i var hw = k (i) g k (i) if k (i) ; if k (i) = : X k p r (x ) N m k (x ) ; ; k + r=! dx k= f (x ) dx : (33) Assume we consider (x ) = N x ; ;, then to ensure that the variance (33) is nite, it requires < and > + r : For more complex problems, it could be impossible to determine analytically what are the necessary conditions on to ensure the variance is nite. In this case, the optimal importance function ;n admits a closed-form and the variance of the associated estimate is zero. IV. Discussion A similar methodology can be developed in any situation where we face a convergent series of integrals/sums; each term being possibly of di erent dimension. In particular, this could be applied to the solving of large linear/nonlinear systems of equations and would be an alternative to [7]. V. Note The rst version of this report was written in. After having updated this report in 4, we realized that a related idea had previously appeared in computer graphics []. However, it
14 is di cult to follow the developments in this paper for people who are not familiar with this eld, including the authors of this report. We believe that this report provides a more general methodology and might also be more accessible. In the 7 update, we have added the reference [7], corrected a couple of typos and updated (once more) the a liations of both authors. References [] D.P. Bertsekas, Dynamic Programming and Optimal Control, Athena Publishers, 99. [] P. Del Moral and L. Miclo, Branching and interacting particle systems approximations of Feynman-Kac formulae with applications to non-linear ltering. In Séminaire de Probabilités XXXIV, ds. J. Azéma, M. mery, M. Ledoux, and M. Yor, Lecture Notes in Mathematics, Berlin: Springer-Verlag, vol. 79, -45,. [3] A. Doucet, J.F.G. de Freitas and N.J. Gordon (eds.), Sequential Monte Carlo Methods in Practice. New York: Springer-Verlag,. [4] P.J. Green, Reversible jump MCMC computation and Bayesian model determination, Biometrika, vol., pp. 7-73, 995. [5] P.J. Green, Trans-dimensional Markov chain Monte Carlo, in Highly Structured Stochastic Systems, Oxford University Press, 3. [6] R.C. Gri ths and S. Tavaré, Simulating probability distributions in the coalescent, Theoretical Population Biology, vol. 46, pp. 3-59, 994. [7] J.H. Halton, Sequential Monte Carlo techniques for solving non-linear systems, Monte Carlo Methods and Applications, vol., no., pp. 3-4, 6. [] R. Rubinstein, Simulation and the Monte Carlo Method, New York: Wiley, 9. [9] J. Spanier and.m. Gelbard, Monte Carlo Principles and Neutron Transportation Problems. Reading, Massachusetts: Addison-Wesley, 969. []. Veach and L.J. Guibas, Metropolis light transport, SIGGRAPH 97 Proceedings, Addison- Wesley, pp ,
Introduction 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 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 comparison of optimal and dynamic control strategies for continuous-time pension plan models
A comparison of optimal and dynamic control strategies for continuous-time pension plan models Andrew J.G. Cairns Department of Actuarial Mathematics and Statistics, Heriot-Watt University, Riccarton,
More informationHandout 4: Deterministic Systems and the Shortest Path Problem
SEEM 3470: Dynamic Optimization and Applications 2013 14 Second Term Handout 4: Deterministic Systems and the Shortest Path Problem Instructor: Shiqian Ma January 27, 2014 Suggested Reading: Bertsekas
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationMartingale Approach to Pricing and Hedging
Introduction and echniques Lecture 9 in Financial Mathematics UiO-SK451 Autumn 15 eacher:s. Ortiz-Latorre Martingale Approach to Pricing and Hedging 1 Risk Neutral Pricing Assume that we are in the basic
More informationLikelihood-based Optimization of Threat Operation Timeline Estimation
12th International Conference on Information Fusion Seattle, WA, USA, July 6-9, 2009 Likelihood-based Optimization of Threat Operation Timeline Estimation Gregory A. Godfrey Advanced Mathematics Applications
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 6 Sequential Monte Carlo methods II February
More information6.231 DYNAMIC PROGRAMMING LECTURE 10 LECTURE OUTLINE
6.231 DYNAMIC PROGRAMMING LECTURE 10 LECTURE OUTLINE Rollout algorithms Cost improvement property Discrete deterministic problems Approximations of rollout algorithms Discretization of continuous time
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationEquilibrium Asset Returns
Equilibrium Asset Returns Equilibrium Asset Returns 1/ 38 Introduction We analyze the Intertemporal Capital Asset Pricing Model (ICAPM) of Robert Merton (1973). The standard single-period CAPM holds when
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 1 Introduction January 16, 2018 M. Wiktorsson
More informationChapter 2 Uncertainty Analysis and Sampling Techniques
Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying
More informationLecture Notes 1
4.45 Lecture Notes Guido Lorenzoni Fall 2009 A portfolio problem To set the stage, consider a simple nite horizon problem. A risk averse agent can invest in two assets: riskless asset (bond) pays gross
More informationApproximating a multifactor di usion on a tree.
Approximating a multifactor di usion on a tree. September 2004 Abstract A new method of approximating a multifactor Brownian di usion on a tree is presented. The method is based on local coupling of the
More informationElif Özge Özdamar T Reinforcement Learning - Theory and Applications February 14, 2006
On the convergence of Q-learning Elif Özge Özdamar elif.ozdamar@helsinki.fi T-61.6020 Reinforcement Learning - Theory and Applications February 14, 2006 the covergence of stochastic iterative algorithms
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More informationEconS Micro Theory I Recitation #8b - Uncertainty II
EconS 50 - Micro Theory I Recitation #8b - Uncertainty II. Exercise 6.E.: The purpose of this exercise is to show that preferences may not be transitive in the presence of regret. Let there be S states
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 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 informationEE266 Homework 5 Solutions
EE, Spring 15-1 Professor S. Lall EE Homework 5 Solutions 1. A refined inventory model. In this problem we consider an inventory model that is more refined than the one you ve seen in the lectures. The
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x
More informationDynamic Programming: An overview. 1 Preliminaries: The basic principle underlying dynamic programming
Dynamic Programming: An overview These notes summarize some key properties of the Dynamic Programming principle to optimize a function or cost that depends on an interval or stages. This plays a key role
More informationRelevant parameter changes in structural break models
Relevant parameter changes in structural break models A. Dufays J. Rombouts Forecasting from Complexity April 27 th, 2018 1 Outline Sparse Change-Point models 1. Motivation 2. Model specification Shrinkage
More informationMonte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
More informationKERNEL PROBABILITY DENSITY ESTIMATION METHODS
5.- KERNEL PROBABILITY DENSITY ESTIMATION METHODS S. Towers State University of New York at Stony Brook Abstract Kernel Probability Density Estimation techniques are fast growing in popularity in the particle
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 3 Importance sampling January 27, 2015 M. Wiktorsson
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 informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}
More informationAccelerated Option Pricing Multiple Scenarios
Accelerated Option Pricing in Multiple Scenarios 04.07.2008 Stefan Dirnstorfer (stefan@thetaris.com) Andreas J. Grau (grau@thetaris.com) 1 Abstract This paper covers a massive acceleration of Monte-Carlo
More informationHandout 8: Introduction to Stochastic Dynamic Programming. 2 Examples of Stochastic Dynamic Programming Problems
SEEM 3470: Dynamic Optimization and Applications 2013 14 Second Term Handout 8: Introduction to Stochastic Dynamic Programming Instructor: Shiqian Ma March 10, 2014 Suggested Reading: Chapter 1 of Bertsekas,
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 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 informationOn solving multistage stochastic programs with coherent risk measures
On solving multistage stochastic programs with coherent risk measures Andy Philpott Vitor de Matos y Erlon Finardi z August 13, 2012 Abstract We consider a class of multistage stochastic linear programs
More informationDistortion operator of uncertainty claim pricing using weibull distortion operator
ISSN: 2455-216X Impact Factor: RJIF 5.12 www.allnationaljournal.com Volume 4; Issue 3; September 2018; Page No. 25-30 Distortion operator of uncertainty claim pricing using weibull distortion operator
More informationThe Azema Yor embedding in non-singular diusions
Stochastic Processes and their Applications 96 2001 305 312 www.elsevier.com/locate/spa The Azema Yor embedding in non-singular diusions J.L. Pedersen a;, G. Peskir b a Department of Mathematics, ETH-Zentrum,
More informationCOS 513: Gibbs Sampling
COS 513: Gibbs Sampling Matthew Salesi December 6, 2010 1 Overview Concluding the coverage of Markov chain Monte Carlo (MCMC) sampling methods, we look today at Gibbs sampling. Gibbs sampling is a simple
More informationA Markov Chain Monte Carlo Approach to Estimate the Risks of Extremely Large Insurance Claims
International Journal of Business and Economics, 007, Vol. 6, No. 3, 5-36 A Markov Chain Monte Carlo Approach to Estimate the Risks of Extremely Large Insurance Claims Wan-Kai Pang * Department of Applied
More informationContinuous-Time Pension-Fund Modelling
. Continuous-Time Pension-Fund Modelling Andrew J.G. Cairns Department of Actuarial Mathematics and Statistics, Heriot-Watt University, Riccarton, Edinburgh, EH4 4AS, United Kingdom Abstract This paper
More informationSlides for DN2281, KTH 1
Slides for DN2281, KTH 1 January 28, 2014 1 Based on the lecture notes Stochastic and Partial Differential Equations with Adapted Numerics, by J. Carlsson, K.-S. Moon, A. Szepessy, R. Tempone, G. Zouraris.
More informationBehavioral Finance and Asset Pricing
Behavioral Finance and Asset Pricing Behavioral Finance and Asset Pricing /49 Introduction We present models of asset pricing where investors preferences are subject to psychological biases or where investors
More informationOption Pricing Using Bayesian Neural Networks
Option Pricing Using Bayesian Neural Networks Michael Maio Pires, Tshilidzi Marwala School of Electrical and Information Engineering, University of the Witwatersrand, 2050, South Africa m.pires@ee.wits.ac.za,
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 informationA note on the existence of unique equivalent martingale measures in a Markovian setting
Finance Stochast. 1, 251 257 1997 c Springer-Verlag 1997 A note on the existence of unique equivalent martingale measures in a Markovian setting Tina Hviid Rydberg University of Aarhus, Department of Theoretical
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationA Practical Implementation of the Gibbs Sampler for Mixture of Distributions: Application to the Determination of Specifications in Food Industry
A Practical Implementation of the for Mixture of Distributions: Application to the Determination of Specifications in Food Industry Julien Cornebise 1 Myriam Maumy 2 Philippe Girard 3 1 Ecole Supérieure
More informationHedging Derivative Securities with VIX Derivatives: A Discrete-Time -Arbitrage Approach
Hedging Derivative Securities with VIX Derivatives: A Discrete-Time -Arbitrage Approach Nelson Kian Leong Yap a, Kian Guan Lim b, Yibao Zhao c,* a Department of Mathematics, National University of Singapore
More informationMeasuring the Wealth of Nations: Income, Welfare and Sustainability in Representative-Agent Economies
Measuring the Wealth of Nations: Income, Welfare and Sustainability in Representative-Agent Economies Geo rey Heal and Bengt Kristrom May 24, 2004 Abstract In a nite-horizon general equilibrium model national
More informationA Comparative Analysis of Crossover Variants in Differential Evolution
Proceedings of the International Multiconference on Computer Science and Information Technology pp. 171 181 ISSN 1896-7094 c 2007 PIPS A Comparative Analysis of Crossover Variants in Differential Evolution
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 SIMPLE DERIVATION OF AND IMPROVEMENTS TO JAMSHIDIAN S AND ROGERS UPPER BOUND METHODS FOR BERMUDAN OPTIONS
A SIMPLE DERIVATION OF AND IMPROVEMENTS TO JAMSHIDIAN S AND ROGERS UPPER BOUND METHODS FOR BERMUDAN OPTIONS MARK S. JOSHI Abstract. The additive method for upper bounds for Bermudan options is rephrased
More informationA Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 129 136 A Numerical Approach to the Estimation of Search Effort in a Search for a Moving
More informationOption Pricing Formula for Fuzzy Financial Market
Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationMarkowitz portfolio theory
Markowitz portfolio theory Farhad Amu, Marcus Millegård February 9, 2009 1 Introduction Optimizing a portfolio is a major area in nance. The objective is to maximize the yield and simultaneously minimize
More informationAn Experimental Study of the Behaviour of the Proxel-Based Simulation Algorithm
An Experimental Study of the Behaviour of the Proxel-Based Simulation Algorithm Sanja Lazarova-Molnar, Graham Horton Otto-von-Guericke-Universität Magdeburg Abstract The paradigm of the proxel ("probability
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationContinuous-Time Consumption and Portfolio Choice
Continuous-Time Consumption and Portfolio Choice Continuous-Time Consumption and Portfolio Choice 1/ 57 Introduction Assuming that asset prices follow di usion processes, we derive an individual s continuous
More informationGENERATION OF STANDARD NORMAL RANDOM NUMBERS. Naveen Kumar Boiroju and M. Krishna Reddy
GENERATION OF STANDARD NORMAL RANDOM NUMBERS Naveen Kumar Boiroju and M. Krishna Reddy Department of Statistics, Osmania University, Hyderabad- 500 007, INDIA Email: nanibyrozu@gmail.com, reddymk54@gmail.com
More informationFaster solutions for Black zero lower bound term structure models
Crawford School of Public Policy CAMA Centre for Applied Macroeconomic Analysis Faster solutions for Black zero lower bound term structure models CAMA Working Paper 66/2013 September 2013 Leo Krippner
More informationarxiv: v1 [q-fin.pr] 1 Nov 2013
arxiv:1311.036v1 [q-fin.pr 1 Nov 013 iance matters (in stochastic dividend discount models Arianna Agosto nrico Moretto Abstract Stochastic dividend discount models (Hurley and Johnson, 1994 and 1998,
More informationOne-Factor Models { 1 Key features of one-factor (equilibrium) models: { All bond prices are a function of a single state variable, the short rate. {
Fixed Income Analysis Term-Structure Models in Continuous Time Multi-factor equilibrium models (general theory) The Brennan and Schwartz model Exponential-ane models Jesper Lund April 14, 1998 1 Outline
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationControl Improvement for Jump-Diffusion Processes with Applications to Finance
Control Improvement for Jump-Diffusion Processes with Applications to Finance Nicole Bäuerle joint work with Ulrich Rieder Toronto, June 2010 Outline Motivation: MDPs Controlled Jump-Diffusion Processes
More informationComputational Finance Least Squares Monte Carlo
Computational Finance Least Squares Monte Carlo School of Mathematics 2019 Monte Carlo and Binomial Methods In the last two lectures we discussed the binomial tree method and convergence problems. One
More informationArbitrage, Martingales, and Pricing Kernels
Arbitrage, Martingales, and Pricing Kernels Arbitrage, Martingales, and Pricing Kernels 1/ 36 Introduction A contingent claim s price process can be transformed into a martingale process by 1 Adjusting
More informationSTOCHASTIC VOLATILITY AND OPTION PRICING
STOCHASTIC VOLATILITY AND OPTION PRICING Daniel Dufresne Centre for Actuarial Studies University of Melbourne November 29 (To appear in Risks and Rewards, the Society of Actuaries Investment Section Newsletter)
More informationYao s Minimax Principle
Complexity of algorithms The complexity of an algorithm is usually measured with respect to the size of the input, where size may for example refer to the length of a binary word describing the input,
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More information36106 Managerial Decision Modeling Monte Carlo Simulation in Excel: Part IV
36106 Managerial Decision Modeling Monte Carlo Simulation in Excel: Part IV Kipp Martin University of Chicago Booth School of Business November 29, 2017 Reading and Excel Files 2 Reading: Handout: Optimal
More informationA Skewed Truncated Cauchy Logistic. Distribution and its Moments
International Mathematical Forum, Vol. 11, 2016, no. 20, 975-988 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2016.6791 A Skewed Truncated Cauchy Logistic Distribution and its Moments Zahra
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 informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationSequential Monte Carlo Samplers
Sequential Monte Carlo Samplers Pierre Del Moral Université Nice Sophia Antipolis, France Arnaud Doucet University of British Columbia, Canada Ajay Jasra University of Oxford, UK Summary. In this paper,
More informationSTA 532: Theory of Statistical Inference
STA 532: Theory of Statistical Inference Robert L. Wolpert Department of Statistical Science Duke University, Durham, NC, USA 2 Estimating CDFs and Statistical Functionals Empirical CDFs Let {X i : i n}
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More informationAnnual risk measures and related statistics
Annual risk measures and related statistics Arno E. Weber, CIPM Applied paper No. 2017-01 August 2017 Annual risk measures and related statistics Arno E. Weber, CIPM 1,2 Applied paper No. 2017-01 August
More informationHIGHER ORDER BINARY OPTIONS AND MULTIPLE-EXPIRY EXOTICS
Electronic Journal of Mathematical Analysis and Applications Vol. (2) July 203, pp. 247-259. ISSN: 2090-792X (online) http://ejmaa.6te.net/ HIGHER ORDER BINARY OPTIONS AND MULTIPLE-EXPIRY EXOTICS HYONG-CHOL
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationFiscal policy and minimum wage for redistribution: an equivalence result. Abstract
Fiscal policy and minimum wage for redistribution: an equivalence result Arantza Gorostiaga Rubio-Ramírez Juan F. Universidad del País Vasco Duke University and Federal Reserve Bank of Atlanta Abstract
More informationStratified Sampling in Monte Carlo Simulation: Motivation, Design, and Sampling Error
South Texas Project Risk- Informed GSI- 191 Evaluation Stratified Sampling in Monte Carlo Simulation: Motivation, Design, and Sampling Error Document: STP- RIGSI191- ARAI.03 Revision: 1 Date: September
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 information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
More informationValuation of Asian Option. Qi An Jingjing Guo
Valuation of Asian Option Qi An Jingjing Guo CONTENT Asian option Pricing Monte Carlo simulation Conclusion ASIAN OPTION Definition of Asian option always emphasizes the gist that the payoff depends on
More informationA Correlated Sampling Method for Multivariate Normal and Log-normal Distributions
A Correlated Sampling Method for Multivariate Normal and Log-normal Distributions Gašper Žerovni, Andrej Trov, Ivan A. Kodeli Jožef Stefan Institute Jamova cesta 39, SI-000 Ljubljana, Slovenia gasper.zerovni@ijs.si,
More information6.231 DYNAMIC PROGRAMMING LECTURE 8 LECTURE OUTLINE
6.231 DYNAMIC PROGRAMMING LECTURE 8 LECTURE OUTLINE Suboptimal control Cost approximation methods: Classification Certainty equivalent control: An example Limited lookahead policies Performance bounds
More informationExtracting Information from the Markets: A Bayesian Approach
Extracting Information from the Markets: A Bayesian Approach Daniel Waggoner The Federal Reserve Bank of Atlanta Florida State University, February 29, 2008 Disclaimer: The views expressed are the author
More informationMarkov Processes and Applications
Markov Processes and Applications Algorithms, Networks, Genome and Finance Etienne Pardoux Laboratoire d'analyse, Topologie, Probabilites Centre de Mathematiques et d'injormatique Universite de Provence,
More informationAppendix for The Long-Run Risks Model and Aggregate Asset Prices: An Empirical Assessment
Appendix for The Long-Run Risks Model and Aggregate Asset Prices: An Empirical Assessment Jason Beeler and John Y. Campbell October 0 Beeler: Department of Economics, Littauer Center, Harvard University,
More informationMEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL
MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL Isariya Suttakulpiboon MSc in Risk Management and Insurance Georgia State University, 30303 Atlanta, Georgia Email: suttakul.i@gmail.com,
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More informationOn Pricing of Discrete Barrier Options
On Pricing of Discrete Barrier Options S. G. Kou Department of IEOR 312 Mudd Building Columbia University New York, NY 10027 kou@ieor.columbia.edu This version: April 2001 Abstract A barrier option is
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 1.1287/opre.11.864ec e-companion ONLY AVAILABLE IN ELECTRONIC FORM informs 21 INFORMS Electronic Companion Risk Analysis of Collateralized Debt Obligations by Kay Giesecke and Baeho
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
More informationFX Smile Modelling. 9 September September 9, 2008
FX Smile Modelling 9 September 008 September 9, 008 Contents 1 FX Implied Volatility 1 Interpolation.1 Parametrisation............................. Pure Interpolation.......................... Abstract
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Further Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Outline
More informationForward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions
Forward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions Pierre Henry-Labordère 1 1 Global markets Quantitative Research, SOCIÉTÉ GÉNÉRALE Outline 1 Introduction 2 Semi-linear PDEs 3
More informationBias Reduction Using the Bootstrap
Bias Reduction Using the Bootstrap Find f t (i.e., t) so that or E(f t (P, P n ) P) = 0 E(T(P n ) θ(p) + t P) = 0. Change the problem to the sample: whose solution is so the bias-reduced estimate is E(T(P
More informationDynamic Hedging and PDE Valuation
Dynamic Hedging and PDE Valuation Dynamic Hedging and PDE Valuation 1/ 36 Introduction Asset prices are modeled as following di usion processes, permitting the possibility of continuous trading. This environment
More informationHomework Assignments
Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)
More information