prostate <- read.delim("../../data/prostate2.txt")
F <- ecdf(prostate$cavol)
plot(F,verticals=TRUE,ylab=expression(hat(F)(x)),xlab="Cancer volume",main="",pch="")

require(boot)
var.boot <- function(x,ind){var(x[ind])}
out <- boot(prostate$cavol,var.boot,1000)
boot.ci(out,type="perc")
hist(out$t,col="gray",xlab=expression(hat(theta)*'*'))

source("http://web.as.uky.edu/statistics/users/pbreheny/764-F11/notes/fun.R")
sqrt(var.interval(out$t))