Как обновить код для 1000 случайных выборок размером от 25 до 1000 повторных выборок для одной случайной выборки размером 25

Question

Код, который у меня есть в настоящее время, демонстрирует взятие 1000 случайных выборок размером 25 для распределения Gamma (5,1) в R. Моя проблема заключается в обновлении этого кода для представления взятия 1000 повторных выборок (с заменой) из взятия 1 одной случайной выборки. образец размером 25 того же распределения.

Как я уже сказал, у меня уже есть код для выполнения 1000 случайных выборок размером 25 для Gamma (5,1)

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

Это сгенерирует 1000 дро из одного распределения rgamma размера 25.

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)
}