несколько симуляций в г

Question

Я пытаюсь запустить 1000 симуляций из экспоненциального распределения в R. Я надеюсь получить список симуляций и вектор их средних значений.Но когда я запускаю код, первые 999 элементов списка равны NULL, и, таким образом, первые 999 средств являются NA.что я делаю не так?

nosim <- 1000
n <- 40
r<- 0.2
means<- vector()
sims <- list()
set.seed(4993)
for(i in nosim){
  sims[[i]] = rexp(n, r)
}
means = sapply(sims,mean)

Answer 1

Мне кажется, проблема в том, что ваша переменная nosim является скаляром, а не последовательностью / массивом.

Вместо этого следует использовать seq (или более быстрый seq.int) (как в * 1006).*):

nosim <- 1000
n <- 40
r<- 0.2
means<- vector()   # not really necessary
sims <- list()
set.seed(4993)
for(i in seq(1,nosim)){    # or use `for(i in seq.int(nosim) ){`
  sims[[i]] = rexp(n, r)
}
means = sapply(sims,mean)

Надеюсь, это поможет