绘制两个ecdf()之间的差额

Question

我有两组来自模拟的100.000次观测结果。由于这两种情况中的一种是“基线”案例，另一种是“治疗”案例，我想要创建一个图来突出这两种模拟的分布差异。

我从这两个种群的ecdf()开始。结果在图中。

​

我想要做的是绘制两个ecdf曲线之间的差分图。

一个简单的ecdf (基线)- Ecdf (处理)不起作用，因为Ecdf返回一个函数；即使使用Hmisc包中的Ecdf也不起作用，因为Ecdf返回一个列表，并且在这种情况下差异'-‘操作符也是定义错误的。

通过运行这段代码，您可以看到上面所描述的场景。

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)