2016-11-23 5 views
0

,R, "replicate"의 결과를 얻기 위해 "for"함수를 사용하여 동일한 결과를 얻을 수 있습니까? <code>for</code> 기능을 사용하여 내가 <code>replicate</code> 기능을 사용하여이 코드 같은 결과를하고 싶은

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

위의 코드는 제대로 작동하지 않습니다.

+0

벡터화 된 (따라서 일반적으로 더 나은) 방법이 완벽하게 작동 할 때 왜 'for'를 사용하려고하는지 궁금해합니다 ...? – LAP

+1

'z <-vector ("numeric", n)'을 루프 밖으로 가져옵니다. 현재의 경우 반복 될 때마다 z가 초기화됩니다 (다시) – etienne

답변

0
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 
+0

너무 늦게 선택합니다. 고마워요 !! –