r作图 如何将多个变量的累积经验分布曲线画在一张图里面？O(∩_∩)O谢谢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如下图，想实现多个累积经验分布曲线画在一张图里且颜色不同，带有图例的效果。
现在只能画出一张图一条曲线。
有什么好的办法吗？
求大神指点！O(∩_∩)O谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：经验分布 多个变量 经验分 张图里 r语言 作图 累积经验分布图

露泫低珠佩3 发表于 2017-6-15 09:50
如下图，想实现多个累积经验分布曲线画在一张图里且颜色不同，带有图例的效果。
现在只能画出一张图一条曲 ...

http://mp.weixin.qq.com/s/0Zm1QBoLTs-1LRUztvR2UA

练手，按原图用基本绘图函数尽量复制

1. legend.text <- c(expression(paste(mu, "=0", ", ", sigma^2, "=0.2")),
   
2.                  expression(paste(mu, "=0", ", ", sigma^2, "=1.0")),
   
3.                  expression(paste(mu, "=0", ", ", sigma^2, "=5.0")),
   
4.                  expression(paste(mu, "=-2", ", ", sigma^2, "=0.5")))
   
5. ylab.text <- expression(paste(Phi[mu], " ", NULL[sigma^2], "(x)"))
   
6. plot(x = c(-5, 5), y = c(0, 1),   # 确定范围并留出与坐标之间的空隙
   
7.      ylab = NA, xlab = NA,    # 不要坐标说明
   
8.      type = "n",                     # 坐标图上留白
   
9.      xaxt = "n", yaxt = "n")    # 不要刻度等信息
   
10. # 画坐标轴数字和刻度
    
11. ## 下x轴
    
12. axis(side = 1, at = -5:5, tck = 1, col = "grey70", mgp = c(0, 0.2, 0)) # 画主网格线
    
13. axis(side = 1, at = c(-5, 5), labels = FALSE, tcl = 0) # 把灰化的坐标轴重新黑化
    
14. axis(side = 1, at = -4.5:4.5,  labels = FALSE, tcl = 0.5) # 画次级刻度
    
15. axis(side = 1, at = seq(-4.75, 4.75, 0.5),  labels = FALSE, tcl = 0.25) # 画再次级刻度
    
16. ## 左y轴
    
17. axis(side = 2, at = seq(0, 1, 0.2), tck = 1, col = "grey70", mgp = c(0, 0.2, 0), las = 1)
    
18. axis(side = 2, at = c(0, 1), labels = FALSE, tcl = 0)
    
19. axis(side = 2, at = seq(0.1, 0.9, 0.1),  labels = FALSE, tcl = 0.5)
    
20. axis(side = 2, at = seq(0.05, 0.95, 0.1),  labels = FALSE, tcl = 0.25)
    
21. ## 上x轴
    
22. axis(side = 3, at = -4.5:4.5,  labels = FALSE, tcl = 0.5)
    
23. axis(side = 3, at = seq(-4.75, 4.75, 0.5),  labels = FALSE, tcl = 0.25)
    
24. ## 右y轴
    
25. axis(side = 4, at = seq(0.1, 0.9, 0.1),  labels = FALSE, tcl = 0.5)
    
26. axis(side = 4, at = seq(0.05, 0.95, 0.1),  labels = FALSE, tcl = 0.25)
    
27. # 画曲线
    
28. curve(pnorm(x, 0, sqrt(0.2)), from = -5, to = 5, n = 1001, add = TRUE, lwd = 3, col = "blue")
    
29. curve(pnorm(x, 0, sqrt(1)), from = -5, to = 5, n = 1001, add = TRUE, lwd = 3, col = "red")
    
30. curve(pnorm(x, 0, sqrt(5)), from = -5, to = 5, n = 1001, add = TRUE, lwd = 3, col = "orange")
    
31. curve(pnorm(x, -2, sqrt(0.5)), from = -5, to = 5, n = 1001, add = TRUE, lwd = 3, col = "green")
    
32. # 写坐标名称
    
33. title(xlab = "x", ylab = ylab.text, cex.lab = 1.8, mgp = c(1.7, 0, 0), font.lab = 3)
    
34. # 写图例
    
35. legend(x = -4.99, y = 0.99, legend = legend.text, border = NULL, lwd=2, col=c("blue", "red", "orange", "green"), bty="0", bg = "white", box.col = "white", cex = 1.5)

复制代码

都是高手呀 又学到了新知识