使用StepReg::stepwise()无常数项建模报错：'BREAK'对象未找到

问题场景 #

使用StepReg::stepwise()进行无常数项的逐步回归分析（构建含一阶差分变量的动态模型），已通过dplyr::mutate()生成数据集df，按照文档要求在公式中加入0或-1移除常数项，代码如下：

formula <- dy   ~ 0 + dx1 + dx2 + dx3 + dx4 + ... + xK
res <- stepwise(formula = formula,
                data = df,
                sle = 0.15,
                sls = 0.1,
                type = "linear",
                #include = c("Prova1", "Prova2", "Prova3", "Prova4", "Prova5", "Prova6"),
                strategy = c('bidirection'),
                metric = c("AICc", "adjRsq"))

运行后触发错误：

Error in checkEnterOrRemove(add_or_remove, best_candidate_model, type, :
'BREAK' object not found.

已尝试0和-1两种移除常数项的写法，问题仍存在。

可能的解决方案 #

• 更新StepReg包到最新版本：这个错误大概率是包的旧版本bug，建议运行install.packages("StepReg")更新到CRAN上的最新版，很多同类问题在新版本中已修复。
• 尝试传入预拟合的无常数项模型：先手动拟合一个不含常数项的线性模型，再将该模型对象传入stepwise()，而非直接传入公式：

  base_model <- lm(dy ~ 0 + dx1 + dx2 + dx3 + dx4 + ... + xK, data = df)
  res <- stepwise(model = base_model,
                  sle = 0.15,
                  sls = 0.1,
                  type = "linear",
                  strategy = c('bidirection'),
                  metric = c("AICc", "adjRsq"))
  

• 替代方案：使用其他逐步回归函数：如果更新包后问题依旧，可以考虑用MASS::stepAIC()替代，它对无常数项模型的支持更成熟，示例代码：

  library(MASS)
  base_model <- lm(dy ~ 0 + dx1 + dx2 + dx3 + dx4 + ... + xK, data = df)
  res <- stepAIC(base_model, direction = "both", trace = TRUE)
  

• 检查数据集的完整性：确保df中没有缺失值、异常值，或者变量类型不匹配的情况，这些也可能触发函数内部的逻辑错误。

关于是否是包的bug #

这个错误属于StepReg包内部的变量作用域问题，在旧版本中确实存在相关bug，并非用户代码写法的问题。更新到最新版后如果问题消失，就可以确认是旧版本的bug导致。

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