Optimal Importance Sampling for Diffusion Processes

Total Page:16

File Type:pdf, Size:1020Kb

Optimal Importance Sampling for Diffusion Processes U.U.D.M. Project Report 2019:38 Optimal Importance Sampling for Diffusion Processes Malvina Fröberg Examensarbete i matematik, 30 hp Handledare: Erik Ekström Examinator: Denis Gaidashev Juni 2019 Department of Mathematics Uppsala University Acknowledgements I would like to express my sincere gratitude to my supervisor, Professor Erik Ekström, for all his help and support. I am thankful for his encouragement and for introducing me to the topic, as well as for the hours spent guiding me. 1 Abstract Variance reduction techniques are used to increase the precision of estimates in numerical calculations and more specifically in Monte Carlo simulations. This thesis focuses on a particular variance reduction technique, namely importance sampling, and applies it to diffusion processes with applications within the field of financial mathematics. Importance sampling attempts to reduce variance by changing the probability measure. The Girsanov theorem is used when chang- ing measure for stochastic processes. However, a change of the probability measure gives a new drift coefficient for the underlying diffusion, which may lead to an increased computational cost. This issue is discussed and formu- lated as a stochastic optimal control problem and studied further by using the Hamilton-Jacobi-Bellman equation with a penalty term to account for compu- tational costs. The objective of this thesis is to examine whether there is an optimal change of measure or an optimal new drift for a diffusion process. This thesis provides examples of optimal measure changes in cases where the set of possible measure changes is restricted but not penalized, as well as examples for unrestricted measure changes but with penalization. 2 Contents 1 Introduction 4 2 Monte Carlo Methods 5 2.1 Monte Carlo Integration and Convergence of Error . .5 2.2 Importance Sampling . .7 2.3 Time Discretization Error . 10 2.3.1 Smooth Coefficients . 14 3 Change of Measure 17 4 Stochastic Optimal Control Problem 20 5 Importance Sampling for Diffusions 24 5.1 Introductory Problem . 24 5.1.1 Constant Coefficients . 26 5.2 Stochastic Process with Controlled Drift . 28 5.2.1 Constant Push in Specified Interval . 30 6 Optimal Importance Sampling for Diffusions 33 6.1 Other Penalty Terms . 36 6.2 The Finite Horizon Version . 37 References 39 3 1 Introduction Variance reduction techniques are used to increase the precision of estimates in nu- merical calculations and more specifically in Monte Carlo simulations. Some of the most commonly used techniques are antithetic variates, control variates and impor- tance sampling. This thesis focuses on importance sampling and applies it to diffusion processes with applications primarily in the field of financial mathematics. Importance sampling attempts to reduce variance by changing the probability mea- sure. The Girsanov theorem is used when changing measure for stochastic processes. After applying importance sampling, the resulting diffusion process that will generate a smaller variance of estimation might have, for example, an exploding drift causing large fluctuations. To handle this numerically, smaller time steps in the Monte Carlo simulation might be needed to compensate the possible loss of precision. We arrive to a problem of finding a balance between the number of sample paths, the size of the time steps and the magnitude of the drift added to the diffusion process when using the Girsanov theorem to change measure in the importance sampling method. To generalize the situation, there are two types of estimation errors that occur in the method of Monte Carlo and importance sampling, namely non-negative variance inherent in the Monte Carlo approach and an error in connection with discretization of the stochastic differential equation. To add a constant drift in the method of importance sampling is presumably not opti- mal. For the technique to be as efficient as possible, we need to allow for an exploding drift somewhere in time and space. However, non-constant (possibly exploding) co- efficients lead to numerical difficulties when simulating trajectories. Indeed, such a situation calls for adaptive methods to distribute mesh points which are more com- plex than standard methods with time steps of a fixed size. We do not wish for the new and improved stochastic process with smaller variance to be difficult to simulate. The central dilemma of this thesis is consequently the trade-off between the improve- ments due to importance sampling and the numerical efficiency of the problem and its simulation. Our next approach is to penalize these circumstances where additional computational costs arise after reducing the variance. This issue is discussed and formulated as a stochastic optimal control problem and studied further using the Hamilton-Jacobi- Bellman equation with a penalty term to account for computational cost. The objec- tive of this thesis is thus to examine whether there is an optimal change of measure or an optimal new drift for a diffusion process. Further follows a discussion on the trade-off between the benefits of importance sampling and its computational costs. In the following section, we introduce importance sampling and touch upon some of the numerical issues in Monte Carlo methods. In Section 3 and 4, we go through needed background material such as change of measure and the Girsanov theorem, but also stochastic optimal control problem and the Hamilton-Jacobi-Bellman equation. In Section 5 and 6, we apply these theorems to carefully selected examples. 4 2 Monte Carlo Methods Monte Carlo methods are a class of computational algorithms that can be applied to vast ranges of problems. They provide approximate solutions and are used in cases where analytical or numerical solutions do not exist or are too difficult to implement. It is a computational algorithm that makes numerical estimations by taking the empirical mean of repeated random sampling. It is an easy way of modeling complex situations which allows for applications in a wide range of fields such as finance and engineering. When simulating Monte Carlo methods there are two main factors that affect the cost-effectiveness: the number of sample paths and the size of the time steps. Let N be the number of sample paths and h = ∆t be the size of the time steps in the general Monte Carlo integration method. Then, according to Seydel (2009) and Hirsa (2012), the numerical errors or rates of convergence depending on h and N are p h = O( h); respectively, p N = O(1= N): Further explanations of these rates of convergence are found in Section 2.1 and 2.3. 2.1 Monte Carlo Integration and Convergence of Error Assume a probability distribution with density f, then the expectation of a function h is Z E[h(X)] = h(x)f(x)dx: R In the one-dimensional case, for a definite integral on some interval I = [a; b], we use the uniform distribution with density 1 1 f = 1 = 1 ; b − a I d(I) I where d(I) denotes the length of the interval I. Let Z b α := d(I)E[h(X)] = h(x)dx: a For independent samples Xi ∼ U[a; b], the law of large numbers implies that the approximation N 1 X α^ := d(I) h(X ) ! α a.s. as N ! 1: N N i i=1 To generalize for the higher dimensional case, let I ⊂ Rm. We want to calculate the integral Z αm := h(x)dx: I 5 Again, we draw independent and uniformly distributed samples X1; :::; XN 2 I, then we get the approximation N 1 X α^m := d (I) h(X ); N m N i i=1 where dm(I) < 1 now is the volume, or the m-dimensional Lebesgue measure, m m of I. Following the law of large numbers, α^N converges almost surely to α = R dm(I)E[h(X)] = I h(x)dx as N ! 1. Let Z m δN := h(x)dx − α^N I be the error. Before deriving the variance of the error, let us examine the zero mean and correlation properties. We have N Z 1 X δ = h(x)dx − d (I) h(X ) N m N i I i=1 N 1 X Z = h(x)dx − d (I)h(X ) N m i i=1 I N 1 X Z 1 = h(x) dx − h(X ) d (I) N d (I) i m i=1 I m N Z dm(I) X 1 = h(x) dx − h(X ) : N d (I) i i=1 I m It is easy to show that R h(x) 1 dx − h(X ) has zero mean and considering X and I dm(I) i i X are independent for i 6= j, then R h(x) 1 dx−h(X ) and R h(x) 1 dx−h(X ) j I dm(I) i I dm(I) j are uncorrelated. We have Z 1 E h(x) dx − h(Xi) = 0 I dm(I) and Z 1 Z 1 E h(x) dx − h(Xi) h(x) dx − h(Xj) = I dm(I) I dm(I) Z 1 Z 1 E h(x) dx − h(Xi) E h(x) dx − h(Xj) = 0: I dm(I) I dm(I) Further, we can have a look at the variance of the error 2 2 Var(δN ) = E[δN ] − (E[δN ]) 2 = E[δN ] 2 2 N "Z 2# (dm(I)) X 1 = E h(x) dx − h(X ) N 2 d (I) i i=1 I m (d (I))2 = m Var(h); N 6 where the variance of h is Z 1 Z 1 2 Var(h) := h2(x) dx − h(x) dx : I dm(I) I dm(I) Thus, the standard deviation of the error δN tends to zero with the order p p N := Var(δN ) = O(1= N): Square integrability of h suffices (h 2 L2), the integrands h need not be smooth (Seydel, 2009).
Recommended publications
  • Black-Scholes Model
    Chapter 5 Black-Scholes Model Copyright c 2008{2012 Hyeong In Choi, All rights reserved. 5.1 Modeling the Stock Price Process One of the basic building blocks of the Black-Scholes model is the stock price process. In particular, in the Black-Scholes world, the stock price process, denoted by St, is modeled as a geometric Brow- nian motion satisfying the following stochastic differential equation: dSt = St(µdt + σdWt); where µ and σ are constants called the drift and volatility, respec- tively. Also, σ is always assumed to be a positive constant. In order to see if such model is \reasonable," let us look at the price data of KOSPI 200 index. In particular Figure 5.1 depicts the KOSPI 200 index from January 2, 2001 till May 18, 2011. If it is reasonable to model it as a geometric Brownian motion, we must have ∆St ≈ µ∆t + σ∆Wt: St In particular if we let ∆t = h (so that ∆St = St+h − St and ∆Wt = Wt+h − Wt), the above approximate equality becomes St+h − St ≈ µh + σ(Wt+h − Wt): St 2 Since σ(Wt+h − Wt) ∼ N(0; σ h), the histogram of (St+h − St)p=St should resemble a normal distribution with standard deviation σ h. (It is called the h-day volatility.) To see it is indeed the case we have plotted the histograms for h = 1; 2; 3; 5; 7 (days) in Figures 5.2 up to 5.1. MODELING THE STOCK PRICE PROCESS 156 Figure 5.1: KOSPI 200 from 2001.
    [Show full text]
  • Lecture 4: Risk Neutral Pricing 1 Part I: the Girsanov Theorem
    SEEM 5670 { Advanced Models in Financial Engineering Professor: Nan Chen Oct. 17 &18, 2017 Lecture 4: Risk Neutral Pricing 1 Part I: The Girsanov Theorem 1.1 Change of Measure and Girsanov Theorem • Change of measure for a single random variable: Theorem 1. Let (Ω; F;P ) be a sample space and Z be an almost surely nonnegative random variable with E[Z] = 1. For A 2 F, define Z P~(A) = Z(!)dP (!): A Then P~ is a probability measure. Furthermore, if X is a nonnegative random variable, then E~[X] = E[XZ]: If Z is almost surely strictly positive, we also have Y E[Y ] = E~ : Z Definition 1. Let (Ω; F) be a sample space. Two probability measures P and P~ on (Ω; F) are said to be equivalent if P (A) = 0 , P~(A) = 0 for all such A. Theorem 2 (Radon-Nikodym). Let P and P~ be equivalent probability measures defined on (Ω; F). Then there exists an almost surely positive random variable Z such that E[Z] = 1 and Z P~(A) = Z(!)dP (!) A for A 2 F. We say Z is the Radon-Nikodym derivative of P~ with respect to P , and we write dP~ Z = : dP Example 1. Change of measure on Ω = [0; 1]. Example 2. Change of measure for normal random variables. • We can also perform change of measure for a whole process rather than for a single random variable. Suppose that there is a probability space (Ω; F) and a filtration fFt; 0 ≤ t ≤ T g. Suppose further that FT -measurable Z is an almost surely positive random variable satisfying E[Z] = 1.
    [Show full text]
  • Risk Neutral Measures
    We will denote expectations and conditional expectations with respect to the new measure by E˜. That is, given a random variable X, Z E˜X = EZ(T )X = Z(T )X dP . CHAPTER 4 Also, given a σ-algebra F, E˜(X | F) is the unique F-measurable random variable such that Risk Neutral Measures Z Z (1.2) E˜(X | F) dP˜ = X dP˜ , F F Our aim in this section is to show how risk neutral measures can be used to holds for all F measurable events F . Of course, equation (1.2) is equivalent to price derivative securities. The key advantage is that under a risk neutral measure requiring the discounted hedging portfolio becomes a martingale. Thus the price of any Z Z 0 ˜ ˜ derivative security can be computed by conditioning the payoff at maturity. We will (1.2 ) Z(T )E(X | F) dP = Z(T )X dP , F F use this to provide an elegant derivation of the Black-Scholes formula, and discuss for all F measurable events F . the fundamental theorems of asset pricing. The main goal of this section is to prove the Girsanov theorem. 1. The Girsanov Theorem. Theorem 1.6 (Cameron, Martin, Girsanov). Let b(t) = (b1(t), b2(t), . , bd(t)) be a d-dimensional adapted process, W be a d-dimensional Brownian motion, and ˜ Definition 1.1. Two probability measures P and P are said to be equivalent define if for every event A, P (A) = 0 if and only if P˜(A) = 0. Z t W˜ (t) = W (t) + b(s) ds .
    [Show full text]
  • The Girsanov Theorem Without (So Much) Stochastic Analysis Antoine Lejay
    The Girsanov theorem without (so much) stochastic analysis Antoine Lejay To cite this version: Antoine Lejay. The Girsanov theorem without (so much) stochastic analysis. 2017. hal-01498129v1 HAL Id: hal-01498129 https://hal.inria.fr/hal-01498129v1 Preprint submitted on 29 Mar 2017 (v1), last revised 22 Sep 2018 (v3) HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. The Girsanov theorem without (so much) stochastic analysis Antoine Lejay March 29, 2017 Abstract In this pedagogical note, we construct the semi-group associated to a stochastic differential equation with a constant diffusion and a Lipschitz drift by composing over small times the semi-groups generated respectively by the Brownian motion and the drift part. Similarly to the interpreta- tion of the Feynman-Kac formula through the Trotter-Kato-Lie formula in which the exponential term appears naturally, we construct by doing so an approximation of the exponential weight of the Girsanov theorem. As this approach only relies on the basic properties of the Gaussian distribution, it provides an alternative explanation of the form of the Girsanov weights without referring to a change of measure nor on stochastic calculus.
    [Show full text]
  • Comparison of Option Pricing Between ARMA-GARCH and GARCH-M Models
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Scholarship@Western Western University Scholarship@Western Electronic Thesis and Dissertation Repository 4-22-2013 12:00 AM Comparison of option pricing between ARMA-GARCH and GARCH-M models Yi Xi The University of Western Ontario Supervisor Dr. Reginald J. Kulperger The University of Western Ontario Graduate Program in Statistics and Actuarial Sciences A thesis submitted in partial fulfillment of the equirr ements for the degree in Master of Science © Yi Xi 2013 Follow this and additional works at: https://ir.lib.uwo.ca/etd Part of the Other Statistics and Probability Commons Recommended Citation Xi, Yi, "Comparison of option pricing between ARMA-GARCH and GARCH-M models" (2013). Electronic Thesis and Dissertation Repository. 1215. https://ir.lib.uwo.ca/etd/1215 This Dissertation/Thesis is brought to you for free and open access by Scholarship@Western. It has been accepted for inclusion in Electronic Thesis and Dissertation Repository by an authorized administrator of Scholarship@Western. For more information, please contact [email protected]. COMPARISON OF OPTION PRICING BETWEEN ARMA-GARCH AND GARCH-M MODELS (Thesis format: Monograph) by Yi Xi Graduate Program in Statistics and Actuarial Science A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science The School of Graduate and Postdoctoral Studies The University of Western Ontario London, Ontario, Canada ⃝c Yi Xi 2013 Abstract Option pricing is a major area in financial modeling. Option pricing is sometimes based on normal GARCH models. Normal GARCH models fail to capture the skewness and the leptokurtosis in financial data.
    [Show full text]
  • Equivalent and Absolutely Continuous Measure Changes for Jump-Diffusion
    The Annals of Applied Probability 2005, Vol. 15, No. 3, 1713–1732 DOI: 10.1214/105051605000000197 c Institute of Mathematical Statistics, 2005 EQUIVALENT AND ABSOLUTELY CONTINUOUS MEASURE CHANGES FOR JUMP-DIFFUSION PROCESSES By Patrick Cheridito1, Damir Filipovic´ and Marc Yor Princeton University, University of Munich and Universit´e Pierre et Marie Curie, Paris VI We provide explicit sufficient conditions for absolute continuity and equivalence between the distributions of two jump-diffusion pro- cesses that can explode and be killed by a potential. 1. Introduction. The purpose of this paper is to give explicit, easy-to- check sufficient conditions for the distributions of two jump-diffusion pro- cesses to be equivalent or absolutely continuous. We consider jump-diffusions that can explode and be killed by a potential. These processes are, in general, not semimartingales. We characterize them by their infinitesimal generators. The structure of the paper is as follows. In Section 2 we introduce the notation and state the paper’s main result, which gives sufficient conditions for the distributions of two jump-diffusions to be equivalent or absolutely continuous. The conditions consist of local bounds on the transformation of one generator into the other one and the assumption that the martingale problem for the second generator has for all initial distributions a unique solution. The formulation of the main theorem involves two sequences of stopping times. Stopping times of the first sequence stop the process before it explodes. The second sequence consists of exit times of the process from regions in the state space where the transformation of the first generator into the second one can be controlled.
    [Show full text]
  • Week 10 Change of Measure, Girsanov
    Week 10 Change of measure, Girsanov Jonathan Goodman November 25, 2013 1 Reweighting Suppose X is a random variable with probability density u(x). Then the ex- pected value of f(X) is Z Eu[ f(X)] = f(x)u(x) dx : (1) Suppose v(x) is a different probability density. The likelihood ratio is u(x) L(x) = : (2) v(x) The expectation (1) can be expressed as Z Eu[ f(X)] = f(x)u(x) dx Z u(x) = f(x) v(x) dx v(x) Eu[ f(X)] = Ev[ f(X)L(X)] : (3) This formula represents reweighting. The likelihood ratio, L(x) turn u into v. The expectation value is unchanged because L is included in the expectation with respect to v. This section is about Girsanov theory, which is reweighting one diffusion process to get a different diffusion process. The specific topics are: 1. When are two diffusions related by a reweighting? Answer: if they have the same noise. You can adjust the drift by reweighting, but not the noise. 2. If two diffusions are related by a reweighting, what likelihood ratio L(x[1;T ])? The path x[0;T ] is the random variable and L is a function of the random variable. Answer: L(x[0;T ]) is given in terms of a stochastic integral. We give a direct derivation, which is straightforward but might be considered complicated. Then we show how to verify that the stochastic integral for- mula is correct if you are presented with it. This is considered (by me) a less useful approach, albeit more common.
    [Show full text]
  • Lecture Notes for Chapter 6 1 Prototype Model: a One-Step
    Lecture Notes for Chapter 6 This is the chapter that brings together the mathematical tools (Brownian motion, It^ocalculus) and the financial justifications (no-arbitrage pricing) to produce the derivative pricing methodology that dominated derivative trading in the last 40 years. We have seen that the no-arbitrage condition is the ultimate check on a pricing model and we have seen how this condition is verified on a binomial tree. Extending this to a model that is based on a stochastic process that can be made into a martingale under certain probability measure, a formal connection is made with so-called risk-neutral pricing, and the probability measure involved leads to a world that is called risk-neutral world. We will examine three aspects involved in this theory: no-arbitrage condition, martingale, and risk-neutral world, and show how these concepts lead to the pricing methodology as we know today. 1 Prototype model: a one-step binomial tree First we take another look at the one-step binomial model and inspect the impli- cations of no-arbitrage condition imposed on this model, and the emergence of the so-called risk-neutral probability measure. To focus on the main issue, let us assume that the risk-free interest rate r = 0. We have this simplest model for a stock price over one period of time: it starts with a price S0, and can move to one of two states after this time period: S+ or S− (assuming S+ > S−). In order that there is no-arbitrage opportunity for any investor, we must have S− < S0 < S+.
    [Show full text]
  • Stochastic Calculus, Week 8 Equivalent Martingale Measures
    Stochastic Calculus, Week 8 Intuition via Binomial Approximation n Equivalent Martingale Measures Let Wt (ω),t∈ [0,T] be defined by √ Xi T Outline W n(ω)= δt Y (ω),t= iδt, δt= , t j n j=1 1. Intuition via binomial approximation where ω =(ω1,...,ωn) represents the combined outcome of n consecutive binary draws, where each draw may be u 2. Cameron-Martin-Girsanov Theorem (up) or d (down); further, 3. Application to derivative pricing: ( 1 if ωi = u, equivalent martingale measures/risk neutral measures, Yi(ω)= −1 if ωi = d. martingale representation theorem 4. Obtaining the Black-Scholes formula once more Let Ω be the set of all possible ω’s, let F be the appropri- ate σ-field, and let P and Q be two probability measures on (Ω, F), such that the draws are independent under both P and Q, and 1 P(ω = u)=p = , i 2 1 √ Q(ω = u)=q = 1 − γ δt , i 2 for some constant γ. 1 n d ˜ n d ˜ As n →∞, Wt (ω) −→ Wt(ω) (convergence in distribu- Then Wt (ω) −→ Wt(ω), where tion), where Brownian motion with drift γ under P, ˜ Standard Brownian motion under P, Wt(ω)= W (ω)= t Standard Brownian motion under Q. Brownian motion with drift − γ under Q. In summary, we can write down the following table with The measures P and Q are equivalent, which means that drifts: for any even A ∈F,wehaveP(A) > 0 if and only if PQ Q(A) > 0. This requires that γ is such that 0 <q<1.
    [Show full text]
  • Maximizing Expected Logarithmic Utility in a Regime-Switching Model with Inside Information
    Maximizing Expected Logarithmic Utility in a Regime-Switching Model with Inside Information Kenneth Tay Advisor: Professor Ramon van Handel Submitted in partial fulfillment of the requirements for the degree of Bachelor of Arts Department of Mathematics Princeton University May 2010 I hereby declare that I am the sole author of this thesis. I authorize Princeton University to lend this thesis to other institutions or individuals for the purpose of scholarly research. Kenneth Tay I further authorize Princeton University to reproduce this thesis by photocopying or by other means, in total or in part, at the request of other institutions or individuals for the purpose of scholarly research. Kenneth Tay Acknowledgements This page is possibly the most important page in this thesis. Without the the love and support of those around me, this work would not have been completed. I would like to dedicate this thesis to the following people: Mum and Dad, for always being there for me. For pouring your heart and soul into my life, and for encouraging me to pursue my passions. Professor van Handel, my thesis advisor. Thank you for guiding me through this entire process, for your patience in explaining ideas, for challenging me to give my very best. Professor Calderbank, my JP advisor and second reader for my thesis, your dedication to the field of applied mathematics has inspired me to explore this area. Professor Xu Jiagu, the professors at SIMO and SIMO X-Men, for nurturing and honing my skills in mathematics in my earlier years. I am deeply honored to have been taught by so many brilliant minds.
    [Show full text]
  • Local Volatility, Stochastic Volatility and Jump-Diffusion Models
    IEOR E4707: Financial Engineering: Continuous-Time Models Fall 2013 ⃝c 2013 by Martin Haugh Local Volatility, Stochastic Volatility and Jump-Diffusion Models These notes provide a brief introduction to local and stochastic volatility models as well as jump-diffusion models. These models extend the geometric Brownian motion model and are often used in practice to price exotic derivative securities. It is worth emphasizing that the prices of exotics and other non-liquid securities are generally not available in the market-place and so models are needed in order to both price them and calculate their Greeks. This is in contrast to vanilla options where prices are available and easily seen in the market. For these more liquid options, we only need a model, i.e. Black-Scholes, and the volatility surface to calculate the Greeks and perform other risk-management tasks. In addition to describing some of these models, we will also provide an introduction to a commonly used fourier transform method for pricing vanilla options when analytic solutions are not available. This transform method is quite general and can also be used in any model where the characteristic function of the log-stock price is available. Transform methods now play a key role in the numerical pricing of derivative securities. 1 Local Volatility Models The GBM model for stock prices states that dSt = µSt dt + σSt dWt where µ and σ are constants. Moreover, when pricing derivative securities with the cash account as numeraire, we know that µ = r − q where r is the risk-free interest rate and q is the dividend yield.
    [Show full text]
  • On Moment Conditions for the Girsanov Theorem See Keong Lee Louisiana State University and Agricultural and Mechanical College, Sklin [email protected]
    Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2006 On moment conditions for the Girsanov Theorem See Keong Lee Louisiana State University and Agricultural and Mechanical College, [email protected] Follow this and additional works at: https://digitalcommons.lsu.edu/gradschool_dissertations Part of the Applied Mathematics Commons Recommended Citation Lee, See Keong, "On moment conditions for the Girsanov Theorem" (2006). LSU Doctoral Dissertations. 998. https://digitalcommons.lsu.edu/gradschool_dissertations/998 This Dissertation is brought to you for free and open access by the Graduate School at LSU Digital Commons. It has been accepted for inclusion in LSU Doctoral Dissertations by an authorized graduate school editor of LSU Digital Commons. For more information, please [email protected]. ON MOMENT CONDITIONS FOR THE GIRSANOV THEOREM A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements for the degree of Doctor of Philosophy in The Department of Mathematics by See Keong Lee B.Sc., University of Science, Malaysia, 1998 M.Sc., University of Science, Malaysia, 2001 M.S., Louisiana State University, 2003 December 2006 Acknowledgments It is a pleasure to express my sincere appreciation to my advisor, Prof. Hui-Hsiung Kuo. Without his invaluable advice, help and guidance, this dissertation would not be completed. I also like to thank him for introducing me to the field of stochastic integration. He has been extremely generous with his insight and knowledge. I also like to thank the committee members, Prof. Augusto Nobile, Prof. Jorge Morales, Prof. Robert Perlis, Prof.
    [Show full text]