correlation_heatmap

数据准备

两种画法共用同一份相关矩阵和 p 值矩阵，所以先算好。

d <- read.csv("content/tessera/data/csv/mtcars.csv")
vars <- c("mpg", "cyl", "disp", "hp", "drat", "wt", "qsec", "vs", "am", "gear")
df <- d[, vars]

cor_mat <- cor(
  df,
  use = "pairwise.complete.obs",  # 有缺失时逐对使用完整观测，不整行删掉
  method = "pearson"              # 线性相关；变量分布很偏时改 "spearman"
)

星号来自 p 值，得单独算一张同样大小的矩阵。这个函数只依赖 base R：

cor_pmat <- function(x, method = "pearson") {
  mat <- as.matrix(x)
  n <- ncol(mat)
  p <- matrix(NA_real_, n, n)
  colnames(p) <- rownames(p) <- colnames(mat)
  for (i in seq_len(n)) {
    for (j in seq_len(n)) {
      p[i, j] <- cor.test(mat[, i], mat[, j], method = method)$p.value
    }
  }
  p
}

p_mat <- cor_pmat(df)

方法 A · corrplot

#| fig: corrplot
#| fig-width: 8
#| fig-height: 6
library(corrplot)

corrplot(
  corr        = cor_mat,
  method      = "color",      # 用色块填格；也可以是 circle / square / number
  type        = "full",       # 画满；upper / lower 只画一半，省掉重复信息
  order       = "hclust",     # 按层次聚类重排变量 —— 块状结构就是这么来的
  hclust.method = "complete", # 聚类的连接方式，影响分块边界怎么切
  col         = colorRampPalette(c("#1991A9", "#E7E9E4", "#495A2E"))(200),
  col.lim     = c(-1, 1),     # 固定色标范围，否则色深随数据变、两张图没法比
  addgrid.col = "white",      # 格线颜色，白线让色块之间断开
  cl.pos      = "r",          # 色标（colour legend）放右边
  cl.ratio    = 0.15,         # 色标占整图宽度的比例
  tl.pos      = "l",          # 变量名（text label）只放左边，不重复放上面
  tl.col      = "black",
  tl.srt      = 0,            # 标签旋转角度，0 = 水平
  p.mat       = p_mat,        # 传入 p 值矩阵，星号才画得出来
  sig.level   = c(0.001, 0.01, 0.05),
  insig       = "label_sig",  # 不显著的格子留白；改成 "blank" 会整格涂掉
  pch.cex     = 0.9,          # 星号大小
  pch.col     = "black",
  title       = "Correlation Matrix of mtcars",
  mar         = c(0, 0, 1.8, 0)  # 上边留白，否则标题会被裁掉（base graphics 的老问题）
)

方法 B · ggcorrplot

返回的是 ggplot 对象，所以能继续叠图层。星号这里自己叠一层 geom_text 而不是用包内置选项——分级、颜色、大小都能自己控制。

#| fig: ggcorrplot
#| fig-width: 8
#| fig-height: 6
library(ggcorrplot)
library(ggplot2)

gg <- ggcorrplot(
  corr          = cor_mat,
  method        = "square",   # 方块；"circle" 会用面积额外编码强度
  type          = "full",
  colors        = c("#1991A9", "#E7E9E4", "#495A2E"),  # 低 / 中 / 高，三个锚点
  ggtheme       = ggplot2::theme_void,  # 从空主题起步，自己加要的部分
  hc.order      = TRUE,       # 同 corrplot 的 order = "hclust"
  hc.method     = "complete",
  outline.color = "grey85"    # 格子描边；设 NA 可以完全去掉
) +
  labs(title = "Correlation Matrix of mtcars") +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16, face = "bold"),
    # 显式给白底：theme_void 是透明的，导出 PNG 时会变成透明背景
    panel.background = element_rect(fill = "white", colour = NA),
    plot.background  = element_rect(fill = "white", colour = NA)
  )

# 星号单独一层：先把 p 值矩阵摊成长表，再按聚类后的顺序对齐
star_mat <- ifelse(p_mat <= 0.001, "***",
             ifelse(p_mat <= 0.01, "**",
              ifelse(p_mat <= 0.05, "*", "")))
diag(star_mat) <- ""   # 对角线是变量和自己，永远显著，标了没意义

star_df <- as.data.frame(as.table(star_mat), stringsAsFactors = FALSE)
names(star_df) <- c("Var1", "Var2", "stars")
star_df <- subset(star_df, nzchar(stars))

# 关键一步：因子水平必须取 ggcorrplot 聚类后的顺序，否则星号会错位到别的格子
star_df$Var1 <- factor(star_df$Var1, levels = levels(gg$data$Var1))
star_df$Var2 <- factor(star_df$Var2, levels = levels(gg$data$Var2))

p <- gg + geom_text(
  data        = star_df,
  aes(x = Var2, y = Var1, label = stars),
  inherit.aes = FALSE,   # 不继承底图的 fill 映射，否则星号会被上色
  size        = 3.8
)

p

导出

页面上这两张是 144 dpi，够看。要投稿或进幻灯片得自己导一遍，而两种画法导出方式根本不同——这也是 base graphics 和 ggplot2 最实际的一处区别。

base graphics 的图不是对象：corrplot() 一执行，图就画在当前设备上了，没有东西可以先存起来再交给别的设备。所以三种格式只能各开一次设备、把同一段绘图代码重跑一遍。

#| eval: false
library(ragg)

# 包成函数只是为了不把 corrplot(...) 那一长串抄三遍
draw <- function() {
  corrplot(corr = cor_mat, method = "color", order = "hclust", ...)  # 同上，原样
}

agg_png("cor_mtcars.png", width = 8, height = 6, units = "in",
        res = 300, background = "white")
draw()
dev.off()   # 忘了这句，文件是空的或者被下一张图接着往里画

# 期刊投稿通常要 TIFF，300 dpi 起，LZW 是无损压缩
agg_tiff("cor_mtcars.tiff", width = 8, height = 6, units = "in",
         res = 300, compression = "lzw", background = "white")
draw()
dev.off()

ggplot 是对象，导出就是把对象交给设备，不用重画：

#| eval: false
ggsave("ggcor_mtcars.png", p, width = 8, height = 6, dpi = 300,
       device = ragg::agg_png)   # 指定 ragg：默认设备在 Windows 上字体渲染很糙

ggsave("ggcor_mtcars.tiff", p, width = 8, height = 6, dpi = 300,
       device = function(...) ragg::agg_tiff(..., compression = "lzw"))

ggsave() 不需要 background —— 前面在 theme() 里显式给的白底就是它的背景。这两件事是同一个：theme_void 不补白底的话，这里导出来就是透明的。

别用 JPEG 存热图。 热图是大块纯色加硬边，JPEG 的有损压缩正好在这种边缘上出振铃，色块交界会糊出一圈杂色——省下的那点体积不值得。JPEG 适合照片，不适合色块图。

两种方法怎么选

1. 只是想看一眼有没有结构 → A。一个函数出图，参数都在一处，改颜色改图例不用翻文档。
2. 这张图要进一套统一风格的报告 → B。它是 ggplot 对象，字体、主题、留白能和别的图对齐。
3. 要在图上叠自定义标注（不只是星号，比如圈出某几个变量对） → B。A 的图层系统是封闭的，加东西只能靠 base graphics 的 text() 硬画坐标。

两张图表达的是同一组系数，选哪个不影响结论，只影响后面好不好改。