使用Stargazer包集成超5个回归模型时触发报错求助

问题分析 #

报错Error in if (is.na(s)) { : the condition has length > 1本质是Stargazer在处理模型统计量时，某个变量返回了长度大于1的向量，而非预期的标量值。结合你的情况：随机数据无问题、SCF数据特定触发，且仅模型数>5时报错，大概率是新版本Stargazer在批量处理SCF数据生成的模型时，对某些统计量（如R²、F值、AIC等）的提取逻辑出现了bug。

可行解决方案 #

• 逐个检查模型统计量：
  对每个模型手动提取关键统计量，确认是否存在向量形式的输出：

  # 示例：检查每个模型的R平方长度
  sapply(list(model1, model2, ..., model7), function(x) length(summary(x)$r.squared))
  

  若返回值有大于1的情况，说明该模型的统计量异常，需手动提取标量（如取第一个值summary(x)$r.squared[1]），或重新拟合模型确保统计量为标量。

• 降级Stargazer版本：
  由于问题是更新后出现的，暂时回滚到之前正常的版本：

  # 卸载当前版本
  remove.packages("stargazer")
  # 安装指定旧版本（示例为5.2.3，可根据你之前的版本调整）
  devtools::install_version("stargazer", version = "5.2.3")
  

• 替换为modelsummary包：
  若降级无法解决，用modelsummary作为替代，它对多模型的兼容性更好，语法接近Stargazer：

  library(modelsummary)
  # 把7个模型放入列表并命名
  model_list <- list(
    "模型1" = model1,
    "模型2" = model2,
    ...,
    "模型7" = model7
  )
  # 生成LaTeX表格
  modelsummary(model_list, output = "latex")
  

• 手动修正模型对象：
  若发现某个模型的特定统计量是向量，手动修改模型对象的属性：

  # 示例：强制将模型的R平方设为标量
  model$r.squared <- as.numeric(model$r.squared[1])
  

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