核密度的 bootstrap 置信区间

R 自带的数据框 faithful 内保存了美国黄石国家公园 Faithful 火山的 272 次爆发持续时间和间歇时间。为估计爆 发持续时间的密度，可以用核密度估计方法，R 函数 density 可以执行此估计，返回 ? 个格子点上的密度曲线 坐标

x <- faithful$eruptions est0 <- density(x) plot(est0)   set.seed(1) resm <- replicate(10000, { x1 <- sample(x, replace=TRUE) density(x1, from=min(est0$x), to=max(est0$x))$y }) CI <- apply(resm, 1, quantile, c(0.025, 0.97

核密度的 bootstrap 置信区间最先出现在Python成神之路。