I estimate a 2-dimensional function in R and want to compute confidence regions for the estimate using bootstrap. That is, I replicate the procedure m times and compute the