R中add1()设置trace=TRUE无进度输出及LRT值缺失问题咨询

trace 参数无进度输出的原因与解决方式 #

• trace 参数的设计逻辑是：add1 遍历每个候选新增项拟合子模型时，通过控制台独立打印当前计算的项名，这类输出不属于最终返回的结果表格内容。R 4.0及以上版本中，触发的打印效果为逐行输出当前正在拟合的候选项名称，例如trying Agriculture:Examination。
• 无输出的常见诱因：
  • R 版本低于4.0：该版本前的add1.lm方法未实现trace打印逻辑，设置为TRUE无任何效果。
  • 输出缓冲拦截：示例中swiss数据集仅49条观测，10个候选交互项的全量拟合耗时极短，部分IDE控制台会将这类短日志输出与最终表格合并，甚至被缓冲机制吞掉。可在调用add1前执行flush.console()强制即时输出。

设置test="Chisq"后未显示LRT数值的原因 #

• 输出表格中的Pr(>Chi)列就是LRT检验的p值，检验流程已经正常运行。
• 未展示LRT卡方统计量的原因来自add1.lm的固定输出规则：
  • 线性模型场景下，add1默认输出列固定为Df、Sum of Sq（残差平方和减少量）、RSS（基准模型残差平方和）、AIC，设置test参数时仅额外输出对应检验的p值，不会单独输出LRT统计量本身。
  • 源码中LRT卡方值仅用于计算p值，未纳入返回结果列。
• 如需获取LRT统计量，可基于返回结果手动计算，参考代码：

  res <- add1(lm1, scope = .~. + .^2, test="Chisq", trace = TRUE)
  n <- nrow(swiss)
  # 提取基准模型RSS
  rss_base <- res["<none>", "RSS"]
  # 计算每个候选项的LRT值
  res$LRT <- round(n * log(rss_base / res$RSS), 3)
  

add1() 核心运行逻辑说明 #

• 从指定scope范围内筛选出基准模型未包含的候选项，逐一拟合“基准模型+单个候选项”的独立子模型，对比子模型与基准模型的拟合差异。
• 所有计算完成后将指标汇总为表格输出，不自动执行变量筛选操作。

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