r -bootstrap 用于截尾均值

Question

我有一个来自未知分布的向量，我们希望通过修剪均值来估计其均值。

vec=c(5, 4, 9, 6, 21, 17, 11, 20, 7, 10, 21, 15, 13, 16, 8)

如何估计修剪均值的概率密度函数（μ的估计量）使用 Bootstrap?

Answer 1

这是 bootstrap 向量 vec 的修剪平均值的方法。
定义一个函数来计算 bootstrapped 值，并使用参数 trim 使其更灵活。然后使用基础包 boot 中的函数 boot 来调用该函数。

bootstrap 个值的平均值就是您想要的估算值。

stat <- function(x, i, trim){
  mean(x[i], trim = trim)
}

R <- 1e3
set.seed(2022)
bmean <- boot::boot(vec, stat, R = R, trim = 0.1)

mean(bmean$t)
#[1] 12.12523

d <- density(bmean$t)
hist(bmean$t, freq = FALSE)
lines(d)

r -bootstrap 用于截尾均值

r -bootstrap for trimmed mean

statistics

r