ブートストラップ信頼区間を使用していくつかの仮説検定を行おうとしていますが、ブートストラップの設定に問題があります。目標は、シータ = 1 およびレベル 0.05 でヌルをテストすることです。以下に使用されているデータのサンプルを提供しました。
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)
ブートストラップ内で、次の数値最適化も試みています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] のエラー <- function(x) dpois(y, lambda = exp(x)) : サブ代入型の非互換型 (クロージャーから論理型まで) 修正