R:以不同颜色在另一个ECDF上绘制一个ECDF

R:以不同颜色在另一个ECDF上绘制一个ECDF,r,plot,ecdf,R,Plot,Ecdf,我有两个累积的经验密度函数,我想把它们相互叠加,以说明这两条曲线的不同。正如前面所指出的,绘制ECDF的函数只是plot(ECDF()),当我阅读精细手册页面时,我确定我可以使用如下方法绘制多个ECDF: require( Hmisc ) set.seed(3) g <- c(rep(1, 20), rep(2, 20)) Ecdf(c( rnorm(20), rnorm(20)), group=g) require(Hmisc) 种子(3) 如果记忆有用的话,我过去也做过。我记得,

我有两个累积的经验密度函数,我想把它们相互叠加,以说明这两条曲线的不同。正如前面所指出的,绘制ECDF的函数只是
plot(ECDF())
,当我阅读精细手册页面时,我确定我可以使用如下方法绘制多个ECDF:

require( Hmisc )
set.seed(3)
g <- c(rep(1, 20), rep(2, 20))  
Ecdf(c( rnorm(20), rnorm(20)), group=g)
require(Hmisc)
种子(3)

如果记忆有用的话,我过去也做过。我记得,您需要欺骗它,因为
Ecdf()
的参数化非常糟糕。我认为在
help(ecdf)
中,它暗示这只是一个step函数的绘图,因此您可以估计两个或多个ecdf,绘制一个,然后通过
lines()进行注释

编辑结果表明,它与

  R> Ecdf(c(rnorm(20), rnorm(20)), group=g, col=c('blue', 'orange'))

正如帮助页面清楚地说明的那样,
col=
参数。但我也发现了一些scriptlet,其中我明确使用了
plot.stepfun()

您可以一次添加一条曲线(每条曲线都有自己的样式),例如


不使用
Ecdf
(看起来不像
Hmisc
可用):

set.seed(3)

是时候从吉尼斯世界纪录切换回咖啡了,JD。有一个
col=
参数,请参阅我编辑的答案。您也可以通过
col
lty
。。。都是矢量化的。手册中有明确的说明,我想我已经读过了。显然我的吉尼斯杯挡住了屏幕的那一部分。哇,他们现在有64盎司的“大口”杯吗?当我在都柏林着陆时,他们把我吸引到一个静脉注射袋里。是静脉注射袋挡住了我的屏幕。
Ecdf(rnorm(20), lwd = 2)
Ecdf(rnorm(20),add = TRUE, col = 'red', lty = 1)
set.seed(3)
mat <- cbind(rnorm(20), rnorm(20))
matplot(apply(mat, 2, sort), seq(20)/20, type='s')