0

我要让向量 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)
}
4

1 回答 1

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]]

现在我不认为我正在以您最终想要的方式使用它，但也许这会给您足够的东西来获得您真正想要的东西。

于 2013-02-02T16:29:18.287 回答