

Line 33: 
Line 33: 
   
   
  ===Binomial Random Variables===  +  ===Example=== 
  Suppose we conduct an experiment observing an ntrial (fixed) Bernoulli process. If we are interested in the RV '''X={Number of successes in the n trials}''', then X is called a '''Binomial RV''' and its distribution is called '''Binomial Distribution'''.  +  Suppose we study N independent trials with results falling in one of k possible categories labeled <math>1,2, cdots, k</math>. Let <math>p_i</math> be the probability of a trial resulting in the <math>i^{th}</math> category, where <math>p_1+p_2+\cdots++p_k =1</math>. Let <math>N_i</math> be the number of trials resulting in the <math>i^{th}</math> category, where <math>N_1+N_2+\cdots++N_k = N</math>. 
   
  ====Examples====
 +  For instance, suppose we have 9 people arriving at a meeting according to the following information: 
  *Roll a standard die ten times. Let X be the number of times {6} turned up. The distribution of the random variable X is a binomial distribution with n = 10 (number of trials) and p = 1/6 (probability of "success={6}). The distribution of X may be explicitly written as (P(X=x) are rounded of, you can compute these exactly by going to [http://socr.ucla.edu/htmls/SOCR_Distributions.html SOCR Distributions] and selecting Binomial):
 +  : P(by Air) = 0.4, P(by Bus) = 0.2, P(by Automobile) = 0.3, P(by Train) = 0.1 
  <center>
 +  
  { class="wikitable" style="textalign:center; width:75%" border="1"
 +  
  
 +  
   x  0  1  2  3  4  5  6  7  8  9  10
 +  
  
 +  
   P(X=x)  0.162  0.323  0.291  0.155  0.0543  0.013  0.0022  0.00025  0.000019  8.269e7  1.654e8
 +  
  }
 +  
  </center>
 +  
   
  <center>[[Image:SOCR_EBook_Dinov_RV_Binomial_013008_Fig1.jpg400px]]</center>
 +  * Compute the following probabilities 
 +  : P(3 by Air, 3 by Bus, 1 by Auto, 2 by Train) = ? 
 +  : P(2 by air) = ? 
   
  * Suppose 10% of the human population carries the greeneye allele. If we choose 1,000 people randomly and let the RV X be the number of greeneyed people in the sample. Then the distribution of X is binomial distribution with n = 1,000 and p = 0.1 (denoted as <math>X \sim B(1,000, 0.1)</math>. In a sample of 1,000, how many are we expecting to have this allele?
 
   
  ===Binomial Modeling===  +  ===SOCR Multinomial Examples=== 
  ====Exact Binomial Model====
 +  
  The Binomial distribution (i.e., biasedcoin tossing experiment) is an exact physical model for any experiment which can be characterized as a series of trials where:
 +  
  * Each trial has only two outcomes: success or failure;
 +  
  * P(success)=p is the same for every trial; and
 +  
  * Trials are independent.
 +  
   +  
  ====Approximate Binomial Model====
 +  
  Suppose we observe an experiment comprised of n identical trials involving a large population (of size N). Assume the population contains a subpopulation of subjects that have a characteristic of interest and the subpopulation proportion is p (0<p<1). Then the distribution of X={the number of outcomes in the sample with that characteristic of interest}, is approximately Binomial(n, p). This approximation is adequate if the ratio n/N < 0.05.
 +  
   +  
  * Example: Polling the US population to see what proportion is/hasbeen married. Because we sample without replacement (can't ask the same individual twice), the second assumption of the (exact) Binomial model is (slightly) violated. Yet, the small ration of sample to population size yields that the Binomial model is approximately valid (i.e., the proportion of subjects that is/hasbeen married does not change significantly as we poll one subject, and therefore remove him/her from the complete pool of subjects we poll).
 +  
   +  
  ===Binomial Probabilities===
 +  
  If the random variable ''X'' follows the Binomial distribution with (fixed) parameters ''n'' (samplesize) and ''p'' (probability of success at one trial), we write ''X'' ~ B(''n'', ''p''). The probability of getting exactly ''x'' successes is given by the Binomial probability (or mass) function: <math>P(X=x)={n\choose k}p^k(1p)^{nk}</math>, for ''x'' = 0, 1, 2, ..., ''n'', where <math>{n\choose k}=\frac{n!}{k!(nk)!}</math> is the binomial coefficient.
 +  
   +  
  This probability expression has an easy and intuitive interpretation. The probability of the ''x'' successes in the ''n'' trials is (''p''<sup>''x''</sup>). Similarly, the probability of the ''nx'' failures is (1 − ''p'')<sup>''nx''</sup>. However, the ''x'' successes can be [[AP_Statistics_Curriculum_2007_Prob_Count  arranged]] anywhere among the ''n'' trials, and there are <math>{n\choose k}=\frac{n!}{k!(nk)!}</math> different ways of arranging the ''x'' successes in a sequence of ''n'' trials, see the [[AP_Statistics_Curriculum_2007_Prob_Count Counting section]].
 +  
   +  
  ===Binomial Expectation and Variance===
 +  
  If X is the random variable representing the r(random) number of heads in n coin toss trials, where the P(Head) = p, i.e., <math>X\sim B(n, p)</math>. Then we have the following expressions for the [[AP_Statistics_Curriculum_2007#Expectation_.28Mean.29_and_Variance  expected value, variance and the standard deviation]] of X :
 +  
  * Mean (expected value): <math>E[X] = n\times p</math>,
 +  
  * Variance: <math>VAR[X] = n\times p \times(1p)</math>, and
 +  
  * Standard deviation: <math>SD[X] = \sqrt{n\times p \times(1p)}</math>
 +  
  * [http://en.wikipedia.org/wiki/Binomial_distribution The complete description of the Binomial distribution is available here]
 +  
   +  
  ===Examples===
 +  
  ====Binomial Coin Toss Example====
 +  
  Refer to the SOCR [[SOCR_EduMaterials_Activities_BinomialCoinExperiment  Binomial Coin Toss Experiment]] and use the [http://socr.ucla.edu/htmls/SOCR_Experiments.html SOCR Binomial Coin Toss Applet] to perform an experiment of tossing a biased coin, P(Head) = 0.3, 5 times and computing the expectation of the number of Heads in such experiment. You should recognize that there are two distinct ways of computing the expected number of Heads:
 +  
   +  
  * Theoretical calculation, using the Binomial Probabilities;<math>E[X]=\sum_x{xP(X=x)} = </math>
 +  
  <math>={5\choose 0}0.3^0(0.7)^{5}+{5\choose 1}0.3^1(0.7)^{4}+{5\choose 2}0.3^2(0.7)^{3}+{5\choose 3}0.3^3(0.7)^{2}+{5\choose 4}0.3^4(0.7)^{1}+{5\choose 5}0.3^5(0.7)^{0} =</math>
 +  
  <math> \cdots = (n\times p) = 5\times 0.3 = 1.5.</math>
 +  
   +  
  * Empirical calculation, using the outcomes 100 repeated coin tosses of 5 coins. The image below illustrates this approximate calculation of the expectation for the number of heads when <math>X\sim B(5, 0.3)</math>. Notice the slight difference between the theoretical expectation (<math>n\times p = 5 \times 0.3 = 1.5</math>) and its empirical approximation of 1.39!
 +  
  <center>[[Image:SOCR_EBook_Dinov_RV_Binomial_013008_Fig2.jpg400px]]</center>
 +  
   +  
  * Notes: Of course, the theoretical calculation is exact and the empirical calculation is only approximate. However, the power of the empirical approximation to the expected number of Heads becomes clear when we increase the number of coins from 5 to 100, or more. In these cases, the exact calculations become very difficult, and for even larger number of coins become intractable.
 +  
   
 <hr>   <hr> 