2

我编写了一个模拟函数，它输出一个包含 100 个元素的向量，我想使用这些*apply函数多次运行该函数，并在每次运行模拟时将重复的输出存储在一个新向量中。

该函数如下所示：

J <- c(1:100)
species_richness <- function(J){
   a <- table(J)
   return(NROW(a))
}

simulation <- function(J,gens,ploton=FALSE,v=0.1){
    species_richness_output <- rep(NA,gens)
    for(rep in 1:gens){
        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)
    }

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

    abuntable <- species_abundance(J)
    print(abuntable)

    octaves <- function(abuntable){
       oct <- (rep(0,log2(sum(abuntable))))

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

       print(oct)
    }   

    # octaves(c(100,64,63,5,4,3,2,2,1,1,1,1))


    if(ploton==TRUE){
       hist(octaves(abuntable))
    }
    print(species_richness(J))

    return(J)
}

simulation(J, 10000,TRUE,v=0.1)

这就是我的函数，它接受 J 一个我之前定义的向量，对其进行操作，然后返回： 新模拟的 100 个元素的向量 J 一个名为 octave 的函数，用于对新向量进行分类 对应于上述“octave”的直方图

我尝试了许多变体： using lapply, mapply putargs=args_from_original_simulation

simulation_repeated <- c(mapply(list, FUN=simulation(args),times=10000))

但是我不断收到 mapply 函数的 match.fun 部分的错误

Error in match.fun(FUN) : 
  'simulation(J, 10000, FALSE, 0.1)' is not a function, character or symbol

尽管我编写的模拟显示为在工作区中保存为函数。

有谁知道这个错误指向什么？

4

1 回答 1

1

在这一行：

simulation_repeated <- c(mapply(list, FUN=simulation(args),times=10000))

您没有提供 mapply 的功能。您（本质上）传递调用的结果simulation(args)并且simulation不返回函数。

于 2013-01-24T23:50:23.930 回答