我正试图在this post的引导程序模拟中完成代码,但有两个主要障碍。首先,我使用自举的几次,我没有使用{boot}
包,依赖replicate
之类的函数,或者简单地循环。第二,我的问题的主要目标,我不知道函数xor
作为输入所接受的内容,以及它最终的作用:
n <- 1000
dat <- c("A", rep("B", n-1))
indicator <- function(x, ndx) xor("A"%in%x[ndx], TRUE)
我看到它会评估"A"
是否在x
中,如果它实际存在则返回FALSE
。问题伴随输入ndx
。它可能代表“索引”,这很可能与随后的引导有关:
p_hat <- function(dat, m=1e3){
foo <- boot(data=dat, statistic=indicator, R=m)
1/mean(foo$t)
}
reps <- replicate(100, p_hat(dat))
但是怎么样?在最后一块代码中,第一个块中的ndx
是什么?