所以,我有一定的正态分布,它的平均值在一个向量中被复制了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个数字加起来的最佳方法是什么?
我不确定我是否完全得到了这个问题 - 您是否可以将其循环并将其存储在向量中,然后对值进行求和?
#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)