Jiaping Wang Department of Mathematical Science 04/22/2013, Monday

1. Chapter 8. Some Approximations to Probability Distributions: Limit TheoremsSections 8.2 -- 8.3: Convergence in Probability and in Distribution Jiaping Wang Department of Mathematical Science 04/22/2013, Monday

2. Outline Convergence in Probability Convergence in Distribution

3. Part1. Convergence in Probability

4. Introduction Suppose that a coin has probability p, with 0≤p≤1, of coming up heads on a single flip. Suppose that we flip the coin n times, what can we say about the fraction of heads observed in the n flips? For example, if p=0.5, we draw different numbers of trials in a simulation, the result is given in the table From here, we can find when n∞, the ratio is closer to 0.5 and thus the difference is closer to zero.

5. Definition 8.1 In mathematical notations, let X denote the number of heads observed in the n tosses. Then E(X)=np, V(X)=np(1-p). One way to measure the closeness of X/n to p is to ascertain the probability that the distance will be less than a pre-assigned small value ε so that . Definition 8.1: The sequence of random variables X1,X2, .., Xn is said to convergence in probability to the constant c, if for every positive number ε,

6. Theorem 8.1 Weak Law of Large Numbers: Let X1,X2, .., Xnbe independent and identical distributed random variables, with E(Xi)=μ and V(Xi)=σ2<∞ for each i=1,…, n. Let Then, for any positive real number ε, Or . Thus, converges in probability toward μ. The proof can be shown based on the Tchebysheff’s theorem with X replaced by and σ2 by σ2/n, then let

7. Theorem 8.2 Suppose that Xn converges in probability toward μ1 and Yn converges in probability toward μ2. Then the following statements are also true. 1. Xn+Ynconverges in probability toward u1+u2. 2. XnYnconvergesin probability toward u1u2. 3. Xn/Ynconverges in probability toward u1/u2, provided u2≠0. 4. converges in probability toward , provided P(Xn≥0)=1.

8. Example 8.1 Let X be a binomial random variable with probability of success p and number of trials n. Show that X/n converges in probability toward p. Answer: We have seen that we can write X as ∑Yi with Yi=1 if the i-th trial results in Success, and Yi=0 otherwise. Then X/n=1/n ∑Yi . Also E(Yi)=p and V(Yi)=p(1-p). Then the conditions of Theorem 8.1 are fulfilled with μ=p and σ2=p(1-p)< ∞ and thus we can conclude that, for any positive ε, limn∞P(|X/n-p| ≥ε)=0.

9. Example 8.2 Suppose that X1, X2, …, Xn are independent and identically distributed random Variables with and all assumed finite. Let S2 denote the sample variance given by Show that S2 converges in probability to V(Xi). Answer: Notice that where The quantity is the average of n independent and identical distributed variables of the form with E()=and V ()=-which is finite. Thus Theorem 8.1 tell us that converges to in probability. Finally, based on Theorem 8.2, we can have converges in probability to - =V(Xi). This example shows that for large samples, the sample variance has a high probability of being close to the population variance.

10. Part 2. Convergence in Distribution

11. Definition 8.2 In the last section, we only study the convergence of certain random variables Toward constants. In this section, we study the probability distributions of certain type random variables as n tends toward infinity. Definition 8.2: Let Xn be a random variable with distribution function Fn(x). Let X be a random variable with distribution function F(x). If limn∞Fn(x)=F(x) At every point x for which F(x) is continuous, then Xn is said to converge in distribution toward X. F(x) is called the limiting distribution function of Xn.

12. Example 8.3 Let X1, X2, …, Xn be independent uniform random variables over the interval (θ, 0) for a negative constant θ. In addition, let Yn=min(X1, X2, …, Xn). Find the limiting distribution of Yn. Answer: The distribution function for the uniform random variable Xi is We know so we can find =

13. Theorem 8.3 Let Xn and X be random variables with moment-generating functions Mn(t) and M(t), respectively. If limn∞Mn(t)=M(t) For all real t, then Xn converges in distribution toward X.

14. Example 8.4 Let Xn be a binomial random variable with n trials and probability p of success on each trial. If n tends toward infinity and p tends zero with np remaining fixed. Show that Xn converges in distribution toward a Poisson random variable. Answer: We know the moment-generating function for the binomial random variables Xn, Mn(t) is given as = based on np=λ . Recall that Letting k=λ(et-1), we have which is the moment generating function of the Poisson random variable. As an example, when n=10 and p=0.1, we can find the true probability from the binomial Distribution is 0.73609 for X is less than 2 and the approximate value from the Poisson Is 0.73575, they are very close. So we can approximate the probability from binomial Distribution by the Poisson distribution when n is large and p is small.

15. Example 8.5 In monitoring for a pollution, an experiment collects a small volume of water and counts the number of bacteria in the sample. Unlike earlier problems, we have only one observation. For purposes of approximating the probability distribution of counts, we can think of the volume as the quantity that is getting large. Let X denote the bacteria count per cubic centimeter of water and assume that X has a Poisson probability distribution with mean λ, which we do by showing that converges in distribution toward a standard normal random variable as λ tends toward infinity. Specifically, if the allowable pollution in a water supply is a count of 110 bacteria per cubic centimeter, approximate the probability that X will be at most 110, assuming that λ=100.

16. Solution Answer: We know the mgf for Poisson random variable X is , thus We can have the mgf of Y as )-1)]. The term )-1) can be written as + Thus MY(t)=exp[++=exp[ When λ∞, MY(t) exp(t2/2) which is the mgf of the standard normal distribution. So we can approximate the probability of the Poisson random variable by the standard normal distribution when λ is large enough (for example, λ25).