如何使用R中的“bootstrap函数”计算置信区间

我正在尝试计算R中的置信区间。由于一些特殊原因，我必须使用“bootstrap”包中的函数来计算。（这意味着我不能使用“boot”包中的函数。）

这是我的代码。

我正在做的是尝试计算 Pearson 相关系数，然后应用 Bootstrap 方法（B = 100）来获得相关系数的估计。但我不知道如何构建 95% 的置信区间。

library(bootstrap) 
data('law')

set.seed(1)
theta <- function(ind) {
  cor(law[ind, 1], law[ind, 2], method = "pearson")
  }
law.boot <- bootstrap(1:15, 100, theta) 
# sd(law$thetastar)
percent.95 <- function(x) {
  quantile(x,  .95)
  }
law.percent.95 <- bootstrap(1:15, 100, theta, func=percent.95)

抱歉，如果我没有说清楚或标记错误的标签。 两次抱歉没有生成数据集（现在提供了），并感谢罗兰教授指出这一点。非常感谢！

通常，在引导之后，我们使用 2.5% 和 97.5% 百分位数作为 95% 置信区间（因为我们从每一边减去 α/2=.025）。也可以看看@索塔尔的 answer以及答案下的评论。

R <- 1e5 - 1  ## number of bootstrap replications
est <- with(law, cor(lsat, gpa))  ## naïve correlation

theta <- function(ind) cor(law[ind, 1], law[ind, 2], method="pearson")
set.seed(1)
B1 <- bootstrap::bootstrap(seq(nrow(law)), R, theta) 
(ci1 <- c(estimate=est, quantile(B1$thetastar, c(.025, .975))))
#  estimate      2.5%     97.5% 
# 0.7763745 0.4594845 0.9620884 

这是从头开始的替代方法：

theta2 <- function(x) with(x, cor(lsat, gpa))
set.seed(1)
B2 <- replicate(R, theta2(law[sample(nrow(law), nrow(law), replace=TRUE), ]))
(ci2 <- c(estimate=est, quantile(B2, c(.025, .975))))
#  estimate      2.5%     97.5% 
# 0.7763745 0.4607644 0.9617970 

最后一种方法是使用boot包中有一个boot.ci功能：

theta3 <- function(data, k) cor(data[k, ])[1,2]
set.seed(1)
B3 <- boot::boot(law, theta3, R=R)
(ci3 <- c(est, boot::boot.ci(B3, type='perc')$percent[4:5]))
# [1] 0.7763745 0.4593727 0.9620923