MOMENT GENERATING FUNCTION AND STATISTICAL DISTRIBUTIONS 1 MOMENT

MOMENT GENERATING FUNCTION AND STATISTICAL DISTRIBUTIONS 1

MOMENT GENERATING FUNCTION The m. g. f. of random variable X is defined as for t Є (-h, h) for some h>0. 2

Properties of m. g. f. • M(0)=E[1]=1 • If a r. v. X has m. g. f. M(t), then Y=a. X+b has a m. g. f. • • M. g. f does not always exists (e. g. Cauchy distribution) 3

Example • Suppose that X has the following p. d. f. Find the m. g. f; expectation and variance. 4

CHARACTERISTIC FUNCTION The c. h. f. of random variable X is defined as for all real numbers t. C. h. f. always exists. 5

Uniqueness Theorem: 1. If two r. v. s have mg. f. s that exist and are equal, then they have the same distribution. 2. If two r, v, s have the same distribution, then they have the same m. g. f. (if they exist) Similar statements are true for c. h. f. 6

Problem • It is sometimes the case that exact values of random variables (Y 1, Y 2, …) cannot be observed, but we can observe they are greater than some fixed value. Let Y 1, Y 2, … be i. i. d. r. v. s. Let a be a fixed number on real line. For i=1, 2, … define, 7

Problem, cont. For example, if a manufacturing process produces parts with strength Yi that are tested to see if they can withstand stress a, then Xi denotes whether the strength is at least a or it is less than a. In such a case, we cannot directly observe the strength Yi of the ith part, but we can observe whether it breaks in stress test. 8

Problem, cont. • Define p=P(Y 1≥a) and q=1 -p, Sn=X 1+X 2+…+Xn. Note that, Sn is the number of Y 1, …, Yn that exceed a. i) Define the characteristic function, say , of a r. v. X ii) Find iv) Find P(Sn=j) 9

Other generating functions • log. M(t) is called cumulant generating function. • is factorial moment generating function. • Note: there is a simple relation between m. g. f. and f. m. g. f. 10

Other generating functions 11

Example • Suppose X has the following p. m. f. Find the expectation and variance of X. • Solution: Let’s use factorial m. g. f. 12

Example 13

STATISTICAL DISTRIBUTIONS 14

Recall • Random variable: A function defined on the sample space S that associates a real number with each outcome in S. 15

Example • Toss three coins • Sample space S={s 1=HHH, s 2=HHT, …, s 6=THT, s 7=TTH, s 8=TTT} • Define X=number of heads: X(s 1)=3, X(s 6)=1, X(s 8)=0 • Define Y=number of tails before first head: Y(s 1)=0, Y(s 6)=1, Y(s 8)=3 16

Random variables • A random variable is continuous if its CDF, F(x)=P(X≤x), is continuous. • A random variable is discrete if its CDF, F(x)=P(X≤x), is a step function. • It is possible for a CDF to have continuous pieces and steps, but we will mostly concentrate on the previous two bullets in this course. 17

SOME DISCRETE PROBABILITY DISTRIBUTIONS Degenerate, Uniform, Bernoulli, Binomial, Poisson, Negative Binomial, Geometric, Hypergeometric 18

DEGENERATE DISTRIBUTION • An rv X is degenerate at point k if The cdf: 19

UNIFORM DISTRIBUTION • A finite number of equally spaced values are equally likely to be observed. • Example: throw a fair die. P(X=1)=…=P(X=6)=1/6 20

BERNOULLI DISTRIBUTION • A Bernoulli trial is an experiment with only two outcomes. An r. v. X has Bernoulli(p) distribution if 21

BERNOULLI DISTRIBUTION • P(X=0)=1 -p and P(X=1)=p • E(X)=p • 22

BINOMIAL DISTRIBUTION • Define an rv Y by Y = total number of successes in n Bernoulli trials. 1. There are n trials (n is finite and fixed). 2. Each trial can result in a success or a failure. 3. The probability p of success is the same for all the trials. 4. All the trials of the experiment are independent. 23

BINOMIAL DISTRIBUTION • Example: • There are black and white balls in a box. Select and record the color of the ball. Put it back and re-pick (sampling with replacement). • n: number of independent and identical trials • p: probability of success (e. g. probability of picking a black ball) • X: number of successes in n trials 24

BINOMIAL THEOREM • For any real numbers x and y and integer n>0 25

BINOMIAL DISTRIBUTION • If Y~Bin(n, p), then 26

POISSON DISTRIBUTION • The number of occurrences in a given time interval can be modeled by the Poisson distribution. • e. g. number of customers to arrive in a bank between 13: 00 and 13: 30. • Another application is in spatial distributions. • e. g. modeling the distribution of bomb hits in an area or the distribution of fish in a lake. 27

POISSON DISTRIBUTION • If X~ Poi(λ), then • E(X)= Var(X)=λ • 28

Relationship between Binomial and Poisson Let =np. The mgf of Poisson( ) The limiting distribution of Binomial rv is the 29 Poisson distribution.

NEGATIVE BINOMIAL DISTRIBUTION (PASCAL OR WAITING TIME DISTRIBUTION) • X: number of Bernoulli trials required to get a fixed number of failures before the r th success; or, alternatively, • Y: number of Bernoulli trials required to get a fixed number of successes, such as r successes. 30

NEGATIVE BINOMIAL DISTRIBUTION (PASCAL OR WAITING TIME DISTRIBUTION) X~NB(r, p) 31

NEGATIVE BINOMIAL DISTRIBUTION • An alternative form of the pdf: Note: Y=X+r 32

GEOMETRIC DISTRIBUTION • Distribution of the number of Bernoulli trials required to get the first success. • It is the special case of the Negative Binomial Distribution r=1. X~Geometric(p) 33

GEOMETRIC DISTRIBUTION • Example: If probability is 0. 001 that a light bulb will fail on any given day, then what is the probability that it will last at least 30 days? • Solution: 34

HYPERGEOMETRIC DISTRIBUTION • A box contains N marbles. Of these, M are red. Suppose that n marbles are drawn randomly from the box without replacement. The distribution of the number of red marbles, x is X~Hypergeometric(N, M, n) It is dealing with finite population. 35

HYPERGEOMETRIC DISTRIBUTION • As N →∞, hypergeometric → binomial. • In that case, sampling with or without replacement does not make much difference (especially if n/N is small). 36

MULTIVARIATE DISTRIBUTIONS 37

EXTENDED HYPERGEOMETRIC DISTRIBUTION • Suppose that a collection consists of a finite number of items, N and that there are k+1 different types; M 1 of type 1, M 2 of type 2, and so on. Select n items at random without replacement, and let Xi be the number of items of type i that are selected. The vector X=(X 1, X 2, …, Xk) has an extended hypergeometric distribution and the joint pdf is 38

MULTINOMIAL DISTRIBUTION • Let E 1, E 2, . . . , Ek+1 be k+1 mutually exclusive and exhaustive events which can occur on any trial of an experiment with P(Ei)=pi, i=1, 2, …, k+1. On n independent trials of the experiment, let Xi be the number of occurrences of the event Ei. Then, the vector X=(X 1, X 2, …, Xk) has a multinomial distribution with joint pdf 39

MULTINOMIAL DISTRIBUTION • Experiment involves drawing with replacement. • Binomial is a special case of multinomial with k+1=2 40

MULTINOMIAL DISTRIBUTION • Consider trinomial case for simplicity. 41

MULTINOMIAL DISTRIBUTION • M. g. f. of X 1: X 1~Bin(n, p 1) Similarly, X 2~Bin(n, p 2) But, Cov(X 1, X 2)≠ 0! Cov(X 1, X 2)=? 42

MULTINOMIAL DISTRIBUTION • Example: Suppose we have a bowl with 10 marbles - 2 red marbles, 3 green marbles, and 5 blue marbles. We randomly select 4 marbles from the bowl, with replacement. What is the probability of selecting 2 green marbles and 2 blue marbles? 43

MULTINOMIAL DISTRIBUTION • n = 4, k+1=3, nred = 0, ngreen = 2, nblue = 2 • pred = 0. 2, pgreen = 0. 3, pblue = 0. 5 • P = [ n! / ( n 1! * n 2! *. . . nk! ) ] * ( p 1 n 1 * p 2 n 2 *. . . * pknk ) P = [ 4! / ( 0! * 2! ) ] * [ (0. 2)0 * (0. 3)2 * (0. 5)2 ] P = 0. 135 44

Problem 1. a) Does a distribution exist for which the m. g. f. ? If yes, find it. If no, prove it. b) Does a distribution exist for which the m. g. f. ? If yes, find it. If no, prove it. 45

Problem 2. An appliance store receives a shipment of 30 microwave ovens, 5 of which are (unknown to the manager) defective. The store manager selects 4 ovens at random, without replacement, and tests to see if they are defective. Let X=number of defectives found. Calculate the pmf and cdf of X. 46

Problem 3. Let X denote the number of “do loops” in a Fortran program and Y the number of runs needed for a novice to debug the program. Assume that the joint density for (X, Y) is given in the following table. 47

Problem x/y 1 2 3 4 0 0. 059 0. 1 0. 05 0. 001 1 0. 093 0. 12 0. 082 0. 003 2 0. 065 0. 102 0. 1 0. 01 3 0. 05 0. 07 0. 02 48

Problem a) Find the probability that a randomly selected program contains at most one “do loop” and requires at least two runs to debug the program. b) Find E[XY]. c) Find the marginal densities for X and Y. Find the mean and variance for both X and Y. d) Find the probability that a randomly selected program requires at least two runs to debug given that it contains exactly one “do loop”. e) Find Cov(X, Y). Find the correlation between X and Y. Based on the observed value correlation, can you claim that X and Y are not independent? Why? 49