创建alist时，如何在运行时提取函数体

Question

对 R 很陌生，我正在努力寻找这个问题的答案，但也许我不知道如何表达这个问题，如果之前有人问过这个问题，很抱歉。

我想在 for 循环中将函数/模型列表拟合到数据中，而不是手动执行每个。为此，我需要在 alist 中指定模型，列出以下内容：

flist <- alist(y ~ dnorm(mu, sigma),
               mu <- a + b*x,     # need to replace this line with general function
               a ~ dnorm(0, 100),
               b ~ dnorm(0, 10),
               sigma ~ dunif(0, 10)
)

要在 for 循环中执行此操作，我需要能够将 a + b*x 替换为在运行时选择的某个函数的主体。所以我正在考虑制作它：

flist <- alist(y ~ dnorm(mu, sigma),
               mu <- body(function[i]),
               a ~ dnorm(0, 100),
               b ~ dnorm(0, 10),
               sigma ~ dunif(0, 10)
)

然而，alist 似乎并没有真正评估body(function[i])，这似乎是alist 的一部分。有没有办法在运行时取出身体？

基本上，我所追求的是要评估 alist 中的部分表达式（body(...)），但不是全部（保持 ~/

parse 可以做到这一点吗？

Answer 1

我不确定我是否理解，但是你需要这样的东西吗？

funs <- list(f1 = function(x) x + 3, f2 = function(x) 2 * x + 2)


flist <- alist(y ~ dnorm(mu, sigma),
        mu <- FUN_BODY,     # need to replace this line with general function
        a ~ dnorm(0, 100),
        b ~ dnorm(0, 10),
        sigma ~ dunif(0, 10))


flist[[2]] <- eval(substitute(substitute(ARG,
  list(FUN_BODY = body(funs[[2]]))),
  list(ARG = flist[[2]])))

# [[1]]
# y ~ dnorm(mu, sigma)
# 
# [[2]]
# mu <- 2 * x + 2
# 
# [[3]]
# a ~ dnorm(0, 100)
# 
# [[4]]
# b ~ dnorm(0, 10)
# 
# [[5]]
# sigma ~ dunif(0, 10)

Answer 2

1. 使用list 而不是alist。正如您所注意到的，alist 的全部意义不是评估它的论点。
2. 不执行分配。目前还不清楚mu <- 部分应该做什么。您是要使用命名参数，还是要指定因变量？

这可能是你想要的：

flist <- list(
    y ~ dnorm(mu, sigma),
    body(myfunction),
    a ~ dnorm(0, 100),
    b ~ dnorm(0, 10),
    sigma ~ dunif(0, 10)
)

但目前尚不清楚如何使用它。我想你可能想使用这样的东西而不是你的第二行：

mu ~ myfunction(somevar)

也就是说，您需要一个使用myfunction的公式。