加权经验累积分布图与latticeExtra

时间:2015-03-29 14:19:44

标签: r lattice

优秀的latticeExtra包有一个函数ecdfplot(),可以创建cumulated empirical distribution plot

可再现的例子:

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

enter image description here

问题:

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

使用Hmisc可以完成此操作:

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

enter image description here

问题:

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

如何做到这一点?

1 个答案:

答案 0 :(得分:0)

由于W的类型为整数,因此解决方案非常简单。

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

enter image description here