如果可能的话,我想使用启动包产生比例的置信区间 . 我有一个矢量,我想设置一个阈值,然后计算低于指定水平的比例 . 之后我想使用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)
我不能让它工作,帮助赞赏
1 回答
你的问题是你的函数
tb
- 它需要两个参数 . 来自帮助文件?boot