如何将 1000 个大小为 25 的随机样本的代码更新为一个大小为 25 的随机样本的 1000 个重采样

Question

我目前的代码演示了为 R 中的分布 Gamma(5,1) 抽取 1000 个大小为 25 的随机样本。我的问题是更新此代码以表示从随机抽取 1 个样本中抽取 1000 个重新样本（带替换）相同分布的大小为 25 的样本。

正如我所说，我已经有代码可以为 Gamma(5,1) 执行 1000 个大小为 25 的随机样本

set.seed(35)
n <- 25
alpha <- 5
beta <- 1
sample.n <- 0
xbar.n <- 0
vars.n <- 0
k=1000
for (i in 1:k){
  sample.n <- rgamma(n, alpha, scale=beta)
  xbar.n[i] <- mean(sample.n)
  vars.n[i] <- var(sample.n)
}

Answer 1

这将从大小为 25 的单个 rgamma 分布生成 1000 次绘制。

sample(size = 1000, x = rgamma(n, alpha, beta), replace = TRUE)

更新：

根据您的评论，我认为这是想要的：

xbar.n <- vector("numeric", k)
vars.n <- vector("numeric", k)
rgdist <- rgamma(n, alpha, scale=beta) ##This generates the gamma 
for (i in 1:k){
  sample.n <- sample(size = n, x = rgdist, replace=T)
  xbar.n[i] <- mean(sample.n)
  vars.n[i] <- var(sample.n)
}