求可兼容psce包、支持导出pcse类回归表为HTML的R工具包

texreg、stargazer报错是因为没有内置pcse类的结果提取规则，以下是两个可直接实现需求的方案： #

方案1：使用gtsummary包（无需手动整理结果） #

gtsummary支持自定义模型提取规则，适配pcse类后可直接导出HTML格式回归表：

• 先安装加载所需依赖包：

install.packages(c("gtsummary", "pcse", "gt"))
library(gtsummary)
library(pcse)
library(gt)

• 注册pcse类的提取方法后直接导出：

# 注册pcse模型的tidy适配方法
tidy.pcse <- function(x, conf.int = TRUE, conf.level = 0.95, ...) {
  res <- data.frame(
    term = names(x$coefficients),
    estimate = as.numeric(x$coefficients),
    std.error = as.numeric(x$se),
    p.value = as.numeric(x$p)
  )
  if (conf.int) {
    critical <- qnorm((1 + conf.level)/2)
    res$conf.low <- res$estimate - critical * res$std.error
    res$conf.high <- res$estimate + critical * res$std.error
  }
  return(res)
}

# 生成格式化回归表并导出为HTML
tbl_regression(Model6_reg_pcse) %>%
  as_gt() %>%
  gtsave(filename = "Model6_regression.html")

方案2：手动提取结果后用htmlTable包导出 #

如果不想自定义适配规则，可直接提取pcse返回的统计量组装成表格后导出：

install.packages("htmlTable")
library(htmlTable)

# 提取所需统计量组装为数据框
reg_result <- data.frame(
  变量 = names(Model6_reg_pcse$coefficients),
  回归系数 = round(Model6_reg_pcse$coefficients, 3),
  面板校正标准误 = round(Model6_reg_pcse$se, 3),
  P值 = round(Model6_reg_pcse$p, 3)
)

# 导出为HTML文件
htmlTable(reg_result, file = "Model6_regression.html")

注意：你问题中提到的psce为拼写误差，pcse包输出的回归结果类名为pcse。

内容的提问来源于stack exchange，提问作者Vera Lind Hansen