### CONFIDENCE INTERVAL FOR P 

### DESCRIPTION:
###    This program will generate the upper and lower bound
###    of a confidence interval for P.

### REQUIRED ARGUMENTS:
###    x = number of successes.
### 	n = number of trials.
###   conf.level = confidence level

conf.level <- .95

x1 <- 45
n1 <- 100
CI1 <- prop.test(x1, n1, conf.level=conf.level)$conf.int
pi1 <- x1/n1
l1 <- x1/n1 - CI1[1]
u1 <- CI1[2] - x1/n1

x2 <- 10
n2 <- 100
CI2 <- prop.test(x2, n2, conf.level=conf.level)$conf.int
pi2 <- x2/n2
l2 <- x2/n2 - CI2[1]
u2 <- CI2[2] - x2/n2

x3 <- 15
n3 <- 100
CI3 <- prop.test(x3, n3, conf.level=conf.level)$conf.int
pi3 <- x3/n3
l3 <- x3/n3 - CI3[1]
u3 <- CI3[2] - x3/n3

x4 <- 45
n4 <- 100
CI4 <- prop.test(x4, n4, conf.level=conf.level)$conf.int
pi4 <- x4/n4
l4 <- x4/n4 - CI4[1]
u4 <- CI4[2] - x4/n4

x5 <- 90
n5 <- 100
CI5 <- prop.test(x5, n5, conf.level=conf.level)$conf.int
pi5 <- x5/n5
l5 <- x5/n5 - CI5[1]
u5 <- CI5[2] - x5/n5
x6 <- 45
n6 <- 100
CI6 <- prop.test(x6, n6, conf.level=conf.level)$conf.int
pi6 <- x6/n6
l6 <- x6/n6 - CI6[1]
u6 <- CI6[2] - x6/n6

x7 <- 10
n7 <- 100
CI7 <- prop.test(x7, n7, conf.level=conf.level)$conf.int
pi7 <- x7/n7
l7 <- x7/n7 - CI7[1]
u7 <- CI7[2] - x7/n7

x8 <- 15
n8 <- 100
CI8 <- prop.test(x8, n8, conf.level=conf.level)$conf.int
pi8 <- x8/n8
l8 <- x8/n8 - CI8[1]
u8 <- CI8[2] - x8/n8

x9 <- 45
n9 <- 100
CI9 <- prop.test(x9, n9, conf.level=conf.level)$conf.int
pi9 <- x9/n9
l9 <- x9/n9 - CI9[1]
u9 <- CI9[2] - x9/n9

x10 <- 90
n10 <- 100
CI10 <- prop.test(x10, n10, conf.level=conf.level)$conf.int
pi10 <- x10/n10
l10 <- x10/n10 - CI10[1]
u10 <- CI10[2] - x10/n10


### settings for plot
x <- c(1,2,3,4,5,6,7,8,9,10)
y <- c(pi1,pi2,pi3,pi4,pi5,pi6,pi7,pi8,pi9,pi10)
lower <- c(l1,l2,l3,l4,l5,l6,l7,l8,l9,l10)
upper <- c(u1,u2,u3,u4,u5,u6,u7,u8,u9,u10)

error.bar(x, y, lower, upper,main="CONFIDENCE INTERVALS FOR Pi",
		sub="",xlab="Car Company",ylab="percent",pch=15,col=7,col=1)


print(CI1)
print(CI2)
print(CI3)
print(CI4)
print(CI5)