我有一个观察到的样本，该样本已建模为7种不同的分布。我想在一个图中显示经验分布函数（Ecdf）和拟合的Ecdf，以便在不带颜色的日记帐中发布时可以识别该图。 当我们允许使用颜色时，我也很想看到答案。 可以从Data Download
下载数据我使用了以下代码，我希望对其进行修改，以使其在黑白打印中更具区别。
#plot x vs F(x)
fun.ecdf <- ecdf(my.obs1)
my.ecdf <- fun.ecdf(sort(my.obs1))
plot(sort(my.obs1),my.ecdf,type="l",lwd=3,xlab="x=area",ylab="F(x)")
lines(sort(my.obs1),y1,lty=2,lwd=2,col=116)
lines(sort(my.obs1),y2,lty=3,lwd=2,col=142)
lines(sort(my.obs1),y3,lty=4,lwd=2,col=96)
lines(sort(my.obs1),y4,lty=5,lwd=2,col=504)
lines(sort(my.obs1),y5,lty=6,lwd=2,col=259)
lines(sort(my.obs1),y6,lty=7,lwd=2,col=373)
lines(sort(my.obs1),y7,lty=8,lwd=2,col=370)
legend(45,0.85, legend = c("empirical distribution","dist A",
"distB","distC","distD","distE","distF","distG"),
col="black",116,142,96,504,259,373,370),
lty = 1:8, cex = 0.8)
答案 0 :(得分：0)
我想您只需要上面的内容，但要用黑白两种颜色，在这种情况下，ggplot2的{{1}}会解决这个问题。
theme_bw