如何在ggplot2中为频谱分析图添加反向次级x轴

要实现主x轴（µHz频率）下方添加反向周期次级x轴，核心是利用ggplot2的scale_x_continuous()配合sec_axis()参数，结合频率与周期的物理转换逻辑来完成，具体操作如下：

关键转换逻辑 #

主x轴的频率（单位：µHz）与周期的转换公式为：

周期（秒）= 1e6 / 频率（µHz）
因为1µHz = 1e-6 Hz，而周期是频率的倒数，所以需要先将µHz转换为Hz（乘以1e-6），再取倒数，简化后就是1e6 / x。如果你的period列已经是计算好的周期值，直接对应这个公式即可。

完整代码示例 #

library(ggplot2)

# 构造模拟频谱数据（含多面分组）
set.seed(123)
spec_data <- data.frame(
  x = rep(seq(10, 100, 10), 4),  # 主x轴：µHz频率
  y = rnorm(40, 5, 1),           # 频谱强度
  facet_group = rep(c("Group1", "Group2", "Group3", "Group4"), each = 10)
)

# 绘制多面折线散点图+反向次级x轴
ggplot(spec_data, aes(x = x, y = y)) +
  geom_point(size = 2, color = "#2c3e50") +
  geom_line(color = "#3498db") +
  facet_wrap(~facet_group) +
  # 设置主x轴与次级反向周期轴
  scale_x_continuous(
    name = "Frequency (µHz)",
    sec.axis = sec_axis(
      trans = ~1e6 / .,  # 频率转周期的公式
      name = "Period (s)",
      trans = "reverse"  # 开启轴反向
    )
  ) +
  theme_minimal() +
  # 调整轴标签间距，避免重叠
  theme(
    axis.title.x = margin(t = 10),
    axis.title.x.top = margin(t = 10),
    strip.text = element_text(size = 11, face = "bold")
  )

注意事项 #

• 如果周期使用其他单位（如分钟），只需修改转换公式：~(1e6 / .)/60，同时更新轴名称为"Period (min)"。
• 次级轴是视觉转换层，所有绘图数据仍基于主x轴的x列，这能保证双轴的数值对应关系完全准确，避免直接映射period列可能出现的错位问题。
• 若多面图设置scales = "free_x"，次级轴会自动适配每个面板的主轴范围，无需额外配置。

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