来自R中的For循环的有问题的输出

Question

rho<-rep(0,4)
for(i in 1:4){
  p=rep(0,4)
  cval=rep(0,times =4)
  for(j in 1:4){
    if(mean(k[j,])==pool[i])
      for (l in 1:4){
        rho[l]=.5*(ran[i,2])+((.75)^.5)*rnorm(1,0,1)
        if(rho[l]<0)
          p[l]=k[j,l]*.6
      }
    cval[j]=sum(p)
  }
  CVAL=sum(cval)}

在上面的代码中，我有一个名为k的4x4矩阵。我还有一个名为ran的4x2矩阵。 pool是一个随机排序的向量，其条目与k行的平均值相匹配。

无论生成什么随机数，cval输出始终为(0,0,b,b)，其中b表示随机数。我相信情况不应该如此。有人在我的代码中看到会导致此错误，感谢任何帮助。