## Section 3.4 Simulations

In the previous section we saw how to apply the binomial formula to find the probability of exactly \(k\) successes in \(n\) independent trials when a success has probability \(p\text{.}\) Sometimes we have a problem we want to solve but we don't know the appropriate formula, or even worse, a formula may not exist. In this case, one common approach is to estimate the probability using *simulations*.

You may already be familiar with simulations. Want to know the probability of rolling a sum of 7 with a pair of dice? Roll a pair of dice many, many, many times and see what proportion of times the sum was 7. The more times you roll the pair of dice, the better the estimate will tend to be. Of course, such experiments can be time consuming or even infeasible.

In this section, we consider simulations using random numbers. Random numbers (or technically, *psuedo-random numbers*) can be produced using a calculator or computer. Random digits are produced such that each digit, *0-9*, is equally likely to come up in each spot. You'll find that occasionally we may have the same number in a row — sometimes multiple times — but in the long run, each digit should appear 1/10th of the time.

Column | |||||||||

Row | 1-5 | 6-10 | 11-15 | 16-20 | |||||

1 | 43087 | 41864 | 51009 | 39689 | |||||

2 | 63432 | 72132 | 40269 | 56103 | |||||

3 | 19025 | 83056 | 62511 | 52598 | |||||

4 | 85117 | 16706 | 31083 | 24816 | |||||

5 | 16285 | 56280 | 01494 | 90240 | |||||

6 | 94342 | 18473 | 50845 | 77757 | |||||

7 | 61099 | 14136 | 39052 | 50235 | |||||

8 | 37537 | 58839 | 56876 | 02960 | |||||

9 | 04510 | 16172 | 90838 | 15210 | |||||

10 | 27217 | 12151 | 52645 | 96218 | |||||

Mika's favorite brand of cereal is running a special where 20% of the cereal boxes contain a prize. Mika really wants that prize. If her mother buys 6 boxes of the cereal over the next few months, what is the probability Mika will get a prize?

To solve this problem using simulation, we need to be able to assign digits to outcomes. Each box should have a 20% chance of having a prize and an 80% chance of not having a prize. Therefore, a valid assignment would be:

Of the ten possible digits (*0, 1, 2, ..., 8, 9*), two of them, i.e. 20% of them, correspond to winning a prize, which exactly matches the odds that a cereal box contains a prize.

In Mika's simulation, one trial will consist of 6 boxes of cereal, and therefore a trial will require six digits (each digit will correspond to one box of cereal). We will repeat the simulation for 20 trials. Therefore we will need 20 sets of 6 digits. Let's begin on row 1 of the random digit table, shown in Table 3.4.1. If a trial consisted of 5 digits, we could use the first 5 digits going across: *43087*. Because here a trial consists of 6 digits, it may be easier to read down the table, rather than read across. We will let trial 1 consist of the first 6 digits in column 1 (*461819*), trial 2 consist of the first 6 digits in column 2 (*339564*), etc. For this simulation, we will end up using the first 6 rows of each of the 20 columns.

In trial 1, there are two *1*'s, so we record that as a success; in this trial there were actually two prizes. In trial 2 there were no *0*'s or *1*'s, therefore we do not record this as a success. In trial 3 there were three prizes, so we record this as a success. The rest of this exercise is left as a Guided Practice problem for you to complete.

###### Guided Practice 3.4.3

Finish the simulation above and report the estimate for the probability that Mika will get a prize if her mother buys 6 boxes of cereal where each one has a 20% chance of containing a prize.^{ 1 }The trials that contain at least one 0 or 1 and therefore are successes are trials: 1, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, and 20. There were 17 successes among the 20 trials, so our estimate of the probability based on this simulation is 17/20 = 0.85.

###### Guided Practice 3.4.4

In the previous example, the probability that a box of cereal contains a prize is 20%. The question presented is equivalent to asking, what is the probability of getting at least one prize in six randomly selected boxes of cereal. This probability question can be solved explicitly using the method of complements. Find this probability. How does the estimate arrived at by simulation compare to this probability?^{ 2 }The true probability is given by \(1 - P(\text{ no prizes in six boxes } ) = 1- 0.8^6 = 0.74\text{.}\) The estimate arrived at by simulation was 11% too high. Note: We only repeated the simulation 20 times. If we had repeated it 1000 times, we would (very likely) have gotten an estimate closer to the true probability.

We can also use simulations to estimate quantities other than probabilities. Consider the following example.

###### Example 3.4.5

Let's say that instead of buying exactly 6 boxes of cereal, Mika's mother agrees to buy boxes of this cereal *until* she finds one with a prize. On average, how many boxes of cereal would one have to buy until one gets a prize?

For this question, we can use the same digit assignment. However, our stopping rule is different. Each trial may require a different number of digits. For each trial, the stopping rule is: look at digits until we encounter a *0* or a *1*. Then, record how many digits/boxes of cereal it took. Repeat the simulation for 20 trials, and then average the numbers from each trial.

Let's begin again at row 1. We can read across or down, depending upon what is most convenient. Since there are 20 columns and we want 20 trials, we will read down the columns. Starting at column 1, we count how many digits (boxes of cereal) we encounter until we reach a *0* or *1* (which represent a prize). For trial 1 we see *461*, so we record 3. For trial 2 we see *3395641*, so we record 7. For trial 3, we see *0*, so we record 1. The rest of this exercise is left as a Guided Practice problem for you to complete.

###### Guided Practice 3.4.6

Finish the simulation above and report your estimate for the average number of boxes of cereal one would have to buy until encountering a prize, where the probability of a prize in each box is 20%.^{ 3 }For the 20 trials, the number of digits we see until we encounter a *0* or *1* is: 3,7,1,4,9, 4,1,2,4,5, 5,1,1,1,3, 8,5,2,2,6. Now we take the average of these 20 numbers to get 74/20 = 3.7.

###### Example 3.4.7

Now, consider a case where the probability of interest is not 20%, but rather 28%. Which digits should correspond to success and which to failure?

This example is more complicated because with only 10 digits, there is no way to select exactly 28% of them. Therefore, each observation will have to consist of *two* digits. We can use two digits at a time and assign pairs of digits as follows:

###### Guided Practice 3.4.8

Assume the probability of winning a particular casino game is 45%. We want to carry out a simulation to estimate the probability that we will win at least 5 times in 10 plays. We will use 30 trials of the simulation. Assign digits to outcomes. Also, how many total digits will we require to run this simulation?^{ 4 }One possible assignment is: *00-44* \(\rightarrow\) win and *45-99* \(\rightarrow\) lose. Each trial requires 10 pairs of digits, so we will need 30 sets of 10 pairs of digits for a total of \(30 \times 10 \times 2 = 600\) digits.

###### Guided Practice 3.4.9

Assume carnival spinner has 7 slots. We want to carry out a simulation to estimate the probability that we will win at least 10 times in 60 plays. Repeat 100 trials of the simulation. Assign digits to outcomes. Also, how many total digits will we require to run this simulation?^{ 5 }Note that \(1/7 = 0.142857...\) This makes it tricky to assign digits to outcomes. The best approach here would be to exclude some of the digits from the simulation. We can assign *0* to success and *1-6* to failure. This corresponds to a \(1/7\) chance of getting a success. If we encounter a *7*, *8*, or *9*, we will just skip over it. Because we don't know how many *7*, *8*, or *9*'s we will encounter, we do not know how many total digits we will end up using for the simulation. (If you want a challenge, try to estimate the total number of digits you would need.)

Does anyone perform simulations like this? Sort of. Simulations are used a lot in statistics, and these often require the same principles covered in this section to properly set up those simulations. The difference is in implementation after the setup. Rather than use a random number table, a statistician will write a program that uses a pseudo-random number generator in a computer to run the simulations very quickly — often times millions of trials each second, which provides much more accurate estimates than running a couple dozen trials by hand.