在ggplot2下手动着色平滑线

Question

我正在绘制一组带有阴影的数据，以遵循“下边界”和“上边界”值。 我对Est值使用stat_smooth进行了平滑处理，但无法获得UB和LB值指定的阴影区域来遵循平滑线。 数据如下：

 Quantiles   Est    LB    UB
 0.10 -4.39 -4.80 -4.00
 0.25 -3.46 -3.72 -3.22
 0.50 -3.11 -3.29 -2.91
 0.75 -2.89 -3.15 -2.60
 0.90 -1.69 -2.21 -1.09

这是我的ggplot代码：

ggplot(data,aes(y=Est,x=Quantiles)) + stat_smooth() + 
geom_ribbon(aes(ymin=LB,ymax=UB),alpha=0.2)

在此先感谢您的帮助！

Answer 1

我的理解是，您希望色带边界遵循平滑的曲线，而不是简单地连接LB和UB点。 在您的情况下，stat_smooth使用黄土方法来计算平滑曲线。 您没有足够的数据点来使用默认的2阶计算真实的平滑黄土曲线，因此loess函数返回一条曲线，该曲线穿过给定的数据数据点，使用二次平滑曲线将它们连接起来，并在警告消息中报告此情况。 。 忽略这些警告，获得平滑功能区的一种方法是计算数据平滑点和功能区边界，然后绘制结果，如下所示：

smooth <- data.frame(Quantiles= seq( min(data$Quantiles), max(data$Quantiles), length.out=100))
smooth$Est <- predict(loess(Est ~ Quantiles, data), newdata=smooth$Quantiles)
smooth$LB  <-  predict(loess(LB ~ Quantiles, data), newdata=smooth$Quantiles)
smooth$UB  <-  predict(loess(UB ~ Quantiles, data), newdata=smooth$Quantiles)
ggplot(data=smooth,aes(y=Est,x=Quantiles)) +  geom_line() +
    geom_ribbon(aes(ymin=LB,ymax=UB),alpha=0.2)