我在R
的给定矩阵上运行泊松。出于某些原因,我想避免在示例中任何rowSums和colSums等于5。我在我的代码中添加了repeat
来解决这个问题但是运行永远不会结束,即使对于小矩阵也是如此。
有什么想法吗? 感谢
mat <- matrix(c(seq(1,60,1)), nrow = 6, ncol = 10, byrow = TRUE)
poiss <- function(mat) {
repeat{
pmat <- rpois(length(mat), 1:3)
dim(pmat) <- dim(mat)
if(any(rowSums(pmat)) == 5 | any(colSums(pmat)) == 5) break
}
return(pmat)
}
p = poiss(mat)