我有以下脚本：
randomdiv <- function(ncells, ndivs, size, accuracy) { sz <- matrix(nrow = ncells, ncol = ndivs)
for (j in 1:ncells) {
total_subunits <- size
for (i in 1:ndivs)
{
accurate_subunits <- (size * accuracy)
random_subunits <- round(size - accurate_subunits)
random_inh <- rbinom(1, random_subunits, 0.5)
accurate_inh <- (accurate_subunits / 2)
total_inh <- 2 * (random_inh + accurate_inh)
sz[j,i] <- total_inh
total_subunits <- total_inh
}
}
return (do.call(rbind, replicate(100, sz, simplify = FALSE)))
}
这样我认为randomdiv(5, 20, 10, 0)将返回一个包含500行的矩阵，其中原始sz矩阵已被复制100次。事实上，情况就是这样。但是，重复是相同的，而不是每个重复是新一代数据，这是我需要的。
我有什么想法可以确保每个复制都是一个新的矩阵，而不是第一个要生成的复制？