Topic 10: The Law of Large Numbers
|
|
- Chrystal Walton
- 5 years ago
- Views:
Transcription
1 Topic : October 6, 2 If we choose adult Europea males idepedetly ad measure their heights, keepig a ruig average, the at the begiig we might see some larger fluctuatios but as we cotiue to make measuremets, we expect to see this ruig average settle ad coverge to the true mea height of Europea males. This pheomea is iformally kow as the law of averages. I probability theory, we call this the law of large umbers. Example. We ca simulate me s heights with idepedet ormal radom variables, mea 68 ad stadard deviatio 3.5. The folowig R commads perform this simulatio ad computes a ruig average of the heights. The results are displayed i Figure. > <-c(:) > x<-rorm(,68,3.5) > s<-cumsum(x) > plot(s/,xlab="",ylim=c(6,7),type="l") Here, we begi with a sequece X, X 2,... of radom variables havig a commo distributio. Their average, S = (X + X X ), is itself a radom variable. If the commo mea for the X i s is µ, the by the liearity property of expectatio, the mea of the average E[ S ] = (EX + EX EX ) = (µ + µ + + µ) = µ = µ () If, i additio, the X i s are idepedet with commo variace σ 2, the by the Pythagorea theorem for the variace of idepedet radom variables, Var( S ) = 2 (Var(X ) + Var(X 2 ) + + Var(X )) = 2 (σ2 + σ σ 2 ) = 2 σ2 = σ2. (2) So the mea of these ruig averages remais at µ but the variace is decreasig to at a rate iversely proportioal to the umber of terms i the sum. For example, the mea of the average height of Europea males is 68 iches, the stadard deviatio is σ/ = 3.5/ =.35 iches. For, males, the mea remais 68 iches, he stadard deviatio is σ/ = 3.5/ =.35 iches The mathematical result, the law of large umbers, tells us that the results of these simulatio could have bee aticipated. Theorem 2. For a sequece of idepedet radom variables X, X 2,... havig a commo distributio, their ruig average S = (X + + X ) has a limit as if ad oly if this sequece of radom variables has a commo mea µ. I this case the limit is µ. c 2 Joseph C. Watkis 29
2 Itroductio to Statistical Methodology s/ s/ s/ s/ Figure : Four simulatios of the ruig average S /, =, 2,..., for idepedet ormal radom variables, mea 68 ad stadard deviatio 3.5. The theorem also states that if the radom variables do ot have a mea, the as the ext example shows, the limit will fail to exist. We shall show with the followig example. Example 3. The stadard Cauchy radom variable X has desity fuctio Let Y = X. I a attempt to compute EY, ote that b b x f X (x) dx = 2 as b. Thus, Y has ifiite mea. b π f X (x) = π + x 2. x + x 2 dx = π l( + x2 ) b = l( + b 2 ) 3
3 Itroductio to Statistical Methodology We ow simulate idepedet Cauchy radom variables. > <-c(:) > y<-abs(rcauchy()) > s<-cumsum(y) > plot(s/,xlab="",ylim=c(-6,6),type="l") These radom variables do ot have a fiite mea. As you ca see i Figure 2 that their ruig averages do ot seem to be covergig. Thus, if we are usig a simulatio strategy that depeds o the law of large umbers, we eed to check that the radom variables have a mea. Exercise 4. Check the failure of the large umber of Cauchy radom variables via simulatios. Note that the shocks ca jump either up or dow i this case. Mote Carlo Itegratio Mote Carlo methods use stochastic simulatios to approximate solutios to questios too difficult to solve aalytically. This approach has see widespread use i fields as diverse as statistical physics, astroomy, populatio geetics, protei chemistry, ad fiace. Our itroductio will focus o examples havig relatively rapid computatios. However, may research groups routiely use Mote Carlo simulatios that ca take weeks of computer time to perform. For example, let X, X 2,... be idepedet radom variables uiformally distributed o the iterval [a, b] ad write f X for the desity of a uiform radom variable o the iterval [a, b]. The, by the law of large umbers, g(x) = g(x i ) Eg(X ) = i= b lsrge. Let s deote by I(g) the value of this itegral. Recall that i calculus, we defied the average of g to be b a b a a g(x)f X (x) dx = b a g(x) dx. We ca also iterpret this itegral as the expected value of g(x ). Thus, Mote Carlo itegratio leads to a procedure for estimatig itegrals. simulatig uiform radom variables X, X 2,..., X o the iterval [a, b], computig the average of g(x ), g(x 2 ),..., g(x ) to estimate the average of g, ad multiplyig by b a to estimate the itegral. b a g(x) dx Example 5. Let g(x) = + cos 3 (x) for x [, π], to fid π g(x) dx. The three steps above become the followig R code. > x<-ruif(,,pi) > g<-sqrt(+cos(x)ˆ3) > pi*mea(g) []
4 Itroductio to Statistical Methodology s/ Figure 2: Four simulatios of the ruig average S /, =, 2,..., for the absolute value of idepedet Cauchy radom variables. Note that the ruig averate does ot seem to be settlig dow ad is subject to shocks. The error i the estimate of the itegral ca be estimated by the variace as give i equatio (2). Var(g(X) ) = Var(g(X )). where σ 2 = Var(g(X )) = E(g(X ) µ g(x)) 2 = b a (g(x) µ g(x )) 2 f X (x) dx. Typically this itegral is more difficult to estimate tha b g(x) dx, our origial itegral of iterest. However, we ca see that the variace of the a estimator is iversely proportioal to, the umber of radom umbers i the simulatio. Thus, the stadard deviatio is iversely proportioal to. Mote Carlo techiques are rarely the best strategy for estimatig oe or eve very low dimesioal itegrals. R 32
5 Itroductio to Statistical Methodology does itegratio umerically usig the fuctio ad the itegrate commads. For example, > g<-fuctio(x){sqrt(+cos(x)ˆ3)} > itegrate(g,,pi) with absolute error < 3.8e-6 However, with oly a small chage i the algorithm, we ca also use this to evaluate high dimesioal multivariate itegrals. For example, i three dimesios, the itegral I(g) = g(x, y, z) dx dy dz ca be estimated usig Mote Carlo itegratio by geeratig three sequeces of uiform radom variables, The, X, X 2,..., X, Y, Y 2,..., Y, ad Z, Z 2,... Z I(g) g(x i, Y i, Z i ). (3) i= Example 6. To obtai a sese of the distributio of the approximatios to I(g), we perform simulatios usig uiform radom variable for each of the three coordiates to perform the Mote Carlo itegratio. The commad Ig<-rep(,) creates a vector of zeros. This is added so that R creates a place ahead of the simulatios to store the results. > Ig<-rep(,) > for(i i :){x<-ruif();y<-ruif();z<-ruif(); g<-32*xˆ3/(3*(y+zˆ4+)); Ig[i]<-mea(g)} > hist(ig) > summary(ig) Mi. st Qu. Media Mea 3rd Qu. Max > var(ig) [] > sd(ig) [] Exercise 7. Estimate the variace ad stadard deviatio of the Mote Carlo estimator for the itegral (3) based o = 5 ad radom umbers. Exercise 8. How may observatios are eeded i estimatig (3) so that the stadard deviatio of the average is.5? To modify this techique for a regio [a, b ] [a 2, b 2 ] [a 3, b 3 ] use idepeet uiform radom variables X i, Y i, ad Z i o the respective itervals, the b b2 b3 g(x i, Y i, Z i ) Eg(X, Y, Z ) = g(x, y, z) dz dy dx. b a b 2 a 2 b 3 a 3 a a 2 a 3 i= Thus, the estimate for the itegral is (b a )(b 2 a 2 )(b 3 a 3 ) g(x i, Y i, Z i ). i= 33
6 Itroductio to Statistical Methodology Histogram of Ig Frequecy 5 5 Figure 3: Histogram of Mote Carlo estimates for the itegral R R σ = Importace Samplig Ig R 32x3 /(y + z 4 + ) dx dy dz. The sample stadard deviatio I may of the large simulatios, the dimesio of the itegral ca be i the hudreds ad the fuctio g ca be very close to zero for large regios i the domai of g. Simple Mote Carlo simulatio will the frequetly choose values for g that are close to zero. These values cotribute very little to the average. Due to this iefficiecy, a more sophisticated strategy is employed. Importace samplig methods begi with the observatio that a better strategy may be to cocetrate the radom poits where the itegrad g is large i absolute value. For example, for the itegral e x/2 dx, (4) x( x) the itegrad is much bigger for values ear x = or x =. (See Figure 4) Thus, we ca hope to have a more accurate estimate by cocetratig our sample poits i these places. With this i mid, we perform the Mote Carlo itegratio begiig with Y, Y 2,... idepedet radom variables with commo desityf Y. The goal is to fid a desity f Y that is large whe g is large ad small whe g is small. The desity f Y is called the importace samplig fuctio or the proposal desity. With this choice of desity, we defie the importace samplig weights To justify this choice, ote that, the sample mea w(y ) = w(y i ) i= w(y) = g(y) f Y (y). (5) w(x)f Y (x) dy = g(x) f Y (x) f Y (x) dx = I(g). Thus, the average of the importace samplig weights, by the strog law of large umbers, still approximates the itegral of g. This is a improvemet over simple Mote Carlo itegratio if the variace decreases, i.e., Var(w(Y )) = (w(y) I(g)) 2 f Y (y) dy = σ 2 f << σ 2. 34
7 Itroductio to Statistical Methodology As the formula shows, this ca be best achieved by havig the weight w(y) be close to the itegral I(g). Referrig to equatio (5), we ca ow see that we should edeavor to have f Y proportioal to g. Example 9. For the itegral (4) we ca use Mote Carlo simulatio based o uiform radom variables. > Ig<-rep(,) > for(i i :){x<-ruif();g<-exp(-x/2)*/sqrt(x*(-x));ig[i]<-mea(g)} > summary(ig) Mi. st Qu. Media Mea 3rd Qu. Max > sd(ig) [] Based o a simulatios, we fid a sample mea value of ad a sample stadard deviatio of.394. Because the itegrad is very large ear x = ad x =, we choose look for a desity f Y to cocetrate the radom samples ear the eds of the itervals. Our choice for the proposal desity is a Beta(/2, /2), the f Y (y) = π y/2 ( y) /2 o the iterval [, ]. Thus the weight w(y) = πe y/ f Y (x)=//sqrt(x( x)).5 g(x) = e x/2 /sqrt(x( x)) x Figure 4: Importace samplig usig the desity fuctio f Y to estimate R g(x) dx. The weight w(x) = πe x/2 is the ratio g(x)/f Y (x). Agai, we perform the simulatio multiple times. > IS<-rep(,) > for(i i :){y<-rbeta(,/2,/2);w<-pi*exp(-y/2);is[i]<-mea(w)} > summary(is) Mi. st Qu. Media Mea 3rd Qu. Max. 35
8 Itroductio to Statistical Methodology > var(is) [].2595 > sd(is) [] Based o simulatios, we fid a sample mea value of ad a sample stadard deviatio of.44, about /9th the size of the Mote Carlo weright. Part of the gai is illusory. Beta radom variables take loger to simulate. If they require a factor more tha 8 times loger to simulate, the the extra work eeded to create a good importace sample is ot helpful i producig a more accurate estimate for the itegral. Simpso s rule ca also be used i this case ad we obtai the aswer Exercise. Evaluate the itegral e x 3 x dx times usig = 2 sample poits usig directly Mote Carlo itegratio ad usig importace samplig with radom variables havig desity f X (x) = x o the iterval [, ]. For the secod part, you will eed to use the probability trasform. Compare the meas ad stadard deviatios of the estimates for the itegral. The itegral is approximately Aswers to Selected Exercises 4. Here are the R commads: > par(mfrow=c(2,2)) > x<-rcauchy() > s<-cumsum(x) > plot (,s/,type="l") > x<-rcauchy() > s<-cumsum(x) > plot (,s/,type="l") > x<-rcauchy() > s<-cumsum(x) > plot (,s/,type="l") > x<-rcauchy() > s<-cumsum(x) > plot (,s/,type="l") This produces the plots below. Notice the differeces for the values o the x-axis 7. The stadard deviatio for the average of observatios is σ/ where σ is the stadard deviatio for a sigle observatio. From the output > sd(ig) [] We have that σ/ = σ/. Thus, σ Cosequetly, for 5 observatios, σ/ For observatios σ/ For σ/ =.5, we have that = (σ/.5) So we eed approximately 4 observatios.. For the direct Mote Carlo simulatio, we have 36
9 Itroductio to Statistical Methodology s/ -2 - s/ s/ s/ > Ig<-rep(,) > for (i i :){x<-ruif(2);g<-exp(-x)/xˆ(/3);ig[i]<-mea(g)} > mea(ig) [] > sd(ig) [] For the importace sampler, the itegral is 3 2 e x f X (x) dx. 37
10 Itroductio to Statistical Methodology To simulate idepedet radom variables with desity f X, we first eed the cumulative distributio fuctio for X, F X (x) = The, to fid the probability trasform, ote that x t dt = x t2/3 = x 2/3. u = F X (x) = x 2/3 ad x = F X (u) = u3/2. Thus, to simulate X, we simulate a uiform radom variable U o the iterval [, ] ad evaluate U 3/2. This leads to the followig R commads for both the simple Mote Carlo ad the importace sample simulatio: > Ig<-rep(,) > for (i i :){u<-ruif(2);g<-exp(-u)/uˆ(/3);ig[i]<-mea(g)} > mea(ig) [].4976 > sd(ig) [] > ISg<-rep(,) > for (i i :){u<-ruif(2);x<-uˆ(3/2); w<-3*exp(-x)/2;isg[i]<-mea(w)} > mea(isg) [].4845 > sd(isg) [].232 Thus, the stadard deviatio usig importace samplig is about 2/7-ths the stadard deviatio usig simple Mote Carlo simulatio. Cosequetly, we will ca decrease the umber of samples usig importace samplig by a factor of (2/7)
Topic 8: Expected Values
Topic 8: Jue 6, 20 The simplest summary of quatitative data is the sample mea. Give a radom variable, the correspodig cocept is called the distributioal mea, the epectatio or the epected value. We begi
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 informationChapter 2 The Monte Carlo Method
Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful
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 informationTopic 10: Introduction to Estimation
Topic 0: Itroductio to Estimatio Jue, 0 Itroductio I the simplest possible terms, the goal of estimatio theory is to aswer the questio: What is that umber? What is the legth, the reactio rate, the fractio
More informationMonte Carlo method and application to random processes
Mote Carlo method ad applicatio to radom processes Lecture 3: Variace reductio techiques (8/3/2017) 1 Lecturer: Eresto Mordecki, Facultad de Ciecias, Uiversidad de la República, Motevideo, Uruguay Graduate
More information1 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 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 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 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 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 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 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 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 informationExpectation and Variance of a random variable
Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio
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 informationMonte Carlo Integration
Mote Carlo Itegratio I these otes we first review basic umerical itegratio methods (usig Riema approximatio ad the trapezoidal rule) ad their limitatios for evaluatig multidimesioal itegrals. Next we itroduce
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 informationThis section is optional.
4 Momet Geeratig Fuctios* This sectio is optioal. The momet geeratig fuctio g : R R of a radom variable X is defied as g(t) = E[e tx ]. Propositio 1. We have g () (0) = E[X ] for = 1, 2,... Proof. Therefore
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013. Large Deviations for i.i.d. Random Variables
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013 Large Deviatios for i.i.d. Radom Variables Cotet. Cheroff boud usig expoetial momet geeratig fuctios. Properties of a momet
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 informationLearning Theory: Lecture Notes
Learig Theory: Lecture Notes Kamalika Chaudhuri October 4, 0 Cocetratio of Averages Cocetratio of measure is very useful i showig bouds o the errors of machie-learig algorithms. We will begi with a basic
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 informationTopic 18: Composite Hypotheses
Toc 18: November, 211 Simple hypotheses limit us to a decisio betwee oe of two possible states of ature. This limitatio does ot allow us, uder the procedures of hypothesis testig to address the basic questio:
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 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 informationNANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS
NANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS STRUCTURE OF EXAMINATION PAPER. There will be oe 2-hour paper cosistig of 4 questios.
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 informationAMS570 Lecture Notes #2
AMS570 Lecture Notes # Review of Probability (cotiued) Probability distributios. () Biomial distributio Biomial Experimet: ) It cosists of trials ) Each trial results i of possible outcomes, S or F 3)
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 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 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 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 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 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 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 informationn outcome is (+1,+1, 1,..., 1). Let the r.v. X denote our position (relative to our starting point 0) after n moves. Thus X = X 1 + X 2 + +X n,
CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 9 Variace Questio: At each time step, I flip a fair coi. If it comes up Heads, I walk oe step to the right; if it comes up Tails, I walk oe
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 informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
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 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 informationAAEC/ECON 5126 FINAL EXAM: SOLUTIONS
AAEC/ECON 5126 FINAL EXAM: SOLUTIONS SPRING 2015 / INSTRUCTOR: KLAUS MOELTNER This exam is ope-book, ope-otes, but please work strictly o your ow. Please make sure your ame is o every sheet you re hadig
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
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 informationSurveying the Variance Reduction Methods
Iteratioal Research Joural of Applied ad Basic Scieces 2013 Available olie at www.irjabs.com ISSN 2251-838X / Vol, 7 (7): 427-432 Sciece Explorer Publicatios Surveyig the Variace Reductio Methods Arash
More informationDirection: This test is worth 150 points. You are required to complete this test within 55 minutes.
Term Test 3 (Part A) November 1, 004 Name Math 6 Studet Number Directio: This test is worth 10 poits. You are required to complete this test withi miutes. I order to receive full credit, aswer each problem
More informationAMS 216 Stochastic Differential Equations Lecture 02 Copyright by Hongyun Wang, UCSC ( ( )) 2 = E X 2 ( ( )) 2
AMS 216 Stochastic Differetial Equatios Lecture 02 Copyright by Hogyu Wag, UCSC Review of probability theory (Cotiued) Variace: var X We obtai: = E X E( X ) 2 = E( X 2 ) 2E ( X )E X var( X ) = E X 2 Stadard
More informationJanuary 25, 2017 INTRODUCTION TO MATHEMATICAL STATISTICS
Jauary 25, 207 INTRODUCTION TO MATHEMATICAL STATISTICS Abstract. A basic itroductio to statistics assumig kowledge of probability theory.. Probability I a typical udergraduate problem i probability, we
More information6. Sufficient, Complete, and Ancillary Statistics
Sufficiet, Complete ad Acillary Statistics http://www.math.uah.edu/stat/poit/sufficiet.xhtml 1 of 7 7/16/2009 6:13 AM Virtual Laboratories > 7. Poit Estimatio > 1 2 3 4 5 6 6. Sufficiet, Complete, ad Acillary
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 informationBayesian Methods: Introduction to Multi-parameter Models
Bayesia Methods: Itroductio to Multi-parameter Models Parameter: θ = ( θ, θ) Give Likelihood p(y θ) ad prior p(θ ), the posterior p proportioal to p(y θ) x p(θ ) Margial posterior ( θ, θ y) is Iterested
More informationAn Introduction to Randomized Algorithms
A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis
More informationMath 10A final exam, December 16, 2016
Please put away all books, calculators, cell phoes ad other devices. You may cosult a sigle two-sided sheet of otes. Please write carefully ad clearly, USING WORDS (ot just symbols). Remember that the
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 400: Georgios Fellouris Homework 5 Due: Friday 24 th, 2017
Stat 400: Georgios Fellouris Homework 5 Due: Friday 4 th, 017 1. A exam has multiple choice questios ad each of them has 4 possible aswers, oly oe of which is correct. A studet will aswer all questios
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More informationMathematical Statistics - MS
Paper Specific Istructios. The examiatio is of hours duratio. There are a total of 60 questios carryig 00 marks. The etire paper is divided ito three sectios, A, B ad C. All sectios are compulsory. Questios
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 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 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 informationLecture 2: Poisson Sta*s*cs Probability Density Func*ons Expecta*on and Variance Es*mators
Lecture 2: Poisso Sta*s*cs Probability Desity Fuc*os Expecta*o ad Variace Es*mators Biomial Distribu*o: P (k successes i attempts) =! k!( k)! p k s( p s ) k prob of each success Poisso Distributio Note
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 informationDiscrete Mathematics for CS Spring 2007 Luca Trevisan Lecture 22
CS 70 Discrete Mathematics for CS Sprig 2007 Luca Trevisa Lecture 22 Aother Importat Distributio The Geometric Distributio Questio: A biased coi with Heads probability p is tossed repeatedly util the first
More informationSTAT Homework 2 - Solutions
STAT-36700 Homework - Solutios Fall 08 September 4, 08 This cotais solutios for Homework. Please ote that we have icluded several additioal commets ad approaches to the problems to give you better isight.
More informationChapter 6 Part 5. Confidence Intervals t distribution chi square distribution. October 23, 2008
Chapter 6 Part 5 Cofidece Itervals t distributio chi square distributio October 23, 2008 The will be o help sessio o Moday, October 27. Goal: To clearly uderstad the lik betwee probability ad cofidece
More informationReview Questions, Chapters 8, 9. f(y) = 0, elsewhere. F (y) = f Y(1) = n ( e y/θ) n 1 1 θ e y/θ = n θ e yn
Stat 366 Lab 2 Solutios (September 2, 2006) page TA: Yury Petracheko, CAB 484, yuryp@ualberta.ca, http://www.ualberta.ca/ yuryp/ Review Questios, Chapters 8, 9 8.5 Suppose that Y, Y 2,..., Y deote a radom
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationHOMEWORK I: PREREQUISITES FROM MATH 727
HOMEWORK I: PREREQUISITES FROM MATH 727 Questio. Let X, X 2,... be idepedet expoetial radom variables with mea µ. (a) Show that for Z +, we have EX µ!. (b) Show that almost surely, X + + X (c) Fid the
More informationSTAT 350 Handout 19 Sampling Distribution, Central Limit Theorem (6.6)
STAT 350 Hadout 9 Samplig Distributio, Cetral Limit Theorem (6.6) A radom sample is a sequece of radom variables X, X 2,, X that are idepedet ad idetically distributed. o This property is ofte abbreviated
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationMATH/STAT 352: Lecture 15
MATH/STAT 352: Lecture 15 Sectios 5.2 ad 5.3. Large sample CI for a proportio ad small sample CI for a mea. 1 5.2: Cofidece Iterval for a Proportio Estimatig proportio of successes i a biomial experimet
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 informationEE 4TM4: Digital Communications II Probability Theory
1 EE 4TM4: Digital Commuicatios II Probability Theory I. RANDOM VARIABLES A radom variable is a real-valued fuctio defied o the sample space. Example: Suppose that our experimet cosists of tossig two fair
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 informationSampling Error. Chapter 6 Student Lecture Notes 6-1. Business Statistics: A Decision-Making Approach, 6e. Chapter Goals
Chapter 6 Studet Lecture Notes 6-1 Busiess Statistics: A Decisio-Makig Approach 6 th Editio Chapter 6 Itroductio to Samplig Distributios Chap 6-1 Chapter Goals After completig this chapter, you should
More informationEco411 Lab: Central Limit Theorem, Normal Distribution, and Journey to Girl State
Eco411 Lab: Cetral Limit Theorem, Normal Distributio, ad Jourey to Girl State 1. Some studets may woder why the magic umber 1.96 or 2 (called critical values) is so importat i statistics. Where do they
More informationUnbiased Estimation. February 7-12, 2008
Ubiased Estimatio February 7-2, 2008 We begi with a sample X = (X,..., X ) of radom variables chose accordig to oe of a family of probabilities P θ where θ is elemet from the parameter space Θ. For radom
More informationA quick activity - Central Limit Theorem and Proportions. Lecture 21: Testing Proportions. Results from the GSS. Statistics and the General Population
A quick activity - Cetral Limit Theorem ad Proportios Lecture 21: Testig Proportios Statistics 10 Coli Rudel Flip a coi 30 times this is goig to get loud! Record the umber of heads you obtaied ad calculate
More informationUnderstanding Samples
1 Will Moroe CS 109 Samplig ad Bootstrappig Lecture Notes #17 August 2, 2017 Based o a hadout by Chris Piech I this chapter we are goig to talk about statistics calculated o samples from a populatio. We
More informationProbability 2 - Notes 10. Lemma. If X is a random variable and g(x) 0 for all x in the support of f X, then P(g(X) 1) E[g(X)].
Probability 2 - Notes 0 Some Useful Iequalities. Lemma. If X is a radom variable ad g(x 0 for all x i the support of f X, the P(g(X E[g(X]. Proof. (cotiuous case P(g(X Corollaries x:g(x f X (xdx x:g(x
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 3 9/11/2013. Large deviations Theory. Cramér s Theorem
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/5.070J Fall 203 Lecture 3 9//203 Large deviatios Theory. Cramér s Theorem Cotet.. Cramér s Theorem. 2. Rate fuctio ad properties. 3. Chage of measure techique.
More informationMONTE CARLO VARIANCE REDUCTION METHODS
MONTE CARLO VARIANCE REDUCTION METHODS M. Ragheb /4/3. INTRODUCTION The questio arises of whether oe ca reduce the variace associated with the samplig of a radom variable? Ideed we ca, but we eed to be
More informationSignals & Systems Chapter3
Sigals & Systems Chapter3 1.2 Discrete-Time (D-T) Sigals Electroic systems do most of the processig of a sigal usig a computer. A computer ca t directly process a C-T sigal but istead eeds a stream of
More informationExercise 4.3 Use the Continuity Theorem to prove the Cramér-Wold Theorem, Theorem. (1) φ a X(1).
Assigmet 7 Exercise 4.3 Use the Cotiuity Theorem to prove the Cramér-Wold Theorem, Theorem 4.12. Hit: a X d a X implies that φ a X (1) φ a X(1). Sketch of solutio: As we poited out i class, the oly tricky
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 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 informationSTAT 515 fa 2016 Lec Sampling distribution of the mean, part 2 (central limit theorem)
STAT 515 fa 2016 Lec 15-16 Samplig distributio of the mea, part 2 cetral limit theorem Karl B. Gregory Moday, Sep 26th Cotets 1 The cetral limit theorem 1 1.1 The most importat theorem i statistics.............
More informationJoint Probability Distributions and Random Samples. Jointly Distributed Random Variables. Chapter { }
UCLA STAT A Applied Probability & Statistics for Egieers Istructor: Ivo Diov, Asst. Prof. I Statistics ad Neurology Teachig Assistat: Neda Farziia, UCLA Statistics Uiversity of Califoria, Los Ageles, Sprig
More informationStat 200 -Testing Summary Page 1
Stat 00 -Testig Summary Page 1 Mathematicias are like Frechme; whatever you say to them, they traslate it ito their ow laguage ad forthwith it is somethig etirely differet Goethe 1 Large Sample Cofidece
More informationCHAPTER 8 FUNDAMENTAL SAMPLING DISTRIBUTIONS AND DATA DESCRIPTIONS. 8.1 Random Sampling. 8.2 Some Important Statistics
CHAPTER 8 FUNDAMENTAL SAMPLING DISTRIBUTIONS AND DATA DESCRIPTIONS 8.1 Radom Samplig The basic idea of the statistical iferece is that we are allowed to draw ifereces or coclusios about a populatio based
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 informationConfidence Interval for Standard Deviation of Normal Distribution with Known Coefficients of Variation
Cofidece Iterval for tadard Deviatio of Normal Distributio with Kow Coefficiets of Variatio uparat Niwitpog Departmet of Applied tatistics, Faculty of Applied ciece Kig Mogkut s Uiversity of Techology
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 informationEECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1
EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum
More informationBIOSTATISTICS. Lecture 5 Interval Estimations for Mean and Proportion. dr. Petr Nazarov
Microarray Ceter BIOSTATISTICS Lecture 5 Iterval Estimatios for Mea ad Proportio dr. Petr Nazarov 15-03-013 petr.azarov@crp-sate.lu Lecture 5. Iterval estimatio for mea ad proportio OUTLINE Iterval estimatios
More informationIntroduction to Econometrics (3 rd Updated Edition) Solutions to Odd- Numbered End- of- Chapter Exercises: Chapter 3
Itroductio to Ecoometrics (3 rd Updated Editio) by James H. Stock ad Mark W. Watso Solutios to Odd- Numbered Ed- of- Chapter Exercises: Chapter 3 (This versio August 17, 014) 015 Pearso Educatio, Ic. Stock/Watso
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationSTAT Homework 1 - Solutions
STAT-36700 Homework 1 - Solutios Fall 018 September 11, 018 This cotais solutios for Homework 1. Please ote that we have icluded several additioal commets ad approaches to the problems to give you better
More information