Ошибка в функции boot.ci () в R

Question

Я пытаюсь вычислить доверительные интервалы начальной загрузки.Вот мой код.

library(boot)
nboot <- 10000 # Number of simulations
alpha <- .01 # alpha level
n <- 1000 # sample size
bootThetaQuantile <- function(x,i) {
  quantile(x[i], probs=.5)
}

raw <- rnorm(n,0, 1) # raw data
( theta.boot.median <- boot(raw, bootThetaQuantile, R=nboot) )
boot.ci(theta.boot.median, conf=(1-alpha)) #this causes no error
boot.ci(theta.boot.median, conf=(1-alpha), type = "percent") #this causes an error

Сообщение об ошибке гласит «Ошибка в ci.out [[4L]]: индекс за пределами».Я очень смущен этим, потому что я не уверен, почему вызов boot.ci вызовет ошибку, если предыдущая строка не вызвала ошибку.

Answer 1

Это потому, что вы должны использовать type = 'perc'.

boot.ci(theta.boot.median, conf=(1-alpha), type = "perc")