尝试使用 for 循环进行人口模拟(2)

Question

很抱歉，我又遇到了 for 循环的问题。我正在尝试将循环的人口估计中的最终数字保存到新矩阵中，但我只能让人口估计显示在第 100 行。我知道它与育种Pop2 有关，但我无法弄清楚。任何帮助将非常感激。请在下面找到代码：

finalPop=matrix(nrow=102, ncol=1)
for(i in 1:100){

SWWAyears=data.frame(iteration=rep(NA,101),pop=NA)
breedingPop<-90000
fallMig<-.825
springMig<-.825
winterSurvival<-rbeta(100,.95,.05)
npFecund<-rbinom(100, 3.0, .9)
pFecund<-rbeta(100, .85,.25)
breedingSurvival<-rbeta(100,.95,.05)

# Set initial starting condition
SWWAyears[1,2]=breedingPop
for(years in 2:101) {
  fallPop<-(SWWAyears[years-1,2]*fallMig)
  for (i in 1:100){
  winterPop<-(fallPop*winterSurvival[i])}
  springPop<-(winterPop*springMig)
        for (i in 1:100){
        summerPop<-(springPop*breedingSurvival[i])
        }
            for(i in 1:100){
            breedingPop2<-((summerPop*.26)*npFecund[i])+((summerPop*.14)*pFecund[i])+(summerPop*.60)
            }
  SWWAyears[years,1]=years
  SWWAyears[years,2]<-breedingPop2             
}
finalPop[i,1]<-breedingPop2
}

Answer 1

我认为您的循环结构存在更多基本问题，并且您没有得到预期的正确结果。但是，您关于仅更新第 100 行的具体问题的原因是：

您的变量i 正在您的“内部”for() 循环中更新，因此当您到达finalPop[i, 1] <- breedingPop2 时，i 始终等于100。

您需要在内部 for() 循环中使用不同的变量，例如 j。

finalPop=matrix(nrow=102, ncol=1)

for(i in 1:100){

  SWWAyears = data.frame(iteration=rep(NA,101),pop=NA)
  breedingPop <- 90000
  fallMig <- .825
  springMig <- .825
  winterSurvival <- rbeta(100,.95,.05)
  npFecund <- rbinom(100, 3.0, .9)
  pFecund <- rbeta(100, .85,.25)
  breedingSurvival <- rbeta(100,.95,.05)

  # Set initial starting condition
  SWWAyears[1,2] = breedingPop

  for(years in 2:101) {
    fallPop <- (SWWAyears[years-1,2]*fallMig)
    for (j in 1:100){
      winterPop <- (fallPop*winterSurvival[j])
    }
    springPop <- (winterPop*springMig)
    for (j in 1:100){
      summerPop <- (springPop*breedingSurvival[j])
    }
    for(j in 1:100){
      breedingPop2 <- ((summerPop*.26)*npFecund[j])+((summerPop*.14)*pFecund[j])+(summerPop*.60)
    }
    SWWAyears[years,1] = years
    SWWAyears[years,2] <- breedingPop2
    }
  finalPop[i,1] <- breedingPop2
}

话虽如此，在R中一般不推荐使用多个嵌套的for()循环；您应该能够使用矩阵乘法/矢量化来获得相同的结果。

其他问题

winterPop 和 summerPop 的值将永远分别为 fallPop * winterSurvival[100] 和 springPop * breedingSurvival[100]。这是你想要的吗？