2

简而言之，我有一个表达式列表，我想应用于数据帧的每一行。这与这个问题非常相似，但有一个细微的区别，因为我没有函数列表，但有表达式列表。这是我正在尝试的：

我有一个方程式项的数据框，我传递给它paste以形成表达式的文本字符串：

# build a test dataframe
a <- c(10, 11, 12)
b <- c(0.5, 06, 0.7)
c <- c(300, 400, 500)
d <- c(1000, 2000, NA)
p <- c(1.32, 1.37, 1.43)
m <- c(0.45, .55, .65)
params <- data.frame(a, b, c, d, p, m)

# create vector of expressions
ff.txt <- paste(params$a, "*2 + ",
    params$b, "*3 + ",
    params$c, "*4 + ",
    params$d, "*5 * ",
    "x * ", 
    params$p, " * ",
    params$m, 
    sep = "")

这将返回存储为文本的表达式向量，例如“10*2 + 0.5*3 + 300*4 + 1000*5 * x * 1.32 * 0.45”。

然后我编写了一个函数来将文本表达式转换为 R 表达式，使用parse：

funs.list <- lapply(funs.txt, function(x) parse(file = "", text = x))

...funs.txt存储为文本的表达式向量在哪里。

现在，这就是我遇到困难的地方。Imaginefuns.list包含两个表达式。如果我这样做：

lapply(funs.list, eval, envir = list(x=c(7, 60)))

...我得到一个 2x2 矩阵，我只希望答案在对角线上。（最终我将拥有大约 1,000 个表达式，作为优化的一部分，它们将被重复调用，因此不需要额外的计算工作。）

我尝试过的另一种方法是：

mapply(do.call, funs.list, lapply(params, list))

...但我从do.call声明中收到一条错误消息'what' must be a character string or a function。

关于从这里去哪里有什么建议吗？

4

2 回答 2

2

如果我理解正确，您想用 x 的第一个值评估第一个表达式，用第二个值评估第二个表达式，等等。

你可以这样做：

mapply(function(ex, x) eval(ex, envir = list(x = x)), funs.list[1:2], c(7, 60))
于 2015-06-23T10:44:33.537 回答
2

我觉得定义一个与您的代数表达式相对应的自定义 R 函数并apply()在行模式下使用该函数将是一个更好的设计决策（您似乎在原始问题中忽略了这一点）。这个函数也x作为输入，所以我也是原始数据帧cbind的向量。x

exp <- function(a, b, c, d, p, m, x) {
           result <- 2*a + 3*b + 4*c + 5*d*x*p*m
           return(result)
       }

x <- c(1, 2, 3)
params <- cbind(params, x)

output <- apply(params, 1, function(x) {
                               exp(x['a']], x['b'], x['c'], x[["d"]],
                                   x['p'], x['m'], x['x'])
                           })
于 2015-06-23T10:49:22.023 回答