如何在stargazer生成的相关矩阵表格中插入均值和标准差

stargazer本身未提供直接合并均值、标准差与相关矩阵的内置参数，你可以通过手动拼接目标矩阵再传入stargazer的方式实现需求，全程不需要引入其他工具包，完全在stargazer框架内完成。具体操作步骤如下：

• 第一步：计算所需的所有统计量，包括变量均值、标准差、相关系数矩阵

library(stargazer)
# 以mtcars数据集的4个变量为示例，可替换为你的实际数据
my_data <- mtcars[, c("mpg", "hp", "wt", "qsec")]

# 计算均值、标准差，保留两位小数
mean_val <- round(sapply(my_data, mean, na.rm = TRUE), 2)
sd_val <- round(sapply(my_data, sd, na.rm = TRUE), 2)
mean_sd <- data.frame(均值 = mean_val, 标准差 = sd_val)

# 计算相关系数矩阵，保留两位小数
cor_matrix <- round(cor(my_data, use = "pairwise.complete.obs"), 2)
# 可选：将上三角（含对角线）置空，符合学术论文相关表的常规展示规范
cor_matrix[upper.tri(cor_matrix, diag = TRUE)] <- ""
cor_matrix <- as.data.frame(cor_matrix)

• 第二步：将描述统计结果与相关矩阵拼接为完整表格

# 按列拼接，前两列为均值、标准差，后续列为相关系数
final_matrix <- cbind(mean_sd, cor_matrix)

• 第三步：传入stargazer输出，注意关闭默认的摘要统计功能

stargazer(final_matrix,
          type = "text", # 需导出LaTeX/HTML格式可分别改为"latex"/"html"
          summary = FALSE, # 关闭stargazer默认的摘要统计计算，使用我们手动拼接的内容
          rownames = TRUE, # 保留变量名作为行标签
          title = "变量描述统计与相关系数矩阵",
          notes = "对角线下方为Pearson相关系数",
          digits = 2)

如果需要将均值、标准差放在表格底部而非前两列，调整拼接逻辑为行绑定rbind()即可，输出格式可根据你的发表需求灵活调整stargazer的相关参数。

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