R语言密度图图例缺失项目名称问题排查与修复求助

问题现象 #

直接用density(mydata$item1)生成的密度图，图例能正常显示对应项目名（比如item1）。但处理带NA值的变量（像item2、item3）时，用density(mydata$item2, na.rm=TRUE)生成的图，图例只会显示x = .，没法正确标注具体是哪个项目。

问题原因 #

核心是清理NA后丢失了变量的名称属性：mydata$item2原本自带item2这个名称，但用na.rm=TRUE生成密度对象时，返回结果里的x向量名称变成了.，导致基础绘图函数识别不了原始变量名，只能显示默认的x = .。

修复方法 #

方法1：手动给密度对象补全名称 #

生成密度对象后，直接给d$x赋值变量名，再绘图：

# 处理含NA的item2
d <- density(mydata$item2, na.rm = TRUE)
# 手动设置x的名称为item2
names(d$x) <- "item2"
# 绘图
plot(d)

方法2：提前清理NA并保留变量名 #

先把清理好NA值的向量存成带名称的变量，再生成密度图：

# 清理item2的NA值，保留变量名
clean_item2 <- mydata$item2[!is.na(mydata$item2)]
# 生成密度对象
d <- density(clean_item2)
# 绘图
plot(d)

方法3：用ggplot2自动处理（更省心） #

如果换用ggplot2包，不仅能自动处理NA和图例，还能更灵活调整图的样式：

library(ggplot2)

ggplot(mydata, aes(x = item2)) +
  geom_density(na.rm = TRUE) +
  labs(x = "item2", title = "item2的密度分布")

测试用数据集 #

mydata <- structure(list(item1 = c(5, 5, 5, 5, 4, 4, 2, 1, 3, 4, 4, 3, 
2, 5, 2, 4, 4, 3, 6, 5, 3, 2, 5, 3, 3, 1, 3, 5, 1, 3, 2, 6, 3, 
5, 4, 4, 3, 5, 6, 3, 2, 6, 6, 5, 2, 2, 2, 3, 3, 3), item2 = c(5, 
4, 5, 1, 2, 2, 3, 2, 2, 2, 2, 3, 2, 5, 1, 4, 4, 3, 3, 5, 3, 2, 
4, 4, 3, 4, 4, 3, 7, NA, 2, 4, 2, 4, 2, 3, 5, 3, 5, 3, 2, 6, 
6, 7, 2, 3, 2, 3, 1, 4), item3 = c(5, 5, 6, 7, 3, 4, 5, 2, 2, 
6, 4, 2, 5, 7, 1, 2, 4, 5, 6, 6, 5, 2, 6, 5, 6, 4, 6, 4, 6, 4, 
6, 5, 5, 6, 6, 6, 5, 6, 7, 5, 5, 7, 7, 6, 2, 6, 6, 6, 5, 3)), row.names = c(NA, 
-50L), class = c("tbl_df", "tbl", "data.frame"))

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