绘制箱线图添加均值时出现xy.coords长度不匹配错误求助

错误原因 #

报错Error in xy.coords(x, y) : 'x' and 'y' lengths differ的核心原因是：
你的siteD$Sampling_Period是包含3个水平的因子，但其中某个水平在数据中没有对应的观测值。此时tapply只会返回存在数据的组的均值（长度为2），而seq(1, nlevels(siteD$Sampling_Period))生成了3个x坐标，两者长度不匹配，引发报错。

排查验证 #

运行以下代码可确认问题：

# 查看因子总水平数
nlevels(siteD$Sampling_Period)
# 计算均值并查看其长度
meanz <- tapply(siteD$Average_Ground, siteD$Sampling_Period, mean, na.rm=T)
length(meanz)
# 查看均值对应的组名
names(meanz)

若length(meanz)小于nlevels(siteD$Sampling_Period)，且names(meanz)缺少某个因子水平，即可验证上述原因。

解决方案 #

提供两种可靠修正方法：

方法1：强制保留所有因子水平的均值结果

通过重新指定因子水平，让tapply为每个因子水平生成结果（无数据组返回NA），确保meanz长度与因子水平数一致：

boxplot(siteD$Average_Ground ~ siteD$Sampling_Period, notch=F)

# 强制保留所有因子水平，无数据组返回NA
meanz <- tapply(siteD$Average_Ground, 
                factor(siteD$Sampling_Period, levels = levels(siteD$Sampling_Period)), 
                mean, na.rm = T)
# 绘制均值，自动跳过NA值
points(seq(1, nlevels(siteD$Sampling_Period)), meanz, 
       col = "red3", pch = 3, lwd = 2, cex = 1.2)

方法2：匹配均值对应的箱线图位置

直接根据均值的组名匹配箱线图的x轴位置，只在有箱线的位置绘制均值：

boxplot(siteD$Average_Ground ~ siteD$Sampling_Period, notch=F)

meanz <- tapply(siteD$Average_Ground, siteD$Sampling_Period, mean, na.rm=T)
# 获取每个均值对应的箱线图x坐标
x_pos <- match(names(meanz), levels(siteD$Sampling_Period))
# 绘制均值
points(x_pos, meanz, col = "red3", pch = 3, lwd = 2, cex = 1.2)

额外建议 #

若箱线图确实显示3组数据（说明3个因子水平均有数据），则可能是代码运行过程中数据集被意外修改。建议在绘制箱线图后立刻计算均值，避免数据篡改。

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