Question: Consider the following bootstrap:

```
library(MASS)
library(boot)

# c)
set.seed(1)
boot.fn= function(data, index) mean(data[index])
output=boot(Boston$medv, boot.fn, 1000)

```

If we run `print(output)`, we get 

```
Call:
boot(data = Boston$medv, statistic = boot.fn, R = 1000)

Bootstrap Statistics :
    original      bias    std. error
t1* 22.53281 0.008517589   0.4119374

```

However, when I examine the `output` object, I cannot find values representing the bootstrap statistics. Where is `original`, `bias` and `std. error` in the actual `output` object returned by boot?

Answer: 
This command seems to work:

```
apply(output$t,2,sd)[1]

```
```
output_tab <- t(rbind(
  # this grabs the original coefficients
  output$t0, 
  # bootstrapped coefficients' standard errors
  apply(output$t, 2, function(x) sd(x, ra.rm = TRUE))  
))  

output_tab
#         [,1]      [,2]
#[1,] 22.53281 0.4106622

```