如何解决R语言绘图时x与y变量长度不匹配的错误？

我正在使用以下R代码绘制按物种分类的散点图：

(dput([winglength v weight melitaeini sex use.csv][1]))
plot(chlosyne.janis$total.weight..mg.[chlosyne.janis$sex == "female"], 
     chlosyne.janis$forewing.length..mm.[chlosyne.janis$sex == "female"], 
     col = "red", 
     pch = 16, 
     xlim = c(0, 30), 
     ylim = c(0, 30), 
     xlab = "weight (mg)", 
     ylab = "forewing length (mm)", 
     main = "Reltaionship of weight to forewing length in the genus Chlosyne")
points(chlosyne.janis$total.weight..mg.[chlosyne.janis$sex == "male"],
       chlosyne.janis$forewing.length..mm.[chlosyne.janis$sex == "male"],
       col = "red", 
       pch = 17) 
points(chlosyne.erodyle$total.weight..mg.[chlosyne.erodyle$sex == "female"],
       chlosyne.erodyle$forewing.length..mm.[chlosyne.erodyle$sex == "female"],
       col = "red", 
       pch = 16) 
points(chlosyne.hippochrome$total.weight..mg.[chlosyne.hippochrome$sex == "male"],
       chlosyne.hippochrome$forewing.length..mm.[chlosyne.hippochrome$sex == "male"],
       col = "red", 
       pch = 17) 
points(c.acastus.neumoegenio$total.weight..mg.[c.acastus.neumoegenio$sex == "male"],
       c.acastus.neumoegenio$forewing.length..mm.[c.acastus.neumoegenio$sex == "male"],
       col = "red", 
       pch = 17)

运行后出现错误：

Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differ

x和y变量均为数值型，我已将xlim和ylim都设置为c(0, 30)，以为能让轴长度一致，但即使只运行plot部分也无法生成图像。尝试修改甚至取消xlim和ylim设置，错误依然存在，请问该如何修改代码以正常生成散点图？

核心原因 #

报错的本质是筛选后x轴（体重）和y轴（翅长）的有效观测数量不匹配。这通常是因为其中一列存在NA值，当你用sex == "female"筛选行时，体重列和翅长列的NA值分布不一致，导致两者筛选后的非NA行数不同，进而触发长度不匹配错误。

分步修复 #

1. 统一筛选逻辑，同步移除NA值
   不要单独对每一列做筛选，先对整个数据框进行行筛选，确保每一行的体重和翅长都有有效值：

   # 处理chlosyne.janis雌性数据：筛选性别为female且体重、翅长均非NA的行
   janis_female <- chlosyne.janis[chlosyne.janis$sex == "female" & 
                                   !is.na(chlosyne.janis$total.weight..mg.) & 
                                   !is.na(chlosyne.janis$forewing.length..mm.), ]
   # 用处理后的数据初始化绘图
   plot(janis_female$total.weight..mg., 
        janis_female$forewing.length..mm., 
        col = "red", 
        pch = 16, 
        xlim = c(0, 30), 
        ylim = c(0, 30), 
        xlab = "weight (mg)", 
        ylab = "forewing length (mm)", 
        main = "Relationship of weight to forewing length in the genus Chlosyne")
   

2. 对所有points调用重复上述处理
   以chlosyne.janis雄性数据为例：

   janis_male <- chlosyne.janis[chlosyne.janis$sex == "male" & 
                                 !is.na(chlosyne.janis$total.weight..mg.) & 
                                 !is.na(chlosyne.janis$forewing.length..mm.), ]
   points(janis_male$total.weight..mg., janis_male$forewing.length..mm., col = "red", pch = 17)
   

   其他物种的points调用都按照这个逻辑修改，确保每组x和y的长度完全一致。

3. 更高效的优化方案：合并数据+ggplot2绘图
   原生plot代码重复度高，容易出错，建议合并所有物种数据后用ggplot2批量处理：

   # 为每个数据框添加物种标识列
   chlosyne.janis$species <- "chlosyne.janis"
   chlosyne.erodyle$species <- "chlosyne.erodyle"
   chlosyne.hippochrome$species <- "chlosyne.hippochrome"
   c.acastus.neumoegenio$species <- "c.acastus.neumoegenio"
   
   # 合并所有数据框
   all_data <- rbind(chlosyne.janis, chlosyne.erodyle, chlosyne.hippochrome, c.acastus.neumoegenio)
   # 移除包含NA值的行，只保留需要的列
   all_data_clean <- na.omit(all_data[, c("species", "sex", "total.weight..mg.", "forewing.length..mm.")])
   
   # 加载ggplot2绘图
   library(ggplot2)
   ggplot(all_data_clean, aes(x = total.weight..mg., y = forewing.length..mm., 
                              shape = sex, color = species)) +
     geom_point(size = 3) +
     xlim(0, 30) +
     ylim(0, 30) +
     labs(x = "weight (mg)", y = "forewing length (mm)", 
          title = "Relationship of weight to forewing length in the genus Chlosyne") +
     theme_bw()
   

   这种方式自动处理NA值，还能通过颜色区分物种、形状区分性别，可视化效果更清晰，也减少了重复代码。

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