R：动态更新公式

如何动态更新公式？

Example:

myvar <- "x"

update(y ~ 1 + x, ~ . -x)
# y ~ 1 (works as intended)

update(y ~ 1 + x, ~ . -myvar)
# y ~ x (doesn't work as intended)

update(y ~ 1 + x, ~ . -eval(myvar))
# y ~ x (doesn't work as intended)

您可以使用paste()内update()call.

myvar <- "x"
update(y ~ 1 + x, paste(" ~ . -", myvar))
# y ~ 1

Edit

正如 @A.Fischer 在评论中指出的，如果myvar是长度 > 1 的向量

myvar <- c("k", "l")
update(y ~ 1 + k + l + m, paste(" ~ . -", myvar))
# y ~ l + m
# Warning message:
# Using formula(x) is deprecated when x is a character vector of length > 1.
#   Consider formula(paste(x, collapse = " ")) instead.

只是“k”被删除，但“l”仍然保留在公式中。

在这种情况下，我们可以将公式转换为字符串，添加/删除我们想要更改的内容并使用以下命令重建公式reformulate， 就像是：

FUN <- function(fo, x, negate=FALSE) {
  foc <- as.character(fo)
  s <- el(strsplit(foc[3], " + ", fixed=T))
  if (negate) {
    reformulate(s[!s %in% x], foc[2], env=.GlobalEnv)
  } else {
    reformulate(c(s, x), foc[2], env=.GlobalEnv)
  }
}

fo <- y ~ 1 + k + l + m

FUN(fo, c("n", "o"))  ## add variables
# y ~ 1 + k + l + m + n + o
FUN(fo, c("k", "l"), negate=TRUE))  ## remove variables
# y ~ 1 + m