Assignment 5

Exercise 1 (15 points)

We will now derive the probability that a given observation is part of a bootstrap sample. Suppose that we obtain a bootstrap sample from a set of \(n\) observations.

1. What is the probability that the first bootstrap observation is not the \(j\)th observation from the original sample? Justify your answer.
2. What is the probability that the second bootstrap observation is not the \(j\)th observation from the original sample? Justify your answer.
3. Argue that the probability that the \(j\)th observation is not in the bootstrap sample is \((1-1/n)^n\).
4. When \(n = 5\) what is the probability that the \(j\)th observation is in the bootstrap sample?
5. When \(n = 100\) what is the probability that the \(j\)th observation is in the bootstrap sample?
6. When \(n = 10,000\) what is the probability that the \(j\)th observation is in the bootstrap sample?
7. Create a plot that displays, for each integer value of \(n\) from 1 to 100,000 the probability that the \(j\)th observation is in the bootstrap sample. Comment on what you observe.
8. We will now investigate numerically that a bootstrap sample of size \(n = 100\) contains the \(j\)th observation. Here \(j = 4\). We repeatedly create bootstrap samples, and each time we record whether or not the fourth observation is contained in the bootstrap sample.

set.seed(  ) # set a seed here

store <- integer(10000)

for (i in seq_along(store)) {
  store[i] <- sum(sample(seq_len(100), replace = TRUE) == 4) > 0
}

mean(store)

Comment on the results obtained.

Exercise 2 (10 points)

Suppose that we use some statistical learning method to make a prediction for the response \(Y\) for a particular value of the predictor \(X\).

1. Carefully describe how we might estimate the standard deviation of our prediction.
2. Is this procedure depends on what statistical learning method we are using?

Exercise 3 (10 points)

This exercise should be answered using the Default data set, which is part of the LSLR package. If you don’t have it installed already you can install it with

install.packages("ISLR")

To load the data set run the following code

library(ISLR)
data("Default")

1. Use the parsnip package to fit a logistic regression on the default data set. default is the response and income and balance are the predictors. Then use summary() on the fitted model to determine the estimated standard errors for the coefficients associated with income and balance. Comment on the estimated standard errors.
2. Use the bootstraps() function from the rsample package to generate 25 bootstraps of Default.
3. Run the following code. Change boots to the name of the bootstrapping object created in the previous question. This will take a minute or two to run. Comment on

# This function takes a `bootstrapped` split object, and fits a logistic model
fit_lr_on_bootstrap <- function(split) {
    logistic_reg() %>%
    set_engine("glm") %>%
    set_mode("classification") %>%
    fit(default ~ income + balance, analysis(split))
}

# This code uses `map()` to run the model inside each split. Then it used
# `tidy()` to extract the model estimates the parameter
boot_models <-
  boots %>% 
  mutate(model = map(splits, fit_lr_on_bootstrap),
         coef_info = map(model, tidy))

# This code extract the estimates for each model that was fit
boot_coefs <-
  boot_models %>% 
  unnest(coef_info)

# This code calculates the standard deviation of the estimate
sd_estimates <- boot_coefs %>%
  group_by(term) %>%
  summarise(std.error_est = sd(estimate))
sd_estimates

1. Comment on the estimated standard errors obtained using the summary() function on the first model and the estimated standard errors you found using the above code.