R，我如何使用函數“ for”從“復制”中獲得結果？

Question

我希望使用replicate功能獲得與此代碼相同的結果，

set.seed(1)
n <- 100
sides <- 6
p <- 1/sides
zs <- replicate(10000,{
x <- sample(1:sides,n,replace=TRUE)
(mean(x==6) - p) / sqrt(p*(1-p)/n)
}) 

通過使用for函數。 我嘗試過的代碼就是這樣。

n<-10000
side<-6
p<-1/side
set.seed(1)
for(i in 1:n){
z<-vector("numeric", n)
x<-sample(1:side, 100, replace=TRUE)
z[i]<-mean(x==6)
}

但是上面的代碼不能很好地工作。

Answer 1

Try this:

set.seed(1)
n <- 100
sides <- 6
p <- 1/sides
zs <- replicate(10000,{
  x <- sample(1:sides,n,replace=TRUE)
  (mean(x==6) - p) / sqrt(p*(1-p)/n)
})

n<-10000
sides<-6
p<-1/sides
set.seed(1)
z<-vector("numeric", n)
for(i in 1:n){
  x<-sample(1:sides, 100, replace=TRUE)
  z[i]<-(mean(x==6) - p) / sqrt(p*(1-p)/100)
}

all(zs == z)
# TRUE