R、 "replicate"の結果を得るために "for"関数を使って同じ結果を得る方法はありますか？

Question

私はfor関数を使用して

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

、replicate機能を使用して、このコードの同じ結果を持っていると思います。私が試みたコードはこれのようなものです。

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