R中公式等价的灵活比较？

Question

考虑这些具有交互作用的各种 R 公式：

x ~ a + b + c + d + c:d + a:b
x ~ c + a + b + d + a:b + d:c
x ~ a + b + c + d + c * d + a * b
x ~ a * b + c * d
x ~ b * a + c * d

出于线性模型之类的目的，这些都是等价的。假设我有一大套公式，我想比较是否有任何重复，但可能有像上面这样不明显的重复。有没有一种简单的方法来进行这种比较？

共有三个挑战：

• 必须删除冗余（d + c * d 等价于 c * d）
• 必须能够以不同的顺序匹配元素（a + b 与 b + a 相同）
• 必须能够匹配通勤交互（c:d 与 d:c 相同）

只是一个带有一些排序的 terms() 调用似乎没有得到它，主要是因为最后一个。

这是我到目前为止的工作方式（写成functional sequence 以便于阅读）：

# uses tidyverse

get.terms <- {
  . %>%              
    terms %>%                  # use terms to get the parts
    attr("term.labels") %>%    # character vector of elements
    str_split(":") %>%         # separate interaction terms (makes list)                        
    map_chr(                   # go through each list item
      ~.x %>% 
      sort %>%                 # if multiples (interaction), sort
      paste0(collapse = ":")   # combine back 
    ) %>%                      # output (now standardized) term list
    sort                       # sort the term list for comparison
  }

# Which gives:
get.terms(x ~ a + b + c + d + c:d + a:b)
get.terms(x ~ c + a + b + d + a:b + d:c)
get.terms(x ~ a + b + c + d + c * d + a * b)
get.terms(x ~ a * b + c * d)
get.terms(x ~ b * a + c * d)

# so you can test:
all.equal(get.terms(x ~ b * a + c * d), get.terms(x ~ c + a + b + d + a:b + d:c))

# would have to add more for this, though:
all.equal(get.terms(foo ~ b * a + c * d), get.terms(bar ~ c + a + b + d + a:b + d:c))

但是对于 R 的这样一个基本部分来说，这似乎是 hacky。

我意识到你可能会在接近尾声时通过列表元素比较来缩短这一点，但额外的步骤是有意的，因为这个想法也是为了能够构建一个标准化的人类可读的公式表示法。更何况整个过程，尤其是交互术语的翻转，似乎没有必要

有人知道更简单或更规范的方法吗？

如果它还可以包含潜在的左侧差异，则可以加分。

如果它可以输出标准化的公式格式（或等效的字符串），则加倍奖励。

Answer 1

我想知道这个想法是否可以扩展到更一般的东西。

frm1 <- x ~ a + b + c + d + c:d + a:b
frm2 <- x ~ c + a + b + d + a:b + d:c

identical(
  sort(attr(terms.formula(frm1), "term.labels")),
  sort(attr(terms.formula(frm2), "term.labels"))
)

[1] TRUE