#Read the data:
a <- read.table("http://www.stat.ucla.edu/~nchristo/statistics13/beeswax.txt", header=TRUE)

#Question (a):
plot(ecdf(a$Hydrocarbon), verticals=TRUE, do.points=FALSE)

hist(a$Hydrocarbon)

boxplot(a$Hydrocarbon)


#Question (b):
quantile(a$Hydrocarbon, probs=c(0.10,0.25,0.50,0.75,0.90))

quantile(a$Hydrocarbon, probs=c(0.10,0.25,0.50,0.75,0.90))
   10%    25%    50%    75%    90% 
13.676 14.070 14.570 15.115 15.470 


#Question ©:
#1% dilution will almost never be detected.  Here are the calculations.  Please get the final answer in the calculations below!
0.99*0.13676 + 0.01*0.85
0.99*0.14070 + 0.01*0.85
0.99*0.14570 + 0.01*0.85
0.99*0.15115 + 0.01*0.85
0.99*0.15470 + 0.01*0.85

#3% dilution will be detected about half the time.  Here are the calculations.
0.97*0.13676 + 0.03*0.85
0.97*0.14070 + 0.03*0.85
0.97*0.14570 + 0.03*0.85
0.97*0.15115 + 0.03*0.85
0.97*0.15470 + 0.03*0.85

#5% dilution will be detected almost all the time.  Here are the calculations.
0.95*0.13676 + 0.05*0.85
0.95*0.14070 + 0.05*0.85
0.95*0.14570 + 0.05*0.85
0.95*0.15115 + 0.05*0.85
0.95*0.15470 + 0.05*0.85