小编alb*_*ren的帖子

我正在尝试使用facet_wrap 在ggplot2 中制作多边形地图。我的变量“作物”中有两个因子水平（大豆、玉米），但是，我得到了三个图：大豆、玉米和一个具有 NA 值的图。此外，前两个方面不显示 NA 值 -

这是我制作地图的代码：

ggplot(study_area.map, aes(x=long, y=lat, group=group)) + 
  geom_polygon(aes(fill=brazil_loss_new2)) + 
  geom_path(colour="black") + 
  facet_wrap(~crop, ncol=2, drop=T) + 
  scale_fill_brewer(na.value="grey", palette="Blues", 
    name="Average production lossess\n per municipality", 
    breaks = levels(study_area.map$brazil_loss_new2), 
    labels = levels(study_area.map$brazil_loss_new2)) + 
  theme() + 
  coord_fixed()

这就是我得到的：

如果我使用 na.omit 我会得到下图（更好，但前两张图中仍然缺少 NA 值）

在此输入图像描述

无论感兴趣的变量是否为 NA，都包括每个变量和城市的行，最终解决问题。这是我一直在寻找的：

具有 NA 值的城市的产量损失