我正在尝试使用bootstrap置信区间进行一些假设测试,但是在设置引导程序时遇到了一些问题。目标是在theta = 1和0.05级测试null。我提供了以下使用的数据样本:
y <- c(3,5,3,1,2,1,2,1,0,2,4,3,1,4,1,2,2,0,4,2,2,4,0,2,1)
在bootstrap中,我还尝试使用以下代码进行数值优化:dpois(data,lambda=exp(x)
。我已经提供了我在下面生成的代码:
nboot<-1000
bootxbar <- rep(NA,nboot)
for (i in 1:nboot) {
data <- y[sample(1:n,replace=TRUE)]
bootxbar[i]<-function(x) dpois(data,lambda=exp(x))
}
quantile(bootxbar,c(0.025,0.975))
错误如下:
bootxbar [i]&lt; - function(x)dpois(y,lambda = exp(x))出错: 子分配类型修复中的不兼容类型(从闭包到逻辑)