Output Analysis and Run-Length Control
|
|
- Logan Goodman
- 6 years ago
- Views:
Transcription
1 IEOR E4703: Mote Carlo Simulatio Columbia Uiversity c 2017 by Marti Haugh Output Aalysis ad Ru-Legth Cotrol I these otes we describe how the Cetral Limit Theorem ca be used to costruct approximate (1 α% cofidece itervals for the quatity, θ, we are tryig to estimate. We also describe methods to estimate the umber of samples that are required to achieve a give cofidece level ad we ed with a discussio of the bootstrap method for performig output aalysis. 1 Output Aalysis Recall the simulatio framework that we use whe we wat to estimate θ := E[h(X where X R. We first simulate X 1,..., X IID ad the set θ = h(x h(x The Strog Law of Large Numbers (SLLN the implies θ θ as w.p. 1. But at this poit we do t kow how large should be so that we ca have cofidece i θ as a estimator of θ. Put aother way, for a fixed value of, what ca we say about the quality of θ? We will ow aswer this questio ad to simplify our otatio we will take Y i := h(x i. 1.1 Cofidece Itervals Oe way to aswer this questio is to use a cofidece iterval. Suppose the that we wat to estimate θ ad we have a radom vector Y = (Y 1,..., Y whose distributio depeds o θ. The we seek L(Y ad U(Y such that P (L(Y θ U(Y = 1 α where 0 α 1 is a pre-specified umber. We the say that [L(Y, U(Y is a 100(1 α% cofidece iterval for θ. Note that [L(Y, U(Y is a radom iterval. However, oce we replace Y with a sample vector, y, the [L(y, U(y becomes a real iterval. We ow discuss the Chebyshev Iequality ad the Cetral Limit Theorem, both of which ca be used to costruct cofidece itervals. The Chebyshev Iequality Sice the Y i s are assumed to be IID we kow the variace of θ is give by Var( θ = 2 where 2 := Var(Y. Clearly a small value of Var( θ implies a more accurate estimate of θ ad this is ideed cofirmed by Chebyshev s Iequality which for ay k > 0 states that ( P θ θ k Var( θ k 2. (1 We ca see from (1 that a smaller value of Var( θ therefore improves our cofidece i θ. We could easily use Chebyshev s Iequality to costruct (how? cofidece itervals for θ but it is geerally very coservative. Exercise 1 Why does Chebyshev s Iequality geerally lead to coservative cofidece itervals? ( Istead, we will use the Cetral Limit Theorem to obtai better estimates of P θ θ k ad as a result, arrower cofidece itervals for θ.
2 Output Aalysis ad Ru-Legth Cotrol 2 The Cetral Limit Theorem The Cetral Limit Theorem is amog the most importat theorems i probability theory ad we state it here for coveiece with the symbol d deotig covergece i distributio. Theorem 1 (Cetral Limit Theorem Suppose Y 1,..., Y are IID ad E[Yi 2 <. The θ θ / d N(0, 1 as where θ = i=1 Y i/, θ := E[Y i ad 2 := Var(Y i. Note that we assume othig about the distributio of the Y i s other tha that E[Yi 2 <. If is sufficietly large i our simulatios, the we ca use the CLT to costruct cofidece itervals for θ := E[Y. We ow describe how to do this. 1.2 A Approximate 100(1 α% Cofidece Iterval for θ Let z 1 α/2 be the the (1 α/2 percetile poit of the N(0, 1 distributio so that P( z 1 α/2 Z z 1 α/2 = 1 α whe Z N(0, 1. Suppose ow that we have simulated IID samples, Y i, for i = 1,...,, ad that we wat to costruct a 100(1 α% CI for θ = E[Y. That is, we wat L(Y ad U(Y such that P (L(Y θ U(Y = 1 α. The CLT implies ( θ θ / is approximately N(0, 1 for large so we have P ( ( θ θ P z 1 α/2 z 1 α/2 ( z 1 α/2 θ θ z 1 α/2 1 α 1 α P ( θ z 1 α/2 θ θ + z 1 α/2 1 α. Our approximate 100(1 α% CI for θ is therefore give by [L(Y, U(Y = [ θ z 1 α/2, θ + z 1 α/2. (2 Recall that θ = (Y Y /, so L ad U are ideed fuctios of Y. There is still a problem, however, as we do ot usually kow 2. We resolve this issue by estimatig 2 with 2 i=1 = (Y i θ 2. 1 It is easy to show that 2 is a ubiased estimator of 2 ad that 2 2 w.p. 1 as. So ow we replace with i (2 to obtai [ [L(Y, U(Y = θ z 1 α/2, θ + z 1 α/2 (3 as our approximate 100(1 α% CI for θ whe is large.
3 Output Aalysis ad Ru-Legth Cotrol 3 Remark 1 Note that whe we obtai sample values of y = (y 1,..., y, the [L(y, U(y becomes a real iterval. The we ca o loger say (why ot? that P (θ [L(y, U(y = 1 α. Istead, we say that we are 100(1 α% cofidet that [L(y, U(y cotais θ. Furthermore, the smaller the value of U(y L(y, the more cofidece we will have i our estimate of θ. Example 1 (Pricig a Europea Call Optio Suppose we wat to estimate the price, C 0, of a call optio o a stock whose price process, S t, is a GBM(µ,. The relevat parameters are r =.05, T = 0.5 years, S 0 = $100, = 0.2 ad strike K = $110. The we kow that C 0 = E Q 0 [e rt max(s T K, 0 where we ca assume that S t GBM(r, uder the risk-eutral probability measure, Q. That is, we assume S T = S 0 exp ( (r 2 /2T + Z where Z N(0, T. Though we ca of course compute C 0 exactly, we ca also estimate C 0 usig Mote Carlo with (3 yieldig a approximate 100(1 α% CI for C 0 with Y i := e rt max(s (i T K, 0 deotig the ith discouted sample payoff of the optio. Based o = 100k samples, we obtai [15.16, as our approximate 95% CI for C 0. Properties of the Cofidece Iterval The width of the cofidece iterval is give by U L = 2 z 1 α/2 ad so the half-width the is (U L/2. The width clearly depeds o α, ad. However, almost surely as, ad is a costat. Therefore, for a fixed α, we eed to icrease if we are to decrease the width of the cofidece iterval. Ideed, sice U L 1, we ca see for example that we would eed to icrease by a factor of four i order to decrease the width of the cofidece iterval by oly a factor of two. 2 Ru-Legth Cotrol Up to this poit we have selected i advace ad the computed the approximate CI. The width of the CI is the a measure of the error i our estimator. Now we will do the reverse by first choosig some error criterio that we wat our estimator to satisfy, ad the choosig so that this criterio is satisfied. There are two types of error that we will cosider: 1. Absolute error, which is give by E a := θ θ ad 2. Relative error, which is give by E r := θ θ. θ Now we kow that θ θ w.p. 1 as so that E a ad E r both 0 as. (If θ = 0 the E r is ot defied. However, i practice ad so the errors will be o-zero. We specify the followig error criterio: Error Criterio: Give 0 α 1 ad ɛ 0, we wat P(E ɛ = 1 α. E is the error type we have specified, i.e., relative or absolute. The goal the is to choose so that the error criterio is approximately satisfied ad this is easily doe. Suppose, for example, that we wat to cotrol absolute error, E a. The, as we saw earlier, P ( θ z 1 α/2 θ θ + z 1 α/2 1 α.
4 Output Aalysis ad Ru-Legth Cotrol 4 ( This the implies P θ θ z 1 α/2 1 α, so i terms of E a we have P (E a z 1 α/2 1 α. If we the wat P (E a ɛ 1 α, it clearly suffices to choose such that = 2 z 2 1 α/2 ɛ 2. If we are workig with relative error, the a similar argumet implies that P (E r ɛ 1 α if There are still some problems, however: = 2 z 2 1 α/2 θ 2 ɛ Whe we are cotrollig E r, we eed to kow ad θ i advace. 2. Whe we are cotrollig E a, we eed to kow i advace. Of course we do ot usually kow or θ i advace. I fact, θ is what we are tryig to estimate! There are two methods we ca use to overcome this problem: the two-stage method ad the sequetial method, both of which we will ow describe. 2.1 The Two-Stage Procedure Suppose we wat to satisfy the coditio P(E a ɛ = 1 α so that we are tryig to cotrol the absolute error. The we saw earlier that we would like to set = 2 z 2 1 α/2 ɛ 2. Ufortuately, we do t kow 2 but we ca solve this problem by first doig a pilot simulatio to estimate it. The idea is to do a small umber, p, of iitial rus to estimate 2. We the use our estimate, 2, to compute a estimate,, of. Fially, we repeat the simulatio, but ow we use rus. We have the followig algorithm. Two-Stage Mote Carlo Simulatio for Estimatig E[h(X / Do pilot simulatio first / for i = 1 to p geerate X i ed for set θ = h(x i /p set 2 = (h(x i θ 2 /(p 1 set = 2 z 2 1 α/2 ɛ 2 / Now do mai simulatio / for i = 1 to geerate X i ed for set θ = h(x i / set 2 = (h(x i θ 2 /( 1 set 100(1 α % CI = [ θ z 1 α/2, θ + z 1 α/2
5 Output Aalysis ad Ru-Legth Cotrol 5 For this method to work, it is importat that θ ad 2 be sufficietly good estimates of θ ad 2. Therefore, it is importat to make p sufficietly large. I practice, we usually take p 50. We ca use a aalogous two-stage procedure if we wat to cotrol the relative error ad have P(E r ɛ = 1 α. 2.2 The Sequetial Procedure Suppose agai that we wish to satisfy the coditio P(E a ɛ = 1 α. The we saw earlier that we would like to set = 2 z1 α/2 2 ɛ 2. I cotrast to the pilot procedure, we do ot precompute durig the sequetial procedure. Istead, we cotiue to geerate samples util z 1 α/2 ɛ where is agai the estimate of based upo the first samples. It is importat that be sufficietly large so that θ ad, 2 are sufficietly good estimates of θ ad 2, respectively. As a result, we typically isist that 50 before we stop. Approximate cofidece itervals are the computed as usual. Questio: Have we allowed ay biases to creep i here? We have the followig algorithm: Sequetial Mote Carlo Simulatio for Estimatig E[h(X set check = 0, = 1 while (check = 0 geerate X set θ = h(x i / set 2 = (h(x i θ 2 /( 1 if ( p ad ɛ else check = 1 ( z 1 α/2 = + 1 ed if ed while set 100(1 α % CI = [ θ z 1 α/2, θ + z 1 α/2 I practice we do ot eed to store every value, h(x i for i = 1,...,, i order to update θ ad. Ideed, we ca update θ ad efficietly by observig that θ = θ 1 + h(x θ 1 ad ( 2 2 = ( θ 1 θ 2 1. If we wat to cotrol the relative error ad have P(E r ɛ = 1 α, the we would simulate samples util (z 1 α/2 θ ɛ.
6 Output Aalysis ad Ru-Legth Cotrol 6 3 Output Aalysis Usig the Bootstrap We ca view 1 our output aalysis problem as oe of estimatig MSE(F := E F [(g(y 1,..., Y θ(f 2 (4 where θ(f = E F [X, g(y 1,..., Y := Ȳ ad F deotes the CDF of Y. I that case, we saw i Sectio 1 how we could use the CLT to costruct approximate cofidece itervals for θ. While this is certaily the most commo cotext i which we ecouter (4, other situatios arise where the CLT caot be easily used to obtai a cofidece iterval for θ(f. For example, if θ(f = Var(Y or θ(f = E [Y Y α, the a alterative method of costructig a cofidece iterval for θ will be required. The bootstrap method provides such a alterative ad i order to describe the method we will assume our problem is to estimate MSE(F as i (4. To begi with, recall that the empirical distributio, F e, is defied to be the CDF of the distributio that places a weight of 1/ o each of the simulated values Y 1,..., Y. The empirical CDF therefore satisfies i=1 F e (y = 1 {Y i y} ad for large it ca be show (ad should be ituitively clear that F e should be 2 a good approximatio to F. Therefore, as log as θ is sufficietly well-behaved, i.e. a cotiuous fuctio of F, the for sufficietly large we should have MSE(F MSE(F e = E Fe [(g(y 1,..., Y θ(f e 2. (5 The quatity MSE(F e is kow as the bootstrap approximatio to MSE(F ad is easy to estimate via simulatio as we shall see below. But first, however, we will cosider a example where MSE(F e ca be computed exactly. Ideed the bootstrap is ot required i this case but it is oetheless istructive to see the calculatios writte out explicitly. Example 2 (Applyig the Bootstrap to the Sample Mea Suppose we wish to estimate θ(f = E F [Y via the estimator θ = g(y 1,..., Y := Ȳ. As oted above, the bootstrap is ot ecessary i this case as we ca apply the CLT directly as i Sectio 1 to obtai cofidece itervals for θ [ (Ȳ 2 or equivaletly, we ca estimate the mea-squared error E θ = 2 / with / 2 = i=1 (y i ȳ 2 /(( 1. Lettig ȳ deote the mea of the observed, i.e. simulated, data-poits y 1,..., y, we obtai that the bootstrap estimator is give by MSE(F e = E Fe [ ( i=1 Y i ( i=1 = Var Y i Fe = Var F e (Y i=1 = (y i ȳ 2 2 ȳ where (6 follows sice E Fe [Y = ȳ, ad (7 follows sice the Y i s are IID F e. We therefore see that the bootstrap approximatio to the MSE is almost idetical to our usual estimator, 2 /. 2 (6 (7 1 We follow Sheldo M. Ross s Simulatio i our developmet of the bootstrap here. 2 Ideed it ca be show that F e(y coverges to F (y uiformly i y w.p. 1 as.
7 Output Aalysis ad Ru-Legth Cotrol 7 I cotrast to Example 2, we caot usually compute MSE(F e explicitly, but as it s a expectatio we ca easily use Mote-Carlo to estimate it. I this case we eed to simulate from F e which is easy to do ad so we obtai the followig bootstrap algorithm for estimatig MSE(F. Bootstrap Simulatio Algorithm for Estimatig MSE(F for i = 1 to B geerate Y 1,..., Y IID from F e set θ i b = g(y 1,..., Y set Zi [ θb b = i θ(f e ed for set MSE(F B = b=1 Z(b /B 2 The Zi b s (or equivaletly the θ i b s are the bootstrap samples ad a value of B = 100 is ofte sufficiet to obtai a sufficietly accurate estimate. I the ext example we apply the bootstrap approach i a historical simulatio cotext where we have real data observatios as opposed to simulated data. (The disadvatage with historical simulatio is that we typically have o cotrol over. Example 3 (Estimatig the Miimum Variace Portfolio Suppose we wish to ivest a fixed sum of moey i two fiacial assets, X ad Z say, that yield radom returs of R x ad R z, respectively. We ivest a fractio θ of our wealth i X, ad the remaiig 1 θ i Z. The goal is to choose θ to miimize the total variace, Var(θR x + (1 θr z, of our ivestmet retur. It is easy to see that the miimizig θ is give by z 2 xz θ = x 2 + z 2 (8 2 xz where 2 x = Var(R x, 2 z = Var(R z ad xz = Cov(R x, R z. I practice, we do ot kow these quatities ad therefore have to estimate them from historical data. We therefore obtai θ = 2 z xz 2 x + 2 z 2 xz. (9 as our estimator of the( miimum variace portfolio with x, 2 z 2 ad xz estimated from historical retur data Y 1,..., Y with Y i := R x (i, R z (i the joit retur i period i. We would like to kow how good a estimator θ is. More specifically, what is the (mea-squared error whe we use θ? We ca aswer this questio usig the bootstrap with θ(f := θ ad g(y 1,..., Y = θ the estimator give by (9. Exercise 2 Provide pseudo-code for estimatig MSE( θ := MSE(F, i Example 3. Exercise 3 Cosider the problem of estimatig θ(f = E [Y Y β for some fixed costat, β. Explai how you would use the bootstrap to estimate MSE(F i this case give Mote-Carlo samples Y 1,..., Y. 3.1 Costructig Bootstrap Cofidece Itervals The bootstrap method is also widely used to costruct cofidece itervals ad here we will cosider the so-called basic bootstrap iterval. Cosider our bootstrap samples θ 1, b..., θ B b ad suppose we wat a 1 α cofidece iterval for θ = θ(f. Let q l ad q u be the α/2 lower- ad upper-sample quatiles, respectively, of the bootstrap samples. The the fractio of bootstrap samples satisfyig q l θ b q u (10
8 Output Aalysis ad Ru-Legth Cotrol 8 is 1 α. But (10 is equivalet to θ q u θ θ b θ q l (11 where θ = g(y 1,..., y is our estimate of θ computed usig the origial data-set. This implies θ q u ad θ q l are the lower ad upper quatiles for θ θ b. The basic bootstrap assumes they are also the quatiles for θ θ. This makes sese ituitively ad ca be justified mathematically as ad if θ is a cotiuous fuctio of F. It therefore follows that θ q u θ θ θ q l (12 will occur i approximately i a fractio 1 α of samples. Addig θ across (12 yields a approximate (1 α% CI for θ of (2 θ q u, 2 θ q l.
1 Introduction to reducing variance in Monte Carlo simulations
Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Output Analysis for Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Output Analysis
More informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More informationDiscrete Mathematics for CS Spring 2008 David Wagner Note 22
CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 22 I.I.D. Radom Variables Estimatig the bias of a coi Questio: We wat to estimate the proportio p of Democrats i the US populatio, by takig
More informationSimulation. Two Rule For Inverting A Distribution Function
Simulatio Two Rule For Ivertig A Distributio Fuctio Rule 1. If F(x) = u is costat o a iterval [x 1, x 2 ), the the uiform value u is mapped oto x 2 through the iversio process. Rule 2. If there is a jump
More informationStatistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.
Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized
More informationLecture 33: Bootstrap
Lecture 33: ootstrap Motivatio To evaluate ad compare differet estimators, we eed cosistet estimators of variaces or asymptotic variaces of estimators. This is also importat for hypothesis testig ad cofidece
More informationEstimation for Complete Data
Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit Theorems Throughout this sectio we will assume a probability space (, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationDS 100: Principles and Techniques of Data Science Date: April 13, Discussion #10
DS 00: Priciples ad Techiques of Data Sciece Date: April 3, 208 Name: Hypothesis Testig Discussio #0. Defie these terms below as they relate to hypothesis testig. a) Data Geeratio Model: Solutio: A set
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit theorems Throughout this sectio we will assume a probability space (Ω, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationSequences and Series of Functions
Chapter 6 Sequeces ad Series of Fuctios 6.1. Covergece of a Sequece of Fuctios Poitwise Covergece. Defiitio 6.1. Let, for each N, fuctio f : A R be defied. If, for each x A, the sequece (f (x)) coverges
More informationFall 2013 MTH431/531 Real analysis Section Notes
Fall 013 MTH431/531 Real aalysis Sectio 8.1-8. Notes Yi Su 013.11.1 1. Defiitio of uiform covergece. We look at a sequece of fuctios f (x) ad study the coverget property. Notice we have two parameters
More informationConvergence of random variables. (telegram style notes) P.J.C. Spreij
Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space
More informationChapter 3. Strong convergence. 3.1 Definition of almost sure convergence
Chapter 3 Strog covergece As poited out i the Chapter 2, there are multiple ways to defie the otio of covergece of a sequece of radom variables. That chapter defied covergece i probability, covergece i
More informationLecture 12: September 27
36-705: Itermediate Statistics Fall 207 Lecturer: Siva Balakrisha Lecture 2: September 27 Today we will discuss sufficiecy i more detail ad the begi to discuss some geeral strategies for costructig estimators.
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationFrequentist Inference
Frequetist Iferece The topics of the ext three sectios are useful applicatios of the Cetral Limit Theorem. Without kowig aythig about the uderlyig distributio of a sequece of radom variables {X i }, for
More informationLecture 19: Convergence
Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may
More informationStatistics 511 Additional Materials
Cofidece Itervals o mu Statistics 511 Additioal Materials This topic officially moves us from probability to statistics. We begi to discuss makig ifereces about the populatio. Oe way to differetiate probability
More informationStat 421-SP2012 Interval Estimation Section
Stat 41-SP01 Iterval Estimatio Sectio 11.1-11. We ow uderstad (Chapter 10) how to fid poit estimators of a ukow parameter. o However, a poit estimate does ot provide ay iformatio about the ucertaity (possible
More informationStochastic Simulation
Stochastic Simulatio 1 Itroductio Readig Assigmet: Read Chapter 1 of text. We shall itroduce may of the key issues to be discussed i this course via a couple of model problems. Model Problem 1 (Jackso
More informationEcon 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1.
Eco 325/327 Notes o Sample Mea, Sample Proportio, Cetral Limit Theorem, Chi-square Distributio, Studet s t distributio 1 Sample Mea By Hiro Kasahara We cosider a radom sample from a populatio. Defiitio
More information32 estimating the cumulative distribution function
32 estimatig the cumulative distributio fuctio 4.6 types of cofidece itervals/bads Let F be a class of distributio fuctios F ad let θ be some quatity of iterest, such as the mea of F or the whole fuctio
More informationSince X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain
Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the
More information1.010 Uncertainty in Engineering Fall 2008
MIT OpeCourseWare http://ocw.mit.edu.00 Ucertaity i Egieerig Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu.terms. .00 - Brief Notes # 9 Poit ad Iterval
More informationLecture 8: Convergence of transformations and law of large numbers
Lecture 8: Covergece of trasformatios ad law of large umbers Trasformatio ad covergece Trasformatio is a importat tool i statistics. If X coverges to X i some sese, we ofte eed to check whether g(x ) coverges
More informationMachine Learning Theory Tübingen University, WS 2016/2017 Lecture 12
Machie Learig Theory Tübige Uiversity, WS 06/07 Lecture Tolstikhi Ilya Abstract I this lecture we derive risk bouds for kerel methods. We will start by showig that Soft Margi kerel SVM correspods to miimizig
More information1 of 7 7/16/2009 6:06 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 6. Order Statistics Defiitios Suppose agai that we have a basic radom experimet, ad that X is a real-valued radom variable
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationA Question. Output Analysis. Example. What Are We Doing Wrong? Result from throwing a die. Let X be the random variable
A Questio Output Aalysis Let X be the radom variable Result from throwig a die 5.. Questio: What is E (X? Would you throw just oce ad take the result as your aswer? Itroductio to Simulatio WS/ - L 7 /
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationMachine Learning Brett Bernstein
Machie Learig Brett Berstei Week Lecture: Cocept Check Exercises Starred problems are optioal. Statistical Learig Theory. Suppose A = Y = R ad X is some other set. Furthermore, assume P X Y is a discrete
More informationThis exam contains 19 pages (including this cover page) and 10 questions. A Formulae sheet is provided with the exam.
Probability ad Statistics FS 07 Secod Sessio Exam 09.0.08 Time Limit: 80 Miutes Name: Studet ID: This exam cotais 9 pages (icludig this cover page) ad 0 questios. A Formulae sheet is provided with the
More informationDistribution of Random Samples & Limit theorems
STAT/MATH 395 A - PROBABILITY II UW Witer Quarter 2017 Néhémy Lim Distributio of Radom Samples & Limit theorems 1 Distributio of i.i.d. Samples Motivatig example. Assume that the goal of a study is to
More informationMachine Learning Brett Bernstein
Machie Learig Brett Berstei Week 2 Lecture: Cocept Check Exercises Starred problems are optioal. Excess Risk Decompositio 1. Let X = Y = {1, 2,..., 10}, A = {1,..., 10, 11} ad suppose the data distributio
More informationARIMA Models. Dan Saunders. y t = φy t 1 + ɛ t
ARIMA Models Da Sauders I will discuss models with a depedet variable y t, a potetially edogeous error term ɛ t, ad a exogeous error term η t, each with a subscript t deotig time. With just these three
More informationProduct measures, Tonelli s and Fubini s theorems For use in MAT3400/4400, autumn 2014 Nadia S. Larsen. Version of 13 October 2014.
Product measures, Toelli s ad Fubii s theorems For use i MAT3400/4400, autum 2014 Nadia S. Larse Versio of 13 October 2014. 1. Costructio of the product measure The purpose of these otes is to preset the
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationLecture Notes 15 Hypothesis Testing (Chapter 10)
1 Itroductio Lecture Notes 15 Hypothesis Testig Chapter 10) Let X 1,..., X p θ x). Suppose we we wat to kow if θ = θ 0 or ot, where θ 0 is a specific value of θ. For example, if we are flippig a coi, we
More informationNotes for Lecture 11
U.C. Berkeley CS78: Computatioal Complexity Hadout N Professor Luca Trevisa 3/4/008 Notes for Lecture Eigevalues, Expasio, ad Radom Walks As usual by ow, let G = (V, E) be a udirected d-regular graph with
More informationLecture 10 October Minimaxity and least favorable prior sequences
STATS 300A: Theory of Statistics Fall 205 Lecture 0 October 22 Lecturer: Lester Mackey Scribe: Brya He, Rahul Makhijai Warig: These otes may cotai factual ad/or typographic errors. 0. Miimaxity ad least
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationMATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4
MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.
More informationLecture 6 Simple alternatives and the Neyman-Pearson lemma
STATS 00: Itroductio to Statistical Iferece Autum 06 Lecture 6 Simple alteratives ad the Neyma-Pearso lemma Last lecture, we discussed a umber of ways to costruct test statistics for testig a simple ull
More informationNotes 5 : More on the a.s. convergence of sums
Notes 5 : More o the a.s. covergece of sums Math 733-734: Theory of Probability Lecturer: Sebastie Roch Refereces: Dur0, Sectios.5; Wil9, Sectio 4.7, Shi96, Sectio IV.4, Dur0, Sectio.. Radom series. Three-series
More informationLECTURE 14 NOTES. A sequence of α-level tests {ϕ n (x)} is consistent if
LECTURE 14 NOTES 1. Asymptotic power of tests. Defiitio 1.1. A sequece of -level tests {ϕ x)} is cosistet if β θ) := E θ [ ϕ x) ] 1 as, for ay θ Θ 1. Just like cosistecy of a sequece of estimators, Defiitio
More informationLecture 27. Capacity of additive Gaussian noise channel and the sphere packing bound
Lecture 7 Ageda for the lecture Gaussia chael with average power costraits Capacity of additive Gaussia oise chael ad the sphere packig boud 7. Additive Gaussia oise chael Up to this poit, we have bee
More information5. Likelihood Ratio Tests
1 of 5 7/29/2009 3:16 PM Virtual Laboratories > 9. Hy pothesis Testig > 1 2 3 4 5 6 7 5. Likelihood Ratio Tests Prelimiaries As usual, our startig poit is a radom experimet with a uderlyig sample space,
More information5. INEQUALITIES, LIMIT THEOREMS AND GEOMETRIC PROBABILITY
IA Probability Let Term 5 INEQUALITIES, LIMIT THEOREMS AND GEOMETRIC PROBABILITY 51 Iequalities Suppose that X 0 is a radom variable takig o-egative values ad that c > 0 is a costat The P X c E X, c is
More informationSOME THEORY AND PRACTICE OF STATISTICS by Howard G. Tucker
SOME THEORY AND PRACTICE OF STATISTICS by Howard G. Tucker CHAPTER 9. POINT ESTIMATION 9. Covergece i Probability. The bases of poit estimatio have already bee laid out i previous chapters. I chapter 5
More informationBig Picture. 5. Data, Estimates, and Models: quantifying the accuracy of estimates.
5. Data, Estimates, ad Models: quatifyig the accuracy of estimates. 5. Estimatig a Normal Mea 5.2 The Distributio of the Normal Sample Mea 5.3 Normal data, cofidece iterval for, kow 5.4 Normal data, cofidece
More informationLet us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.
Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,
More informationOutput Analysis (2, Chapters 10 &11 Law)
B. Maddah ENMG 6 Simulatio Output Aalysis (, Chapters 10 &11 Law) Comparig alterative system cofiguratio Sice the output of a simulatio is radom, the comparig differet systems via simulatio should be doe
More informationClases 7-8: Métodos de reducción de varianza en Monte Carlo *
Clases 7-8: Métodos de reducció de variaza e Mote Carlo * 9 de septiembre de 27 Ídice. Variace reductio 2. Atithetic variates 2 2.. Example: Uiform radom variables................ 3 2.2. Example: Tail
More informationEcon 325 Notes on Point Estimator and Confidence Interval 1 By Hiro Kasahara
Poit Estimator Eco 325 Notes o Poit Estimator ad Cofidece Iterval 1 By Hiro Kasahara Parameter, Estimator, ad Estimate The ormal probability desity fuctio is fully characterized by two costats: populatio
More informationENGI 4421 Confidence Intervals (Two Samples) Page 12-01
ENGI 44 Cofidece Itervals (Two Samples) Page -0 Two Sample Cofidece Iterval for a Differece i Populatio Meas [Navidi sectios 5.4-5.7; Devore chapter 9] From the cetral limit theorem, we kow that, for sufficietly
More informationECE 901 Lecture 12: Complexity Regularization and the Squared Loss
ECE 90 Lecture : Complexity Regularizatio ad the Squared Loss R. Nowak 5/7/009 I the previous lectures we made use of the Cheroff/Hoeffdig bouds for our aalysis of classifier errors. Hoeffdig s iequality
More informationAdvanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology
Advaced Aalysis Mi Ya Departmet of Mathematics Hog Kog Uiversity of Sciece ad Techology September 3, 009 Cotets Limit ad Cotiuity 7 Limit of Sequece 8 Defiitio 8 Property 3 3 Ifiity ad Ifiitesimal 8 4
More informationAdvanced Stochastic Processes.
Advaced Stochastic Processes. David Gamarik LECTURE 2 Radom variables ad measurable fuctios. Strog Law of Large Numbers (SLLN). Scary stuff cotiued... Outlie of Lecture Radom variables ad measurable fuctios.
More informationDiscrete Mathematics and Probability Theory Fall 2016 Walrand Probability: An Overview
CS 70 Discrete Mathematics ad Probability Theory Fall 2016 Walrad Probability: A Overview Probability is a fasciatig theory. It provides a precise, clea, ad useful model of ucertaity. The successes of
More informationProblem Set 4 Due Oct, 12
EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios
More informationLecture Chapter 6: Convergence of Random Sequences
ECE5: Aalysis of Radom Sigals Fall 6 Lecture Chapter 6: Covergece of Radom Sequeces Dr Salim El Rouayheb Scribe: Abhay Ashutosh Doel, Qibo Zhag, Peiwe Tia, Pegzhe Wag, Lu Liu Radom sequece Defiitio A ifiite
More information4. Partial Sums and the Central Limit Theorem
1 of 10 7/16/2009 6:05 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 4. Partial Sums ad the Cetral Limit Theorem The cetral limit theorem ad the law of large umbers are the two fudametal theorems
More informationSDS 321: Introduction to Probability and Statistics
SDS 321: Itroductio to Probability ad Statistics Lecture 23: Cotiuous radom variables- Iequalities, CLT Puramrita Sarkar Departmet of Statistics ad Data Sciece The Uiversity of Texas at Austi www.cs.cmu.edu/
More informationChapter 6 Principles of Data Reduction
Chapter 6 for BST 695: Special Topics i Statistical Theory. Kui Zhag, 0 Chapter 6 Priciples of Data Reductio Sectio 6. Itroductio Goal: To summarize or reduce the data X, X,, X to get iformatio about a
More informationDepartment of Mathematics
Departmet of Mathematics Ma 3/103 KC Border Itroductio to Probability ad Statistics Witer 2017 Lecture 19: Estimatio II Relevat textbook passages: Larse Marx [1]: Sectios 5.2 5.7 19.1 The method of momets
More informationMath 113, Calculus II Winter 2007 Final Exam Solutions
Math, Calculus II Witer 7 Fial Exam Solutios (5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute x x + dx The check your aswer usig the Evaluatio Theorem Solutio: I this
More informationMATH 472 / SPRING 2013 ASSIGNMENT 2: DUE FEBRUARY 4 FINALIZED
MATH 47 / SPRING 013 ASSIGNMENT : DUE FEBRUARY 4 FINALIZED Please iclude a cover sheet that provides a complete setece aswer to each the followig three questios: (a) I your opiio, what were the mai ideas
More informationIntroductory statistics
CM9S: Machie Learig for Bioiformatics Lecture - 03/3/06 Itroductory statistics Lecturer: Sriram Sakararama Scribe: Sriram Sakararama We will provide a overview of statistical iferece focussig o the key
More informationA statistical method to determine sample size to estimate characteristic value of soil parameters
A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig
More information17. Joint distributions of extreme order statistics Lehmann 5.1; Ferguson 15
17. Joit distributios of extreme order statistics Lehma 5.1; Ferguso 15 I Example 10., we derived the asymptotic distributio of the maximum from a radom sample from a uiform distributio. We did this usig
More informationParameter, Statistic and Random Samples
Parameter, Statistic ad Radom Samples A parameter is a umber that describes the populatio. It is a fixed umber, but i practice we do ot kow its value. A statistic is a fuctio of the sample data, i.e.,
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 6 9/23/2013. Brownian motion. Introduction
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/5.070J Fall 203 Lecture 6 9/23/203 Browia motio. Itroductio Cotet.. A heuristic costructio of a Browia motio from a radom walk. 2. Defiitio ad basic properties
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
More informationConfidence intervals summary Conservative and approximate confidence intervals for a binomial p Examples. MATH1005 Statistics. Lecture 24. M.
MATH1005 Statistics Lecture 24 M. Stewart School of Mathematics ad Statistics Uiversity of Sydey Outlie Cofidece itervals summary Coservative ad approximate cofidece itervals for a biomial p The aïve iterval
More informationECE 901 Lecture 14: Maximum Likelihood Estimation and Complexity Regularization
ECE 90 Lecture 4: Maximum Likelihood Estimatio ad Complexity Regularizatio R Nowak 5/7/009 Review : Maximum Likelihood Estimatio We have iid observatios draw from a ukow distributio Y i iid p θ, i,, where
More informationDirection: This test is worth 250 points. You are required to complete this test within 50 minutes.
Term Test October 3, 003 Name Math 56 Studet Number Directio: This test is worth 50 poits. You are required to complete this test withi 50 miutes. I order to receive full credit, aswer each problem completely
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 21 11/27/2013
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 21 11/27/2013 Fuctioal Law of Large Numbers. Costructio of the Wieer Measure Cotet. 1. Additioal techical results o weak covergece
More information7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals
7-1 Chapter 4 Part I. Samplig Distributios ad Cofidece Itervals 1 7- Sectio 1. Samplig Distributio 7-3 Usig Statistics Statistical Iferece: Predict ad forecast values of populatio parameters... Test hypotheses
More informationIntroduction to Machine Learning DIS10
CS 189 Fall 017 Itroductio to Machie Learig DIS10 1 Fu with Lagrage Multipliers (a) Miimize the fuctio such that f (x,y) = x + y x + y = 3. Solutio: The Lagragia is: L(x,y,λ) = x + y + λ(x + y 3) Takig
More informationLECTURE 8: ASYMPTOTICS I
LECTURE 8: ASYMPTOTICS I We are iterested i the properties of estimators as. Cosider a sequece of radom variables {, X 1}. N. M. Kiefer, Corell Uiversity, Ecoomics 60 1 Defiitio: (Weak covergece) A sequece
More informationThe standard deviation of the mean
Physics 6C Fall 20 The stadard deviatio of the mea These otes provide some clarificatio o the distictio betwee the stadard deviatio ad the stadard deviatio of the mea.. The sample mea ad variace Cosider
More informationIntroduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT
Itroductio to Extreme Value Theory Laures de Haa, ISM Japa, 202 Itroductio to Extreme Value Theory Laures de Haa Erasmus Uiversity Rotterdam, NL Uiversity of Lisbo, PT Itroductio to Extreme Value Theory
More informationEstimation of a population proportion March 23,
1 Social Studies 201 Notes for March 23, 2005 Estimatio of a populatio proportio Sectio 8.5, p. 521. For the most part, we have dealt with meas ad stadard deviatios this semester. This sectio of the otes
More informationAda Boost, Risk Bounds, Concentration Inequalities. 1 AdaBoost and Estimates of Conditional Probabilities
CS8B/Stat4B Sprig 008) Statistical Learig Theory Lecture: Ada Boost, Risk Bouds, Cocetratio Iequalities Lecturer: Peter Bartlett Scribe: Subhrasu Maji AdaBoost ad Estimates of Coditioal Probabilities We
More informationAn Introduction to Asymptotic Theory
A Itroductio to Asymptotic Theory Pig Yu School of Ecoomics ad Fiace The Uiversity of Hog Kog Pig Yu (HKU) Asymptotic Theory 1 / 20 Five Weapos i Asymptotic Theory Five Weapos i Asymptotic Theory Pig Yu
More informationLecture 11 and 12: Basic estimation theory
Lecture ad 2: Basic estimatio theory Sprig 202 - EE 94 Networked estimatio ad cotrol Prof. Kha March 2 202 I. MAXIMUM-LIKELIHOOD ESTIMATORS The maximum likelihood priciple is deceptively simple. Louis
More informationBinomial Distribution
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 1 2 3 4 5 6 7 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Overview Example: coi tossed three times Defiitio Formula Recall that a r.v. is discrete if there are either a fiite umber of possible
More informationRecall the study where we estimated the difference between mean systolic blood pressure levels of users of oral contraceptives and non-users, x - y.
Testig Statistical Hypotheses Recall the study where we estimated the differece betwee mea systolic blood pressure levels of users of oral cotraceptives ad o-users, x - y. Such studies are sometimes viewed
More informationRates of Convergence by Moduli of Continuity
Rates of Covergece by Moduli of Cotiuity Joh Duchi: Notes for Statistics 300b March, 017 1 Itroductio I this ote, we give a presetatio showig the importace, ad relatioship betwee, the modulis of cotiuity
More informationIntegrable Functions. { f n } is called a determining sequence for f. If f is integrable with respect to, then f d does exist as a finite real number
MATH 532 Itegrable Fuctios Dr. Neal, WKU We ow shall defie what it meas for a measurable fuctio to be itegrable, show that all itegral properties of simple fuctios still hold, ad the give some coditios
More informationElement sampling: Part 2
Chapter 4 Elemet samplig: Part 2 4.1 Itroductio We ow cosider uequal probability samplig desigs which is very popular i practice. I the uequal probability samplig, we ca improve the efficiecy of the resultig
More informationIntro to Learning Theory
Lecture 1, October 18, 2016 Itro to Learig Theory Ruth Urer 1 Machie Learig ad Learig Theory Comig soo 2 Formal Framework 21 Basic otios I our formal model for machie learig, the istaces to be classified
More informationApplication to Random Graphs
A Applicatio to Radom Graphs Brachig processes have a umber of iterestig ad importat applicatios. We shall cosider oe of the most famous of them, the Erdős-Réyi radom graph theory. 1 Defiitio A.1. Let
More informationST5215: Advanced Statistical Theory
ST525: Advaced Statistical Theory Departmet of Statistics & Applied Probability Tuesday, September 7, 2 ST525: Advaced Statistical Theory Lecture : The law of large umbers The Law of Large Numbers The
More informationNYU Center for Data Science: DS-GA 1003 Machine Learning and Computational Statistics (Spring 2018)
NYU Ceter for Data Sciece: DS-GA 003 Machie Learig ad Computatioal Statistics (Sprig 208) Brett Berstei, David Roseberg, Be Jakubowski Jauary 20, 208 Istructios: Followig most lab ad lecture sectios, we
More informationKLMED8004 Medical statistics. Part I, autumn Estimation. We have previously learned: Population and sample. New questions
We have previously leared: KLMED8004 Medical statistics Part I, autum 00 How kow probability distributios (e.g. biomial distributio, ormal distributio) with kow populatio parameters (mea, variace) ca give
More informationMath 220B Final Exam Solutions March 18, 2002
Math 0B Fial Exam Solutios March 18, 00 1. (1 poits) (a) (6 poits) Fid the Gree s fuctio for the tilted half-plae {(x 1, x ) R : x 1 + x > 0}. For x (x 1, x ), y (y 1, y ), express your Gree s fuctio G(x,
More information