我正在尝试使用for循环构建两人口模型。我尝试了这些步骤:首先,我使用名为rmax的随机数生成函数runif创建了一个向量,该代码用于代码的逻辑增长方程。我手动将rmax中的第一个元素设置为NA,用作占一个人口数量的占位符。我将输入参数设置为K和时间步数(30个时间步)。然后,将第一时间(先前设置为NA)的初始总体大小的占位符设置为15。我还希望循环开始时间2。但是,当我运行循环时,我不断收到错误消息“ “替换不是替换长度的倍数”“ 29倍。这是第二次使用R,我是一个初学者,任何建议都将非常有用
谢谢
rmax<-runif(30,min=-0.2,max=0.5)
rmax[1]<-NA
rmax
k<-500
tsp<-30
N<-rep(NA, tsp)
N[1]<-15
for(i in 2:tsp){N[i]<-N[i-1] + N[i-1] * rmax * (1-N[i-1]/k)}