Как сгенерировать вектор 10000 раз и каждый раз сохранять определенную позицию? - PullRequest
0 голосов
/ 24 апреля 2019

Итак, у меня есть определенное нормальное распределение, и его среднее значение повторяется 10000 раз в векторе. Мне нужно найти 9-й элемент в векторе, который находится за пределами контроля и сохранить его положение, что я смог сделать. Теперь, как повторить это действие 10000 раз и сохранить 10000 позиций как целые числа? Это то, что я до сих пор:

LCI = 1000 - 3*4^(1/2)

LCS = 1000 + 3*4/4^(1/2)

k = replicate(10000, (mean(rnorm(4, mean=1000, sd=4))))

k[k < LCI | k > LCS]

n = which(k %in% c(k[k < LCI | k > LCS]))

n[9]

Бонусный вопрос: тогда, как лучше всего сложить все эти 10000 номеров, которые я получил?

1 Ответ

0 голосов
/ 24 апреля 2019

Я не уверен, что у меня возник вопрос полностью - не могли бы вы просто повторить это и сохранить его в векторе, а затем суммировать значения?

#empty vector  
l <- c()

#making constant
LCI = 1000 - 3*4^(1/2)

LCS = 1000 + 3*4/4^(1/2)

#for loop with you code
for(i in 1:1000){


  k = replicate(10000, (mean(rnorm(4, mean=1000, sd=4))))

  #k[k < LCI | k > LCS] (don't need for the loop)

  n = which(k %in% c(k[k < LCI | k > LCS]))

  #store in empty vector as integer
  l[i] <- as.integer(n[9])
}

#check output
l

#sum of all values
sum(l)
...