Доверительные интервалы Bootstrap в R?

Как я могу вычислить 95% доверительные интервалы бутстрапа для T.hat с помощью базового доверительного интервала бутстрапа и студенческого интервала бутстрапа?

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)
}

Спасибо!


person user3347124    schedule 08.04.2014    source источник
comment
Я бы посмотрел на пакет boot (cran.r-project.org/ web/packages/boot/boot.pdf) -- в частности, функция boot.ci вычисляет их за вас.   -  person Jealie    schedule 08.04.2014
comment
@Jealie Хорошо, спасибо, но что будет с объектом boot.out для моего бутстрапа? Не знаете, как его вычислить?   -  person user3347124    schedule 09.04.2014
comment
посмотрите, поможет ли вам этот вопрос: stackoverflow.com/questions/13726315/bootstrap -coverage-in-r, так как включает пример boot.out   -  person Jealie    schedule 10.04.2014