使用Viridis绘制小提琴图并去除NA值时遇美学映射错误

核心原因 #

报错的本质是：你的ggplot美学映射中，某一个映射参数引用了原数据框的列（长度与去NA后的withoutNAs不一致），而非直接使用withoutNAs内的列名。原数据框行数多于withoutNAs的61行，导致映射向量长度不匹配。

解决步骤 #

1. 检查去NA逻辑
   确保你是对目标列（或全数据框）正确过滤NA，比如：

   # 精准过滤目标列的NA（推荐）
   withoutNAs <- 原数据框 %>% filter(!is.na(你的数值列))
   # 或全数据框去NA
   withoutNAs <- na.omit(原数据框)
   

   确认nrow(withoutNAs)返回61，与报错提示的长度一致。

2. 修正美学映射写法
   绝对不要用原数据框$列名的方式在映射中引用数据，要直接写列名——ggplot会自动从你指定的data = withoutNAs中取值。

   错误示例（会触发报错）：

   ggplot(withoutNAs, aes(x = 原数据框$分组列, y = 原数据框$数值列, fill = 原数据框$分组列)) +
     geom_violin() +
     scale_fill_viridis_d()
   

   正确示例：

   ggplot(withoutNAs, aes(x = 分组列, y = 数值列, fill = 分组列)) +
     geom_violin() +
     scale_fill_viridis_d() # 离散变量用_d，连续变量用_c
   

3. Viridis配色适配
   根据你的fill变量类型选择对应配色函数：

   • 离散分组：scale_fill_viridis_d()
   • 连续数值：scale_fill_viridis_c()

验证逻辑 #

你制作reprex时问题消失，是因为reprex中你大概率是规范地使用了新数据框的列名，而非跨数据框引用——这也反向验证了原代码的引用方式错误。

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