如何获得R中每行几乎有一个数字1的随机矩阵

时间:2016-01-29 13:36:38

标签: r

我已经完成了捕获 - 重新捕获历史的随机矩阵,数字1和0到“for循环”。我需要每次生成矩阵时,每行几乎有一个数字1。但是,如果我尝试“while循环”,总会有一行没有。

M = 10
oc=5
p = runif(oc, 1 , 1)
phi = runif(oc, 0.2, 0.2) 
hc.mat = matrix(0, M, oc)
z.mat = matrix(0, M, oc)

#For the occasion 1
z.mat[,1] = rbinom(M, 1, phi[1])
hc.mat[,1] = rbinom(M, 1, p[1]*z.mat[,1])

#For occasions >1
for (i in 2:oc) { 
  z.mat[,i] = rbinom(M, 1, phi[i])
  hc.mat[,i] = rbinom(M, 1, 1*z.mat[,i]) 
  j <- sum(hc.mat[1:length(M),])
  while (j >= 1) {
    break; print(hc.mat)
  } 
}
hc.mat

0 个答案:

没有答案