如何使用stargazer包生成按作物分组的可发表级汇总统计表格？

核心实现步骤 #

先按作物分组计算所需统计量，再用stargazer将结果格式化为符合发表要求的表格，以下是完整R代码示例：

示例代码 #

假设你的数据框名为crop_data，包含crop（作物类型）和yield（产量）两个变量：

# 加载依赖包
library(stargazer)
library(dplyr)

# 按作物分组计算汇总统计（以均值、样本量、标准差为例，可按需扩展）
summary_stats <- crop_data %>%
  group_by(crop) %>%
  summarise(
    平均产量 = mean(yield, na.rm = TRUE),
    样本量 = n(),
    标准差 = sd(yield, na.rm = TRUE)
    # 可根据需求添加中位数、最小值、最大值等统计项
  )

# 转换为stargazer兼容的矩阵格式
stats_matrix <- as.matrix(summary_stats[, -1])
rownames(stats_matrix) <- summary_stats$crop

# 生成可发表级表格
stargazer(stats_matrix,
          type = "text", # 需LaTeX格式改"latex"，HTML格式改"html"
          title = "各作物产量汇总统计",
          column.labels = c("平均产量", "样本量", "标准差"),
          rownames = TRUE,
          digits = 2, # 设置小数保留位数
          header = FALSE)

自定义调整说明 #

• 若需更多统计量，直接在summarise()函数内添加对应计算逻辑即可
• type参数可根据最终使用场景灵活切换，适配不同文档格式需求
• 表格标题、列标签、小数位数等均可通过stargazer参数自定义调整

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