使用ggplot绘制NMDS图时遇数据框行数不匹配错误求助

错误原因 #

arguments imply differing number of rows: 12, 26 本质是绘图用的两个数据源行数不匹配：NMDS站点坐标数据（12行，对应12个站点）与生境数据data_2（26行）无法一一对应，导致ggplot无法完成变量映射。

解决步骤（无需修改原始数据结构） #

1. 提取每个站点的唯一生境信息 #

假设data_2包含站点ID列（例：Site）和生境列（例：Habitat），先从data_2中提取每个站点对应的唯一生境（确保单个站点仅归属一个生境）：

library(dplyr)
# 去重，保留每个站点的唯一生境记录
site_habitat <- data_2 %>%
  distinct(Site, Habitat)

2. 提取NMDS站点坐标并关联站点ID #

假设你的NMDS分析对象名为nmds，提取站点得分并添加站点ID（默认与原始物种矩阵的行名对应）：

# 提取NMDS站点得分
nmds_scores <- as.data.frame(scores(nmds, display = "sites"))
# 添加站点ID，行名对应原始站点顺序
nmds_scores$Site <- rownames(nmds_scores)

3. 合并数据框 #

将NMDS坐标和生境数据按站点ID合并，得到行数匹配的绘图数据集：

nmds_combined <- merge(nmds_scores, site_habitat, by = "Site")

4. 绘制带颜色和标签的NMDS图 #

用合并后的数据集绘图，即可同时显示站点点、生境区分色和站点标签：

library(ggplot2)
ggplot(nmds_combined, aes(x = NMDS1, y = NMDS2, color = Habitat)) +
  geom_point(size = 3) +  # 显示站点点
  geom_text(aes(label = Site), hjust = 0.5, vjust = -1) +  # 显示站点标签
  theme_bw()

排查要点 #

• 确认nmds_scores$Site和site_habitat$Site的拼写、大小写完全一致，否则会合并失败
• 如果data_2中存在单个站点对应多个生境的情况，需先核实数据准确性（单个站点仅能归属一个生境才能实现正确着色）

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