library(boot)

boot_returns <- function(data, indices){
  d <- data[indices]
  Return.calculate(d, method="log")
}

bret <- boot(end_eq, boot_returns, R=1000)

dim(bret$t)
bretMeans <- colMeans(bret$t, na.rm = TRUE)