ggplot2绘制偏RDA双图：var_fit箭头长度不匹配问题排查

核心排查方向 #

• 坐标缩放参数不统一
  Base R的plot.cca默认会给排序轴加缩放（比如scaling=1或2），而ggplot2用fortify提取数据时，默认缩放方式可能和Base R不一样。必须确保两者的scaling参数完全一致：

  • Base R里要明确写plot(rda_result, scaling=1)
  • ggplot2中用fortify(rda_result, scaling=1)，或者手动提取分数时指定scaling值

• envfit计算的基准对象不一致
  要保证var_fit（也就是envfit的结果）是基于同一个RDA对象计算的，不能在Base R里用原始rda对象，到ggplot2里却误用转换后的对象。统一用这段代码计算：

  var_fit <- envfit(rda_result$CCA, env_data, permutations = 999)
  

• 箭头绘制的缩放逻辑差异
  Base R的plot(var_fit)会自动匹配当前绘图的轴范围调整箭头长度，但ggplot2手动画箭头时如果没做同样的缩放，就会出现长度偏差。可以手动计算缩放比例：

  # 提取排序轴的最大范围
  axis_max <- max(abs(rda_scores$RDA1), abs(rda_scores$RDA2))
  # 提取箭头的最大长度
  arrow_max <- max(abs(var_fit$vectors$arrows))
  # 计算缩放系数
  arrow_scale <- axis_max / arrow_max
  # 对箭头数据做缩放
  var_arrows <- var_fit$vectors$arrows * arrow_scale
  

• fortify的额外数据转换
  用ggvegan的autoplot时，内部可能对排序分数做了额外转换，和手动fortify提取的数据不一致。建议直接从原始RDA对象提取分数，别依赖fortify默认输出：

  rda_scores <- scores(rda_result, display = "sites", scaling = 1)
  rda_scores <- as.data.frame(rda_scores)
  

验证步骤 #

1. 打印两种方式的轴范围：Base R用par("usr")，ggplot2用ggplot_build(plt)$layout$panel_params[[1]]$x.range，确认轴的上下限是否一致。
2. 输出var_fit$vectors$arrows的数值，对比Base R和ggplot2中有没有做相同的缩放处理。
3. 把所有环节的scaling参数统一成1或2，再重新绘图对比。

修正后示例代码 #

Base R版本 #

# 运行偏RDA
rda_result <- rda(comm ~ env1 + env2 + Condition(conf1 + conf2), data = dat)
# 计算环境拟合
var_fit <- envfit(rda_result$CCA, dat[, c("env1", "env2")], perm = 999)
# 指定scaling绘图
plot(rda_result, scaling = 1, display = "sites")
plot(var_fit, scaling = 1, add = TRUE, col = "red")

ggplot2版本 #

# 统一scaling提取站点分数
site_scores <- as.data.frame(scores(rda_result, display = "sites", scaling = 1))
# 统一scaling提取环境箭头数据
env_arrows <- as.data.frame(scores(var_fit, display = "vectors", scaling = 1))
env_arrows$env <- rownames(env_arrows)

# 绘图时固定坐标轴比例
ggplot(site_scores, aes(x = RDA1, y = RDA2)) +
  geom_point() +
  geom_segment(data = env_arrows, 
               aes(x = 0, xend = RDA1, y = 0, yend = RDA2),
               arrow = arrow(length = unit(0.2, "cm")),
               color = "red") +
  coord_fixed()

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