Normal Distribution. Lecture 6: More Binomial Distribution. Properties of the Unit Normal Distribution. Unit Normal Distribution

Similar documents
Lecture 3 Gaussian Probability Distribution

8 Laplace s Method and Local Limit Theorems

Continuous Random Variables

Solution for Assignment 1 : Intro to Probability and Statistics, PAC learning

1 Probability Density Functions

Lecture 21: Order statistics

Chapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses

Chapter 5 : Continuous Random Variables

For the percentage of full time students at RCC the symbols would be:

Expectation and Variance

7 - Continuous random variables

CS 109 Lecture 11 April 20th, 2016

CS667 Lecture 6: Monte Carlo Integration 02/10/05

Section 11.5 Estimation of difference of two proportions

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

Tests for the Ratio of Two Poisson Rates

The steps of the hypothesis test

Read section 3.3, 3.4 Announcements:

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

38.2. The Uniform Distribution. Introduction. Prerequisites. Learning Outcomes

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

Best Approximation. Chapter The General Case

3.4 Numerical integration

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Orthogonal Polynomials

Section 17.2 Line Integrals

12 TRANSFORMING BIVARIATE DENSITY FUNCTIONS

p(t) dt + i 1 re it ireit dt =

AP Calculus Multiple Choice: BC Edition Solutions

Monte Carlo method in solving numerical integration and differential equation

Math 426: Probability Final Exam Practice

Quantum Physics II (8.05) Fall 2013 Assignment 2

Taylor Polynomial Inequalities

Tutorial 4. b a. h(f) = a b a ln 1. b a dx = ln(b a) nats = log(b a) bits. = ln λ + 1 nats. = log e λ bits. = ln 1 2 ln λ + 1. nats. = ln 2e. bits.

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.

Summary: Method of Separation of Variables

4.1. Probability Density Functions

The Regulated and Riemann Integrals

Riemann Sums and Riemann Integrals

AQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system

Math 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8

Recitation 3: More Applications of the Derivative

Problem Set 9. Figure 1: Diagram. This picture is a rough sketch of the 4 parabolas that give us the area that we need to find. The equations are:

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

Math 360: A primitive integral and elementary functions

X Z Y Table 1: Possibles values for Y = XZ. 1, p

Riemann Sums and Riemann Integrals

CHM Physical Chemistry I Chapter 1 - Supplementary Material

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

The final exam will take place on Friday May 11th from 8am 11am in Evans room 60.

Administrivia CSE 190: Reinforcement Learning: An Introduction

Math 113 Exam 2 Practice

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

Review of Gaussian Quadrature method

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

4181H Problem Set 11 Selected Solutions. Chapter 19. n(log x) n 1 1 x x dx,

Improper Integrals, and Differential Equations

4.4 Areas, Integrals and Antiderivatives

Chapter 0. What is the Lebesgue integral about?

Review of Probability Distributions. CS1538: Introduction to Simulations

Lecture 14: Quadrature

Lecture 12: Numerical Quadrature

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

Fundamental Theorem of Calculus

RATE OF CONVERGENCE OF POLYA S URN TO THE BETA DISTRIBUTION JOHN DRINANE

MAA 4212 Improper Integrals

10. AREAS BETWEEN CURVES

Math 8 Winter 2015 Applications of Integration

MATH 144: Business Calculus Final Review

Math 31S. Rumbos Fall Solutions to Assignment #16

Lecture 1. Functional series. Pointwise and uniform convergence.

DISCRETE MATHEMATICS HOMEWORK 3 SOLUTIONS

Numerical Integration

Mathacle. PSet Stats, Concepts In Statistics Level Number Name: Date: , the DeMoivre-Laplace Theorem states: ( ) n x.

Overview of Calculus I

7.2 The Definite Integral

Line Integrals. Partitioning the Curve. Estimating the Mass

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

5.2 Exponent Properties Involving Quotients

A sequence is a list of numbers in a specific order. A series is a sum of the terms of a sequence.

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

How can we approximate the area of a region in the plane? What is an interpretation of the area under the graph of a velocity function?

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Population bottleneck : dramatic reduction of population size followed by rapid expansion,

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

1 2-D Second Order Equations: Separation of Variables

Section Areas and Distances. Example 1: Suppose a car travels at a constant 50 miles per hour for 2 hours. What is the total distance traveled?

different methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).

Chapter 3 Solving Nonlinear Equations

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Math 1B, lecture 4: Error bounds for numerical methods

Definite integral. Mathematics FRDIS MENDELU

Chapter 1. Basic Concepts

Transcription:

Norml Distribution Lecture 6: More Binomil Distribution If X is rndom vrible with norml distribution with men µ nd vrince σ 2, X N (µ, σ 2, then P(X = x = f (x = 1 e 1 (x µ 2 2 σ 2 σ Sttistics 104 Colin Rundel Februry 1, 2012 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 1 / 21 Unit Norml Distribution The unit norml distribution is specil cse of the norml distribution where µ = 0 nd σ = 1, Z N (0, 1. P(Z = z = φ(z = 1 e 1 2 z2 The re under the unit norml curve from to is given by P(z = 1 e t2 /2 dt = Φ( 3 2 1 0 1 2 3 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 2 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 3 / 21

The re under the unit norml curve from to is given by P(z = 1 e t2 /2 dt = 1 Φ(x The re under the unit norml curve from to b where b is given by P( z b = 1 b e t2 /2 dt = Φ(b Φ( 3 2 1 0 1 2 3 3 2 1 0 b 1 2 3 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 4 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 5 / 21 Evluting Φ The re under the unit norml curve outside of to b where b is given by P( z b = 1 1 b e t2 /2 dt = Φ( + (1 Φ(b = 1 (Φ(b Φ( The function Φ(x hs no simple closed form solution, sometimes written in terms of the error function erf Φ(x = 1 x e x2 /2 = 1 [ ( ] x 1 + erf 2 2 where erf(x = 2 π x 0 e t2 This doesn t seem to get us very fr, but we cn tke the Tylor expnsion of e t2 nd evlute the integrl of ech term to get series for erf(x. 3 2 1 0 b 1 2 3 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 6 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 7 / 21

Evluting Φ - In prctice Evluting Φ - Some prctice 362 APPENDIX C. DISTRIBUTION TABLES In prctice you will never hve to evlute Φ explicitly, we either look Φ(x up in tble or use computer to clculte it (pnorm in R. Working with your neighbor(s see if you cn work out the following probbilities for rndom vrible Z N(0, 1. negtive Z Second deciml plce of Z 0.09 0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0.00 Z 0.0002 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 3.4 0.0003 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0005 0.0005 0.0005 3.3 0.0005 0.0005 0.0005 0.0006 0.0006 0.0006 0.0006 0.0006 0.0007 0.0007 3.2 0.0007 0.0007 0.0008 0.0008 0.0008 0.0008 0.0009 0.0009 0.0009 0.0010 3.1 0.0010 0.0010 0.0011 0.0011 0.0011 0.0012 0.0012 0.0013 0.0013 0.0013 3.0 0.0014 0.0014 0.0015 0.0015 0.0016 0.0016 0.0017 0.0018 0.0018 0.0019 2.9 0.0019 0.0020 0.0021 0.0021 0.0022 0.0023 0.0023 0.0024 0.0025 0.0026 2.8 0.0026 0.0027 0.0028 0.0029 0.0030 0.0031 0.0032 0.0033 0.0034 0.0035 2.7 0.0036 0.0037 0.0038 0.0039 0.0040 0.0041 0.0043 0.0044 0.0045 0.0047 2.6 0.0048 0.0049 0.0051 0.0052 0.0054 0.0055 0.0057 0.0059 0.0060 0.0062 2.5 0.0064 0.0066 0.0068 0.0069 0.0071 0.0073 0.0075 0.0078 0.0080 0.0082 2.4 0.0084 0.0087 0.0089 0.0091 0.0094 0.0096 0.0099 0.0102 0.0104 0.0107 2.3 0.0110 0.0113 0.0116 0.0119 0.0122 0.0125 0.0129 0.0132 0.0136 0.0139 2.2 0.0143 0.0146 0.0150 0.0154 0.0158 0.0162 0.0166 0.0170 0.0174 0.0179 2.1 0.0183 0.0188 0.0192 0.0197 0.0202 0.0207 0.0212 0.0217 0.0222 0.0228 2.0 0.0233 0.0239 0.0244 0.0250 0.0256 0.0262 0.0268 0.0274 0.0281 0.0287 1.9 0.0294 0.0301 0.0307 0.0314 0.0322 0.0329 0.0336 0.0344 0.0351 0.0359 1.8 0.0367 0.0375 0.0384 0.0392 0.0401 0.0409 0.0418 0.0427 0.0436 0.0446 1.7 0.0455 0.0465 0.0475 0.0485 0.0495 0.0505 0.0516 0.0526 0.0537 0.0548 1.6 0.0559 0.0571 0.0582 0.0594 0.0606 0.0618 0.0630 0.0643 0.0655 0.0668 1.5 0.0681 0.0694 0.0708 0.0721 0.0735 0.0749 0.0764 0.0778 0.0793 0.0808 1.4 0.0823 0.0838 0.0853 0.0869 0.0885 0.0901 0.0918 0.0934 0.0951 0.0968 1.3 Sttistics 104 (Colin Rundel 0.0985 0.1003 0.1020 Lecture 0.10386: 0.1056 More Binomil 0.1075 0.1093 Distribution 0.1112 0.1131 0.1151 1.2 Februry 1, 2012 8 / 21 0.1170 0.1190 0.1210 0.1230 0.1251 0.1271 0.1292 0.1314 0.1335 0.1357 1.1 0.1379 0.1401 0.1423 0.1446 0.1469 0.1492 0.1515 0.1539 0.1562 0.1587 1.0 0.1611 0.1635 0.1660 0.1685 0.1711 0.1736 0.1762 0.1788 0.1814 0.1841 0.9 0.1867 0.1894 0.1922 0.1949 0.1977 0.2005 0.2033 0.2061 0.2090 0.2119 0.8 0.2148 0.2177 0.2206 0.2236 0.2266 0.2296 0.2327 0.2358 0.2389 0.2420 0.7 0.2451 0.2483 0.2514Chpter 0.2546 2.1-2.3 0.2578 0.2611 0.2643 0.2676 0.2709 0.2743 0.6 0.2776 0.2810 0.2843 0.2877 0.2912 0.2946 0.2981 0.3015 0.3050 0.3085 0.5 0.3121 0.3156 0.3192 0.3228 0.3264 0.3300 0.3336 0.3372 0.3409 0.3446 0.4 0.3483 0.3520 0.3557 0.3594 0.3632 0.3669 0.3707 0.3745 0.3783 0.3821 0.3 0.3859 0.3897 0.3936 0.3974 0.4013 0.4052 0.4090 0.4129 0.4168 0.4207 0.2 0.4247 0.4286 0.4325 0.4364 0.4404 0.4443 0.4483 0.4522 0.4562 0.4602 0.1 0.4641 0.4681 0.4721 0.4761 0.4801 0.4840 0.4880 0.4920 0.4960 0.5000 0.0 For Z 3.50, the probbility is less thn or equl to 0.0002. Empiricl Rule P(z < 1 P(z > 2.22 P( 1.53 z 2.75 P(0.75 z 1.43 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 9 / 21 Stndrdizing Norml Distributions Everything we just discussed pplies only to the unit norml distribution, but this doesn t come up very often in useful problems. P( 1 z 1 = Φ(1 Φ( 1 = 0.683 P( 2 z 2 = Φ(2 Φ( 2 = 0.954 Let X be normlly distributed rndom vrible with men µ nd vrince σ 2 then we define the rndom vrible Z such tht ( x µ Z = N(0, 1 σ In previous clsses you hve probbly referred to this s z-score. We will see why this works soon when we get to functions of rndom vribles. P( x b = b ( 1 exp 1 σ 2 (x µ 2 b µ µ σ 2 dx = Φ Φ σ σ P( 3 z 3 = Φ(3 Φ( 3 = 0.997 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 10 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 11 / 21

Another wy to see DeMoivre-Lplce We sw nother useful pproximtion lst week - Stirling s pproximtion to the fctoril function which is prticulrly good for lrge n. n! n n n e n Stirling s pproximtion is bsed on the Stirling Series n! = ( n ( n 1 n 1 + e 2 1 (6n 1 + 1 2 3 (6n 2 139 2 3 (2 3 5(6n 3 571 2 6 (2 3 5(6n 4 + We strt with the definition of the Binomil distribution nd use Sterling s pproximtion for the fctorils n! P(X = k = (n k! k! pk (1 p n k n n n e n (n k(n k n k e n+k ( k k k e k pk q (n k n k(n k ( np ( k nq n k k n k Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 12 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 13 / 21 Set x = k np, so k = np + x nd n k = nq x then ( 1 np np+x ( nq P(X = np + x npq np + x nq x ( np + x np x ( nq x c np nq nq+x c (1 + x/np np x (1 x/nq nq+x nq x Like lst time, we will use n pproximtion for log(1 ± x which is bsed on the Tylor expnsions log(1 + x = x 1 2 x2 + 1 3 x3 + ( 1 n 1 n xn log(1 x = x 1 2 x2 1 3 x3 1 n xn Tke the log of both sides of the eqution nd use the 1st nd 2nd term of the Tylor expnsion for the pproximtion to get log P(X = np + x c (np + x log(1 + x/np (nq x log(1 x/nq ( ( x c (np + x np x2 2(np 2 (nq x x nq x2 2(nq 2 P(X = k = c x x2 np + x2 2np + x3 2(np 2 + x x2 nq + x2 2nq c 2qx2 2npq + qx2 2npq + 2px2 2npq + px2 2npq + x2 c 2npq ( 1 exp 1 npq 2 (k np 2 npq x3 2(nq 2 ( x 3 2(np 2 x3 2(nq 2 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 14 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 15 / 21

Improving the pproximtion Improving the pproximtion, cont. Tke for exmple Binomil distribution where n = 20 nd p = 10, we should be ble to pproximte the distribution with X N (10, 5. Binomil probbility: 13 ( 20 P(7 x 13 = 0.5 x (1 0.5 20 x x x=7 Nive pproximtion: P(7 x 13 Φ 13 10 7 10 Φ 5 5 2 3 4 5 6 7 8 9 10 12 14 16 18 It is cler tht our pproximtion is missing 1/2 of P(X = 7 nd P(X = 13, s n this error is very smll. In this cse P(X = 7 = P(X = 13 = 0.073 so our pproximtion is off by 7%. Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 16 / 21 Continuity corrected pproximtion: 13 + 1/2 10 7 1/2 10 P(7 x 13 Φ Φ 5 5 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 17 / 21 Improving the pproximtion, cont. Improving the pproximtion, cont. When we scle up n to 200 nd our rnge interest to P(70 x 130: Binomil probbility: Nive pproximtion: P(70 x 130 = P(70 x 130 Φ 130 x=70 ( 200 x 0.5 x (1 0.5 200 x 130 100 70 100 Φ 50 50 Continuity corrected pproximtion: 130 + 1/2 100 70 1/2 100 P(70 x 130 Φ Φ 500 50 This correction lso lets us do, modertely useless, things like clculte the probbility for prticulr vlue of k. Such s, wht is the chnce of 50 Heds in 100 tosses of slightly unfir coin (p = 0.55? Binomil probbility: P(x = 50 = ( 100 0.55 50 (1 0.55 50 50 Nive pproximtion: 50 55 50 55 P(x = 50 Φ Φ 4.97 4.97 Continuity corrected pproximtion: 50 + 1/2 55 50 1/2 55 P(x = 50 Φ Φ 4.97 4.97 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 18 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 19 / 21

Smple Problem (1 Smple Problem (2-2.2.13 Roll fir die 500 times, wht s the probbility of rolling t lest 100 ones? A pollster wishes to know the percentge p of people in popultion who intend to vote for prticulr cndidte. How lrge must rndom smple with replcement be in order to be t lest 95% sure tht the smple percentge is within one percentge point of p? Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 20 / 21 Sttistics 104 (Colin Rundel Lecture 6: More Binomil Distribution Februry 1, 2012 21 / 21