使用bootstrap置信区间进行假设检验

时间:2015-12-16 05:47:19

标签: r bootstrapping hypothesis-test

我正在尝试使用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))出错:     子分配类型修复中的不兼容类型(从闭包到逻辑)

0 个答案:

没有答案