#set document(title: "5.2 Binomial Probability Distribution", author: "Kathryn Kozak") #set page(width: 8.5in, height: auto, margin: 1in) #import "@preview/cetz:0.5.2" #set text(font: ("STIX Two Text", "Libertinus Serif", "New Computer Modern"), size: 10.5pt, lang: "en") #show math.equation: set text(font: ("STIX Two Math", "New Computer Modern Math")) #set par(justify: true, leading: 0.62em, spacing: 0.9em) #set enum(spacing: 1.1em) // room between list items so tall inline fractions don't collide #set list(spacing: 1.1em) #set table(stroke: 0.5pt + rgb("#c7ccd3")) #let BLUE = rgb("#183B6F") // brand navy — section bars + example/solution labels (white on navy 11.09:1) #let ORANGE = rgb("#A94509") // brand primary-700 — AA-safe deep orange for TEXT (5.93:1 on white; raw brand #F37021 is 2.94:1 and must never carry text) #let RED = rgb("#DC2626") // brand error-600 #let GREEN = rgb("#059669") // brand success-600 (decoration only; small green text uses green-text #007942) #show heading.where(level: 1): it => block(width: 100%, above: 0pt, below: 16pt, fill: gradient.linear(BLUE, rgb("#2C5AA0")), inset: (x: 14pt, y: 12pt), radius: 3pt, text(fill: white, weight: "bold", size: 19pt, it.body)) #show heading.where(level: 2): it => block(width: 100%, above: 18pt, below: 10pt, fill: BLUE, inset: (x: 10pt, y: 6pt), radius: 2pt, text(fill: white, weight: "bold", size: 12pt, it.body)) #show heading.where(level: 3): it => text(fill: ORANGE, weight: "bold", size: 12.5pt, it.body) #show heading.where(level: 4): it => text(fill: BLUE, weight: "bold", size: 10.5pt, it.body) #let examplebox(label, title, body) = block(width: 100%, breakable: true, fill: rgb("#EFF1F5"), stroke: 0.5pt + rgb("#CFDDF0"), radius: 4pt, inset: 10pt, above: 12pt, below: 12pt)[ #block(below: 6pt)[#box(fill: BLUE, inset: (x: 6pt, y: 2pt), radius: 2pt, text(fill: white, weight: "bold", size: 8.5pt, label)) #h(0.4em) #strong[#title]] #body] // rail = decorative left rule (raw brand token); labelcolor = AA-safe label text shade #let notebox(label, rail, labelcolor, tint, body) = block(width: 100%, breakable: true, fill: tint, stroke: (left: 3pt + rail), inset: (left: 10pt, rest: 8pt), radius: (right: 4pt), above: 11pt, below: 11pt)[ #text(fill: labelcolor, weight: "bold", size: 7.5pt, tracking: 0.5pt)[#upper(label)] #linebreak() #body] #let solutionbox(body) = block(above: 4pt, below: 8pt)[ #text(fill: BLUE, weight: "bold", size: 8.5pt)[Solution] #linebreak() #body] #let figph(msg) = block(width: 100%, height: 60pt, fill: rgb("#f6f7f9"), stroke: (paint: rgb("#c7ccd3"), dash: "dashed"), radius: 4pt, inset: 10pt)[ #align(center + horizon, text(fill: rgb("#889"), style: "italic", size: 9pt, msg))] // Standardize inlined figure sizes: measure the natural CeTZ canvas, then scale to a // consistent envelope (aspect-aware; see build_typst.py FIG_* constants). Unlike the // print preamble, dimensions are FLOORED: in an editor a user can trim a figure to a // degenerate 1-D shape (a bare line), and w/h or tw/w would then divide by zero. #let _STD_W = 3.5 #let _WIDE_W = 5.6 #let _MAX_H = 3.4 #let _ASPECT_WIDE = 2.2 #let _UPSCALE_MAX = 1.15 #let stdfig(body) = context { let m = measure(body) let w = calc.max(m.width / 1in, 0.01) let h = calc.max(m.height / 1in, 0.01) let tw = if w / h > _ASPECT_WIDE { _WIDE_W } else { _STD_W } let s = calc.min(tw / w, _MAX_H / h, _UPSCALE_MAX) align(center, box(scale(x: s * 100%, y: s * 100%, reflow: true, body))) } #show figure: set block(breakable: false) #set figure(gap: 8pt) #show figure.caption: set text(size: 8.5pt, fill: rgb("#555")) == 5.2#h(0.6em)Binomial Probability Distribution Section 5.1 introduced the concept of a probability distribution. The focus of the section was on discrete #strong[probability distributions] (pdf). To find the pdf for a situation, you usually needed to actually conduct the experiment and collect data. Then you can calculate the experimental probabilities. Normally you cannot calculate the theoretical probabilities instead. However, there are certain types of experiment that allow you to calculate the theoretical probability. One of those types is called a #strong[Binomial Experiment]. === Properties of a binomial experiment (or Bernoulli trial) + Fixed number of trials, #emph[n], which means that the experiment is repeated a specific number of times. + The #emph[n] trials are independent, which means that what happens on one trial does not influence the outcomes of other trials. + There are only two outcomes, which are called a success and a failure. + The probability of a success doesn’t change from trial to trial, where p = probability of success and #emph[q] = probability of failure, #emph[q =]1-#emph[p]. If you know you have a binomial experiment, then you can calculate binomial probabilities. This is important because binomial probabilities come up often in real life. Examples of binomial experiments are: - Toss a fair coin ten times, and find the probability of getting two heads. - Question twenty people in class, and look for the probability of more than half being women? - Shoot five arrows at a target, and find the probability of hitting it five times? To develop the process for calculating the probabilities in a binomial experiment, consider Example #math.equation(block: false, alt: "1")[$1$]. #examplebox("Example 1")[Deriving the Binomial Probability Formula][ Suppose you are given a 3 question multiple-choice test. Each question has 4 responses and only one is correct. Suppose you want to find the probability that you can just guess at the answers and get 2 questions right. (Teachers do this all the time when they make up a multiple-choice test to see if students can still pass without studying. In most cases the students can’t.) To help with the idea that you are going to guess, suppose the test is in Martian. + What is the random variable? + Is this a binomial experiment? + What is the probability of getting 2 questions right? + What is the probability of getting zero right, one right, and all three right? #solutionbox[ a. x = number of correct answers b. + There are 3 questions, and each question is a trial, so there are a fixed number of trials. In this case, #emph[n] = 3. + Getting the first question right has no affect on getting the second or third question right, thus the trials are independent. + Either you get the question right or you get it wrong, so there are only two outcomes. In this case, the success is getting the question right. + The probability of getting a question right is one out of four. This is the same for every trial since each question has 4 responses. In this case, #math.equation(block: false, alt: "p equals the fraction 1 over 4 and q equals 1 minus the fraction 1 over 4 equals the fraction 3 over 4")[$p = frac(1, 4) " and " q = 1 − frac(1, 4) = frac(3, 4)$] This is a binomial experiment, since all of the properties are met. c. To answer this question, start with the sample space. SS = {RRR, RRW, RWR, WRR, WWR, WRW, RWW, WWW}, where RRW means you get the first question right, the second question right, and the third question wrong. The same is similar for the other outcomes. Now the event space for getting 2 right is {RRW, RWR, WRR}. What you did in chapter four was just to find three divided by eight. However, this would not be right in this case. That is because the probability of getting a question right is different from getting a question wrong. What else can you do? Look at just P(RRW) for the moment. Again, that means P(RRW) = P(R on 1st, R on 2nd, and W on 3rd) Since the trials are independent, then P(RRW) = P(R on 1st, R on 2nd, and W on 3rd) = P(R on 1st) \* P(R on 2nd) \* P(W on 3rd) Just multiply p \* p \* q #math.equation(block: true, alt: "P open parenthesis R R W close parenthesis equals the fraction 1 over 4 * the fraction 1 over 4 * the fraction 3 over 4 equals open parenthesis the fraction 1 over 4 close parenthesis squared open parenthesis the fraction 3 over 4 close parenthesis to the power 1")[$P ( upright(R) upright(R) upright(W) ) = frac(1, 4) * frac(1, 4) * frac(3, 4) = attach(( frac(1, 4) ), t: 2) attach(( frac(3, 4) ), t: 1)$] The same is true for P(RWR) and P(WRR). To find the probability of 2 correct answers, just add these three probabilities together. You get #math.equation(block: true, alt: "P open parenthesis 2 correct answers close parenthesis equals P open parenthesis R R W close parenthesis plus P open parenthesis R W R close parenthesis plus P open parenthesis W R R close parenthesis; equals open parenthesis the fraction 1 over 4 close parenthesis squared open parenthesis the fraction 3 over 4 close parenthesis to the power 1 plus open parenthesis the fraction 1 over 4 close parenthesis squared open parenthesis the fraction 3 over 4 close parenthesis to the power 1 plus open parenthesis the fraction 1 over 4 close parenthesis squared open parenthesis the fraction 3 over 4 close parenthesis to the power 1; equals 3 open parenthesis the fraction 1 over 4 close parenthesis squared open parenthesis the fraction 3 over 4 close parenthesis to the power 1")[$P ( 2 " correct answers " ) = P ( upright(R) upright(R) upright(W) ) + P ( upright(R) upright(W) upright(R) ) + P ( upright(W) upright(R) upright(R) ) \ = attach(( frac(1, 4) ), t: 2) attach(( frac(3, 4) ), t: 1) + attach(( frac(1, 4) ), t: 2) attach(( frac(3, 4) ), t: 1) + attach(( frac(1, 4) ), t: 2) attach(( frac(3, 4) ), t: 1) \ = 3 attach(( frac(1, 4) ), t: 2) attach(( frac(3, 4) ), t: 1)$] d. You could go through the same argument that you did above and come up with the following: #figure(table( columns: 2, align: left, inset: 6pt, table.header([r right], [P(r right)]), [0 right], [#math.equation(block: false, alt: "1 to the power * open parenthesis the fraction 1 over 4 close parenthesis to the power 0 open parenthesis the fraction 3 over 4 close parenthesis cubed")[$1^(*) attach(( frac(1, 4) ), t: 0) attach(( frac(3, 4) ), t: 3)$]], [1 right], [#math.equation(block: false, alt: "3 to the power * open parenthesis the fraction 1 over 4 close parenthesis to the power 1 open parenthesis the fraction 3 over 4 close parenthesis squared")[$3^(*) attach(( frac(1, 4) ), t: 1) attach(( frac(3, 4) ), t: 2)$]], [2 right], [#math.equation(block: false, alt: "3 * open parenthesis the fraction 1 over 4 close parenthesis squared open parenthesis the fraction 3 over 4 close parenthesis to the power 1")[$3 * attach(( frac(1, 4) ), t: 2) attach(( frac(3, 4) ), t: 1)$]], [3 right], [#math.equation(block: false, alt: "1 to the power * open parenthesis the fraction 1 over 4 close parenthesis cubed open parenthesis the fraction 3 over 4 close parenthesis to the power 0")[$1^(*) attach(( frac(1, 4) ), t: 3) attach(( frac(3, 4) ), t: 0)$]], )) Hopefully you see the pattern that results. You can now write the general formula for the probabilities for a Binomial experiment ] ] First, the random variable in a binomial experiment is #emph[x] = number of successes. Be careful, a success is not always a good thing. Sometimes a success is something that is bad, like finding a defect. A success just means you observed the outcome you wanted to see happen. #notebox("Note", rgb("#8a94a6"), rgb("#556666"), rgb("#f7f8fa"))[ Binomial Formula for the probability of #emph[r] successes in #emph[n] trials is #math.equation(block: true, alt: "P open parenthesis x equals r close parenthesis equals sub n C sub r p to the power r q to the power n times r where sub n C sub r equals the fraction n ! over r ! open parenthesis n minus r close parenthesis !")[$P ( x = r ) = n C_(r) p^(r) q^(n · r) attach(" where ", b: n) C_(r) = frac(n !, r ! ( n − r ) !)$] ] The #math.equation(block: false, alt: "sub n C sub r")[$n C_(r)$] is the number of combinations of #emph[n] things taking #emph[r] at a time. It is read “#emph[n] choose #emph[r]”. Some other common notations for #emph[n] choose #emph[r] are #math.equation(block: false, alt: "C sub n , r")[$C_(n , r)$], and #math.equation(block: false, alt: "n; r")[$n \ r$]. #emph[n]! means you are multiplying #math.equation(block: false, alt: "n to the power * open parenthesis n minus 1 close parenthesis to the power * open parenthesis n minus 2 close parenthesis to the power * and so on to the power * 2 to the power * 1")[$n^(*) ( n − 1 )^(*) ( n − 2 )^(*) …^(*) 2^(*) 1$]. As an example, #math.equation(block: false, alt: "5 ! equals 5 to the power * 4 to the power * 3 to the power * 2 to the power * 1 equals 120")[$5 ! = 5^(*) 4^(*) 3^(*) 2^(*) 1 = 120$]. When solving problems, make sure you define your random variable and state what #emph[n, p, q,] and#emph[r] are. Without doing this, the problems are a great deal harder. #examplebox("Example 2")[Calculating Binomial Probabilities][ When looking at a person’s eye color, it turns out that 1% of people in the world has green eyes ("What percentage of," 2013). Consider a group of 20 people. + State the random variable. + Argue that this is a binomial experiment. + Find the probability that none have green eyes. + Find the probability that nine have green eyes. + Find the probability that at most three have green eyes. + Find the probability that at most two have green eyes. + Find the probability that at least four have green eyes. + In Europe, four people out of twenty have green eyes. Is this unusual? What does that tell you? #solutionbox[ a. #emph[x] = number of people with green eyes b. + There are 20 people, and each person is a trial, so there are a fixed number of trials. In this case, #math.equation(block: false, alt: "n equals 20")[$n = 20$]. + If you assume that each person in the group is chosen at random the eye color of one person doesn’t affect the eye color of the next person, thus the trials are independent. + Either a person has green eyes or they do not have green eyes, so there are only two outcomes. In this case, the success is a person has green eyes. + The probability of a person having green eyes is 0.01. This is the same for every trial since each person has the same chance of having green eyes. #emph[p] = 0.01 and #emph[q] = 1 - 0.01 = 0.99 c. #math.equation(block: false, alt: "P open parenthesis x equals 0 close parenthesis equals sub 20 C sub 0 open parenthesis 0.01 close parenthesis to the power 0 open parenthesis 0.99 close parenthesis to the power 20 minus 0 approximately equals 0.818")[$P ( x = 0 ) = 20 C_(0) ( 0.01 )^(0) ( 0.99 )^(20 − 0) ≈ 0.818$] d. #math.equation(block: false, alt: "P open parenthesis x equals 9 close parenthesis equals sub 20 C sub 9 open parenthesis 0.01 close parenthesis to the power 9 open parenthesis 0.99 close parenthesis to the power 20 minus 9 approximately equals 1.50 times 10 to the power minus 13 approximately equals 0.000")[$P ( x = 9 ) = 20 C_(9) ( 0.01 )^(9) ( 0.99 )^(20 − 9) ≈ 1.50 × 10^(− 13) ≈ 0.000$] e. At most three means that three is the highest value you will have. Find the probability of #emph[x] is less than or equal to three. #math.equation(block: true, alt: "P open parenthesis x less than or equal to 3 close parenthesis equals P open parenthesis x equals 0 close parenthesis plus P open parenthesis x equals 1 close parenthesis plus P open parenthesis x equals 2 close parenthesis plus P open parenthesis x equals 3 close parenthesis; equals sub 20 C sub 0 open parenthesis 0.01 close parenthesis to the power 0 open parenthesis 0.99 close parenthesis to the power 20 plus sub 20 C sub 1 open parenthesis 0.01 close parenthesis to the power 1 open parenthesis 0.99 close parenthesis to the power 19; plus sub 20 C sub 2 open parenthesis 0.01 close parenthesis squared open parenthesis 0.99 close parenthesis to the power 18 plus sub 20 C sub 3 open parenthesis 0.01 close parenthesis cubed open parenthesis 0.99 close parenthesis to the power 17; approximat")[$P ( x ≤ 3 ) = P ( x = 0 ) + P ( x = 1 ) + P ( x = 2 ) + P ( x = 3 ) \ = 20 C_(0) ( 0.01 )^(0) ( 0.99 )^(20) + 20 C_(1) ( 0.01 )^(1) ( 0.99 )^(19) \ + 20 C_(2) ( 0.01 )^(2) ( 0.99 )^(18) + 20 C_(3) ( 0.01 )^(3) ( 0.99 )^(17) \ ≈ 0.818 + 0.165 + 0.016 + 0.001 > 0.999$] The reason the answer is written as being greater than 0.999 is because the answer is actually 0.9999573791, and when that is rounded to three decimal places you get 1. But 1 means that the event will happen, when in reality there is a slight chance that it won’t happen. It is best to write the answer as greater than 0.999 to represent that the number is very close to 1, but isn’t 1. f. #math.equation(block: true, alt: "P open parenthesis x less than or equal to 2 close parenthesis equals P open parenthesis x equals 0 close parenthesis plus P open parenthesis x equals 1 close parenthesis plus P open parenthesis x equals 2 close parenthesis; equals sub 20 C sub 0 open parenthesis 0.01 close parenthesis to the power 0 open parenthesis 0.99 close parenthesis to the power 20 plus sub 20 C sub 1 open parenthesis 0.01 close parenthesis to the power 1 open parenthesis 0.99 close parenthesis to the power 19 plus sub 20 C sub 2 open parenthesis 0.01 close parenthesis squared open parenthesis 0.99 close parenthesis to the power 18; approximately equals 0.818 plus 0.165 plus 0.016 approximately equals 0.999")[$P ( x ≤ 2 ) = P ( x = 0 ) + P ( x = 1 ) + P ( x = 2 ) \ = 20 C_(0) ( 0.01 )^(0) ( 0.99 )^(20) + 20 C_(1) ( 0.01 )^(1) ( 0.99 )^(19) + 20 C_(2) ( 0.01 )^(2) ( 0.99 )^(18) \ ≈ 0.818 + 0.165 + 0.016 ≈ 0.999$] g. At least four means four or more. Find the probability of x being greater than or equal to four. That would mean adding up all the probabilities from four to twenty. This would take a long time, so it is better to use the idea of complement. The complement of being greater than or equal to four is being less than four. That would mean being less than or equal to three. Part (e) has the answer for the probability of being less than or equal to three. Just subtract that number from 1. #math.equation(block: true, alt: "P open parenthesis x greater than or equal to 4 close parenthesis equals 1 minus P open parenthesis x less than or equal to 3 close parenthesis equals 1 minus 0.999 equals 0.001")[$P ( x ≥ 4 ) = 1 − P ( x ≤ 3 ) = 1 − 0.999 = 0.001$] Actually the answer is less than 0.001, but it is fine to write it this way. h. Since the probability of finding four or more people with green eyes is much less than 0.05, it is unusual to find four people out of twenty with green eyes. That should make you wonder if the proportion of people in Europe with green eyes is more than the 1% for the general population. If this is true, then you may want to ask why Europeans have a higher proportion of green-eyed people. That of course could lead to more questions. ] ] The binomial formula is cumbersome to use, so you can find the probabilities by using technology. On the TI-83/84 calculator, the commands on the TI-83/84 calculators when the number of trials is equal to #emph[n] and the probability of a success is equal to #emph[p] are #math.equation(block: false, alt: "binompdf open parenthesis n , p , r close parenthesis")[$"binompdf" ( n , p , r )$] when you want to find #emph[P] (x=r) and #math.equation(block: false, alt: "binomcdf open parenthesis n , p , r close parenthesis")[$"binomcdf" ( n , p , r )$] when you want to find #math.equation(block: false, alt: "P open parenthesis x less than or equal to r close parenthesis")[$P ( x ≤ r )$]. If you want to find #math.equation(block: false, alt: "P open parenthesis x greater than or equal to r close parenthesis")[$P ( x ≥ r )$], then you use the property that #math.equation(block: false, alt: "P open parenthesis x greater than or equal to r close parenthesis equals 1 minus P open parenthesis x less than or equal to r minus 1 close parenthesis")[$P ( x ≥ r ) = 1 − P ( x ≤ r − 1 )$], since #math.equation(block: false, alt: "x greater than or equal to r")[$x ≥ r$] and #math.equation(block: false, alt: "x less than r")[$x < r$] or #math.equation(block: false, alt: "x less than or equal to r minus 1")[$x ≤ r − 1$] are complementary events. Both binompdf and binomcdf commands are found in the DISTR menu. Using R, the commands are #math.equation(block: false, alt: "P open parenthesis x equals r close parenthesis equals dbinom open parenthesis r , n , p close parenthesis and P open parenthesis x less than or equal to r close parenthesis equals pbinom open parenthesis r , n , p close parenthesis")[$P ( x = r ) = " dbinom " ( r , n , p ) " and " P ( x ≤ r ) = " pbinom " ( r , n , p )$]. #examplebox("Example 3")[using the binomial command on the ti-83/84][ When looking at a person’s eye color, it turns out that 1% of people in the world has green eyes ("What percentage of," 2013). Consider a group of 20 people. + State the random variable. + Find the probability that none have green eyes. + Find the probability that nine have green eyes. + Find the probability that at most three have green eyes. + Find the probability that at most two have green eyes. + Find the probability that at least four have green eyes. #solutionbox[ a. #emph[x] = number of people with green eyes b. You are looking for #emph[P] (x=0). Since this problem is x=0, you use the binompdf command on the TI-83/84 or dbinom command on R. On the TI83/84, you go to the DISTR menu, select the binompdf, and then type into the parenthesis your #emph[n], #emph[p], and #emph[r] values into your calculator, making sure you use the comma to separate the values. The command will look like #math.equation(block: false, alt: "binompdf open parenthesis 20 , .01 , 0 close parenthesis")[$"binompdf" ( 20 , ".01" , 0 )$] and when you press ENTER you will be given the answer. (If you have the new software on the TI-84, the screen looks a bit different.) #figure(figph[TI-84 calculator screen showing the command binompdf(20,.01,0) and the result 0.8179069376.], alt: "TI-84 calculator screen showing the command binompdf(20,.01,0) and the result 0.8179069376.", caption: [Figure #math.equation(block: false, alt: "1")[$1$]: Calculator Results for binompdf]) On R, the command would look like dbinom(0, 20, 0.01) dbinom(0, 20, 0.01) #emph[P] (x=0) = 0.8179. Thus there is an 81.8% chance that in a group of 20 people none of them will have green eyes. #notebox("Note", rgb("#8a94a6"), rgb("#556666"), rgb("#f7f8fa"))[ #emph[binompdf and binomcdf without the DISTR menu] The Binomial panel computes the cumulative P(x ≤ k) — the calculator's binomcdf. Cumulative up to k = 0 is the same number as binompdf(20, .01, 0); for P(x = 9) subtract P(x ≤ 8) from P(x ≤ 9); part e is the complement 1 − P(x ≤ 3). - P(x = 0) = 0.8179 — the binompdf(20, .01, 0) screen - P(x ≤ 3) = 0.99996 — the binomcdf(20, .01, 3) screen ] c. In this case you want to find the #emph[P] (x=9). Again, you will use the binompdf command or the dbinom command. Following the procedure above, you will have binompdf(20, .01, 9) on the TI-83/84 or dbinom(9,20,0.01) on R. Your answer is #math.equation(block: false, alt: "P open parenthesis x equals 9 close parenthesis equals 1.50 times 10 to the power minus 13")[$P ( x = 9 ) = 1.50 × 10^(− 13)$]. (Remember when the calculator gives you #math.equation(block: false, alt: "1.50 E minus 13")[$1.50 E − 13$] and R give you #math.equation(block: false, alt: "1.50 e minus 13")[$1.50 e − 13$], this is how they display scientific notation.) The probability that out of twenty people, nine of them have green eyes is a very small chance. dbinom(9, 20, 0.01) d. At most three means that three is the highest value you will have. Find the probability of #emph[x] being less than or equal to three, which is #math.equation(block: false, alt: "P open parenthesis x less than or equal to 3 close parenthesis")[$P ( x ≤ 3 )$]. This uses the binomcdf command on the TI-83/84 and pbinom command in R. You use the command on the TI-83/84 of binomcdf(20, .01, 3) and the command on R of pbinom(3,20,0.01) pbinom(3, 20, 0.01) #figure(figph[TI-83/84 calculator screen showing binomcdf(20,.01,3) with result 0.9999573791.], alt: "TI-83/84 calculator screen showing binomcdf(20,.01,3) with result 0.9999573791.", caption: [Figure #math.equation(block: false, alt: "2")[$2$]: Calculator Results for binomcdf]) Your answer is 0.99996. Thus there is a really good chance that in a group of 20 people at most three will have green eyes. (Note: don’t round this to one, since one means that the event will happen, when in reality there is a slight chance that it won’t happen. It is best to write the answer out to enough decimal points so it doesn’t round off to one. e. You are looking for #math.equation(block: false, alt: "P open parenthesis x less than or equal to 2 close parenthesis")[$P ( x ≤ 2 )$]. Again use binomcdf or pbinom. Following the procedure above you will have #math.equation(block: false, alt: "binomcdf open parenthesis 20 , .01 , 2 close parenthesis")[$"binomcdf" ( 20 , ".01" , 2 )$] on the TI-83/84 and pbinom(2,20,0.01), with #math.equation(block: false, alt: "P open parenthesis x less than or equal to 2 close parenthesis equals 0.998996")[$P ( x ≤ 2 ) = 0.998996$]. Again there is a really good chance that at most two people in the room will have green eyes. pbinom(2, 20, 0.01) f. At least four means four or more. Find the probability of #emph[x] being greater than or equal to four. That would mean adding up all the probabilities from four to twenty. This would take a long time, so it is better to use the idea of complement. The complement of being greater than or equal to four is being less than four. That would mean being less than or equal to three. Part (e) has the answer for the probability of being less than or equal to three. Just subtract that number from 1. #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 4 close parenthesis equals 1 minus P open parenthesis x less than or equal to 3 close parenthesis equals 1 minus 0.99996 equals 0.00004")[$P ( x ≥ 4 ) = 1 − P ( x ≤ 3 ) = 1 − 0.99996 = 0.00004$] You can also find this answer by doing the following on TI-83/84: #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 4 close parenthesis equals 1 minus P open parenthesis x less than or equal to 3 close parenthesis equals 1 minus binomcdf open parenthesis 20 , .01 , 3 close parenthesis equals 1 minus 0.99996 equals 0.00004")[$P ( x ≥ 4 ) = 1 − P ( x ≤ 3 ) = 1 − " binomcdf " ( 20 , ".01" , 3 ) = 1 − 0.99996 = 0.00004$] on R: #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 4 close parenthesis equals 1 minus P open parenthesis x less than or equal to 3 close parenthesis equals 1 minus pbinom open parenthesis 3 , 20 , .01 close parenthesis equals 1 minus 0.99996 equals 0.0004")[$P ( x ≥ 4 ) = 1 − P ( x ≤ 3 ) = 1 − " pbinom " ( 3 , 20 , ".01" ) = 1 − 0.99996 = 0.0004$] Again, this is very unlikely to happen. ] ] There are other technologies that will compute binomial probabilities. #examplebox("Example 4")[calculating binomial probabilities][ According to the Center for Disease Control (CDC), about 1 in 88 children in the U.S. have been diagnosed with autism ("CDC-data and statistics,," 2013). Suppose you consider a group of 10 children. + State the random variable. + Argue that this is a binomial experiment. + Find the probability that none have autism. + Find the probability that seven have autism. + Find the probability that at least five have autism. + Find the probability that at most two have autism. + Suppose five children out of ten have autism. Is this unusual? What does that tell you? #solutionbox[ a. #emph[x] = number of children with autism b. + There are 10 children, and each child is a trial, so there are a fixed number of trials. In this case, #emph[n] = 10. + If you assume that each child in the group is chosen at random, then whether a child has autism does not affect the chance that the next child has autism. Thus the trials are independent. + Either a child has autism or they do not have autism, so there are two outcomes. In this case, the success is a child has autism. + The probability of a child having autism is 1/88. This is the same for every trial since each child has the same chance of having autism. #math.equation(block: false, alt: "p equals the fraction 1 over 88")[$p = frac(1, 88)$] and #math.equation(block: false, alt: "q equals 1 minus the fraction 1 over 88 equals the fraction 87 over 88")[$q = 1 − frac(1, 88) = frac(87, 88)$]. c. Using the formula: #math.equation(block: true, alt: "P open parenthesis x equals 0 close parenthesis equals sub 10 C sub 0 open parenthesis the fraction 1 over 88 close parenthesis to the power 0 open parenthesis the fraction 87 over 88 close parenthesis to the power 10 minus 0 approximately equals 0.892")[$P ( x = 0 ) = 10 C_(0) attach(( frac(1, 88) ), t: 0) attach(( frac(87, 88) ), t: 10 − 0) ≈ 0.892$] Using the TI-83/84 Calculator: #math.equation(block: true, alt: "P open parenthesis x equals 0 close parenthesis equals binompdf open parenthesis 10 , 1 divided by 88 , 0 close parenthesis approximately equals 0.892")[$P ( x = 0 ) = " binompdf " ( 10 , 1 ÷ 88 , 0 ) ≈ 0.892$] Using R: #math.equation(block: true, alt: "P open parenthesis x equals 0 close parenthesis equals pbinom open parenthesis 0 , 10 , 1 / 88 close parenthesis approximately equals 0.892")[$P ( x = 0 ) = " pbinom " ( 0 , 10 , 1 / 88 ) ≈ 0.892$] pbinom(0, 10, 1/88) d. Using the formula: #math.equation(block: true, alt: "P open parenthesis x equals 7 close parenthesis equals sub 10 C sub 7 open parenthesis the fraction 1 over 88 close parenthesis to the power 7 open parenthesis the fraction 87 over 88 close parenthesis to the power 10 minus 7 approximately equals 0.000")[$P ( x = 7 ) = 10 C_(7) attach(( frac(1, 88) ), t: 7) attach(( frac(87, 88) ), t: 10 − 7) ≈ 0.000$] Using the TI-83/84 Calculator: #math.equation(block: true, alt: "P open parenthesis x equals 7 close parenthesis equals binompdf open parenthesis 10 , 1 divided by 88 , 7 close parenthesis approximately equals 2.84 times 10 to the power minus 12")[$P ( x = 7 ) = " binompdf " ( 10 , 1 ÷ 88 , 7 ) ≈ 2.84 × 10^(− 12)$] Using R: #math.equation(block: true, alt: "P open parenthesis x equals 7 close parenthesis equals dbinom open parenthesis 7 , 10 , 1 / 88 close parenthesis approximately equals 2.84 times 10 to the power minus 12")[$P ( x = 7 ) = op("dbinom") ( 7 , 10 , 1 / 88 ) ≈ 2.84 × 10^(− 12)$] dbinom(7, 10, 1/88) e. Using the formula: #math.equation(block: true, alt: "P open parenthesis x greater than or equal to 5 close parenthesis equals P open parenthesis x equals 5 close parenthesis plus P open parenthesis x equals 6 close parenthesis plus P open parenthesis x equals 7 close parenthesis; plus P open parenthesis x equals 8 close parenthesis plus P open parenthesis x equals 9 close parenthesis plus P open parenthesis x equals 10 close parenthesis; equals sub 10 C sub 5 open parenthesis the fraction 1 over 88 close parenthesis to the power 5 open parenthesis the fraction 78 over 88 close parenthesis to the power 10 minus 5 plus sub 10 C sub 6 open parenthesis the fraction 1 over 88 close parenthesis to the power 6 open parenthesis the fraction 78 over 88 close parenthesis to the power 10 minus 6; plus sub 10 C sub 7 open parenthesis the fraction 1 over")[$P ( x ≥ 5 ) = P ( x = 5 ) + P ( x = 6 ) + P ( x = 7 ) \ + P ( x = 8 ) + P ( x = 9 ) + P ( x = 10 ) \ = 10 C_(5) attach(( frac(1, 88) ), t: 5) attach(( frac(78, 88) ), t: 10 − 5) + 10 C_(6) attach(( frac(1, 88) ), t: 6) attach(( frac(78, 88) ), t: 10 − 6) \ + 10 C_(7) attach(( frac(1, 88) ), t: 7) attach(( frac(78, 88) ), t: 10 − 7) + 10 C_(8) attach(( frac(1, 88) ), t: 8) attach(( frac(78, 88) ), t: 10 − 8) \ + 10 C_(9) attach(( frac(1, 88) ), t: 9) attach(( frac(78, 88) ), t: 10 − 9) + 10 C_(10) attach(( frac(1, 88) ), t: 10) attach(( frac(78, 88) ), t: 10 − 10) \ = 0.000 + 0.000 + 0.000 + 0.000 + 0.000 + 0.000 \ = 0.000$] Using the TI-83/84 Calculator: To use the calculator you need to use the complement. #math.equation(block: true, alt: "P open parenthesis x greater than or equal to 5 close parenthesis equals 1 minus P open parenthesis x less than 5 close parenthesis; equals 1 minus P open parenthesis x less than or equal to 4 close parenthesis; equals 1 minus binomcdf open parenthesis 10 , 1 divided by 88 , 4 close parenthesis; approximately equals 1 minus 0.9999999 equals 0.000")[$P ( x ≥ 5 ) = 1 − P ( x < 5 ) \ = 1 − P ( x ≤ 4 ) \ = 1 − " binomcdf " ( 10 , 1 ÷ 88 , 4 ) \ ≈ 1 − 0.9999999 = 0.000$] Using R: To use R you need to use the complement. #math.equation(block: true, alt: "P open parenthesis x greater than or equal to 5 close parenthesis equals 1 minus P open parenthesis x less than 5 close parenthesis; equals 1 minus P open parenthesis x less than or equal to 4 close parenthesis; equals 1 minus pbinom open parenthesis 4 , 10 , 1 / 88 close parenthesis; approximately equals 1 minus 0.9999999 equals 0.000")[$P ( x ≥ 5 ) = 1 − P ( x < 5 ) \ = 1 − P ( x ≤ 4 ) \ = 1 − " pbinom " ( 4 , 10 , 1 / 88 ) \ ≈ 1 − 0.9999999 = 0.000$] pbinom(4, 10, 1/88) Notice, the answer is given as 0.000, since the answer is less than 0.000. Don’t write 0, since 0 means that the event is impossible to happen. The event of five or more is improbable, but not impossible. f. Using the formula: #math.equation(block: true, alt: "P open parenthesis x less than or equal to 2 close parenthesis equals P open parenthesis x equals 0 close parenthesis plus P open parenthesis x equals 1 close parenthesis plus P open parenthesis x equals 2 close parenthesis; equals sub 10 C sub 0 open parenthesis the fraction 1 over 88 close parenthesis to the power 0 open parenthesis the fraction 78 over 88 close parenthesis to the power 10 minus 0 plus sub 10 C sub 1 open parenthesis the fraction 1 over 88 close parenthesis to the power 1 open parenthesis the fraction 78 over 88 close parenthesis to the power 10 minus 1; plus sub 10 C sub 2 open parenthesis the fraction 1 over 88 close parenthesis squared open parenthesis the fraction 78 over 88 close parenthesis to the power 10 minus 2; equals 0.892 plus 0.103 plus 0.005 greater than 0.")[$P ( x ≤ 2 ) = P ( x = 0 ) + P ( x = 1 ) + P ( x = 2 ) \ = 10 C_(0) attach(( frac(1, 88) ), t: 0) attach(( frac(78, 88) ), t: 10 − 0) + 10 C_(1) attach(( frac(1, 88) ), t: 1) attach(( frac(78, 88) ), t: 10 − 1) \ + 10 C_(2) attach(( frac(1, 88) ), t: 2) attach(( frac(78, 88) ), t: 10 − 2) \ = 0.892 + 0.103 + 0.005 > 0.999$] Using the TI-83/84 Calculator: #math.equation(block: true, alt: "P open parenthesis x less than or equal to 2 close parenthesis equals binomcdf open parenthesis 10 , 1 divided by 88 , 2 close parenthesis approximately equals 0.9998")[$P ( x ≤ 2 ) = " binomcdf " ( 10 , 1 ÷ 88 , 2 ) ≈ 0.9998$] Using R: #math.equation(block: true, alt: "P open parenthesis x less than or equal to 2 close parenthesis equals pbinom open parenthesis 2 , 10 , 1 / 88 close parenthesis approximately equals 0.9998")[$P ( x ≤ 2 ) = " pbinom " ( 2 , 10 , 1 / 88 ) ≈ 0.9998$] pbinom(2, 10, 1/88) g. Since the probability of five or more children in a group of ten having autism is much less than 5%, it is unusual to happen. If this does happen, then one may think that the proportion of children diagnosed with autism is actually more than 1/88. ] ] === Homework #notebox("Your Turn", rgb("#8a94a6"), rgb("#556666"), rgb("#f7f8fa"))[ + Suppose a random variable, #emph[x], arises from a binomial experiment. If #emph[n] = 14, and #emph[p] = 0.13, find the following probabilities using the binomial formula. + #emph[P] (x=5) + #emph[P] (x=8) + #emph[P] (x=12) + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 4 close parenthesis")[$P ( x ≤ 4 )$] + #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 8 close parenthesis")[$P ( x ≥ 8 )$] + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 12 close parenthesis")[$P ( x ≤ 12 )$] + Suppose a random variable, #emph[x], arises from a binomial experiment. If #emph[n] = 22, and #emph[p] = 0.85, find the following probabilities using the binomial formula. + #emph[P] (x=18) + #emph[P] (x=5) + #emph[P] (x=20) + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 3 close parenthesis")[$P ( x ≤ 3 )$] + #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 18 close parenthesis")[$P ( x ≥ 18 )$] + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 20 close parenthesis")[$P ( x ≤ 20 )$] + Suppose a random variable, #emph[x], arises from a binomial experiment. If #emph[n] = 10, and #emph[p] = 0.70, find the following probabilities using the binomial formula. + #emph[P] (x=2) + #emph[P] (x=8) + #emph[P] (x=7) + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 3 close parenthesis")[$P ( x ≤ 3 )$] + #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 7 close parenthesis")[$P ( x ≥ 7 )$] + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 4 close parenthesis")[$P ( x ≤ 4 )$] + Suppose a random variable, #emph[x], arises from a binomial experiment. If #emph[n] = 6, and #emph[p] = 0.30, find the following probabilities using the binomial formula. + #emph[P] (x=1) + #emph[P] (x=5) + #emph[P] (x=3) + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 3 close parenthesis")[$P ( x ≤ 3 )$] + #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 5 close parenthesis")[$P ( x ≥ 5 )$] + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 4 close parenthesis")[$P ( x ≤ 4 )$] + Suppose a random variable, #emph[x], arises from a binomial experiment. If #emph[n] = 17, and #emph[p] = 0.63, find the following probabilities using the binomial formula. + #emph[P] (x=8) + #emph[P] (x=15) + #emph[P] (x=14) + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 12 close parenthesis")[$P ( x ≤ 12 )$] + #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 10 close parenthesis")[$P ( x ≥ 10 )$] + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 7 close parenthesis")[$P ( x ≤ 7 )$] + Suppose a random variable, #emph[x], arises from a binomial experiment. If #emph[n] = 23, and #emph[p] = 0.22, find the following probabilities using the binomial formula. + #emph[P] (x=21) + #emph[P] (x=6) + #emph[P] (x=12) + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 14 close parenthesis")[$P ( x ≤ 14 )$] + #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 17 close parenthesis")[$P ( x ≥ 17 )$] + #math.equation(block: false, alt: "P open parenthesis x less than or equal to 9 close parenthesis")[$P ( x ≤ 9 )$] + Approximately 10% of all people are left-handed ("11 little-known facts," 2013). Consider a grouping of fifteen people. + State the random variable. + Argue that this is a binomial experiment Find the probability that + None are left-handed. + Seven are left-handed. + At least two are left-handed. + At most three are left-handed. + At least seven are left-handed. + Seven of the last 15 U.S. Presidents were left-handed. Is this unusual? What does that tell you? + According to an article in the American Heart Association’s publication Circulation, 24% of patients who had been hospitalized for an acute myocardial infarction did not fill their cardiac medication by the seventh day of being discharged (Ho, Bryson & Rumsfeld, 2009). Suppose there are twelve people who have been hospitalized for an acute myocardial infarction. + State the random variable. + Argue that this is a binomial experiment Find the probability that + All filled their cardiac medication. + Seven did not fill their cardiac medication. + None filled their cardiac medication. + At most two did not fill their cardiac medication. + At least three did not fill their cardiac medication. + At least ten did not fill their cardiac medication. + Suppose of the next twelve patients discharged, ten did not fill their cardiac medication, would this be unusual? What does this tell you? + Eyeglassomatic manufactures eyeglasses for different retailers. In March 2010, they tested to see how many defective lenses they made, and there were 16.9% defective lenses due to scratches. Suppose Eyeglassomatic examined twenty eyeglasses. + State the random variable. + Argue that this is a binomial experiment Find the probability that + None are scratched. + All are scratched. + At least three are scratched. + At most five are scratched. + At least ten are scratched. + Is it unusual for ten lenses to be scratched? If it turns out that ten lenses out of twenty are scratched, what might that tell you about the manufacturing process? + The proportion of brown M&M’s in a milk chocolate packet is approximately 14% (Madison, 2013). Suppose a package of M&M’s typically contains 52 M&M’s. + State the random variable. + Argue that this is a binomial experiment Find the probability that + Six M&M’s are brown. + Twenty-five M&M’s are brown. + All of the M&M’s are brown. + Would it be unusual for a package to have only brown M&M’s? If this were to happen, what would you think is the reason? Answer 1. a. P(x=5) = 0.0212, b. P(x=8) = #math.equation(block: false, alt: "1.062 times 10 to the power minus 4")[$1.062 × 10^(− 4)$], c. P(x=12) = #math.equation(block: false, alt: "1.605 times 10 to the power minus 9")[$1.605 × 10^(− 9)$], d. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 4 close parenthesis equals 0.973")[$P ( x ≤ 4 ) = 0.973$], e. #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 8 close parenthesis equals 1.18 times 10 to the power minus 4")[$P ( x ≥ 8 ) = 1.18 × 10^(− 4)$], f. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 12 close parenthesis equals 0.99999")[$P ( x ≤ 12 ) = 0.99999$] 3. a. #math.equation(block: false, alt: "P open parenthesis x equals 2 close parenthesis equals 0.0014")[$P ( x = 2 ) = 0.0014$], b. #math.equation(block: false, alt: "P open parenthesis x equals 8 close parenthesis equals 0.2335")[$P ( x = 8 ) = 0.2335$], c. #math.equation(block: false, alt: "P open parenthesis x equals 7 close parenthesis equals 0.2668")[$P ( x = 7 ) = 0.2668$], d. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 3 close parenthesis equals 0.0106")[$P ( x ≤ 3 ) = 0.0106$], e. #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 7 close parenthesis equals 0.6496")[$P ( x ≥ 7 ) = 0.6496$], f. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 4 close parenthesis equals 0.0473")[$P ( x ≤ 4 ) = 0.0473$] 5. a. #math.equation(block: false, alt: "P open parenthesis x equals 8 close parenthesis equals 0.0784")[$P ( x = 8 ) = 0.0784$], b. #math.equation(block: false, alt: "P open parenthesis x equals 15 close parenthesis equals 0.0182")[$P ( x = 15 ) = 0.0182$], c. #math.equation(block: false, alt: "P open parenthesis x equals 14 close parenthesis equals 0.0534")[$P ( x = 14 ) = 0.0534$], d. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 12 close parenthesis equals 0.8142")[$P ( x ≤ 12 ) = 0.8142$], e. #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 10 close parenthesis equals 0.7324")[$P ( x ≥ 10 ) = 0.7324$], f. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 7 close parenthesis equals 0.0557")[$P ( x ≤ 7 ) = 0.0557$] 7. a. See solutions, b. See solutions, c. P(x=0) = 0.2059, d. #math.equation(block: false, alt: "P open parenthesis x equals 7 close parenthesis equals 2.770 times 10 to the power minus 4")[$P ( x = 7 ) = 2.770 × 10^(− 4)$], e. #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 2 close parenthesis equals 0.4510")[$P ( x ≥ 2 ) = 0.4510$], f. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 3 close parenthesis equals 0.944")[$P ( x ≤ 3 ) = 0.944$], g.#math.equation(block: false, alt: "P open parenthesis x greater than or equal to 7 close parenthesis equals 3.106 times 10 to the power minus 4")[$P ( x ≥ 7 ) = 3.106 × 10^(− 4)$], h. See solutions 9. a. See solutions, b. See solutions, c. #math.equation(block: false, alt: "P open parenthesis x equals 0 close parenthesis equals 0.0247")[$P ( x = 0 ) = 0.0247$], d. #math.equation(block: false, alt: "P open parenthesis x equals 20 close parenthesis equals 3.612 times 10 to the power minus 16")[$P ( x = 20 ) = 3.612 × 10^(− 16)$], e. #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 3 close parenthesis equals 0.6812")[$P ( x ≥ 3 ) = 0.6812$], f. #math.equation(block: false, alt: "P open parenthesis x less than or equal to 5 close parenthesis equals 0.8926")[$P ( x ≤ 5 ) = 0.8926$], g. #math.equation(block: false, alt: "P open parenthesis x greater than or equal to 10 close parenthesis equals 6.711 times 10 to the power minus 4")[$P ( x ≥ 10 ) = 6.711 × 10^(− 4)$], h. See solutions ]