函数不填充 R 中的空向量

Question

我要让向量 simualted_results 获取“模拟”返回的值，该值会根据迭代产生一个不同长度的向量。

最初我有这段代码可以运行，但速度很慢：

simulated_results<-NULL
while(as.numeric(Sys.time())-start<duration){
  simulated_results <- cbind(simulated_results,simulation(J,4* (length(J)^2),0.0007,duration,start))

 }

但它很慢所以我修改了它：

start<-as.numeric(Sys.time())
duration<-10
simulated_results<-NULL
simulated_results <- cbind(simulated_results,
                       replicate(n=10000,expr=(while(as.numeric(Sys.time())-start<duration)
                         {simulation(J,4*(length(J)^2),0.0007,duration,start)})))

现在使用新代码，我的问题是，尽管一切都在运行，但我无法将模拟结果传递给 simualted_results，而是 simualted_results jsut 采用 NULL 值的列向量 我没有收到任何错误消息

我将非常感谢任何帮助！

仿真代码供参考：

iter<-as.numeric(Sys.getenv("PBS_ARRAY_INDEX"))

if(iter <= 40){J<-1:500
}else if(iter <= 80){J<-1:1500
}else if(iter <= 120){J<-1:2500
}else if(iter <= 160){J<-1:5000}

set.seed(iter)
simulation <- function(J,gens,v=0.1,duration,start){

  species_richness <- function(J){
    a <- table(J)
    return(NROW(a))
  }

  start<-as.numeric(Sys.time())

  species_richness_output <- rep(NA,gens)
  for(rep in 1:gens){

    if (as.numeric(Sys.time())-start<duration){

      index1 <- sample(1:length(J),1)

      if(runif(1,0,1) < v){
        J[index1] <- (rep+100)
      } 
      else{
        index2 <- sample(1:length(J),1)
        while(index1==index2) {
          index2 <- sample(1:length(J),1)
        }
        J[index1] <- J[index2]
      }
      species_richness_output[rep] <- species_richness(J)} else break
  }

  species_abundance <- function(J){
    a <- table(J)
    return(a)
  }

  abuntable <- species_abundance(J)

  octaves <- function(abuntable)
  {
    oct<-rep(0,floor(log2(length(J))+1))

    for(i in 1:length(abuntable)){
      oct2 <- floor(log2(abuntable[i])+1)
      oct[oct2] <- oct[oct2]+1
    }

    return(oct)
  }   

  octaves(abuntable)
}

Answer 1

我同意@Nathan G 的观点，但确实引起了我的注意：您正在尝试cbind 两个不能绑定在一起的东西，因为它们具有不同的维度。我们不知道你的simulation 函数返回什么样的数据类型，但它显然不是NULL。考虑一下：

df1 <- NULL
df2 <- data.frame(x = 1:10, y = 11:20)
cbind(df1, df2)
cbind(df2, df1)

两个cbind 语句都会出错。你有错误吗？如果发生这种情况，您应该将simulated_results 初始化为NULL，而不是作为simulation 返回的任何函数的空版本。

编辑

iter = 10
set.seed(iter)
J <- 1:1500
# critical to preallocate the list size for speed
res <- vector("list", iter)
for (i in 1: iter) {
    res[[i]] <- simulation(J,4* (length(J)^2),0.0007,duration = 10,start)
}
str(res)
res[[1]]

现在我不认为我使用它的方式完全符合您的最终意图，但也许这足以让您获得您真正想要的东西。