如何传递公式至线性模型构造器并解决boxcox调用报错问题？

以下代码可以正常打印线性模型摘要并生成Box-Cox图：

library(MASS)
library(monomvn)
formula_string <- "y ~ x1 + x2 + x3 + x4"
formula <- as.formula(formula_string)
#linear_model <- lm(formula, data = cement)
linear_model <- lm(y ~ x1 + x2 + x3 + x4, data = cement)
summary_of_linear_model <- summary(linear_model)
print(summary_of_linear_model)
results_of_Box_Cox_Method <- boxcox(linear_model)

但切换注释，改用linear_model <- lm(formula, data = cement)创建模型时，虽然能正常打印线性模型摘要，但执行boxcox(linear_model)会报错：

Error in x$formula : object of type 'closure' is not subsettable

报错原因 #

问题出在你加载的monomvn包——这个包里有一个和R基础包重名的formula函数。当你自定义了名为formula的对象后，R会优先调用monomvn包的formula函数，而非你创建的公式对象。boxcox函数尝试读取模型中的公式时，误将其识别为闭包类型的函数，因此无法执行子集操作，导致报错。

解决方法 #

方法1：修改自定义对象名，避免冲突 #

把你创建的公式对象改名，彻底避开命名冲突：

library(MASS)
library(monomvn)
formula_string <- "y ~ x1 + x2 + x3 + x4"
my_formula <- as.formula(formula_string)
linear_model <- lm(my_formula, data = cement)
summary_of_linear_model <- summary(linear_model)
print(summary_of_linear_model)
results_of_Box_Cox_Method <- boxcox(linear_model)

方法2：明确调用基础包的as.formula函数 #

创建公式时，指定使用stats基础包的as.formula函数，防止被monomvn的同名函数覆盖：

library(MASS)
library(monomvn)
formula_string <- "y ~ x1 + x2 + x3 + x4"
formula <- stats::as.formula(formula_string)
linear_model <- lm(formula, data = cement)
summary_of_linear_model <- summary(linear_model)
print(summary_of_linear_model)
results_of_Box_Cox_Method <- boxcox(linear_model)

方法3：卸载或不加载monomvn包 #

如果你的代码不需要用到monomvn包的功能，直接不加载它即可：

library(MASS)
formula_string <- "y ~ x1 + x2 + x3 + x4"
formula <- as.formula(formula_string)
linear_model <- lm(formula, data = cement)
summary_of_linear_model <- summary(linear_model)
print(summary_of_linear_model)
results_of_Box_Cox_Method <- boxcox(linear_model)

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