2

我想知道是否可以使用 lapply() 函数来改变输入的值，类似于：

a1<-runif(100)
a2<-function(i){
a1[i]<-a1[i-1]*a1[i];a1[i]
}
a3<-lapply(2:100,a2)

我正在寻找类似于 for() 循环的东西，但使用 lapply() 基础设施。我无法让 rapply() 做到这一点。

原因是“真正的”a2 函数是一个困难的函数，只有当 a1[i-1] 的值满足某些条件时才需要进行评估。

重新措辞：所以我试图用 lapply() 类型的东西替换下面代码中的 for()：

    a1<-runif(100)
    a2<-function(i, a1){
        a1[i]<-a1[i-1]*2
        a1[i]
    }
    a3<-as.numeric(lapply(2:100, a2, a1=a1))
#compare the output of a3 with that of a1 after the recursive loop
    a2<-a1 #saved for comparison
    for(i in 2:length(a1)){
        a1[i]<-a1[i-1]*2
    }
cbind(a1[2:100],a3)
#actually this is would be like writting a lapply() version of the cumprod() function
cbind(a1,cumprod(a2))

R 邮件列表建议查看 Reduce() 函数....如：

a1<-runif(100)
cadd<-function(x) Reduce("*", x, accumulate = TRUE)
cadd(a1)

这给出了与 cumprod(a1) 相同的结果......但比循环还要慢：

a1<-runif(100000)
cadd<-function(x) Reduce("*", x, accumulate = TRUE)
looop<-function(a1){
j<-length(a1)
    for(i in 2:j){
        a1[i]<-a1[i-1]*a1[i]
    }
a1
}

> system.time(cadd(a1))
   user  system elapsed 
  1.344   0.004   1.353 
> system.time(cumprod(a1))
   user  system elapsed 
  0.004   0.000   0.002 
> system.time(loop(a1))
   user  system elapsed 
  0.772   0.000   0.775 
>

任何的想法 ？

4

1 回答 1

2

编辑：在您澄清之后：不，我不相信您可以使用 apply 函数来递归地做这样的事情。应用函数的全部意义在于它同时应用于向量/矩阵。

您可能还想在 stackoverflow 上查看这个相关问题。

我的旧答案：

试试这个：

a1<-runif(100)
a2<-function(i, a1){
    a1[i]<-a1[i-1]*a1[i]
    a1[i]
}
a3 <- as.numeric(lapply(2:100, a2, a1=a1))

与for循环不同，您需要在lapply. return 也是一个列表，因此您需要将其转换回您想要的任何形式。

您可能还想查看 plyr 包，了解执行此类操作的简单方法。

除此之外，您可以在没有循环的情况下进行操作：

a3 <- a1[-length(a1)] * a1[-1]

换句话说，这些语句是完全等价的：

> all((a1[-length(a1)] * a1[-1]) == as.numeric(lapply(2:100, a2, a1=a1)))
[1] TRUE

但是第一个版本更可取，因为它没有迭代。

于 2009-10-13T13:25:22.427 回答