优秀的latticeExtra包有一个函数ecdfplot()
,可以创建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图形系统。
如何做到这一点?
答案 0 :(得分:0)
由于W
的类型为整数,因此解决方案非常简单。
ecdfplot(~ rep(x= X, times= W), data=df, xlab= "X")