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)