绘制两个ecdf（）之间的差异

Question

我有两组来自模拟的100.000个观测值。 由于这两个案例中的一个是基线＆＃39;案件和另一个是一个＆＃39;治疗＆＃39;例如，我想创建一个图表，突出显示两个模拟的分布差异。

我从两个人群的ecdf（）开始。结果如图所示。

我想要做的是在两条ecdf曲线之间绘制差异的图。

简单的ecdf（基线） - ecdf（处理）不起作用，因为ecdf返回一个函数;即使使用Hmisc软件包中的Ecdf也行不通，因为Ecdf会返回一个列表并再次输入差异＆＃39; - ＆＃39;在这种情况下，运算符定义不明确。

通过运行此代码，您可以进入上图中描述的场景

a <- runif(10000)
b <- rnorm(10000,0.5,0.5)
plot(ecdf(a))
lines(ecdf(b), col='red')

任何提示都非常受欢迎。

Answer 1

那么评估一下这些功能？

decdf <- function(x, baseline, treatment)  ecdf(baseline)(x) - ecdf(treatment)(x)