R中vegan包RDA/PCA绘图站点标签改为纯数字的问题

核心原因 #

vegan的text.cca（调用text()处理排序对象时的底层函数）默认会给纯数字的站点行名添加row/sit前缀，防止标签被误识别为绘图参数。你只需手动指定标签即可覆盖这个默认行为。

方法1：绘图时直接指定站点标签 #

修改绘图代码中站点文本的部分，通过labels参数传入你的纯数字行名，替换默认生成的带前缀标签：

# 原站点文本代码
# text(tbRDA, display = "sites", col = "dark grey", cex = 1, family = "A")

# 修改后：手动传入纯数字站点ID
text(tbRDA, display = "sites", 
     labels = rownames(spe.hell),  # 用物种数据的行名作为标签
     col = "dark grey", cex = 1, family = "A")

如果环境数据行名和物种数据完全一致，也可以用rownames(env_data)替代rownames(spe.hell)。

方法2：提前修改RDA对象的站点标签 #

如果希望所有后续绘图调用都默认使用纯数字标签，可以直接修改RDA对象中的站点得分行名：

# 修改RDA对象内的站点得分行名
rownames(tbRDA$CA$u) <- rownames(spe.hell)

# 之后正常绘图即可，text()会自动使用纯数字标签
ordiplot(tbRDA, type ="n", family = "A")
points(tbRDA, display = "species", pch = 3, cex = 0.5, col = "dark grey", family = "A")
text(tbRDA, display = "species", col = "dark grey", cex = 1, family = "A")
text(tbRDA, display="bp", col=1, family = "A")
points(tbRDA, display ="sites", family = "A")
text(tbRDA, display = "sites", col = "dark grey", cex = 1, family = "A")

额外验证：确保数据行名正确 #

虽然你已确认行名，可再做一次验证避免读取数据时丢失：

# 检查物种数据行名
print(rownames(spe.hell))
# 检查环境数据行名
print(rownames(env_data))
# 确认两者完全匹配
all(rownames(spe.hell) == rownames(env_data))  # 应返回TRUE

如果读取数据时行名丢失，重新设置：

# 若csv第一列是站点ID，读取时保留为行名
species_data = read.csv("species_data.csv", sep = ";", dec = ",", row.names = 1)
# 或手动设置纯数字行名
rownames(species_data) <- as.character(1:nrow(species_data))
# 重新做Hellinger转换
spe.hell <- decostand(species_data, "hellinger")

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