Stochastic Simulation
|
|
- Julianna Rich
- 5 years ago
- Views:
Transcription
1 Stochastic Simulation APPM 7400 Lesson 5: Generating (Some) Continuous Random Variables September 12, 2018 esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
2 Generating Continuous Random Variables Suppose that X is a continuous random variable with cdf F (x). esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
3 Generating Continuous Random Variables Suppose that X is a continuous random variable with cdf F (x). Furthermore, suppose that you can invert F. esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
4 Generating Continuous Random Variables Suppose that X is a continuous random variable with cdf F (x). Furthermore, suppose that you can invert F. If U unif (0, 1), the random variable F 1 (U) has the same distribution as X! esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
5 Generating Continuous Random Variables Proof: P(F 1 (U) x) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
6 Generating Continuous Random Variables Proof: P(F 1 (U) x) = P(U F (x)) F non-decreasing F 1 non-decreasing esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
7 Generating Continuous Random Variables Proof: P(F 1 (U) x) = P(U F (x)) = F (x) F non-decreasing cdf of unif(0,1) F 1 non-decreasing esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
8 Generating Continuous Random Variables Proof: P(F 1 (U) x) = P(U F (x)) = F (x) = P(X x) F non-decreasing F 1 non-decreasing cdf of unif(0,1) definition of F esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
9 Generating Continuous Random Variables Example: Generating exponential rate λ r.v. s: pdf: f (x) = λe λx, x 0 esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
10 Generating Continuous Random Variables Example: Generating exponential rate λ r.v. s: pdf: f (x) = λe λx, x 0 cdf: F (x) = x 0 λe λu du = 1 e λx, x 0 esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
11 Generating Continuous Random Variables Example: Generating exponential rate λ r.v. s: pdf: f (x) = λe λx, x 0 cdf: F (x) = x 0 λe λu du = 1 e λx, x 0 F 1 (x) = 1 λ ln(1 x) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
12 Generating Continuous Random Variables Example: Generating exponential rate λ r.v. s: pdf: f (x) = λe λx, x 0 cdf: F (x) = x 0 λe λu du = 1 e λx, x 0 F 1 (x) = 1 λ ln(1 x) Now plug in a uniform! esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
13 Generating Continuous Random Variables Results: (100,000 values, λ = 3.5) Density sample Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
14 Generating Continuous Random Variables The Standard Normal Distribution: f (x) = 1 2π e 1 2 x2, < x < Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
15 Generating Continuous Random Variables The Standard Normal Distribution: f (x) = 1 2π e 1 2 x2, < x < There is no closed form expression for the cdf! Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
16 Generating Continuous Random Variables The Box-Muller Transformation: Let U 1 and U 2 be independent unif (0, 1) r.v. s. Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
17 Generating Continuous Random Variables The Box-Muller Transformation: Let U 1 and U 2 be independent unif (0, 1) r.v. s. Then X 1 = 2 ln U 1 cos(2πu 2 ) is normally distributed with mean 0 and variance 1. Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
18 Generating Continuous Random Variables The Box-Muller Transformation: Let U 1 and U 2 be independent unif (0, 1) r.v. s. Then X 1 = 2 ln U 1 cos(2πu 2 ) is normally distributed with mean 0 and variance 1. So is X 2 = 2 ln U 1 sin(2πu 2 ) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
19 Generating Continuous Random Variables The Box-Muller Transformation: Let U 1 and U 2 be independent unif (0, 1) r.v. s. Then X 1 = 2 ln U 1 cos(2πu 2 ) is normally distributed with mean 0 and variance 1. So is X 2 = 2 ln U 1 sin(2πu 2 ) and X 1 and X 2 are independent! Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
20 Generating Continuous Random Variables So far, we have ways to simulate normal random variables random variables for which we have a closed-form invertible cdf Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
21 Generating Continuous Random Variables So far, we have ways to simulate normal random variables random variables for which we have a closed-form invertible cdf If neither of these apply to your simulation problem, move on to the next algorithm. Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
22 Generating Continuous Random Variables So far, we have ways to simulate normal random variables random variables for which we have a closed-form invertible cdf If neither of these apply to your simulation problem, move on to the next algorithm. It requires that you be able to uniformly bound the target pdf by an integrable function. Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
23 Generating Continuous Random Variables So far, we have ways to simulate normal random variables random variables for which we have a closed-form invertible cdf If neither of these apply to your simulation problem, move on to the next algorithm. It requires that you be able to uniformly bound the target pdf by an integrable function. If you can t do that, move on to the next algorithm... there is no catch-all procedure so we d better build up a big arsenal of simulation techniques! Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
24 Generating Continuous Random Variables : (due to von Neumann) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
25 Generating Continuous Random Variables : want to simulate a r.v. X with pdf f (due to von Neumann) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
26 Generating Continuous Random Variables : (due to von Neumann) want to simulate a r.v. X with pdf f need to find another function g so that g(x) f (x) x Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
27 Generating Continuous Random Variables : (due to von Neumann) want to simulate a r.v. X with pdf f need to find another function g so that g(x) f (x) x normalize g to a pdf h(x) = 1 c g(x) where c = g(x) dx Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
28 The Algorithm: Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
29 The Algorithm: generate Y with pdf h Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
30 The Algorithm: generate Y with pdf h generate U unif (0, 1) (indep of Y ) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
31 The Algorithm: generate Y with pdf h generate U unif (0, 1) (indep of Y ) if U f (Y )/g(y ), accept Y, set X = Y (Done!) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
32 The Algorithm: generate Y with pdf h generate U unif (0, 1) (indep of Y ) if U f (Y )/g(y ), accept Y, set X = Y (Done!) otherwise, reject Y, return to Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
33 Proof: (discrete case) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
34 Proof: (discrete case) Let X be the value produced by this algorithm esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
35 Proof: (discrete case) Let X be the value produced by this algorithm We will show that P(X = x) = f (x) where f is the target density. esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
36 Proof: (discrete case) Let X be the value produced by this algorithm We will show that P(X = x) = f (x) where f is the target density. Let N be the number of trials until a success? Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
37 Proof: (discrete case) Let X be the value produced by this algorithm We will show that P(X = x) = f (x) where f is the target density. Let N be the number of trials until a success? Q: Is N a geometric random variable? Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
38 Proof: (discrete case) Let X be the value produced by this algorithm We will show that P(X = x) = f (x) where f is the target density. Let N be the number of trials until a success? Q: Is N a geometric random variable? A: Yes! Success probability? Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
39 Proof: (discrete case) Consider the joint probability P(X = x, N = n) We will write down an expression for this and then marginalize out the N. Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
40 P(X = x, N = n) = P (n 1 failures, ) U n 1 > f (y n 1) g(y n 1 ), Y n = x, U n f (x) g(x) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
41 P(X = x, N = n) = P (n 1 failures, ) U n 1 > f (y n 1) g(y n 1 ), Y n = x, U n f (x) g(x) ( ) indep = P(n 1 failures) P Y n = x, U n f (x) g(x) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
42 P(X = x, N = n) = P (n 1 failures, ) U n 1 > f (y n 1) g(y n 1 ), Y n = x, U n f (x) g(x) ( ) indep = P(n 1 failures) P Y n = x, U n f (x) g(x) ( ) ident = P(n 1 failures) P Y = x, U f (x) g(x) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
43 So, P(X = x) = n=1 P(X = x, N = n) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
44 So, P(X = x) = n=1 P(X = x, N = n) ( = P Y = x, U f (x) g(x) ) n=1 P(n 1 failures) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
45 So, P(X = x) = n=1 P(X = x, N = n) ( = P ( = P Y = x, U f (x) g(x) Y = x, U f (x) g(x) ) n=1 P(n 1 failures) ) [blah (no x s)] esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
46 So, P(X = x) = n=1 P(X = x, N = n) ( = P ( = P Y = x, U f (x) g(x) Y = x, U f (x) g(x) ( indep = P(Y = x) P U f (x) g(x) ) n=1 P(n 1 failures) ) [blah (no x s)] ) [blah (no x s)] esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
47 So, P(X = x) = n=1 P(X = x, N = n) ( = P ( = P Y = x, U f (x) g(x) Y = x, U f (x) g(x) ( indep = P(Y = x) P U f (x) g(x) = h(x) f (x) g(x) [blah (no x s)] ) n=1 P(n 1 failures) ) [blah (no x s)] ) [blah (no x s)] (since f (x) g(x)) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
48 So, P(X = x) = h(x) f (x) g(x) [blah (no x s)] esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
49 So, P(X = x) = h(x) f (x) g(x) [blah (no x s)] = 1 f (x) c g(x) g(x) [blah (no x s)] esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
50 So, P(X = x) = h(x) f (x) g(x) [blah (no x s)] = 1 f (x) c g(x) g(x) [blah (no x s)] = 1 c f (x) [blah (no x s)] esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
51 So, P(X = x) = h(x) f (x) g(x) [blah (no x s)] = 1 f (x) c g(x) g(x) [blah (no x s)] = 1 c f (x) [blah (no x s)] Sum both sides over x: P(X = x) = x x { } 1 f (x) [blah (no x s)] c esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
52 P(X = x) = x x { } 1 f (x) [blah (no x s)] c Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
53 P(X = x) = x x { } 1 f (x) [blah (no x s)] c { } 1 1 = 1 [blah (no x s)] c Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
54 P(X = x) = x x { } 1 f (x) [blah (no x s)] c { } 1 1 = 1 [blah (no x s)] c { 1 c } [blah (no x s)] = 1 Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
55 P(X = x) = x x { } 1 f (x) [blah (no x s)] c { } 1 1 = 1 [blah (no x s)] c { 1 c } [blah (no x s)] = 1 P(X = x) = f (x) Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
56 Example: Simulate X Γ(2, 3) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
57 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
58 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 Will try to bound it by a function of the form g(x) = k e 2x. esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
59 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 Will try to bound it by a function of the form g(x) = k e 2x. ke 2x 9xe 3x = e 3x (ke x 9x) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
60 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 Will try to bound it by a function of the form g(x) = k e 2x. ke 2x 9xe 3x = e 3x (ke x 9x) Hmmm... k = 9 will work... f esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1 x
61 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 Will try to bound it by a function of the form g(x) = k e 2x. ke 2x 9xe 3x = e 3x (ke x 9x) Hmmm... k = 9 will work but we might try to do better since the accept probabilities are of the form f /g. f esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1 x
62 Example: Simulate X Γ(2, 3) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
63 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
64 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 Will try to bound it by a function of the form g(x) = 9k e 2x. esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
65 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 Will try to bound it by a function of the form g(x) = 9k e 2x. 9ke 2x 9xe 3x = 9e 3x (ke x x) esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
66 Example: Simulate X Γ(2, 3) Target Density: f (x) = 9xe 3x, x > 0 Will try to bound it by a function of the form g(x) = 9k e 2x. 9ke 2x 9xe 3x = 9e 3x (ke x x) k = e 1 will work esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1 x
67 Example: Simulate X Γ(2, 3) In fact, it will be best (in terms of maximizing the acceptance rate) if we can, with our choice of g, closely mimic the target density f x esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
68 Example: Simulate X Γ(2, 3) In fact, it will be best (in terms of maximizing the acceptance rate) if we can, with our choice of g, closely mimic the target density f. But, choosing something of the form g(x) = x c 1 e c 2x that will be 0 for x = 0 will result in us having to draw values (for possible acceptance) from a density that is probably at least as difficult to draw from as the original f! x esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
69 Example: Simulate X Γ(2, 3) In fact, it will be best (in terms of maximizing the acceptance rate) if we can, with our choice of g, closely mimic the target density f. We could use a piecewise uniform/exponential x esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
70 Example: Simulate X Γ(2, 3) In fact, it will be best (in terms of maximizing the acceptance rate) if we can, with our choice of g, closely mimic the target density f. We could use a piecewise uniform/exponential but this is probably not worth the trouble. Let s see how fast the algorithm is with g(x) = 9e 1 e 2x x esson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
71 Results: (100,000 values) Density mean accept time = s.d. = Density Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
72 Odds and Ends... Incidentally, how would one draw from this (unnormalized) density: f (x) { 9e 2, 0 < x < 0.5 9e 1 e 2x, x Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
73 Odds and Ends... Incidentally, how would one draw from this (unnormalized) density: could normalize then compute and and invert the piecewise continuous cdf Lesson 5: Generating (Some) Continuous Randomx Variables Stochastic Simulation September 12, / 1
74 Odds and Ends... Incidentally, how would one draw from this (unnormalized) density: could normalize then compute and and invert the piecewise continuous cdf could normalize with c, then w.p. 9e 1 /c draw a unif (0, 1) w.p. 1 9e 1 /c draw an exp(rate = 2) and shift it by adding Lesson 5: Generating (Some) Continuous Randomx Variables Stochastic Simulation September 12, / 1
75 Odds and Ends... Results: (second method candidate draws) Density Lesson 5: Generating (Some) Continuous Random Variables Stochastic Simulation September 12, / 1
IEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationUQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.
UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.
More informationcontinuous rv Note for a legitimate pdf, we have f (x) 0 and f (x)dx = 1. For a continuous rv, P(X = c) = c f (x)dx = 0, hence
continuous rv Let X be a continuous rv. Then a probability distribution or probability density function (pdf) of X is a function f(x) such that for any two numbers a and b with a b, P(a X b) = b a f (x)dx.
More informationReliability and Risk Analysis. Survival and Reliability Function
Reliability and Risk Analysis Survival function We consider a non-negative random variable X which indicates the waiting time for the risk event (eg failure of the monitored equipment, etc.). The probability
More informationImportance Sampling for Option Pricing. Steven R. Dunbar. Put Options. Monte Carlo Method. Importance. Sampling. Examples.
for for January 25, 2016 1 / 26 Outline for 1 2 3 4 2 / 26 Put Option for A put option is the right to sell an asset at an established price at a certain time. The established price is the strike price,
More informationRandom Variables Handout. Xavier Vilà
Random Variables Handout Xavier Vilà Course 2004-2005 1 Discrete Random Variables. 1.1 Introduction 1.1.1 Definition of Random Variable A random variable X is a function that maps each possible outcome
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Moments of a distribubon Measures of
More informationThe Normal Distribution
The Normal Distribution The normal distribution plays a central role in probability theory and in statistics. It is often used as a model for the distribution of continuous random variables. Like all models,
More information10. Monte Carlo Methods
10. Monte Carlo Methods 1. Introduction. Monte Carlo simulation is an important tool in computational finance. It may be used to evaluate portfolio management rules, to price options, to simulate hedging
More informationHand and Spreadsheet Simulations
1 / 34 Hand and Spreadsheet Simulations Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 9/8/16 2 / 34 Outline 1 Stepping Through a Differential Equation 2 Monte
More informationSTAT 825 Notes Random Number Generation
STAT 825 Notes Random Number Generation What if R/Splus/SAS doesn t have a function to randomly generate data from a particular distribution? Although R, Splus, SAS and other packages can generate data
More informationChapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as
Lecture 0 on BST 63: Statistical Theory I Kui Zhang, 09/9/008 Review for the previous lecture Definition: Several continuous distributions, including uniform, gamma, normal, Beta, Cauchy, double exponential
More informationIEOR 165 Lecture 1 Probability Review
IEOR 165 Lecture 1 Probability Review 1 Definitions in Probability and Their Consequences 1.1 Defining Probability A probability space (Ω, F, P) consists of three elements: A sample space Ω is the set
More information6. Continous Distributions
6. Continous Distributions Chris Piech and Mehran Sahami May 17 So far, all random variables we have seen have been discrete. In all the cases we have seen in CS19 this meant that our RVs could only take
More informationGenerating Random Variables and Stochastic Processes
IEOR E4703: Monte Carlo Simulation Columbia University c 2017 by Martin Haugh Generating Random Variables and Stochastic Processes In these lecture notes we describe the principal methods that are used
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 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 informationChapter 2: Random Variables (Cont d)
Chapter : Random Variables (Cont d) Section.4: The Variance of a Random Variable Problem (1): Suppose that the random variable X takes the values, 1, 4, and 6 with probability values 1/, 1/6, 1/, and 1/6,
More informationComputer labs. May 10, A list of matlab tutorials can be found under
Computer labs May 10, 2018 A list of matlab tutorials can be found under http://snovit.math.umu.se/personal/cohen_david/teachlinks.html Task 1: The following MATLAB code generates (pseudo) uniform random
More informationAuction theory. Filip An. U.U.D.M. Project Report 2018:35. Department of Mathematics Uppsala University
U.U.D.M. Project Report 28:35 Auction theory Filip An Examensarbete i matematik, 5 hp Handledare: Erik Ekström Examinator: Veronica Crispin Quinonez Augusti 28 Department of Mathematics Uppsala University
More informationLesson Plan for Simulation with Spreadsheets (8/31/11 & 9/7/11)
Jeremy Tejada ISE 441 - Introduction to Simulation Learning Outcomes: Lesson Plan for Simulation with Spreadsheets (8/31/11 & 9/7/11) 1. Students will be able to list and define the different components
More information4-2 Probability Distributions and Probability Density Functions. Figure 4-2 Probability determined from the area under f(x).
4-2 Probability Distributions and Probability Density Functions Figure 4-2 Probability determined from the area under f(x). 4-2 Probability Distributions and Probability Density Functions Definition 4-2
More informationReview for Final Exam Spring 2014 Jeremy Orloff and Jonathan Bloom
Review for Final Exam 18.05 Spring 2014 Jeremy Orloff and Jonathan Bloom THANK YOU!!!! JON!! PETER!! RUTHI!! ERIKA!! ALL OF YOU!!!! Probability Counting Sets Inclusion-exclusion principle Rule of product
More information2. The sum of all the probabilities in the sample space must add up to 1
Continuous Random Variables and Continuous Probability Distributions Continuous Random Variable: A variable X that can take values on an interval; key feature remember is that the values of the variable
More informationSTOR Lecture 15. Jointly distributed Random Variables - III
STOR 435.001 Lecture 15 Jointly distributed Random Variables - III Jan Hannig UNC Chapel Hill 1 / 17 Before we dive in Contents of this lecture 1. Conditional pmf/pdf: definition and simple properties.
More informationContinuous random variables
Continuous random variables probability density function (f(x)) the probability distribution function of a continuous random variable (analogous to the probability mass function for a discrete random variable),
More informationExam M Fall 2005 PRELIMINARY ANSWER KEY
Exam M Fall 005 PRELIMINARY ANSWER KEY Question # Answer Question # Answer 1 C 1 E C B 3 C 3 E 4 D 4 E 5 C 5 C 6 B 6 E 7 A 7 E 8 D 8 D 9 B 9 A 10 A 30 D 11 A 31 A 1 A 3 A 13 D 33 B 14 C 34 C 15 A 35 A
More informationNormal Distribution. Definition A continuous rv X is said to have a normal distribution with. the pdf of X is
Normal Distribution Normal Distribution Definition A continuous rv X is said to have a normal distribution with parameter µ and σ (µ and σ 2 ), where < µ < and σ > 0, if the pdf of X is f (x; µ, σ) = 1
More informationAdvanced Financial Economics Homework 2 Due on April 14th before class
Advanced Financial Economics Homework 2 Due on April 14th before class March 30, 2015 1. (20 points) An agent has Y 0 = 1 to invest. On the market two financial assets exist. The first one is riskless.
More informationStatistical Computing (36-350)
Statistical Computing (36-350) Lecture 16: Simulation III: Monte Carlo Cosma Shalizi 21 October 2013 Agenda Monte Carlo Monte Carlo approximation of integrals and expectations The rejection method and
More informationThe Normal Distribution
Will Monroe CS 09 The Normal Distribution Lecture Notes # July 9, 207 Based on a chapter by Chris Piech The single most important random variable type is the normal a.k.a. Gaussian) random variable, parametrized
More informationECE 340 Probabilistic Methods in Engineering M/W 3-4:15. Lecture 10: Continuous RV Families. Prof. Vince Calhoun
ECE 340 Probabilistic Methods in Engineering M/W 3-4:15 Lecture 10: Continuous RV Families Prof. Vince Calhoun 1 Reading This class: Section 4.4-4.5 Next class: Section 4.6-4.7 2 Homework 3.9, 3.49, 4.5,
More informationWeek 1 Quantitative Analysis of Financial Markets Probabilities
Week 1 Quantitative Analysis of Financial Markets Probabilities Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
More information4-1. Chapter 4. Commonly Used Distributions by The McGraw-Hill Companies, Inc. All rights reserved.
4-1 Chapter 4 Commonly Used Distributions 2014 by The Companies, Inc. All rights reserved. Section 4.1: The Bernoulli Distribution 4-2 We use the Bernoulli distribution when we have an experiment which
More informationChapter 2. Random variables. 2.3 Expectation
Random processes - Chapter 2. Random variables 1 Random processes Chapter 2. Random variables 2.3 Expectation 2.3 Expectation Random processes - Chapter 2. Random variables 2 Among the parameters representing
More informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
More informationCh4. Variance Reduction Techniques
Ch4. Zhang Jin-Ting Department of Statistics and Applied Probability July 17, 2012 Ch4. Outline Ch4. This chapter aims to improve the Monte Carlo Integration estimator via reducing its variance using some
More informationReview. Binomial random variable
Review Discrete RV s: prob y fctn: p(x) = Pr(X = x) cdf: F(x) = Pr(X x) E(X) = x x p(x) SD(X) = E { (X - E X) 2 } Binomial(n,p): no. successes in n indep. trials where Pr(success) = p in each trial If
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -26 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -26 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Hydrologic data series for frequency
More informationWelcome to Stat 410!
Welcome to Stat 410! Personnel Instructor: Liang, Feng TA: Gan, Gary (Lingrui) Instructors/TAs from two other sessions Websites: Piazza and Compass Homework When, where and how to submit your homework
More informationSlides for Risk Management
Slides for Risk Management Introduction to the modeling of assets Groll Seminar für Finanzökonometrie Prof. Mittnik, PhD Groll (Seminar für Finanzökonometrie) Slides for Risk Management Prof. Mittnik,
More informationChapter 6. Importance sampling. 6.1 The basics
Chapter 6 Importance sampling 6.1 The basics To movtivate our discussion consider the following situation. We want to use Monte Carlo to compute µ E[X]. There is an event E such that P(E) is small but
More informationChapter 4: Commonly Used Distributions. Statistics for Engineers and Scientists Fourth Edition William Navidi
Chapter 4: Commonly Used Distributions Statistics for Engineers and Scientists Fourth Edition William Navidi 2014 by Education. This is proprietary material solely for authorized instructor use. Not authorized
More informationLecture 23. STAT 225 Introduction to Probability Models April 4, Whitney Huang Purdue University. Normal approximation to Binomial
Lecture 23 STAT 225 Introduction to Probability Models April 4, 2014 approximation Whitney Huang Purdue University 23.1 Agenda 1 approximation 2 approximation 23.2 Characteristics of the random variable:
More informationFor more information about how to cite these materials visit
Author(s): Kerby Shedden, Ph.D., 2010 License: Unless otherwise noted, this material is made available under the terms of the Creative Commons Attribution Share Alike 3.0 License: http://creativecommons.org/licenses/by-sa/3.0/
More informationA GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies
A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL Stephen Chin and Daniel Dufresne Centre for Actuarial Studies University of Melbourne Paper: http://mercury.ecom.unimelb.edu.au/site/actwww/wps2009/no181.pdf
More informationBivariate Birnbaum-Saunders Distribution
Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators
More informationGenerating Random Numbers
Generating Random Numbers Aim: produce random variables for given distribution Inverse Method Let F be the distribution function of an univariate distribution and let F 1 (y) = inf{x F (x) y} (generalized
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 informationOverview. Transformation method Rejection method. Monte Carlo vs ordinary methods. 1 Random numbers. 2 Monte Carlo integration.
Overview 1 Random numbers Transformation method Rejection method 2 Monte Carlo integration Monte Carlo vs ordinary methods 3 Summary Transformation method Suppose X has probability distribution p X (x),
More informationSTAT Chapter 4/6: Random Variables and Probability Distributions
STAT 251 - Chapter 4/6: Random Variables and Probability Distributions We use random variables (RV) to represent the numerical features of a random experiment. In chapter 3, we defined a random experiment
More informationResults for option pricing
Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform
More informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
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 informationCentral Limit Theorem (CLT) RLS
Central Limit Theorem (CLT) RLS Central Limit Theorem (CLT) Definition The sampling distribution of the sample mean is approximately normal with mean µ and standard deviation (of the sampling distribution
More informationExpected Utility And Risk Aversion
Expected Utility And Risk Aversion Econ 2100 Fall 2017 Lecture 12, October 4 Outline 1 Risk Aversion 2 Certainty Equivalent 3 Risk Premium 4 Relative Risk Aversion 5 Stochastic Dominance Notation From
More informationLecture 10: Point Estimation
Lecture 10: Point Estimation MSU-STT-351-Sum-17B (P. Vellaisamy: MSU-STT-351-Sum-17B) Probability & Statistics for Engineers 1 / 31 Basic Concepts of Point Estimation A point estimate of a parameter θ,
More informationName: CS3130: Probability and Statistics for Engineers Practice Final Exam Instructions: You may use any notes that you like, but no calculators or computers are allowed. Be sure to show all of your work.
More informationTutorial 11: Limit Theorems. Baoxiang Wang & Yihan Zhang bxwang, April 10, 2017
Tutorial 11: Limit Theorems Baoxiang Wang & Yihan Zhang bxwang, yhzhang@cse.cuhk.edu.hk April 10, 2017 1 Outline The Central Limit Theorem (CLT) Normal Approximation Based on CLT De Moivre-Laplace Approximation
More informationPROBABILITY. Wiley. With Applications and R ROBERT P. DOBROW. Department of Mathematics. Carleton College Northfield, MN
PROBABILITY With Applications and R ROBERT P. DOBROW Department of Mathematics Carleton College Northfield, MN Wiley CONTENTS Preface Acknowledgments Introduction xi xiv xv 1 First Principles 1 1.1 Random
More informationMath489/889 Stochastic Processes and Advanced Mathematical Finance Homework 5
Math489/889 Stochastic Processes and Advanced Mathematical Finance Homework 5 Steve Dunbar Due Fri, October 9, 7. Calculate the m.g.f. of the random variable with uniform distribution on [, ] and then
More information14.30 Introduction to Statistical Methods in Economics Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 4.30 Introduction to Statistical Methods in Economics Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More information18.05 Problem Set 3, Spring 2014 Solutions
8.05 Problem Set 3, Spring 04 Solutions Problem. (0 pts.) (a) We have P (A) = P (B) = P (C) =/. Writing the outcome of die first, we can easily list all outcomes in the following intersections. A B = {(,
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 informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 13 Lecture 13 November 15, 217 Derivation of the Black-Scholes-Merton
More informationChapter 6: Point Estimation
Chapter 6: Point Estimation Professor Sharabati Purdue University March 10, 2014 Professor Sharabati (Purdue University) Point Estimation Spring 2014 1 / 37 Chapter Overview Point estimator and point estimate
More informationRandom Samples. Mathematics 47: Lecture 6. Dan Sloughter. Furman University. March 13, 2006
Random Samples Mathematics 47: Lecture 6 Dan Sloughter Furman University March 13, 2006 Dan Sloughter (Furman University) Random Samples March 13, 2006 1 / 9 Random sampling Definition We call a sequence
More information1 Auctions. 1.1 Notation (Symmetric IPV) Independent private values setting with symmetric riskneutral buyers, no budget constraints.
1 Auctions 1.1 Notation (Symmetric IPV) Ancient market mechanisms. use. A lot of varieties. Widespread in Independent private values setting with symmetric riskneutral buyers, no budget constraints. Simple
More informationCHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION
CHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION Szabolcs Sebestyén szabolcs.sebestyen@iscte.pt Master in Finance INVESTMENTS Sebestyén (ISCTE-IUL) Choice Theory Investments 1 / 65 Outline 1 An Introduction
More information1 Rare event simulation and importance sampling
Copyright c 2007 by Karl Sigman 1 Rare event simulation and importance sampling Suppose we wish to use Monte Carlo simulation to estimate a probability p = P (A) when the event A is rare (e.g., when p
More informationExpected Value and Variance
Expected Value and Variance MATH 472 Financial Mathematics J Robert Buchanan 2018 Objectives In this lesson we will learn: the definition of expected value, how to calculate the expected value of a random
More informationCS 237: Probability in Computing
CS 237: Probability in Computing Wayne Snyder Computer Science Department Boston University Lecture 12: Continuous Distributions Uniform Distribution Normal Distribution (motivation) Discrete vs Continuous
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 2 Random number generation January 18, 2018
More information1 Geometric Brownian motion
Copyright c 05 by Karl Sigman Geometric Brownian motion Note that since BM can take on negative values, using it directly for modeling stock prices is questionable. There are other reasons too why BM is
More informationBasic notions of probability theory: continuous probability distributions. Piero Baraldi
Basic notions of probability theory: continuous probability distributions Piero Baraldi Probability distributions for reliability, safety and risk analysis: discrete probability distributions continuous
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulation Efficiency and an Introduction to Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University
More informationCommonly Used Distributions
Chapter 4: Commonly Used Distributions 1 Introduction Statistical inference involves drawing a sample from a population and analyzing the sample data to learn about the population. We often have some knowledge
More informationFavorite Distributions
Favorite Distributions Binomial, Poisson and Normal Here we consider 3 favorite distributions in statistics: Binomial, discovered by James Bernoulli in 1700 Poisson, a limiting form of the Binomial, found
More informationBuilding Consistent Risk Measures into Stochastic Optimization Models
Building Consistent Risk Measures into Stochastic Optimization Models John R. Birge The University of Chicago Graduate School of Business www.chicagogsb.edu/fac/john.birge JRBirge Fuqua School, Duke University
More informationChapter 4 Continuous Random Variables and Probability Distributions
Chapter 4 Continuous Random Variables and Probability Distributions Part 2: More on Continuous Random Variables Section 4.5 Continuous Uniform Distribution Section 4.6 Normal Distribution 1 / 27 Continuous
More informationMATH 476/567 ACTUARIAL RISK THEORY FALL 2016 PROFESSOR WANG
MATH 476/567 ACTUARIAL RISK THEORY FALL 206 PROFESSOR WANG Homework 5 (max. points = 00) Due at the beginning of class on Tuesday, November 8, 206 You are encouraged to work on these problems in groups
More informationExam STAM Practice Exam #1
!!!! Exam STAM Practice Exam #1 These practice exams should be used during the month prior to your exam. This practice exam contains 20 questions, of equal value, corresponding to about a 2 hour exam.
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 informationDiscrete Random Variables and Probability Distributions. Stat 4570/5570 Based on Devore s book (Ed 8)
3 Discrete Random Variables and Probability Distributions Stat 4570/5570 Based on Devore s book (Ed 8) Random Variables We can associate each single outcome of an experiment with a real number: We refer
More informationSTA258H5. Al Nosedal and Alison Weir. Winter Al Nosedal and Alison Weir STA258H5 Winter / 41
STA258H5 Al Nosedal and Alison Weir Winter 2017 Al Nosedal and Alison Weir STA258H5 Winter 2017 1 / 41 NORMAL APPROXIMATION TO THE BINOMIAL DISTRIBUTION. Al Nosedal and Alison Weir STA258H5 Winter 2017
More informationLecture Notes 6. Assume F belongs to a family of distributions, (e.g. F is Normal), indexed by some parameter θ.
Sufficient Statistics Lecture Notes 6 Sufficiency Data reduction in terms of a particular statistic can be thought of as a partition of the sample space X. Definition T is sufficient for θ if the conditional
More informationLecture 7: Computation of Greeks
Lecture 7: Computation of Greeks Ahmed Kebaier kebaier@math.univ-paris13.fr HEC, Paris Outline 1 The log-likelihood approach Motivation The pathwise method requires some restrictive regularity assumptions
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 informationMATH 142 Business Mathematics II
MATH 142 Business Mathematics II Summer, 2016, WEEK 2 JoungDong Kim Week 2: 4.1, 4.2, 4.3, 4.4, 4.5 Chapter 4 Rules for the Derivative Section 4.1 Derivatives of Powers, Exponents, and Sums Differentiation
More informationChapter 2 ( ) Fall 2012
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 2 (2.1-2.6) Fall 2012 Definitions and Notation There are several equivalent ways to characterize the probability distribution of a survival
More informationRisk Measurement in Credit Portfolio Models
9 th DGVFM Scientific Day 30 April 2010 1 Risk Measurement in Credit Portfolio Models 9 th DGVFM Scientific Day 30 April 2010 9 th DGVFM Scientific Day 30 April 2010 2 Quantitative Risk Management Profit
More informationDefinition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ.
9 Point estimation 9.1 Rationale behind point estimation When sampling from a population described by a pdf f(x θ) or probability function P [X = x θ] knowledge of θ gives knowledge of the entire population.
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 informationMath 361. Day 8 Binomial Random Variables pages 27 and 28 Inv Do you have ESP? Inv. 1.3 Tim or Bob?
Math 361 Day 8 Binomial Random Variables pages 27 and 28 Inv. 1.2 - Do you have ESP? Inv. 1.3 Tim or Bob? Inv. 1.1: Friend or Foe Review Is a particular study result consistent with the null model? Learning
More informationNumerical valuation for option pricing under jump-diffusion models by finite differences
Numerical valuation for option pricing under jump-diffusion models by finite differences YongHoon Kwon Younhee Lee Department of Mathematics Pohang University of Science and Technology June 23, 2010 Table
More informationSolving the Black-Scholes Equation
Solving the Black-Scholes Equation An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2014 Initial Value Problem for the European Call The main objective of this lesson is solving
More information(Practice Version) Midterm Exam 1
EECS 126 Probability and Random Processes University of California, Berkeley: Fall 2014 Kannan Ramchandran September 19, 2014 (Practice Version) Midterm Exam 1 Last name First name SID Rules. DO NOT open
More informationDescribing Uncertain Variables
Describing Uncertain Variables L7 Uncertainty in Variables Uncertainty in concepts and models Uncertainty in variables Lack of precision Lack of knowledge Variability in space/time Describing Uncertainty
More informationApproximating random inequalities with. Edgeworth expansions
Approximating random inequalities with Edgeworth expansions John D. Cook November 3, 2012 Random inequalities of the form Abstract Prob(X > Y + δ often appear as part of Bayesian clinical trial methods.
More information