Thiết kế tháp dân số bằng phần mềm R

Chỉ cần có số liệu dân số theo tuổi và giới tính trong một khoảng thời gian nào đó, ta có thể thiết kế một tháp dân số đẹp mắt bằng phần mềm thống kê R, đây là phần mềm thống kê càng được nhiều người đón nhận hiện nay. Tôi có số liệu dân số Việt Nam năm 1979, từ số liệu này ta dễ dàng tóm lược lại theo bảng dưới đây bằng Excel:

          Để thiết kế được tháp dân số trong R ta cần hỗ trợ gói plotrix và một số mã lệnh sau:

# Gọi gói plotrix vận hành trong R

>library(plotrix)

# Nhập tỷ lệ phần trăm Nam giới theo nhóm tuổi:

>Nam=c(7.48,7.45,6.89,5.60,4.32,3.30,2.23,1.83,1.74,1.89,1.56,1.29,1.03,0.79,0.54,0.35,0.12,0.08)

# Nhập tỷ lệ phần trăm Nam giới theo nhóm tuổi:

>Nu=c(7.14,7.13,6.46,5.80,4.93,3.75,2.49,2.09,2.06,2.11,1.71,1.66,1.26,1.06,0.82,0.59,0.26,0.17)

# Nhập nhóm tuổi

>Tuoi=c("0-4","5-9","10-14","15-19","20-24","25-29","30-34","35-39","40-44","45-49","50-54","55-59","60-64","65-69","70-74","75-79","80-44","85+")

# Tạo hiệu ứng màu theo giới và tuổi cho đẹp mắt.

>mcol=color.gradient(c(0,0,0.5,1),c(0,0,0.5,1),c(1,1,0.5,1),22)

>fcol=color.gradient(c(1,1,0.5,1),c(0.5,0.5,0.5,1),c(0.5,0.5,0.5,1),22)

>par(mar=pyramid.plot(Nam,Nu,labels=Tuoi,labelcex=0.8, ndig=2,main="THÁP TUỔI NĂM 1979",lxcol=mcol,rxcol=fcol, gap=1,show.values=TRUE))

Ghi chú: có thể thay đổi labelcex để thay đổi cỡ chữ, ndig để thay đổi số thập phân của tỷ lệ phần trăm, gap đề thay độ khoảng cách nhóm tuổi.

Trước khi chỉnh sửa bởi Paint

Tháp dân số trên có thể chấp nhận được nhưng còn có tiếng Anh (do gói plotrix mặc định) và lỗi phông tiếng việt. Để khắc phục có lẽ nhờ đến vài thủ thuật của chương trình Paint thì ta có được một tháp dân số đẹp hơn.

Sau khi chỉnh sửa bởi Paint

Bs Bùi Văn Dủ