我已经完成了捕获 - 重新捕获历史的随机矩阵,数字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