复现顶刊中这张图，换数据即可用（进阶）！

pythonic生物人2026年04月08日约 2 分钟

又看到SCI上一张Mantel test相关性网络Heatmap新玩法

Butterfly diagram illustrating the correlation between immune cell infiltration, immune-related functions, and CTRGPS. s12885-024-12952-w fig7 h

这张图展示了CTRGPS分别和immune cell infiltration、immune-related functions之间的Mantel检验结果，揭示了三组数据之间的潜在关联关系，可辅助进一步研究immune cell infiltration和immune-related functions之间的关联关系。

本次简单介绍一下R中如何实现这张图，原理很简单，例如，有varechem、varespec、varechem_subset三组数据：

✅第1步，对varechem、varespec数据进行mantel_test

# 关键代码
mantel <- mantel_test(varespec, varechem,
  spec_select = list(
    Spec01 = 1:7,
    Spec02 = 8:18,
    Spec03 = 19:37,
    Spec04 = 38:44
  )
)

✅第2步，对varechem、varespec数据进行可视化，相关热图（左下角）和Mantel检验结果（右上角）

# 关键代码
qcorrplot(correlate(varechem), type = "lower", diag = FALSE) +
  geom_square() +
  geom_couple(aes(colour = pd, size = rd),
    data = mantel,
    curvature = nice_curvature()
  )

✅第3步，对varechem_subset、varespec数据进行mantel_test

#关键代码
mantel <- mantel_test(varespec, varechem_subset,
  spec_select = list(
    Spec01 = 1:7,
    Spec02 = 8:18,
    Spec03 = 19:37,
    Spec04 = 38:44
  )
)

✅第4步，对varechem_subset、varespec数据进行可视化，相关热图（右下角）和Mantel检验结果（左下角）

# 关键代码
qcorrplot(correlate(varechem_subset), type = "upper", diag = FALSE) +
  geom_square() +
  geom_couple(aes(colour = pd, size = rd),
    data = mantel,
    curvature = nice_curvature()
  )

✅第5步，两张图通过varespec数据连接（共用Spec01、Spec02、Spec03、Spec04）

# 拼接
p_final <- p_left + p_right + 
  plot_layout(widths = c(1, 1)) &
  guides(
    size = guide_legend(title = "Mantel's r", override.aes = list(colour = "grey35"), order = 2),
    colour = guide_legend(title = "Mantel's p", override.aes = list(size = 3), order = 1),
    fill = guide_colorbar(title = "Pearson's r", order = 3)
  )

p_final