GGPLOT2：設置（非線性）用於α-

Question

值我想繪製鏡像95％的密度曲線和映射α位密度：

foo <- function(mw, sd, lower, upper) { 
x <- seq(lower, upper, length=500) 
dens <- dnorm(x, mean=mw, sd=sd, log=TRUE) 
dens0 <- dens -min(dens) 
return(data.frame(dens0, x)) 
} 

df.rain <- foo(0,1,-1,1) 

library(ggplot2) 


drf <- ggplot(df.rain, aes(x=x, y=dens0))+ 
geom_line(aes(alpha=..y..))+ 
geom_line(aes(x=x, y=-dens0, alpha=-..y..))+ 
stat_identity(geom="segment", aes(xend=x, yend=0, alpha=..y..))+ 
stat_identity(geom="segment", aes(x=x, y=-dens0, xend=x, yend=0, alpha=-..y..)) 
drf 

這工作得很好，但我想使邊緣和中間的對比更加突出，即我希望邊緣幾乎是白色的，只有中間部分是黑色的。我一直在篡改scale_alpha()，但沒有運氣。有任何想法嗎？

編輯：最終，我想繪製幾個雨滴，即單個滴劑會很小，但陰影仍應清晰可見。

Answer 1

雖然都琢磨你的答案我居然發現正是我一直在尋找。最簡單的方法是簡單地使用帶有灰色矢量的scale_colour_gradientn。

library(RColorBrewer) 
grey <- brewer.pal(9,"Greys") 

drf <- ggplot(df.rain, aes(x=x, y=dens0, col=dens0))+ 
stat_identity(geom="segment", aes(xend=x, yend=0))+ 
stat_identity(geom="segment", aes(x=x, y=-dens0, xend=x, yend=0))+ 
scale_colour_gradientn(colours=grey) 
drf 

Answer 2

不是映射dens0到alpha的，我想它映射到color：

drf <- ggplot(df.rain, aes(x=x, y=dens0))+ 
    geom_line(aes(color=..y..))+ 
    geom_line(aes(x=x, y=-dens0, color=-..y..))+ 
    stat_identity(geom="segment", aes(xend=x, yend=0, color=..y..))+ 
    stat_identity(geom="segment", aes(x=x, y=-dens0, xend=x, yend=0, color=-..y..)) 

現在我們仍然有顏色的反差主要存在於尾部。使用兩種顏色有助於一點（注意，在顏色開關處於0.25）：

drf + scale_color_gradient2(midpoint = 0.25) 

最後，以包括dens0值的分佈，I類基礎顏色的中點比例上在數據的中間值：

drf + scale_color_gradient2(midpoint = median(df.rain$dens0)) 

注意！：但是，無論您如何調整數據，數據中的大多數對比度都是數據集中更極端的值。試圖通過混淆非線性尺度或通過像我這樣調整顏色尺度來掩蓋這一點，可能會呈現真實數據的錯誤畫面。

Answer 3

下面是一個使用的解決方案geom_ribbon（）而不是geom_line（）

df.rain$group <- seq_along(df.rain$x) 
tmp <- tail(df.rain, -1) 
tmp$group <- tmp$group - 1 
tmp$dens0 <- head(df.rain$dens0, -1) 
dataset <- rbind(head(df.rain, -1), tmp) 
ggplot(dataset, aes(x = x, ymin = -dens0, ymax = dens0, group = group, 
    alpha = dens0)) + geom_ribbon() + scale_alpha(range = c(0, 1)) 

ggplot(dataset, aes(x = x, ymin = -dens0, ymax = dens0, group = group, 
    fill = dens0)) + geom_ribbon() + 
    scale_fill_gradient(low = "white", high = "black") 

看到了保羅的回答改變顏色。

dataset9 <- merge(dataset, data.frame(study = 1:9)) 
ggplot(dataset9, aes(x = x, ymin = -dens0, ymax = dens0, group = group, 
    alpha = dens0)) + geom_ribbon() + scale_alpha(range = c(0, 0.5)) + 
    facet_wrap(~study)