基本的なブートストラップの信頼区間とスチューデント化されたブートストラップの区間によって、T.hat の 95% のブートストラップ信頼区間を計算するにはどうすればよいですか?
failtimes<-c(3, 5, 7, 18, 43, 85, 91, 98, 100, 130, 230, 487)
T.hat=length(failtimes)/sum(failtimes)
#Bootstrap estimate of bias
B <- 999
n <- length(failtimes)
T.b <- numeric(B)
for (b in 1:B) {
y <- sample(failtimes, size = n, replace = TRUE)
T.b[b] <- mean(y)
}
ありがとう!