有没有比字符串操作更好的替代方案来以编程方式构建公式？

其他人的函数似乎都采用公式对象，然后在内心深处对它们进行黑暗魔法，我很嫉妒。

我正在编写一个适合多个模型的函数。这些模型的部分公式保持不变，部分公式从一个模型到下一个模型有所变化。笨拙的方法是让用户将公式部分作为字符串输入，对它们进行一些字符操作，然后使用as.formula.

但在走这条路之前，我只是想确保我没有忽略一些更简洁的方法，该方法允许函数接受标准 R 格式的公式（例如从其他使用公式的对象中提取的公式）。

我想要像...

> LHS <- y~1; RHS <- ~a+b; c(LHS,RHS);
y ~ a + b
> RHS2 <- ~c;
> c(LHS, RHS, RHS2);
y ~ a + b + c

or...

> LHS + RHS;
y ~ a + b
> LHS + RHS + RHS2;
y ~ a + b + c

...但不幸的是这两种语法都不起作用。有人知道是否有什么可以做到的吗？谢谢。

reformulate会做你想做的事。

reformulate(termlabels = c('x','z'), response = 'y')
## y ~ x + z

或者没有拦截

reformulate(termlabels = c('x','z'), response = 'y', intercept = FALSE)
## y ~ x + z - 1

请注意，您不能构造具有多个的公式reponses例如x+y ~z+b

reformulate(termlabels = c('x','y'), response = c('z','b'))
z ~ x + y

从现有的术语中提取术语formula（给出你的例子）

attr(terms(RHS), 'term.labels')
## [1] "a" "b"

获得响应略有不同，这是一种简单的方法（对于单个变量响应）。

as.character(LHS)[2]
## [1] 'y'


combine_formula <- function(LHS, RHS){
  .terms <- lapply(RHS, terms)
  new_terms <- unique(unlist(lapply(.terms, attr, which = 'term.labels')))
  response <- as.character(LHS)[2]

  reformulate(new_terms, response)


}


combine_formula(LHS, list(RHS, RHS2))

## y ~ a + b + c
## <environment: 0x577fb908>

我认为将响应指定为字符向量会更明智，例如

combine_formula2 <- function(response, RHS, intercept = TRUE){
  .terms <- lapply(RHS, terms)
  new_terms <- unique(unlist(lapply(.terms, attr, which = 'term.labels')))
  response <- as.character(LHS)[2]

  reformulate(new_terms, response, intercept)


}
combine_formula2('y', list(RHS, RHS2))

你还可以定义一个+使用公式的运算符（更新为公式对象设置新方法）

`+.formula` <- function(e1,e2){
  .terms <- lapply(c(e1,e2), terms)
  reformulate(unique(unlist(lapply(.terms, attr, which = 'term.labels'))))
}

RHS + RHS2
## ~a + b + c

您还可以使用update.formula using .明智地

 update(~a+b, y ~ .)
 ##  y~a+b