General Principles in Random Variates Generation
|
|
- Lydia Marsh
- 6 years ago
- Views:
Transcription
1 General Principles in Random Variates Generation E. Moulines and G. Fort Telecom ParisTech June 2015 Bibliography : Luc Devroye, Non-Uniform Random Variate Generator, Springer-Verlag (1986) available on the webpage of the author http ://luc.devroye.org/books-luc.html
2 Outline Uniform Random number generators Pesudo-Random sequences Inversion Methods for distributions on R Cumulative distribution function Algorithm Sampling real valued random variables When F 1 (u) is not explicit Simulation of Gaussian variables Box-Müller Marsaglia-Bray Inversion method The rejection method (on R d ) Sampling under a curve Conditional distribution Algorithm Design parameters Example 1 : a normal generator based on the rejection algorithm Example 2 : Gaussian from Cauchy
3 Uniform Random number generators Uniform Random number generators Pesudo-Random sequences Inversion Methods for distributions on R Cumulative distribution function Algorithm Sampling real valued random variables When F 1 (u) is not explicit Simulation of Gaussian variables Box-Müller Marsaglia-Bray Inversion method The rejection method (on R d ) Sampling under a curve Conditional distribution Algorithm Design parameters Example 1 : a normal generator based on the rejection algorithm Example 2 : Gaussian from Cauchy
4 Uniform Random number generators Pesudo-Random sequences Pseudo-random sequences (1/3) The key ingredient for Monte Carlo methods is a generator of random number. We will see in this lesson that the generator of uniform random variable on [0,1] is a central tool for sampling more general distributions. Therefore, how to sample from a uniform distribution on [0, 1]? i.e. how to obtain a sequence of numbers u 1,, u n, that can be considered as a path of the random sequence (U 1,, U n, ) where the r.v. (U n) n are independent and with the same distribution U([0, 1]). Well, we are only able to produce a sequence (u 1,, u n, ) from a machine which is a pseudo random sequence.
5 Uniform Random number generators Pesudo-Random sequences Pseudo-random sequences (2/3) Let us start with the question of sampling a sequence of {0, 1}-valued r.v. with a (fair) coin. Examples (L = 12) All of these sequences of length L occur with probability 1/2 L. But are they random?
6 Uniform Random number generators Pesudo-Random sequences Pseudo-random sequences (2/3) Let us start with the question of sampling a sequence of {0, 1}-valued r.v. with a (fair) coin. Examples (L = 12) All of these sequences of length L occur with probability 1/2 L. But are they random? Champernowne sequence : }{{ 0 } }{{} 3 1 } {{ 0 0 } 4 1 } {{ 0 1 } 5
7 Uniform Random number generators Pesudo-Random sequences Pseudo-random sequences (3/3) For any length L and any binary motif (a 1,, a L), the Law of Large Numbers for i.i.d. r.v. implies 1 n n k=1 1 (a1,,a L )(U (k 1)L+1,, U kl ) n 1 2 L w.p.1 A sequence (u n) n is said -uniform if for any L 1 and any (a 1,, a L) {0, 1} L, 1 n n k=1 1 (a1,,a L )(u (k 1)L+1,, u kl ) 1 2 L It can be proved that the Champernowne sequence is -uniform...
8 Uniform Random number generators Pesudo-Random sequences Pseudo-random generator The definition of a generator consists in a finite sequence of states (x 1,, x M ), a mapping T : x k+1 = T (x k ) a mapping S : u k = S(x k ) an output sequence (u n) n. The output sequence is called the random sequence produced by this generator. Example : The best known and still most widely used generators are the simple linear congruential generators. x 0 (seed) x k = (ax k 1 + b) modm u k = x k M The properties (e.g. cycles) of this generator depends on a, b, M (M = 2 32 )
9 Uniform Random number generators Pesudo-Random sequences Bibliography More complex generators exist, with improved properties. See e.g. Cours d E. Moulines, sur le site pédagogique du cours MDI 345, rubriques : Introduction aux méthodes de Monte Carlo et à la Simulation Méthodes de simulation de v.a. uniformes Cours de B. Ycart méthodes de Monte Carlo, http ://ljk.imag.fr/membres/bernard.ycart/polys/polys.html Livre de S. Asmussen and P.W. Glynn Stochastic Simulation, Springer. Chapter II section 2.1. Livre de P. Glasserman Monte Carlo methods in Financial Engineering, Springer. Chapter II section 2.1.
10 Inversion Methods for distributions on R Uniform Random number generators Pesudo-Random sequences Inversion Methods for distributions on R Cumulative distribution function Algorithm Sampling real valued random variables When F 1 (u) is not explicit Simulation of Gaussian variables Box-Müller Marsaglia-Bray Inversion method The rejection method (on R d ) Sampling under a curve Conditional distribution Algorithm Design parameters Example 1 : a normal generator based on the rejection algorithm Example 2 : Gaussian from Cauchy
11 Inversion Methods for distributions on R Cumulative distribution function Cumulative distribution function Let F be a cumulative distribution function on R, F (x) = P(X x). Properties of F : The function x F (x) is non-decreasing. At every point, F has left limits lim F (y) y x exists F is continuous to the right For any x R, F (x) = lim y x + F (y) F (x + ) F (x ) = P(X = x).
12 General Principles in Random Variates Generation Inversion Methods for distributions on R Cumulative distribution function The quantile function Here are examples of cumulative distribution functions F : F is not necessarily invertible. The generalized inverse of the function F, denoted F 1, is defined on (0, 1) by F 1 (p) = inf{y R, F (y) p}, 0 < p < 1. This function is called the quantile function. On the figures above : [left] F 1 : (0, 1) R ; [center] F 1 : (0, 1) (0, ) ; [right] F 1 : (0, 1) (2, 8)
13 Inversion Methods for distributions on R Cumulative distribution function Properties of the quantile function F 1 (p) = inf{y R, F (y) p}, 0 < p < 1. F 1 is a proper inverse if and only if F is continuous and strictly increasing. Properties : F 1 is nondecreasing. F F 1 (p) p ; equality can fail only if F is discontinuous at F 1 (p). F 1 F (x) x ; equality fails iff x is in the interior or at the right end of a flat of F. F 1 (p) x if and only if p F (x).
14 Inversion Methods for distributions on R Cumulative distribution function 1 p 0.8 p=f(x) F(F 1 (p)) p x F 1 (p) F 1 (p) =F 1 (F(x)) F 1 (p)
15 Inversion Methods for distributions on R Algorithm Simulation using the quantile transformation Theorem (Inverse CDF method) Let F be a cumulative distribution function and U U([0, 1]). Then the cumulative distribution function of the r.v. F 1 (U) is F. This is called the inverse transform sampling method (also known as the inverse probability integral transform). The interest for this method stems from the fact that many programming languages have the ability to generate pseudorandom numbers which are (almost) distributed according to i.i.d. standard uniform r.v. see the first section of these slides, anf refs therein Theorem If the cumulative distribution function F of the r.v. X is continuous, then F (X) U([0, 1]). This is called the probability integral transformation
16 Inversion Methods for distributions on R Sampling real valued random variables Sampling discrete random variables Case 1 : X takes values in {a 1,, a n} with probabilities p 1,, p n 0 if x < a 1 F (x) = p p j 1 if a j 1 x < a j 1 if x a n and for any p (0, 1), { F 1 a1 if 0 < p p 1 (p) = a j ifp p j 1 < p p p j
17 Inversion Methods for distributions on R Sampling real valued random variables Sampling discrete random variables Case 1 : X takes values in {a 1,, a n} with probabilities p 1,, p n 0 if x < a 1 F (x) = p p j 1 if a j 1 x < a j 1 if x a n and for any p (0, 1), { F 1 a1 if 0 < p p 1 (p) = a j ifp p j 1 < p p p j Hence, the algorithm is Draw U U([0, 1]) Find J {1,, n} such that p p J 1 < U p p J. Return a J.
18 Inversion Methods for distributions on R Sampling real valued random variables Case 2 : X takes values in the countable set {a 1,, a j, } with probabilities p 1,, p j,. The algorithm is Draw U U([0, 1]) Find J 1 such that p p J 1 < U p p J. Return a J. Remark : strategies have ben developed to solve efficiently find the index J such that p p J 1 < U p p J (see TD and the book Non-uniform Random Variate Generation by L. Devroye)
19 Inversion Methods for distributions on R Sampling real valued random variables Case 3 : Sampling random variable with continuous (and explicit) cdf support f F (x) X = F 1 (U) Exponential x > 0 λe λx 1 e λx λ 1 log(u) σ 1 Cauchy R + 1 atan ( ) x σtan(πu) π(x 2 +σ ) 2 π σ x Rayleigh R σ e x2 /2σ 2 1 e x2 /2σ 2 σ 2 log(u) Pareto x b > 0 1 ( b x ab a x a+1 ) a b U 1/a The cdf of a Gaussian random variable does not have an explicit expression. We will see later how to sample such a distribution from a generator of uniform random variables.
20 Inversion Methods for distributions on R When F 1 (u) is not explicit When F 1 is not explicit (1/2) The inversion method is exact when an explicit form of F 1 is known. In other cases, we must solve the equation F (x) = p and this requires an infinite amount of time if F is continuous. Any stopping rule that we use with the numerical method leads necessarily to an inexact algorithm...
21 Inversion Methods for distributions on R When F 1 (u) is not explicit When F 1 is not explicit (2/2) : The Bissection Method When F 1 (u) is not explicit, the solution x such that F (x) = u can be numerically approximated by the Bissection method Algorithm : Find an initial interval [a, b] to which the solution belongs repeat X (a + b)/2 if F (X) u then a X else b X end if until b a 2δ Return X Rmk : This algorithm may never work if, for fixed u, there is an interval of solution to F (x) = u. Nevertheless, it can be proved that the set is of null Lebesgue-measure. {u (0, 1) : x < y s.t.f (x) = F (y) = u}
22 Simulation of Gaussian variables Uniform Random number generators Pesudo-Random sequences Inversion Methods for distributions on R Cumulative distribution function Algorithm Sampling real valued random variables When F 1 (u) is not explicit Simulation of Gaussian variables Box-Müller Marsaglia-Bray Inversion method The rejection method (on R d ) Sampling under a curve Conditional distribution Algorithm Design parameters Example 1 : a normal generator based on the rejection algorithm Example 2 : Gaussian from Cauchy
23 Simulation of Gaussian variables Box-Müller Box-Müller algorithm (1/2) Set Z 1, Z 2, R, Θ such that Z 1 = R cos(θ) Z 2 = R sin(θ) The algorithm is based on the following property : (Z 1, Z 2) N (0, I 2) iff (R, Θ) are independent random variables, R has the distribution of the square root of a chi-square with two degrees of freedom or, equivalently, R is a Rayleigh distribution with parameter σ = 1 Θ has the uniform distribution [0, 2π].
24 Simulation of Gaussian variables Box-Müller Box-Müller algorithm (1/2) Set Z 1, Z 2, R, Θ such that Z 1 = R cos(θ) Z 2 = R sin(θ) The algorithm is based on the following property : (Z 1, Z 2) N (0, I 2) iff (R, Θ) are independent random variables, R has the distribution of the square root of a chi-square with two degrees of freedom or, equivalently, R is a Rayleigh distribution with parameter σ = 1 Θ has the uniform distribution [0, 2π]. Algorithm Generate a pair (U 1, U 2) of independent variables uniform on [0, 1]. R 2 log(u 1) V 2πU 2 Z 1 R cos(v ), Z 2 R sin(v ) Return Z 1, Z 2.
25 Simulation of Gaussian variables Box-Müller Box-Müller algorithm (2/2) Show that 1 R is a Rayleigh distribution. 2 V is a uniform distribution on [0, 2π]. 3 R and V are independent. and prove that (Z 1, Z 2) are independent r.v. with distribution N (0, 1). The main drawback of this method is that it requires the use of the sinus and cosinus functions, which can be expensive to evaluate. The following method is a way to overcome this problem.
26 Simulation of Gaussian variables Marsaglia-Bray Marsaglia-Bray algorithm (1/2) repeat Sample two independent uniform [0, 1] random variables U 1 and U 2 U 1 2U 1 1, U 2 2U 2 1 until U1 2 + U2 2 1 Y 2 log(u1 2 + U 2 2) U Z 1 Y 1 U, Z U 2 1 +U2 2 2 Y 2. U 2 1 +U2 2 Return Z 1, Z 2.
27 Simulation of Gaussian variables Marsaglia-Bray Marsaglia-Bray algorithm (2/2) Set T = U 1 U U 2 2 W = U U What is the distribution of (U 1, U 2) at the end of the repeat/until loop? 2 Show that T and W are independent, W is a uniform distribution on [0, 1] and T has the same distribution as cos(θ) when Θ U([0, 1]). 3 By using the result of the Box-Müller algorithm, show that Z 1, Z 2 are independent r.v. with normal distribution N (0, 1). 4 What is the acceptance probability of the repeat/until test? rmk : π/ This method substitutes the call to the sinus, cosinus functions for the following algorithm : instead of computing (cos(θ), sin(θ)) from Θ, a point on the unit ball is drawn, then projected on the unit sphere and its cartesian coordinates are then considered.
28 Simulation of Gaussian variables Inversion method Inversion Method The cumulative distribution function Φ of a standard Gaussian random variable satisfies Φ 1 (1 u) = Φ 1 (u), 0 < u < 1, and it suffices to approximate Φ 1 on [1/2, 1]. We may approximate Φ 1 using a rational function 3 Φ 1 n=0 an(u 1/2)2n+1 (u) 3 n=0 bn(u, 1/2)2n where 0.5 u 0.92 with a precision better than 10e 5 using the value computed in Springer Beasley (1977) For 0.92 u 1, Moro suggest to approximate Φ 1 using Φ 1 (u) 8 c n [log (log(1 u))] n, 0.92 u 1, n=0
29 The rejection method (on R d ) Uniform Random number generators Pesudo-Random sequences Inversion Methods for distributions on R Cumulative distribution function Algorithm Sampling real valued random variables When F 1 (u) is not explicit Simulation of Gaussian variables Box-Müller Marsaglia-Bray Inversion method The rejection method (on R d ) Sampling under a curve Conditional distribution Algorithm Design parameters Example 1 : a normal generator based on the rejection algorithm Example 2 : Gaussian from Cauchy
30 The rejection method (on R d ) Sampling under a curve Tools 1 and 2 : Sampling under the curve Theorem (i) Let X be a random vector with density g on R d, and let U be a uniform r.v. on [0, 1], independent of X. Let c > 0 be an arbitrary constant. Then, (X, cug(x)) is uniformly distributed on the set A = {(x, v) : x R d, 0 v cg(x)}. (ii) Conversely, if (X, V ) is a random vector in R d+1 uniformly distributed on this set A, then X has density g on R d Function x > cg (x) 0.8 Function x > g(x) Point with coordinates: (X, c g(x)) X
31 The rejection method (on R d ) Conditional distribution Tool 3 : First draw hitting a subset (1/2) Theorem Let X 1, X 2,... be a sequence of i.i.d. random vectors taking values in A R d, and let B A be a Borel set such that P(X 1 B) = p > 0. Let Y be the first r.v. X i taking values in A. Then (i) The distribution of Y is given by P(Y A) = P(X 1 A X 1 B) (ii) Let I = inf{i 1, X i B}. Then I is a geometric r.v. with parameter p i.e. P(I = k) = (1 p) k 1 p k 1.
32 The rejection method (on R d ) Conditional distribution Tool 3 : First draw hitting a subset (2/2) repeat Draw i.i.d. samples X k with distribution π until X k B Return X k. 1 If π is the uniform distribution on A, then the algorithm returns a sample with a uniform distribution on B. 2 The number of loops to return a value is a Geometric random variable with parameter p (and expected value 1/p).
33 The rejection method (on R d ) Algorithm Rejection Method for sampling under the distribution f on X R d below : either f, g are densities w.r.t. Lebesgue on R d ; or f, g are probability mass functions on an (at most) countable set X of R d The basic version of the rejection algorithm assumes the existence of a distribution g on R d and the knowledge of a constant c 1 such that f(x) cg(x) for all x X Algorithm : repeat Generate independently X with distribution g and U U([0, 1]) until U f(x) cg(x) Return X. Key result : The distribution of the r.v. X is f.
34 The rejection method (on R d ) Design parameters On the design parameters The three things we need before running the rejection algorithm are 1 a dominating distribution g 2 a simple method for generating random variates with distribution g 3 knowledge of c. Basically g must have heavier tails and sharper peaks than f. The dominating measure should be chosen with care! The number of iterations in order to return a value X is a Geometric r.v. with parameter 1/c : we should keep c as small as possible!
35 The rejection method (on R d ) Design parameters Development of good rejection algorithms Generally speaking, g is chosen from a class of classical densities. This class includes the uniform density, triangular densities and most densities that can be generated quickly by the inversion method. One generally starts with a family of dominating densities g (say a parametric family {g θ, θ Θ}) and chooses the density within this class for which c is the smallest. This approach sometimes leads to some difficult optimization problem.
36 The rejection method (on R d ) Example 1 : a normal generator based on the rejection algorithm Example 1 : A normal generator by rejection from the Laplace density The Laplace density is given by g(x) exp( x ) x R and the Normal density is given by f(x) = 1 2π exp( 0.5 x 2 ) x R Show that the following algorithm is a normal generator by rejection from the Laplace density repeat Generate independently an exponential random variate X (with parameter λ = 1) two r.v. U and V with distribution U([0, 1]). If U < 1/2 set X X. until V e 1/2 X e X2 /2 return X
37 The rejection method (on R d ) Example 2 : Gaussian from Cauchy Example 2 : Gaussian from Cauchy (1/2) We want to sample a standard Gaussian distribution on R. The family of dominating densities is the Cauchy family with scale parameter θ g θ (x) = θ 1 π θ 2 + x. 2 There is no need to consider a translation parameter as well because both f and g θ are unimodal with peak at 0. The optimal rejection constant is defined by It is given by c θ = c θ = sup x { 2π eθ eθ2 /2 f(x) g θ (x). θ < 2 θ π/2 θ 2 The function c θ has only one minimum at θ = 1, and the minimal value is c 1 = 2π/e.
38 The rejection method (on R d ) Example 2 : Gaussian from Cauchy An example : Gaussian from Cauchy (2/2) Show that the following algorithm is a normal generator by rejection from the Cauchy distribution. set α e/2 repeat Sample independently the r.v. U and V with uniform distribution on [0, 1]. Set X tan(πv ) until U α(1 + X 2 )e X2 /2 Return X The rejection constant is near 1.52 ; it is no match for most normal generators developed further on.
2 Random Variable Generation
2 Random Variable Generation Most Monte Carlo computations require, as a starting point, a sequence of i.i.d. random variables with given marginal distribution. We describe here some of the basic methods
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More information( ) ( ) Monte Carlo Methods Interested in. E f X = f x d x. Examples:
Monte Carlo Methods Interested in Examples: µ E f X = f x d x Type I error rate of a hypothesis test Mean width of a confidence interval procedure Evaluating a likelihood Finding posterior mean and variance
More informationMultivariate Simulations
Multivariate Simulations Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Prague, Czech Republic 21.11.2011 Katarína Starinská Multivariate Simulations 1 /
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca October 22nd, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationLecture Notes 1 Probability and Random Variables. Conditional Probability and Independence. Functions of a Random Variable
Lecture Notes 1 Probability and Random Variables Probability Spaces Conditional Probability and Independence Random Variables Functions of a Random Variable Generation of a Random Variable Jointly Distributed
More information2905 Queueing Theory and Simulation PART IV: SIMULATION
2905 Queueing Theory and Simulation PART IV: SIMULATION 22 Random Numbers A fundamental step in a simulation study is the generation of random numbers, where a random number represents the value of a random
More informationIndependent Events. Two events are independent if knowing that one occurs does not change the probability of the other occurring
Independent Events Two events are independent if knowing that one occurs does not change the probability of the other occurring Conditional probability is denoted P(A B), which is defined to be: P(A and
More informationLecture Notes 1 Probability and Random Variables. Conditional Probability and Independence. Functions of a Random Variable
Lecture Notes 1 Probability and Random Variables Probability Spaces Conditional Probability and Independence Random Variables Functions of a Random Variable Generation of a Random Variable Jointly Distributed
More informationStat 451 Lecture Notes Simulating Random Variables
Stat 451 Lecture Notes 05 12 Simulating Random Variables Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 6 in Givens & Hoeting, Chapter 22 in Lange, and Chapter 2 in Robert & Casella 2 Updated:
More informationMath Stochastic Processes & Simulation. Davar Khoshnevisan University of Utah
Math 5040 1 Stochastic Processes & Simulation Davar Khoshnevisan University of Utah Module 1 Generation of Discrete Random Variables Just about every programming language and environment has a randomnumber
More informationNonuniform Random Variate Generation
Nonuniform Random Variate Generation 1 Suppose we have a generator of i.i.d. U(0, 1) r.v. s. We want to generate r.v. s from other distributions, such as normal, Weibull, Poisson, etc. Or other random
More informationRandom variable X is a mapping that maps each outcome s in the sample space to a unique real number x, x. X s. Real Line
Random Variable Random variable is a mapping that maps each outcome s in the sample space to a unique real number,. s s : outcome Sample Space Real Line Eamples Toss a coin. Define the random variable
More informationST5215: Advanced Statistical Theory
Department of Statistics & Applied Probability Wednesday, October 5, 2011 Lecture 13: Basic elements and notions in decision theory Basic elements X : a sample from a population P P Decision: an action
More informationStatistics for scientists and engineers
Statistics for scientists and engineers February 0, 006 Contents Introduction. Motivation - why study statistics?................................... Examples..................................................3
More informationI. ANALYSIS; PROBABILITY
ma414l1.tex Lecture 1. 12.1.2012 I. NLYSIS; PROBBILITY 1. Lebesgue Measure and Integral We recall Lebesgue measure (M411 Probability and Measure) λ: defined on intervals (a, b] by λ((a, b]) := b a (so
More informationS6880 #7. Generate Non-uniform Random Number #1
S6880 #7 Generate Non-uniform Random Number #1 Outline 1 Inversion Method Inversion Method Examples Application to Discrete Distributions Using Inversion Method 2 Composition Method Composition Method
More informationCh3. Generating Random Variates with Non-Uniform Distributions
ST4231, Semester I, 2003-2004 Ch3. Generating Random Variates with Non-Uniform Distributions This chapter mainly focuses on methods for generating non-uniform random numbers based on the built-in standard
More informationRandom processes and probability distributions. Phys 420/580 Lecture 20
Random processes and probability distributions Phys 420/580 Lecture 20 Random processes Many physical processes are random in character: e.g., nuclear decay (Poisson distributed event count) P (k, τ) =
More informationCourse: ESO-209 Home Work: 1 Instructor: Debasis Kundu
Home Work: 1 1. Describe the sample space when a coin is tossed (a) once, (b) three times, (c) n times, (d) an infinite number of times. 2. A coin is tossed until for the first time the same result appear
More informationChapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued
Chapter 3 sections Chapter 3 - continued 3.1 Random Variables and Discrete Distributions 3.2 Continuous Distributions 3.3 The Cumulative Distribution Function 3.4 Bivariate Distributions 3.5 Marginal Distributions
More informationMTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6
MTH739U/P: Topics in Scientific Computing Autumn 16 Week 6 4.5 Generic algorithms for non-uniform variates We have seen that sampling from a uniform distribution in [, 1] is a relatively straightforward
More informationChapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued
Chapter 3 sections 3.1 Random Variables and Discrete Distributions 3.2 Continuous Distributions 3.3 The Cumulative Distribution Function 3.4 Bivariate Distributions 3.5 Marginal Distributions 3.6 Conditional
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationGenerating the Sample
STAT 80: Mathematical Statistics Monte Carlo Suppose you are given random variables X,..., X n whose joint density f (or distribution) is specified and a statistic T (X,..., X n ) whose distribution you
More informationDUBLIN CITY UNIVERSITY
DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr
More informationLecture 3: Random variables, distributions, and transformations
Lecture 3: Random variables, distributions, and transformations Definition 1.4.1. A random variable X is a function from S into a subset of R such that for any Borel set B R {X B} = {ω S : X(ω) B} is an
More informationFundamental Tools - Probability Theory II
Fundamental Tools - Probability Theory II MSc Financial Mathematics The University of Warwick September 29, 2015 MSc Financial Mathematics Fundamental Tools - Probability Theory II 1 / 22 Measurable random
More information2. Variance and Covariance: We will now derive some classic properties of variance and covariance. Assume real-valued random variables X and Y.
CS450 Final Review Problems Fall 08 Solutions or worked answers provided Problems -6 are based on the midterm review Identical problems are marked recap] Please consult previous recitations and textbook
More informationProbability Theory and Statistics. Peter Jochumzen
Probability Theory and Statistics Peter Jochumzen April 18, 2016 Contents 1 Probability Theory And Statistics 3 1.1 Experiment, Outcome and Event................................ 3 1.2 Probability............................................
More informationElementary Probability. Exam Number 38119
Elementary Probability Exam Number 38119 2 1. Introduction Consider any experiment whose result is unknown, for example throwing a coin, the daily number of customers in a supermarket or the duration of
More informationLecture 2 : CS6205 Advanced Modeling and Simulation
Lecture 2 : CS6205 Advanced Modeling and Simulation Lee Hwee Kuan 21 Aug. 2013 For the purpose of learning stochastic simulations for the first time. We shall only consider probabilities on finite discrete
More informationTheorem 2.1 (Caratheodory). A (countably additive) probability measure on a field has an extension. n=1
Chapter 2 Probability measures 1. Existence Theorem 2.1 (Caratheodory). A (countably additive) probability measure on a field has an extension to the generated σ-field Proof of Theorem 2.1. Let F 0 be
More informationGenerating Random Variables
Generating Random Variables Christian Robert Université Paris Dauphine and CREST, INSEE George Casella University of Florida Keywords and Phrases: Random Number Generator, Probability Integral Transform,
More informationSTA205 Probability: Week 8 R. Wolpert
INFINITE COIN-TOSS AND THE LAWS OF LARGE NUMBERS The traditional interpretation of the probability of an event E is its asymptotic frequency: the limit as n of the fraction of n repeated, similar, and
More informationDS-GA 1002 Lecture notes 2 Fall Random variables
DS-GA 12 Lecture notes 2 Fall 216 1 Introduction Random variables Random variables are a fundamental tool in probabilistic modeling. They allow us to model numerical quantities that are uncertain: the
More informationSTAT2201. Analysis of Engineering & Scientific Data. Unit 3
STAT2201 Analysis of Engineering & Scientific Data Unit 3 Slava Vaisman The University of Queensland School of Mathematics and Physics What we learned in Unit 2 (1) We defined a sample space of a random
More informationProbability and Measure
Probability and Measure Robert L. Wolpert Institute of Statistics and Decision Sciences Duke University, Durham, NC, USA Convergence of Random Variables 1. Convergence Concepts 1.1. Convergence of Real
More information1 Sequences of events and their limits
O.H. Probability II (MATH 2647 M15 1 Sequences of events and their limits 1.1 Monotone sequences of events Sequences of events arise naturally when a probabilistic experiment is repeated many times. For
More informationSample Spaces, Random Variables
Sample Spaces, Random Variables Moulinath Banerjee University of Michigan August 3, 22 Probabilities In talking about probabilities, the fundamental object is Ω, the sample space. (elements) in Ω are denoted
More information40.530: Statistics. Professor Chen Zehua. Singapore University of Design and Technology
Singapore University of Design and Technology Lecture 9: Hypothesis testing, uniformly most powerful tests. The Neyman-Pearson framework Let P be the family of distributions of concern. The Neyman-Pearson
More informationSo we will instead use the Jacobian method for inferring the PDF of functionally related random variables; see Bertsekas & Tsitsiklis Sec. 4.1.
2011 Page 1 Simulating Gaussian Random Variables Monday, February 14, 2011 2:00 PM Readings: Kloeden and Platen Sec. 1.3 Why does the Box Muller method work? How was it derived? The basic idea involves
More informationMonte Carlo and cold gases. Lode Pollet.
Monte Carlo and cold gases Lode Pollet lpollet@physics.harvard.edu 1 Outline Classical Monte Carlo The Monte Carlo trick Markov chains Metropolis algorithm Ising model critical slowing down Quantum Monte
More informationGenerating Random Variates 2 (Chapter 8, Law)
B. Maddah ENMG 6 Simulation /5/08 Generating Random Variates (Chapter 8, Law) Generating random variates from U(a, b) Recall that a random X which is uniformly distributed on interval [a, b], X ~ U(a,
More informationBrief Review of Probability
Maura Department of Economics and Finance Università Tor Vergata Outline 1 Distribution Functions Quantiles and Modes of a Distribution 2 Example 3 Example 4 Distributions Outline Distribution Functions
More informationChapter 2 Continuous Distributions
Chapter Continuous Distributions Continuous random variables For a continuous random variable X the probability distribution is described by the probability density function f(x), which has the following
More informationChapter 1. Sets and probability. 1.3 Probability space
Random processes - Chapter 1. Sets and probability 1 Random processes Chapter 1. Sets and probability 1.3 Probability space 1.3 Probability space Random processes - Chapter 1. Sets and probability 2 Probability
More information1 Probability theory. 2 Random variables and probability theory.
Probability theory Here we summarize some of the probability theory we need. If this is totally unfamiliar to you, you should look at one of the sources given in the readings. In essence, for the major
More informationMonte Carlo Integration II & Sampling from PDFs
Monte Carlo Integration II & Sampling from PDFs CS295, Spring 2017 Shuang Zhao Computer Science Department University of California, Irvine CS295, Spring 2017 Shuang Zhao 1 Last Lecture Direct illumination
More informationRandom variable X is a mapping that maps each outcome s in the sample space to a unique real number x, < x <. ( ) X s. Real Line
Random Variable Random variable is a mapping that maps each outcome s in the sample space to a unique real number, <
More information2 Functions of random variables
2 Functions of random variables A basic statistical model for sample data is a collection of random variables X 1,..., X n. The data are summarised in terms of certain sample statistics, calculated as
More informationConditional Distributions
Conditional Distributions The goal is to provide a general definition of the conditional distribution of Y given X, when (X, Y ) are jointly distributed. Let F be a distribution function on R. Let G(,
More information1 Review of Probability
1 Review of Probability Random variables are denoted by X, Y, Z, etc. The cumulative distribution function (c.d.f.) of a random variable X is denoted by F (x) = P (X x), < x
More informationGeneration from simple discrete distributions
S-38.3148 Simulation of data networks / Generation of random variables 1(18) Generation from simple discrete distributions Note! This is just a more clear and readable version of the same slide that was
More informationSystems Simulation Chapter 7: Random-Number Generation
Systems Simulation Chapter 7: Random-Number Generation Fatih Cavdur fatihcavdur@uludag.edu.tr April 22, 2014 Introduction Introduction Random Numbers (RNs) are a necessary basic ingredient in the simulation
More informationfunctions Poisson distribution Normal distribution Arbitrary functions
Physics 433: Computational Physics Lecture 6 Random number distributions Generation of random numbers of various distribuition functions Normal distribution Poisson distribution Arbitrary functions Random
More informationLecture 1: Random number generation, permutation test, and the bootstrap. August 25, 2016
Lecture 1: Random number generation, permutation test, and the bootstrap August 25, 2016 Statistical simulation 1/21 Statistical simulation (Monte Carlo) is an important part of statistical method research.
More informationPart IA Probability. Definitions. Based on lectures by R. Weber Notes taken by Dexter Chua. Lent 2015
Part IA Probability Definitions Based on lectures by R. Weber Notes taken by Dexter Chua Lent 2015 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures.
More information1 Stat 605. Homework I. Due Feb. 1, 2011
The first part is homework which you need to turn in. The second part is exercises that will not be graded, but you need to turn it in together with the take-home final exam. 1 Stat 605. Homework I. Due
More information17. Convergence of Random Variables
7. Convergence of Random Variables In elementary mathematics courses (such as Calculus) one speaks of the convergence of functions: f n : R R, then lim f n = f if lim f n (x) = f(x) for all x in R. This
More informationLecture 21: Convergence of transformations and generating a random variable
Lecture 21: Convergence of transformations and generating a random variable If Z n converges to Z in some sense, we often need to check whether h(z n ) converges to h(z ) in the same sense. Continuous
More informationMathematical Statistics 1 Math A 6330
Mathematical Statistics 1 Math A 6330 Chapter 2 Transformations and Expectations Mohamed I. Riffi Department of Mathematics Islamic University of Gaza September 14, 2015 Outline 1 Distributions of Functions
More informationLecture Quantitative Finance Spring Term 2015
on bivariate Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 07: April 2, 2015 1 / 54 Outline on bivariate 1 2 bivariate 3 Distribution 4 5 6 7 8 Comments and conclusions
More informationSimulation. Where real stuff starts
1 Simulation Where real stuff starts ToC 1. What is a simulation? 2. Accuracy of output 3. Random Number Generators 4. How to sample 5. Monte Carlo 6. Bootstrap 2 1. What is a simulation? 3 What is a simulation?
More informationNon-parametric Inference and Resampling
Non-parametric Inference and Resampling Exercises by David Wozabal (Last update. Juni 010) 1 Basic Facts about Rank and Order Statistics 1.1 10 students were asked about the amount of time they spend surfing
More informationStatistics (1): Estimation
Statistics (1): Estimation Marco Banterlé, Christian Robert and Judith Rousseau Practicals 2014-2015 L3, MIDO, Université Paris Dauphine 1 Table des matières 1 Random variables, probability, expectation
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationLecture 15: MCMC Sanjeev Arora Elad Hazan. COS 402 Machine Learning and Artificial Intelligence Fall 2016
Lecture 15: MCMC Sanjeev Arora Elad Hazan COS 402 Machine Learning and Artificial Intelligence Fall 2016 Course progress Learning from examples Definition + fundamental theorem of statistical learning,
More informationLecture 20. Randomness and Monte Carlo. J. Chaudhry. Department of Mathematics and Statistics University of New Mexico
Lecture 20 Randomness and Monte Carlo J. Chaudhry Department of Mathematics and Statistics University of New Mexico J. Chaudhry (UNM) CS 357 1 / 40 What we ll do: Random number generators Monte-Carlo integration
More informationChapter 4: Monte Carlo Methods. Paisan Nakmahachalasint
Chapter 4: Monte Carlo Methods Paisan Nakmahachalasint Introduction Monte Carlo Methods are a class of computational algorithms that rely on repeated random sampling to compute their results. Monte Carlo
More information2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is
More informationMath 321 Final Examination April 1995 Notation used in this exam: N. (1) S N (f,x) = f(t)e int dt e inx.
Math 321 Final Examination April 1995 Notation used in this exam: N 1 π (1) S N (f,x) = f(t)e int dt e inx. 2π n= N π (2) C(X, R) is the space of bounded real-valued functions on the metric space X, equipped
More information1 Acceptance-Rejection Method
Copyright c 2016 by Karl Sigman 1 Acceptance-Rejection Method As we already know, finding an explicit formula for F 1 (y), y [0, 1], for the cdf of a rv X we wish to generate, F (x) = P (X x), x R, is
More informationThe Arzelà-Ascoli Theorem
John Nachbar Washington University March 27, 2016 The Arzelà-Ascoli Theorem The Arzelà-Ascoli Theorem gives sufficient conditions for compactness in certain function spaces. Among other things, it helps
More informationNumerical Methods I Monte Carlo Methods
Numerical Methods I Monte Carlo Methods Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Dec. 9th, 2010 A. Donev (Courant Institute) Lecture
More informationScientific Computing: Monte Carlo
Scientific Computing: Monte Carlo Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 April 5th and 12th, 2012 A. Donev (Courant Institute)
More informationMinicourse on: Markov Chain Monte Carlo: Simulation Techniques in Statistics
Minicourse on: Markov Chain Monte Carlo: Simulation Techniques in Statistics Eric Slud, Statistics Program Lecture 1: Metropolis-Hastings Algorithm, plus background in Simulation and Markov Chains. Lecture
More information1: PROBABILITY REVIEW
1: PROBABILITY REVIEW Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 1: Probability Review 1 / 56 Outline We will review the following
More informationOrdered Sample Generation
Ordered Sample Generation Xuebo Yu November 20, 2010 1 Introduction There are numerous distributional problems involving order statistics that can not be treated analytically and need to simulated through
More informationNonlife Actuarial Models. Chapter 14 Basic Monte Carlo Methods
Nonlife Actuarial Models Chapter 14 Basic Monte Carlo Methods Learning Objectives 1. Generation of uniform random numbers, mixed congruential method 2. Low discrepancy sequence 3. Inversion transformation
More informationGenerating Random Variables
Generating Random Variables These slides are created by Dr. Yih Huang of George Mason University. Students registered in Dr. Huang's courses at GMU can make a single machine-readable copy and print a single
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 8 10/1/2008 CONTINUOUS RANDOM VARIABLES
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 8 10/1/2008 CONTINUOUS RANDOM VARIABLES Contents 1. Continuous random variables 2. Examples 3. Expected values 4. Joint distributions
More informationP-adic Functions - Part 1
P-adic Functions - Part 1 Nicolae Ciocan 22.11.2011 1 Locally constant functions Motivation: Another big difference between p-adic analysis and real analysis is the existence of nontrivial locally constant
More information15-388/688 - Practical Data Science: Basic probability. J. Zico Kolter Carnegie Mellon University Spring 2018
15-388/688 - Practical Data Science: Basic probability J. Zico Kolter Carnegie Mellon University Spring 2018 1 Announcements Logistics of next few lectures Final project released, proposals/groups due
More informationThree examples of a Practical Exact Markov Chain Sampling
Three examples of a Practical Exact Markov Chain Sampling Zdravko Botev November 2007 Abstract We present three examples of exact sampling from complex multidimensional densities using Markov Chain theory
More informationRandom Variate Generation
Random Variate Generation 28-1 Overview 1. Inverse transformation 2. Rejection 3. Composition 4. Convolution 5. Characterization 28-2 Random-Variate Generation General Techniques Only a few techniques
More informationSimulation. Where real stuff starts
Simulation Where real stuff starts March 2019 1 ToC 1. What is a simulation? 2. Accuracy of output 3. Random Number Generators 4. How to sample 5. Monte Carlo 6. Bootstrap 2 1. What is a simulation? 3
More informationPCMI Introduction to Random Matrix Theory Handout # REVIEW OF PROBABILITY THEORY. Chapter 1 - Events and Their Probabilities
PCMI 207 - Introduction to Random Matrix Theory Handout #2 06.27.207 REVIEW OF PROBABILITY THEORY Chapter - Events and Their Probabilities.. Events as Sets Definition (σ-field). A collection F of subsets
More informationMarkov Chain Monte Carlo (MCMC)
Markov Chain Monte Carlo (MCMC Dependent Sampling Suppose we wish to sample from a density π, and we can evaluate π as a function but have no means to directly generate a sample. Rejection sampling can
More informationLESSON RELATIONS & FUNCTION THEORY
2 Definitions LESSON RELATIONS & FUNCTION THEORY Ordered Pair Ordered pair of elements taken from any two sets P and Q is a pair of elements written in small brackets and grouped together in a particular
More informationMAT 271E Probability and Statistics
MAT 71E Probability and Statistics Spring 013 Instructor : Class Meets : Office Hours : Textbook : Supp. Text : İlker Bayram EEB 1103 ibayram@itu.edu.tr 13.30 1.30, Wednesday EEB 5303 10.00 1.00, Wednesday
More informationMATH/STAT 235A Probability Theory Lecture Notes, Fall 2013
MATH/STAT 235A Probability Theory Lecture Notes, Fall 2013 Dan Romik Department of Mathematics, UC Davis December 30, 2013 Contents Chapter 1: Introduction 6 1.1 What is probability theory?...........................
More informationMFM Practitioner Module: Quantitative Risk Management. John Dodson. September 23, 2015
MFM Practitioner Module: Quantitative Risk Management September 23, 2015 Mixtures Mixtures Mixtures Definitions For our purposes, A random variable is a quantity whose value is not known to us right now
More informationProblem Set 1 Sept, 14
EE6: Random Processes in Systems Lecturer: Jean C. Walrand Problem Set Sept, 4 Fall 06 GSI: Assane Gueye This problem set essentially reviews notions of conditional expectation, conditional distribution,
More informationA GENERAL METHOD FOR PRODUCING RANDOM V ARIABtES IN A COMPUTER. George Marsaglia. Boeing Scientific Research Laboratories Seattle, Washington
A GENERAL METHOD FOR PRODUCING RANDOM V ARIABtES IN A COMPUTER George Marsaglia Boeing Scientific Research Laboratories Seattle, Washington SUMMARY Many random variables can be approximated quite closely
More informationMonte Carlo Methods. Jesús Fernández-Villaverde University of Pennsylvania
Monte Carlo Methods Jesús Fernández-Villaverde University of Pennsylvania 1 Why Monte Carlo? From previous chapter, we want to compute: 1. Posterior distribution: π ³ θ Y T,i = f(y T θ,i)π (θ i) R Θ i
More informationWeak convergence. Amsterdam, 13 November Leiden University. Limit theorems. Shota Gugushvili. Generalities. Criteria
Weak Leiden University Amsterdam, 13 November 2013 Outline 1 2 3 4 5 6 7 Definition Definition Let µ, µ 1, µ 2,... be probability measures on (R, B). It is said that µ n converges weakly to µ, and we then
More informationTopology. Xiaolong Han. Department of Mathematics, California State University, Northridge, CA 91330, USA address:
Topology Xiaolong Han Department of Mathematics, California State University, Northridge, CA 91330, USA E-mail address: Xiaolong.Han@csun.edu Remark. You are entitled to a reward of 1 point toward a homework
More informationStochastic Convergence, Delta Method & Moment Estimators
Stochastic Convergence, Delta Method & Moment Estimators Seminar on Asymptotic Statistics Daniel Hoffmann University of Kaiserslautern Department of Mathematics February 13, 2015 Daniel Hoffmann (TU KL)
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More information