ggplot2疊加正態分佈：Bug？

Question

我想疊加兩個正態分佈。但是，如果我限制XLIM一個狹窄的INTERVALL，情節看起來很怪異：

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')