我有两组来自模拟的 100.000 个观察值。由于这两个案例之一是“基线”案例，另一个是“治疗”案例，我想创建一个突出显示两个模拟分布差异的图。
我从两个群体的 ecdf() 开始。结果在图片中。
我想做的是绘制两条 ecdf 曲线之间的差异图。
一个简单的 ecdf(baseline) - ecdf(treatment) 不起作用，因为 ecdf 返回一个函数；即使使用 Hmisc 包中的 Ecdf 也不起作用，因为 Ecdf 返回一个列表，并且在这种情况下，differene '-' 运算符的定义不明确。
通过运行此代码，您可以到达上图描述的场景
a <- runif(10000)
b <- rnorm(10000,0.5,0.5)
plot(ecdf(a))
lines(ecdf(b), col='red')
任何提示都会受到欢迎。