请问：如何将某个步骤（该步骤中有FOR语句）重复N次？

levyangel 发表于 2011-10-21 17:34
比如将以下过程重复1000次，计算z的均值
m=10;c=5;a=2;b=4
k=rgamma(1,a,b)

1. m=10;c=5;a=2;b=4
   
2. k=rgamma(1,a,b)
   
3. z<-c()
   
4. for(i in 1:1000)
   
5. {
   
6.         w=runif(m)
   
7.         S=sum(-log(w)/k)
   
8.         y=10-2*k
   
9.         y.star=1-(a+m)/(S+b)
   
10.         z[i]=(y-y.star)^2
    
11. }

复制代码

这是循环的形式，另外可以确定的是，这个是可以向量化运算的

我的实际问题比这个要复杂，我再向量化试一下，谢谢你们的热心帮助

shenbaiseshatan 发表于 2011-10-22 19:45
这是循环的形式，另外可以确定的是，这个是可以向量化运算的

您的程序里的S应该是 S=sum(-2*log(w)/k)

qoiqpwqr 发表于 2011-10-22 21:24
您的程序里的S应该是 S=sum(-2*log(w)/k)

还是你慧眼，哈哈