3

我想添加一些函数 f1,f2,...,fn 以便我有一个新函数产生 f(x)=f1(x)+...+fn(x) （称为逐点加法）。所以我有一个功能列表并尝试过

Reduce("funadd",fun.list)

其中 funadd 定义为

funadd <- function(f1,f2){
    retfun <- function(x){
        f1(x)+f2(x)
    }
    retfun
}

在两个函数上测试 funadd 时，它可以完美运行。但是，当我尝试评估 Reduce 命令的结果时，出现错误

Error: evaluation nested too deeply: infinite recursion / options(expressions=)?
4

2 回答 2

5

有趣的是Reduce不起作用......请注意，“手动减少”有效：

f <- function(x) x^2
g <- function(x) x^3
h <- function(x) x^4
x <- runif(3)

f(x)+g(x)+h(x)
#[1] 0.9760703 0.1873004 0.1266966

funadd(funadd(f,g),h)(x)
#[1] 0.9760703 0.1873004 0.1266966

或者，您可以使用这个：

funadd2 <- function(...){
    function(x) Reduce(`+`, lapply(list(...), function(f) f(x)))
}

funadd2(f,g,h)(x)
#[1] 0.9760703 0.1873004 0.1266966

编辑：这是怎么回事：

查看 的源代码Reduce，我们可以看到它（大致）有一个循环来执行此操作：

init <- f
init <- funadd(init, g)

如果有更多元素（init <- funadd(init, h)，...），则继续。

这会导致引用f在第一次循环迭代中丢失：

init(x)
# Error: evaluation nested too deeply: infinite recursion / options(expressions=)?

发生这种情况是因为f1最后一个retfun指向自身：

identical(environment(init)$f1, init, ignore.environment=FALSE)
# [1] TRUE

正如@Vincent 发现的那样，这也可以通过强制参数来解决，即通过制作一个避免对f1and进行惰性求值的本地副本f2：

funadd3 <- function(f1,f2){
    f1.save <- f1
    f2.save <- f2
    retfun <- function(x){
        f1.save(x)+f2.save(x)
    }
    retfun
}

Reduce(funadd3, list(f,g,h))(x)
# [1] 0.9760703 0.1873004 0.1266966
于 2013-09-11T12:48:39.587 回答
5

强制对参数进行评估可以解决问题。

funadd <- function(f1,f2){
    force(f1)
    force(f2)
    retfun <- function(x){
        f1(x)+f2(x)
    }
    retfun
}
r <- Reduce( funadd, list( f, g, h ) )
r(x)  # works
于 2013-09-11T13:13:44.430 回答