我一直在尝试计算更大的函数,其中一部分是带有2个条件的while循环。在一定范围的值(x_min和x_max是在整个函数内计算)中,k的Foreach值,我试图计算具有来自k本身是其一部分的分布的值的矩阵。 while循环评估满足分布的必要条件,而foreach函数应为k的每个元素计算while循环。由于我不知道k中元素的确切数量,我认为问题可能是I的预定。我可以得到的最好的是循环内的无穷计算(或R的简单崩溃)。我很感谢任何建议!
I<-mat.or.vec(n,10000)
k<-x[x_min < x & x < x_max]
foreach(k) %do% {
for(i in 1:100){
check=0
while(check==0){
I<-replicate(n=100,rbinom(n= 250, size=1, prob = k/250))
if(sum(I[,i])==k) check=1
}
}
}
不幸地改变顺序不起作用。 它似乎仍有问题。我试图提取矩阵,但它报告“NULL”。
n是250,x_min和x_max由公式(约2-8)定义,k在上面给出的公式中定义,x是0到大约10之间的值(也在公式内计算)。我会为你提供整个公式,但它很大,我可以缩小这些部分的问题,所以我想让问题尽可能简单。感谢您的帮助和评论!