ggplot2に正規分布をスーパーインポーズする：Bug？

Question

2つの正規分布を重ね合わせたいと思います。私は狭いintervallにXLIM制限する場合は、プロットは非常に奇妙に見える：

library(ggplot2) 
x=seq(1,6,.01) 
dat <- data.frame(
    c1 = dnorm(x,4.95, .72), 
    c2 = dnorm(x,4.85, .91), x = x 
) 
    ggplot(data=dat, aes(x=x)) + 
    geom_polygon(aes(y=dat[,1]), fill="red", alpha=0.6) + 
    geom_polygon(aes(y=dat[,2]), fill="blue", alpha=0.6) + 
    xlab("Scale") + ylab("") + xlim(1,6) 

範囲を大きくすると、それが動作：

x=seq(1,10,.01) 
    dat <- data.frame(
     c1 = dnorm(x,4.95, .72), 
     c2 = dnorm(x,4.85, .91), x = x 
    ) 
     ggplot(data=dat, aes(x=x)) + 
     geom_polygon(aes(y=dat[,1]), fill="red", alpha=0.6) + 
     geom_polygon(aes(y=dat[,2]), fill="blue", alpha=0.6) + 
     xlab("Scale") + ylab("") + xlim(1,10) 

（出力については図を参照）。しかし、x軸を1から6に制限したいと思います。どんな助けでも大歓迎です！ また、プロットの右側にラベルを追加する方法はありますか？

Answer 1

これはバグではありません.y軸がx軸に再結合する前に範囲を切り捨てるので、ポリゴンが最初の点に折り返されます。 X軸にポリゴンを配置する場合は、geom_ribbonを使用できます。

ggplot(dat) + 
    aes(x, ymin = 0) + 
    geom_ribbon(aes(ymax = c1), fill = 'red', alpha = 0.6) + 
    geom_ribbon(aes(ymax = c2), fill = 'blue', alpha = 0.6) + 
    xlim(1, 6) 

Answer 2

geom_areaとposition = "identity"を使用できます。ところで、

library(tidyverse) 

ggplot(dat %>% gather(var, val, -x), aes(x = x, y = val, fill = var)) + 
    geom_area(alpha = 0.6, position = 'identity') + 
    labs(x = "Scale", y = NULL, fill = 'Variable')