自举比例置信区间

Question

如果可能的话，我想使用启动包生成比例的置信区间。 我有一个矢量，我想设置一个阈值，然后计算低于指定水平的比例。 之后我想使用boot包中的bootstrap函数来计算比例的置信区间。 我到目前为止的简单例子：

library(boot)

vec <- abs(rnorm(1000)*10) #generate example vector

data_to_tb <- vec

tb <- function(data) {
  sum(data < 10, na.rm = FALSE)/length(data) #function for generating the proportion
}

tb(data_to_tb)

boot(data = data_to_tb, statistic = tb, R = 999)
quantile(boot.out$t, c(.025,.975))

但是，我收到此错误消息：

> boot(data = data_to_tb, statistic = tb, R = 999)
Error in statistic(data, original, ...) : unused argument (original)

我无法让它工作，帮助赞赏

Answer 1

你的问题是你的函数tb - 它需要两个参数。来自帮助文件?boot

  

统计信息应用于数据时返回向量的函数   包含感兴趣的统计数据。当sim =＆＃34;参数＆＃34;时，   统计的第一个参数必须是数据。对于每个复制品a   将传递由ran.gen返回的模拟数据集。在所有其他   案件统计必须至少有两个论点。