r语言 - 从引导程序获取跨样本的平均值

Question

我想在 20 个采样数据中获取均值和 sds，但不确定如何做到这一点。我当前的代码可以为每个样本提供平均值，而不是跨样本。

## create data
data <- round(rnorm(100, 5, 3))
data[1:10]
## obtain 20 boostrap samples
## display the first of the boostrap samples
resamples <- lapply(1:20, function(i) sample(data, replace = T))
resamples[1]

## calculate the means for each bootstrap sample
r.mean <- sapply(resamples, mean)
r.median
## calculate the sd of the distribution of medians 
sqrt(var(r.median))

从上面的代码中，我从每个采样数据中得到了 20 个均值，以及均值分布的 sd。如何从 20 个样本的分布中获得 100 个均值，每个均值？标准差也一样？

非常感谢！！

Answer 1

尽管@konvas的答案可能是您想要的，但在引导方面，我仍然会看一下基本包boot。

看看下面的示例是否可以让你更接近你正在尝试做的事情。

set.seed(6929)    # Make the results reproducible
data <- round(rnorm(100, 5, 3))
boot_mean <- function(data, indices) mean(data[indices])
boot_sd <- function(data, indices) sd(data[indices])
Runs <- 100
r.mean <- boot::boot(data, boot_mean, Runs)
r.sd <- boot::boot(data, boot_sd, Runs)
r.mean$t
r.sd$t
sqrt(var(r.mean$t))
#          [,1]
#[1,] 0.3152989
sd(r.mean$t)
#[1] 0.3152989

现在，查看自举均值和标准误差的分布。

op <- par(mfrow = c(1, 2))
hist(r.mean$t)
hist(r.sd$t)
par(op)

Answer 2

用你的样本制作一个矩阵

mat <- do.call(rbind, resamples)

然后

rowMeans(mat)

将为您提供"样本内"平均值和

colMeans(mat) 

"跨样本"的意思。对于其他量，例如标准偏差，您可以使用apply，例如matrixStats包中的apply(mat, 1, sd)或功能，例如matrixStats::rowSds(mat).