r - 標準エラー シェーディングで ggplot2 stat_ecdf プロットを作成する

Question

3回の治療からのデータがあり、各回の投与ごとに3つの複製があります。

df <- data.frame(dose=c(rep(1,300),rep(3,300),rep(5,300)),
replicate=rep(c(rep("X1",100),rep("X2",100),rep("X3",100)),3),
value=c(rnorm(300,1,1),rnorm(300,3,1),rnorm(300,5,1)),stringsAsFactors=F)

df$dose <- factor(df$dose,levels=c(1,3,5))

cdf プロットを使用して表示したい。各複製ごとに、次のように 3 つの用量の cdfs を簡単にプロットできます。

for(r in c("X1","X2","X3")){
 ggplot(dplyr::filter(df,replicate==r),aes(x=value,color=dose))+
   stat_ecdf(geom="step")+ 
   theme_bw()+
   theme(panel.border=element_blank(),strip.background=element_blank())
}

しかし、各用量のすべての複製を 1 つの図に表示する方法を探していますstat_smooth。

これは達成できますか？

また、これまたは単一のレプリケートのプロットのいずれか:

r <- "X1"
ggplot(dplyr::filter(df,replicate==r),aes(x=value,color=dose))+
  stat_ecdf(geom="step")+ 
  theme_bw()+
  theme(panel.border=element_blank(),strip.background=element_blank())

各 ecdfs の下の領域を計算する方法はありますか?