Vẽ biểu đồ với ggplot2

Bài này mình giới thiệu cách xác định vùng phân bố tập trung của đối tượng điều tra bằng package ggplot2. Thay vì, xác định phân bố type histogram thông thường, hay phân bố point kết hợp với type histogram hoặc kết hợp với type density trong package ggExtra. Trong phạm vi bài này sẽ kết hợp thêm layer (density2d) để xác định vùng phân bố tập trung của đối tượng điều tra.

#>rungtrong

#>attach(rungtrong)

#>head(rungtrong)

#>library(ggplot2)

> r=ggplot(data=rungtrong, aes(x=d,y=h))+geom_point()+geom_density2d()+ggtitle("Fig 1")

# result Fig. 1

# Hoặc

> r=ggplot(data=rungtrong, aes(x=d,y=h))+geom_point()+stat_density2d(aes(alpha=..level.., fill=..level..), size=2)

# result Fig.2

# Kết quả

 > r=ggplot(data=rungtrong, aes(x=d,y=h))+geom_point()+stat_density2d(aes(alpha=..level.., fill=..level..), size=2,geom="polygon")+ scale_fill_gradient(low="blue", high="red")+ ggtitle("Fig.3")+theme_bw()+theme_classic()+xlab("Diameter, cm")+ ylab("Top height tree, m") + geom_smooth()

# result Fig. 3