我如何要求R按列编译矩阵？

Question

我尝试通过首先应用二项式生成矩阵sz，然后从pombe_new_subs的相应列添加值，并将此组合值输入为size以用于以下列

经过多次挫折之后，下面的代码就是我最终得到的结果，它只是不起作用 - 我遇到的问题是;

# Error in sz[j, i + 1] = sz[, i] + pombe_new_subs[, i] : 
# number of items to replace is not a multiple of replacement length

pombe_new_subs <- rmultinom(3, 15, prob = c(0.3, 0.3, 0.3))

randomdiv <- function(nchrom, ndivs, size) {
    sz <- matrix(nrow = nchrom, ncol = ndivs) 
    for (j in 1:nchrom) {
        n <- size
        for (i in 1:ndivs) {
            n <- rbinom(1, n, 0.5)  
            sz[j,i] <- n
        }
        sz[j,i+1] = sz[ ,i] + pombe_new_subs[ ,i]
        sz[j, i+1] <- n
    }
    return (sz)
}

randomdiv(3, 3, 10)

我知道这可能是一个相当简单的循环练习，但挫折已经完全接管了。