Data <- read.delim("nerve-pulse.txt")
Fhat <- ecdf(Data$time)
Fhat(0.1)
Fhat(0.6)

plot(Fhat,verticals=TRUE,do.points=FALSE,xlab="Time",ylab=expression(hat(F)(x)),main="")

x <- Data$time
n <- length(x)
x <- sort(Data$time)
l1 <- u1 <- l2 <- u2 <- numeric(n)
for (i in 1:n) {
  ci <- binom.test(sum(x <= x[i]),n)$conf.int
  l1[i] <- ci[1]
  u1[i] <- ci[2]
  epsilon <- sqrt((1/(2*n))*log(2/.05))
  l2[i] <- max(Fhat(x[i])-epsilon,0)
  u2[i] <- min(Fhat(x[i])+epsilon,1)
}
plot(Fhat,verticals=TRUE,do.points=FALSE,xlab="Time",ylab=expression(hat(F)(x)),main="")
polygon.step(x,l2[-n],u2[-n],col="gray85",border=FALSE)
polygon.step(x,l1[-n],u1[-n],col="gray70",border=FALSE)
plot(Fhat,verticals=TRUE,do.points=FALSE,add=TRUE)