使用 latticeExtra 的加权经验累积分布图

Question

优秀latticeExtra package has a function ecdfplot(), which creates a cumulated empirical distribution plot。

可重现的例子：

df= data.frame(X= 1:100, W= 100:1)
library(latticeExtra)
ecdfplot(~X, data=df)

问题：

我需要使用 df$W 变量来应用权重，而不是根据观察频率进行计数。

使用 Hmisc 可以做到这一点：

library(Hmisc)
Ecdf(df$X, weights= df$W)

问题：

我想做同样的事情，但使用 lattice 图表系统。

如何做到这一点？

Answer 1

由于W是整数类型，解决方法非常简单。

ecdfplot(~ rep(x= X, times= W), data=df, xlab= "X")

Weighted empirical cumulated distribution plot with latticeExtra