Kết hợp 3 packages (ggplot2, ggthemes và gridExtra) để vẽ biểu đồ

Sự hoàn chỉnh của một biểu đồ phản ảnh những ý nghĩa mà những lời văn khó mà diễn tả hết được. “... một hình ảnh có giá trị bằng 1.000 chữ. Hình ảnh và biểu đồ dễ gây ấn tượng hơn là con số hay chữ, ...” (1). Trong phạm vi bài viết chưa hẳn đã là một biểu đồ hoàn chỉnh, chỉ là ví dụ về sự kết hợp giữa các GÓI khác nhau trong R để tạo nên một biểu đồ tương đối ĐẸP, LOGIC và quan trọng hơn là, ĐỘC ĐÁO, KHOA HỌC có thể dùng cho việc in ấn, công bố hoặc trình diễn kết quả nghiên cứu, ...

Thứ nhất

library(ggplot2)

 > p=ggplot(kll_Quang_Tri, aes(x=dg, y=hvn_m))+ geom_point(aes(color=Age, size=Thanchinh))+ geom_smooth(method="lm")+ theme_bw()+ggtitle("Fig.1")

 # result Fig.1

Hoặc             

> p=ggplot(kll_Quang_Tri, aes(x=dg, y=hvn_m))+ geom_point(aes(color=Age, size=Thanchinh))+ geom_smooth(method="lm")+ theme_bw()+ggtitle("Fig.2")+ xlab("Đường kính gốc, cm")+ ylab("Chiều cao cây, m")+ theme_classic()

# result Fig.2

Nếu sử dụng gói ggplot2 thì kết quả tương đối như hình 1 và hình 2 bên trên. Tuy nhiên, nếu kết hợp thêm gói ggthemes sẽ LÀM ĐẸP hơn cho ggplot2.

Thứ hai

library (ggthemes)

 > p=ggplot(kll_Quang_Tri, aes(x=dg, y=hvn_m))+ geom_point(aes(color=Age, size=Thanchinh))+ geom_smooth(method="lm")+ theme_bw()+ xlab("Đường kính gốc, cm")+ ylab("Chiều cao cây, m")+ theme_classic()

> p1=p+ geom_rangeframe()+ theme_tufte()+ ggtitle("Fig.3")

 # chú ý đến thang trị số ở trục tung và trục hoành để thấy sự khác biệt giữa hình 3, 4 và hình 1, 2.

# result Fig.3

Hoặc

> p1=p+ geom_rangeframe()+ theme_tufte()+scale_x_continuous(breaks=extended_range_breaks()(kll_Quang_Tri$dg))+scale_y_continuous(breaks=extended_range_breaks()(kll_Quang_Tri$hvn_m))+ ggtitle("Fig.4")

# result Fig.4

Kết quả như hình 4 là tương đối đẹp cũng có thể in ấn, công bố được. Tuy nhiên, biểu đồ mới thể hiện phân bố điểm (point) trên trục tung và trục hoành theo độ tuổi (1,25-2,25) và theo số thân chính (từ 2-8). Trên cùng một biểu đồ, vừa thể hiện phân bố điểm (point) vừa thể hiện phân bố histogram (hình 5), phân bố đường (hình 6), thể hiện được tính ĐỘC ĐÁO và KHOA HỌC thì gridExtra là gói bạn cần.

Thứ ba

library (gridExtra)

Với code như sau:

> ggMarginal(p1, type="histogram", color="white", size=2)

# result Fig.5

Hoặc

> ggMarginal(p1, type="density", size=4, color="red")

 # result Fig.6

(1) Source: http://tuanvannguyen.blogspot.com/