可以使用rrapply::rrapply()递归搜索元素并返回索引：

library(rrapply)

find_ind <- \(f, value) {
  rrapply(
    f,
    condition = \(x) x == value,
    f = \(x, .xpos) .xpos,
    how = "flatten"
  )
}

find_ind(f, "x")

[[1]]
[1] 3 2 2 2 3

[[2]]
[1] 3 2 3 2 3

确认这些是公式中x的索引：

lapply(find_ind(f, "x"), \(i) f[[i]])

[[1]]
x

[[2]]
x

下面是另一个在rrapply::rrapply()中使用how = "melt"选项的解决方案。这将返回语法树中每个级别和每个叶子的索引，之后我们可以简单地过滤对应于x值的行：

library(rrapply)

f <- y ~ (a*x) +(b*x) + (c*z)

rrapply(f, how = "melt") |>
  subset(value == "x")

#>    L1 L2 L3 L4 L5 value
#> 8   3  2  2  2  3     x
#> 12  3  2  3  2  3     x

或者更短，使用rrapply::rrapply()中的condition参数：

rrapply(f, condition = \(x) x == "x", how = "melt")
#>   L1 L2 L3 L4 L5 value
#> 1  3  2  2  2  3     x
#> 2  3  2  3  2  3     x