R中基本仿真模型的循环结构

Question

我正在尝试编写一个模拟人口增长的基本模型（其初始大小是从正态分布中随机抽取的），然后每个'年'按用户定义的数量增长（目前代码中有2个人）以下为了论证而且）。产生的输出仅显示一次模拟的结果，并且在此模拟中，人口根本没有增长，即每个'年'人口没有增长/没有增加到前几年的人口。我假设我已经在循环结构中填充了一些东西并且热衷于任何建议！

n.years <- 3
n.sim <- 5
store.growth <- matrix(ncol=3,nrow= (n.years * n.sim))

for (i in 1:n.sim) {
 init.pop.size <- rnorm(1,100,10)
 for (j in 1:n.years){
    #grow population
    grow.pop <- init.pop.size + 5
    store.growth[j,] <- cbind(grow.pop, n.years, n.sim)
 }
}
store.growth