R语言 双坐标轴组合图形可视化实现

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. 1、plotrix包 twoord.plot用法和参数解释
   
2. 
   
3. # 1、用法/Usage：
   
4. twoord.plot(lx,ly,rx,ry,data=NULL,main="",xlim=NULL,lylim=NULL,rylim=NULL,
   
5. mar=c(5,4,4,4),lcol=1,rcol=2,xlab="",lytickpos=NA,ylab="",ylab.at=NA,
   
6. rytickpos=NA,rylab="",rylab.at=NA,lpch=1,rpch=2,
   
7. type="b",xtickpos=NULL,xticklab=NULL,halfwidth=0.4,axislab.cex=1,
   
8. do.first=NULL,xaxt="s",...)
   
9. 
   
10. # 2、参数解释：
    
11. 1、lx,ly,rx,ry：左坐标轴和右坐标轴的值
    
12.    注意：lx和rx值至少相邻，并且可能重叠
    
13. 2、data：数据框
    
14. 3、main：标题
    
15. 4、xlim、lylim,rylim：限制横坐标值、左右纵坐标值的范围
    
16. 5、mar：设置图形边界距，默认值为(5,4,4,4)
    
17. 6、lcol、rcol：左右坐标轴的颜色
    
18. 7、xlab：设置横坐标轴标签
    
19. 8、lytickpos：设置左坐标轴刻度标签的位置
    
20. 9、ylab：设置左坐标轴标签
    
21. 10、ylab.at：设置左坐标轴标签位置
    
22. 11、rytickpos：设置右坐标轴刻度标签的位置
    
23. 12、rylab：设置右坐标轴标签
    
24. 13、rylab.at：设置右坐标轴标签位置
    
25. 14、type：指定图形类型 type
    
26.   "p" for points,"l" for lines,"b" for both,"c" for the lines part alone of "b",
    
27.   "o" for both ‘overplotted’,"h" for ‘histogram’ like (or ‘high-density’) vertical lines,
    
28.   "s" for stair steps,"S" for other steps, see ‘Details’ below,
    
29.   "n" for no plotting.
    
30. 15、xtickpos：设置横坐标轴刻度标签位置
    
31. 16、xticklab：设置横坐标轴刻度标签
    
32. ……
    
33. 详见R-gui帮助文档，代码：help(plotrix)
    
34. 
    
35. 
    
36. ### 声  明：本内容为作者借助R3.6.3和Rstudio及相关包制作而成，仅供学习交流，咨询交流加wx：huyanggs 或Email：huyanggs@hotmail.com
    
37. ### 主  题：双(Y)坐标轴制图
    
38. ### 数据源：随机生成的数据集，未设置set.seed()
    
39. ### author：@拴小林Nobeli
    
40. ### 时  间：2020/3/29
    
41. 
    
42. ## 1、加载plotrix包
    
43. # install.packages("plotrix")
    
44. library(plotrix)
    
45. 
    
46. 
    
47. ## 2、简单双坐标折线图(x周数值包含)
    
48. xval1 <- 2:10
    
49. xval2 <- 1:15
    
50. going_up <- seq(3,7,by=0.5)+rnorm(9)
    
51. going_down <- rev(60:74)+rnorm(15)
    
52. 
    
53. twoord.plot(lx = xval1,ly = going_up, rx = xval2, ry = going_down,xlab="Sequence",
    
54.             ylab="Ascending values",rylab="Descending values",lcol=4,
    
55.             main="Plot with two ordinates - points and lines",
    
56.             do.first="plot_bg();grid(col=\"white\",lty=2)")
    
57. 
    
58. 
    
59. ## 3、简单双坐标折线图(x周数值重叠一致)
    
60. xval <- seq.Date(as.Date("2017-01-01"),
    
61.                  as.Date("2017-01-15"), by="day")
    
62. 
    
63. going_up<-seq(3,7,length.out = 15)+rnorm(15)
    
64. going_down<-rev(60:74)+rnorm(15)
    
65. twoord.plot(xval, going_up, xval, going_down,xlab="Sequence",
    
66.             ylab="Ascending values",rylab="Descending values",lcol=4,
    
67.             main="Plot with two ordinates - points and lines",
    
68.             xtickpos=as.numeric(xval), xticklab=as.character(xval),
    
69.             do.first="plot_bg();grid(col=\"white\",lty=1)")
    
70. # 注释：xtickpos=as.numeric(xval), xticklab=as.character(xval)两个参数是控制x轴显示日期的关键，请知晓~
    
71. 
    
72. 
    
73. ## 4、折线&柱形混合双坐标图形
    
74. twoord.plot(xval1, going_up, xval2, going_down, xlab="Sequence",
    
75.             lylim=c(-3,8), rylim=c(50,100), type=c("l","bar"),
    
76.             ylab="Ascending values", rylab="Descending values",
    
77.             text = list(going_up,going_down),
    
78.             lytickpos=3:7, rytickpos=seq(55,75,by=5), ylab.at=5,rylab.at=65,
    
79.             main="Lines on left axis, bars on right axis(含点数值标签）",
    
80.             lcol="red", rcol=4,do.first="plot_bg(\"white\")")
    
81. 
    
82. ## 4.添加数值标签到图表中（补充：此部分为@拴小林Nobeli编写，因没有找到x相应参数，如有相关参数，望交流，感谢！！）
    
83. #注意：
    
84. # 1.添加left线形图形标签可直接使用text函数，
    
85. # 2.但right柱形图形标签不适合直接使用text标签，需要将右侧ry值等比例转换为ly值
    
86. # 3.转换公式：（ry-min（rylim））/（max（rylim）-min（rylim））* （max（lylim）-min（lylim））+ min（lylim）+位置调整值
    
87. text(xval1, going_up,as.numeric(round(going_up)))
    
88. text(xval2,
    
89.      (going_down-min(c(50,100)))/(max(c(50,100))-min(c(50,100)))*(max(c(-3,8))-min(c(-3,8)))+min(c(-3,8)+0.5),
    
90.      round(going_down))
    
91. 
    
92. 
    
93. 
    
94. ## 5、带概率密度的双坐标图
    
95. xhist<-hist(rnorm(100),plot=FALSE)
    
96. xdens<-dnorm(seq(-3,3,by=0.05))
    
97. twoord.plot(xhist$mids,xhist$counts,seq(-3,3,by=0.05),
    
98.             xdens,type=c("bar","l"),lcol=4,rcol=2,ylab="Counts",
    
99.             rylab="Density",main="Histogram and density curve",
    
100.             halfwidth=0.2,lylim=c(0,max(xhist$counts)+1),rylim=c(0,0.45),lwd=2)

复制代码

扫码加我 拉你入群

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

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

分享0 收藏1 回帖