基准R中缺少值的面积图

我想绘制一个面积图，其中多边形的底面为零，并且在每个数据中断处（即开始、结束和可能的NAs/NaN），数据线通过垂直线段连接到底面

我画了这个：

我不得不强制垂直向下的部分，在那里，系列赛被NAs中断，我在0年代就完成了NAs的转型。但这不会生成垂直线段，而是生成到达以下0的多边形线。我解决了系列赛开始和结束时的问题，在系列赛的两边加了一个（y=0，x=0）点

但如果NAs在系列中，这并不能解决问题

有什么想法吗

下面是一个示例代码（不同的图像）：

花粉我会使用ggplot2：

pollen <- c(45, 257.4, 24.67, 54.6, 89.4, 297, 471.25, 1256.5, 312.25, 969.2, 787.5, 425, NaN, 76.6, 42.67, 38.5, 20.2, 5.67, 15.8, 13.2, 11, 6.25, 6.67, 2.3, 0.5, 30.8, 3.75, 3, 2, 2.2, 3.25, 4.5, 9.6, 15.8, 200.2, NaN)

weeks.vec <- c(5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40)


DF <- data.frame(pollen, weeks.vec)

library(ggplot2)
ggplot(DF, aes(x = weeks.vec, y = pollen)) +
  geom_ribbon(aes(ymin = 0, ymax = pollen), 
              colour = NA, fill = "red", alpha = 0.3) +
  geom_line(colour = "red") + 
  geom_point(colour = "red", size = 3) +
  xlab("Week") + ylab("Pollen") +
  theme_bw()

请创建一个最小的可复制示例。添加了代码，如果满意，请删除下一票。基本r代码有效！ggplot给了我这个错误：“在if（class（trial）！=“try error”）return（trial）：条件的长度大于1，只使用第一个元素”并在NA位置绘制点，通过完整的线段连接。ggplot2代码与您提供的示例数据一起工作。如果它不适用于您的真实数据，那么示例数据就不能代表它！我复制粘贴了代码！无论如何，基本的r代码完全可以工作，我一直在等待其他答案来了解更多，其他方面你的答案是绝对正确的！
pollen <- c(45, 257.4, 24.67, 54.6, 89.4, 297, 471.25, 1256.5, 312.25, 969.2, 787.5, 425, NaN, 76.6, 42.67, 38.5, 20.2, 5.67, 15.8, 13.2, 11, 6.25, 6.67, 2.3, 0.5, 30.8, 3.75, 3, 2, 2.2, 3.25, 4.5, 9.6, 15.8, 200.2, NaN)

weeks.vec <- c(5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40)


DF <- data.frame(pollen, weeks.vec)

library(ggplot2)
ggplot(DF, aes(x = weeks.vec, y = pollen)) +
  geom_ribbon(aes(ymin = 0, ymax = pollen), 
              colour = NA, fill = "red", alpha = 0.3) +
  geom_line(colour = "red") + 
  geom_point(colour = "red", size = 3) +
  xlab("Week") + ylab("Pollen") +
  theme_bw()

plot.ts(y = pollen, x = weeks.vec, col = 'red', 
        ylab = 'Pollen',  xlab = 'Weeks', lwd = 3, 
        xy.labels = F, xy.lines = T)

g <- cumsum(!is.finite(pollen))
for (i in unique(g)) {
  y <- pollen[g == i]
  x <- weeks.vec[g == i]
  x <- x[is.finite(y)]
  y <- y[is.finite(y)]
  x <- c(x, rev(x))
  y <- c(y, y * 0)
  polygon(y = y, x = x, density = NA,border = NA, col = rgb(1,0,0, .3))
}