The Problem
Example # 353: Estimate the success rate for the population using a confidence interval. To estimate this difference, we collect data. The data are a series of “Success” and “Failure” values. For our sample, the data are
“Success”, “Success”, “Failure”, “Success”, “Failure”, “Success”, “Failure”, “Success”, “Failure”, “Failure”, “Success”, “Success”, “Success”, “Failure”, “Failure”, “Failure”, “Failure”, “Success”, “Success”, “Failure”, “Success”, “Success”, “Success”, “Failure”, “Success”, “Failure”, “Success”, “Success”, “Success”, “Success”, “Failure”, “Success”, “Success”, “Failure”, “Failure”, “Success”, “Failure”, “Failure”, “Failure”, “Failure”, “Failure”
With this information, calculate the endpoints of the symmetric 99% confidence interval.
Information given:
To summarize the above, the values of import are:
Summary statistics from the problem
| \( x \)
| = |
21 |
|---|
| \( n \)
| = |
41 |
|---|
| \( \hat{p} \)
| = |
0.5122 |
|---|
| \( \alpha \)
| = |
0.01 |
|---|
Note that there is no value given for p0. This is because confidence intervals are based solely on the data, and not on any hypothesized values.
It may be helpful if you calculate these values yourself. Once you have, you can check your answers by hovering your mouse over the grey spaces to see if you calculated them correctly.
Your Answer
![[Correct!]](images/correct.png "You are correct. Well done!")
You got the correct endpoints to the 99% confidence: (0.3112, 0.7132). Congratulations!
![[Incorrect!]](images/incorrect.png "You are not correct.")
Unfortunately, your answer was not correct. Either try again or click on “Show Solution” below to see how to obtain the correct answer.
Assistance
Hide Solution
$$ \begin{align}
\text{Confidence Limits} &= \hat{p} \pm Z(\alpha/2) \sqrt{\frac{\hat{p}(1-\hat{p})}{n}} \\[3em]
&= 0.5122 \pm Z(0.01/2) \sqrt{\frac{ 0.5122 \left(1- 0.5122\right)}{ 41}} \\[1em]
&= 0.5122 \pm \left(2.575\right) \sqrt{\frac{ 0.2499}{ 41}} \\[1em]
&= 0.5122 \pm \left(2.575\right) \sqrt{ 0.0061} \\[1em]
&= 0.5122 \pm \left(2.575\right) 0.078064 \\[1em]
&= 0.5122 \pm 0.201014 \end{align}
$$
Thus, we are 99% confident that the success rate in the population is between 0.3112 and 0.7132.
Note that 0.201014 is the margin of error, which is usually symbolized as E. So, for a sample like this of size n = 41, polling companies will report the results as “51.2% plus or minus 20.1 points.” As you may expect, larger sample sizes will produce smaller margins of error.
Hide the R Code
This formulation of the confidence interval is pedagogically simple to understand. That is why it is used in introductory textbooks. It is actually based on the Normal approximation to the Binomial distribution. There are several improvements to the test. For these reasons, R does not have a built-in function to calculate this version o fthe confidence interval. The following code echoes the above calculations to provide the endpoints of the confidence interval.
Copy and paste the following code into your R script window, then run it from there.
sample = c("Success", "Success", "Failure", "Success", "Failure", "Success", "Failure", "Success", "Failure", "Failure", "Success", "Success", "Success", "Failure", "Failure", "Failure", "Failure", "Success", "Success", "Failure", "Success", "Success", "Success", "Failure", "Success", "Failure", "Success", "Success", "Success", "Success", "Failure", "Success", "Success", "Failure", "Failure", "Success", "Failure", "Failure", "Failure", "Failure", "Failure")
x = sum(sample=="Success")
n = length(sample)
phat = x/n
alpha = 0.01
se2 = phat*(1-phat)/n
lcl = phat + qnorm(alpha/2)*sqrt(se2)
ucl = phat - qnorm(alpha/2)*sqrt(se2)
lcl; ucl
In the R output, the endpoints of the confidence interval are the numbers output after running the last line. Note that R will give you calculations that are more accurate and more precise than doing the calculations by hand.
Hide the Excel Code
This formulation of the confidence interval is pedagogically simple to understand. That is why it is used in introductory textbooks. It is actually based on the Normal approximation to the Binomial distribution. There are several improvements to the test. For such reasons, Excel does not have a built-in function to perform these calculations. The following code echoes the above calculations to provide the confidence interval.
Copy and paste the following code into your Excel window, making sure the value sample ends up in A1 after pasting.
How to calculate the test statistic in Excel.
| sample | | | |
| Success |
|
alpha: |
0.01 |
| Success |
|
x: |
=COUNTIF(A:A,"Success") |
| Failure |
|
n: |
=COUNTIF(A:A,"Success")+COUNTIF(A:A,"Failure") |
| Success |
|
p-hat: |
=D3/D4 |
| Failure |
|
lcl: |
=D5+NORM.S.INV(D2/2)*sqrt(D5*(1-D5)/D4) |
| Success |
|
ucl: |
=D5-NORM.S.INV(D2/2)*sqrt(D5*(1-D5)/D4) |
| Failure |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Success |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Success |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
| Failure |
|
|
|
The endpoints of the 99% confidence interval are the numbers calculated in cells D5 and D6. Again, when you paste this code into Excel, make sure that you start the pasting in cell A1. To help with that, you may want to also copy this notice. It seems to help.