如何从R语言公式中安全提取竖线右侧的固定效应

你的原方法依赖字符串分割，当公式存在空格变化、复杂表达式（如交互项、变量变换项）时容易失效，推荐使用基于公式语法解析的可靠方法：

方法一：基础R原生实现（无需额外包） #

利用R公式的语法结构，结合all.vars()提取变量名：

formula <- mpg ~ wt | cyl + am
# 定位竖线右侧的表达式部分
right_expr <- formula[[3]][[3]]
# 提取所有变量名
fe <- all.vars(right_expr)
# 结果：[1] "cyl" "am"

这个方法能兼容各种复杂场景，比如右侧是cyl:am + I(cyl^2)时，依然能正确提取cyl和am。

方法二：用formula.tools包简化操作 #

如果经常处理公式，formula.tools包提供了更直观的公式拆分工具：

# 先安装包（首次使用时）
# install.packages("formula.tools")
library(formula.tools)

formula <- mpg ~ wt | cyl + am
# 按竖线拆分公式
split_form <- splitFormula(formula, split = "|")
# 提取右侧部分的变量
fe <- all.vars(split_form$right)
# 结果：[1] "cyl" "am"

内容的提问来源于stack exchange，提问作者pachadotdev