自定义RDA图为何会使箭头变短且重叠？

问题根源 #

默认plot.rda()绘制解释变量箭头时，会自动对bp（解释变量）得分做尺度适配，使其匹配当前绘图的坐标轴范围；而直接通过scores()提取的bp得分是原始回归系数，未经过该适配处理，导致箭头显示过短、重叠。

解决方法 #

方法1：利用vegan内置的箭头绘制逻辑（推荐） #

直接调用vegan包的arrows()函数，它会自动处理缩放，无需手动调整：

# 替换自定义箭头绘制代码段
vegan::arrows(my_rda2, display = "bp", col = "red", lwd = 3)
# 标签坐标适配缩放后的箭头位置
text(scores(my_rda2, display = "bp", scaling = 1) - c(0.05, 0.03), 
     labels = rownames(sc_bp), 
     col = "red", cex = 1, font = 1)

方法2：手动缩放bp得分（适合精细控制场景） #

如果坚持使用基础R的arrows()，需要将bp得分按站点得分的尺度进行缩放：

# 统一提取各元素得分（保持scaling参数一致）
sc_si <- scores(my_rda2, display="sites", choices=c(1,2), scaling=1)
sc_bp <- scores(my_rda2, display="bp", choices=c(1,2), scaling=1)

# 计算缩放因子：让bp得分范围匹配站点得分范围
si_extent <- apply(sc_si, 2, function(x) max(x) - min(x))
bp_extent <- apply(sc_bp, 2, function(x) max(x) - min(x))
scale_fac <- si_extent / bp_extent

# 缩放bp得分
sc_bp_scaled <- sc_bp * scale_fac

# 绘制缩放后的箭头和标签
arrows(0, 0, sc_bp_scaled[,1], sc_bp_scaled[,2], col = "red", lwd = 3)
text(x = sc_bp_scaled[,1] - 0.05, 
     y = sc_bp_scaled[,2] - 0.03, 
     labels = rownames(sc_bp), 
     col = "red", cex = 1, font = 1)

关键注意事项 #

• 避免手动固定xlim/ylim为c(-1,1)这类窄范围，默认plot.rda()会根据所有得分的极值自动调整坐标轴，手动固定会压缩箭头显示。若需自定义范围，需先计算站点、物种、解释变量得分的全局极值，再设置合理边界。
• 确保所有元素的scaling参数完全一致（比如统一用scaling=1），避免尺度不匹配导致的显示异常。

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