Vẽ biểu đồ hộp với package ggplot 2 - b1

Trong bài này mình xin được giới thiệu cách vẽ biểu đồ hộp trong r với package ggplot2 mà mình mới tìm hiểu và thực hành. Biểu đồ hộp cung cấp rất nhiều thông tin về phân bố của một biến số kèm theo những chỉ số thống kê cơ bản. Biểu đồ hộp dưới đây mô tả dữ liệu về một số chỉ tiêu sinh trưởng (tree height, stump diameter...) của cây keo lá liềm trồng trên vùng cát ven biển các tỉnh Bắc Trung Bộ mà mình và nhóm đang thực hiện. Biểu đồ hộp dưới đây cung cấp cho chúng ta một số thông tin như:

• Trung vị thể hiện bằng thanh ngang màu đậm, nằm giữa thanh bách phân vị thứ 75 trên cùng và thanh bách phân vị thứ 25 dưới cùng của biểu đồ hộp.
• Bách phân vị thứ 75, thanh ngang màu nhạt hơn nằm phía trên của thanh trung vị. Kết quả cho thấy bách phân vị thứ 75 là 1,51 m (hình 1, TN) và 3,66 cm (hình 1, DC). Nghĩa là có 75% đối tượng cây keo lá liềm điều tra có sinh trưởng chiều cao cây là 1,51 m ở công thức thí nghiệm (TN) và 75% đối tượng điều tra có sinh trưởng chiều cao cây là 1,04 m ở công thức đối chứng (DC).
• Bách phân vị thứ 25, thanh ngang màu nhạt hơn nằm phía dưới của thanh trung vị. Kết quả cho thấy, có 25% đối tượng điều tra có sinh trưởng chiều cao cây đạt 0,74 m (TN) và 0,70 m (DC).

# library(ggplot2)

> rph1

> attach(rph1)

> head(rph1)

> b=ggplot(data=rph1, aes(x=rph1$CTTN, y=Dt))

> b+geom_boxplot(aes(fill=CTTN))+theme_bw() + theme_classic() + ylab("Hvn")+xlab("CTTN")+theme(legend.position="top")

# result Fig.1

# Thay đổi xlab, ylab

> b=ggplot(data=rph1, aes(x=rph1$CTTN, y=Dt))

> b+geom_boxplot(aes(fill=CTTN))+theme_bw() + theme_classic() + ylab("Duong kinh tan")+xlab("CTTN")+theme(legend.position="top")

# result Fig.2

# Vẽ biểu đồ tương quan

> p =ggplot(data=rph1, aes(x=Dgoc, y=Hvn))

> p+geom_point(aes(col=CTTN, size=Sothan))+ geom_smooth(method="lm")+theme_bw()+ theme_classic()+ xlab("Dgoc") + ylab("Hvn")

# result Fig.3