小编Hal*_*son的帖子

我正在绘制50-100个实验的结果.每个实验都会产生一个时间序列.我可以绘制所有时间序列的意大利面条图,但我想要的是时间序列羽流的密度图.(类似于下图中下方面板中的灰色阴影:http://www.ipcc.ch/graphics/ar4-wg1/jpg/fig-6-14.jpg)

我可以通过2d binning或binhex"做"这样做,但结果可能更漂亮(参见下面的示例).

这是一个代码,用于再现模拟数据的羽流图(使用ggplot2和reshape2).

# mock data: random walk plus a sinus curve.
# two envelopes for added contrast.
tt=10*sin(c(1:100)/(3*pi))
rr=apply(matrix(rnorm(5000),100,50),2,cumsum) +tt
rr2=apply(matrix(rnorm(5000),100,50),2,cumsum)/1.5 +tt

# stuff data into a dataframe and melt it.
df=data.frame(c(1:100),cbind(rr,rr2) )
names(df)=c("step",paste("ser",c(1:100),sep=""))
dfm=melt(df,id.vars = 1)

# ensemble average
ensemble_av=data.frame(step=df[,1],ensav=apply(df[,-1],1,mean))
ensemble_av$variable=as.factor("Mean")


ggplot(dfm,aes(step,value,group=variable))+
  stat_binhex(alpha=0.2) + geom_line(alpha=0.2) +
  geom_line(data=ensemble_av,aes(step,ensav,size=2))+
  theme(legend.position="none") 

有没有人知道一个很好的方法来获得带渐变的阴影信封.我也试过geom_ribbon,但没有给出任何沿羽流密度变化的迹象.binhex做到了这一点,但没有美学上令人愉悦的结果.