Sử dụng package ggplot2, ggExtra và gridExtra để vẽ biểu đồ

Trong phạm vi bài này, mình có sử dụng kết hợp 3 packages gồm: ggplot2, ggExtra và gridExtra để vẽ biểu đồ dạng điểm, tương quan, kết hợp giữa biểu đồ tương quan với type histogram hoặc type density... để tạo nên những biểu đồ sinh động, đẹp, logic và khoa học hơn. Package ggplot2 là một gói rất hữu ích để vẽ biểu đồ trong nghiên cứu để thể hiện, biểu thị các kết quả nghiên cứu, những thông số khoa học... mà không đơn giản chúng ta chỉ dùng những lời văn để có thể diễn tả hết được. Để tạo nên một biểu đồ đẹp hơn, khoa học hơn nữa, việc kết hợp thêm các layer, các packages khác là rất cần thiết.

# library(ggplot2)

# Nếu sử dụng ggplot2 (result Fig.1)

> s=ggplot(data=st_nqb1, aes(x=D1.3,y=Hvn))

>  p=s+ geom_point(size=2, color="blue")+ geom_smooth(color="red")+xlab("diameter, cm") + ylab("tree height, m")+ theme_bw(base_size =13)+ggtitle("Fig.1")

# result Fig.1

# Nếu kết hợp ggplot2 + ggExtra 

# library(ggExtra)

> ggMarginal(p, type="histogram", xparams=list(binwidth=0.3,color="white",fill="green"), yparams=list(binwidth=0.2, color="white"))

# result Fig.2

# hoặc type =” density”

# result Fig.3

# Nếu kết hợp ggplot2 + gridExtra 

# library(gridExtra)

> pMain=ggplot(st_nqb1, aes(x=D1.3, y=Hvn))+geom_point(size=1.5, color="green")+geom_smooth(color="red")+ggtitle("Fig.4")

> pTop=ggplot(st_nqb1, aes(x=D1.3))+geom_histogram()

> pRight=ggplot(st_nqb1, aes(x=Hvn))+ geom_histogram()+ coord_flip()

> pEmpty=ggplot(st_nqb1, aes(x=D1.3, y=Hvn))+ geom_blank()+ theme(axis.text=element_blank(), axis.title=element_blank(),line=element_blank(), panel.background=element_blank())

> grid.arrange(pTop,pEmpty, pMain,pRight, ncol=2, nrow=2, widths=c(3,1), heights=c(1,3))

#result Fig.4