如何在R语言base包barplot中为NA值显示空条形？

在ggplot2中，当数据包含NA值时，geom_col()会自动保留对应分组的标签，同时不绘制该组的条形，呈现空条形效果：

library(ggplot2)
df <- data.frame(group = LETTERS[1:5], value = c(NA, 1, 4, 3, 2))
ggplot(df, aes(x = group, y = value)) +
  geom_col()
#> Warning: Removed 1 rows containing missing values (`position_stack()`).

但直接使用基础包的barplot()函数时，带NA值的分组会被直接忽略，不会在x轴上显示：

df <- data.frame(group = LETTERS[1:5], value = c(NA, 1, 4, 3, 2))
barplot(value ~ group, data = df)

如果将NA替换为0，虽能保留分组标签，但会显示极小的条形，无法达到ggplot2的空条形效果：

df[is.na(df)] <- 0
barplot(value ~ group, data = df)

方法一：将数据转换为矩阵传入barplot #

通过把数据转换为矩阵（保留所有分组作为列名），barplot()会保留所有分组标签，且对NA值对应的位置不绘制条形，完美复刻ggplot2的空条形效果：

df <- data.frame(group = LETTERS[1:5], value = c(NA, 1, 4, 3, 2))
# 构造矩阵，列名为分组标签
bar_matrix <- matrix(df$value, nrow = 1, dimnames = list(NULL, df$group))
# 绘制条形图
barplot(bar_matrix)

方法二：替换NA为0并隐藏对应条形 #

将NA替换为0后，通过设置条形颜色和边框为透明，让0值条形不可见，同时保留分组标签：

df <- data.frame(group = LETTERS[1:5], value = c(NA, 1, 4, 3, 2))
# 把NA替换为0
df$value[is.na(df$value)] <- 0
# 绘制时，给0值条形设置透明样式
barplot(
  value ~ group, 
  data = df,
  col = ifelse(df$value == 0, NA, "steelblue"),
  border = ifelse(df$value == 0, NA, "black")
)

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