R envelope 曲线的置信区间

说明

该函数根据在多个固定点处评估的曲线的引导复制来计算曲线的整体和逐点置信包络。

用法

envelope(boot.out = NULL, mat = NULL, level = 0.95, index = 1:ncol(mat))

参数

boot.out	"boot" 类的对象，其中 boot.out$t 包含多个固定点处的曲线重复项。
mat	自举矩阵复制了多个固定点处的曲线值。如果未提供boot.out，则这是必需的参数，否则设置为boot.out$t。
level	所需信封的置信度。默认情况下是找到 95% 的置信区间。它可以是标量或长度为 2 的向量。如果它是标量，则可以在该级别找到逐点包络和整体包络。如果 是向量，则第一个元素给出逐点包络的级别，第二个元素给出整体包络的级别。
index	包含引导复制的 mat 的列数。这可用于确保在引导程序中计算的其他统计数据不被视为函数的值。

细节

通过简单地查看每个点的重复的分位数即可找到逐点包络。然后使用 Davison 和 Hinkley (1997) 的方程 (4.17) 计算该包络的总体误差。然后找到一系列逐点包络，直到其中一个的总体误差大约等于所需的水平。如果找不到这样的信封，则返回的信封将仅包含 mat 的每列的极值。

值

包含以下组件的列表：

point	计算了一个矩阵，其中两行对应于与引导重复相同的点处的上下逐点置信包络值。
overall	类似于 point 的矩阵，但包含控制总体误差的包络线。
k.pt	用于逐点包络的分位数。
err.pt	具有两个分量的向量，第一个分量给出逐点包络的逐点错误率，第二个分量给出该包络的总体错误率。
k.ov	用于整体包络的分位数。
err.ov	具有两个分量的向量，第一个分量给出整个包络的逐点错误率，第二个分量给出该包络的整体错误率。
err.nom	长度为 2 的向量，给出逐点包络和整体包络的标称错误率。

例子

# Testing whether the final series of measurements of the gravity data
# may come from a normal distribution.  This is done in Examples 4.7 
# and 4.8 of Davison and Hinkley (1997).
grav1 <- gravity$g[gravity$series == 8]
grav.z <- (grav1 - mean(grav1))/sqrt(var(grav1))
grav.gen <- function(dat, mle) rnorm(length(dat))
grav.qqboot <- boot(grav.z, sort, R = 999, sim = "parametric",
                    ran.gen = grav.gen)
grav.qq <- qqnorm(grav.z, plot.it = FALSE)
grav.qq <- lapply(grav.qq, sort)
plot(grav.qq, ylim = c(-3.5, 3.5), ylab = "Studentized Order Statistics",
     xlab = "Normal Quantiles")
grav.env <- envelope(grav.qqboot, level = 0.9)
lines(grav.qq$x, grav.env$point[1, ], lty = 4)
lines(grav.qq$x, grav.env$point[2, ], lty = 4)
lines(grav.qq$x, grav.env$overall[1, ], lty = 1)
lines(grav.qq$x, grav.env$overall[2, ], lty = 1)

参考

Davison, A.C. and Hinkley, D.V. (1997) Bootstrap Methods and Their Application. Cambridge University Press.

也可以看看

boot , boot.ci