Setup Data
susenas_individu <- read.csv("D:\\IPB\\4\\Visdat\\2023\\2023 Maret JABAR - SUSENAS KOR INDIVIDU PART1.csv",sep=",")
str(susenas_individu)
## 'data.frame': 84688 obs. of 183 variables:
## $ X : int 0 1 2 3 4 5 6 7 8 9 ...
## $ URUT : int 500001 500001 500001 500001 500002 500002 500003 500003 500003 500004 ...
## $ PSU : int 12448 12448 12448 12448 31373 31373 12092 12092 12092 31135 ...
## $ SSU : int 123442 123442 123442 123442 311039 311039 119908 119908 119908 308689 ...
## $ WI1 : int 12435 12435 12435 12435 31360 31360 12079 12079 12079 31122 ...
## $ WI2 : int 123427 123427 123427 123427 311024 311024 119893 119893 119893 308674 ...
## $ R101 : int 32 32 32 32 32 32 32 32 32 32 ...
## $ R102 : int 7 7 7 7 72 72 6 6 6 72 ...
## $ R105 : int 2 2 2 2 1 1 2 2 2 1 ...
## $ R401 : int 1 2 3 4 1 2 1 2 3 1 ...
## $ R403 : int 1 3 6 6 1 3 1 2 3 1 ...
## $ R404 : int 4 2 1 1 3 1 2 2 1 2 ...
## $ R405 : int 2 2 2 1 2 1 1 2 2 1 ...
## $ R407 : int 68 46 16 6 62 38 54 41 10 40 ...
## $ R408 : int 0 5 0 0 0 0 1 1 0 1 ...
## $ R409 : int 20 22 0 0 16 0 27 18 0 16 ...
## $ R406A : int 5 2 23 17 6 28 10 12 21 26 ...
## $ R406B : int 4 10 1 8 12 7 11 11 4 3 ...
## $ R406C : int 1954 1976 2007 2016 1960 1984 1968 1981 2012 1982 ...
## $ R410 : int 2 2 2 2 1 1 2 2 2 2 ...
## $ R501 : int 1 1 0 0 2 0 1 1 0 1 ...
## $ R502 : int 1 2 3 4 1 1 1 1 1 1 ...
## $ R503 : int 0 1 2 2 0 1 0 0 2 0 ...
## $ R504 : int 1 1 1 1 1 1 1 1 1 1 ...
## $ R506 : int 2 2 2 1 2 1 1 2 2 1 ...
## $ R507 : int 1 1 1 1 1 1 1 1 1 1 ...
## $ R508 : int 5 2 2 2 2 1 5 5 1 1 ...
## $ R509 : int 1 1 1 1 1 1 1 1 1 1 ...
## $ R601 : int 32 32 32 32 32 32 32 32 32 32 ...
## $ R602 : int 7 7 7 7 72 72 6 6 6 72 ...
## $ R603 : int 32 32 32 32 32 32 32 32 32 32 ...
## $ R604 : int 7 7 7 7 72 72 6 6 6 72 ...
## $ R605 : int 0 0 0 3 0 0 0 0 4 0 ...
## $ R606 : int 0 0 0 1 0 0 0 0 0 0 ...
## $ R607 : int 1 1 1 1 1 1 1 1 1 1 ...
## $ R608 : int 1 1 1 1 5 5 1 1 1 1 ...
## $ R609 : int 5 5 5 5 5 5 5 5 5 5 ...
## $ R610 : int 3 3 3 2 3 3 3 3 2 3 ...
## $ R611 : int 0 0 0 1 0 0 0 0 1 0 ...
## $ R612 : int 3 19 15 3 3 3 13 8 3 8 ...
## $ R613 : int 8 8 1 1 8 8 8 8 5 8 ...
## $ R614 : int 3 19 8 25 3 3 13 8 25 8 ...
## $ R615 : int 0 0 5 5 0 0 0 0 5 0 ...
## $ R616 : int 0 0 5 5 0 0 0 0 5 0 ...
## $ R617 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ R618 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ R619 : int 0 0 3 1 0 0 0 0 2 0 ...
## $ R620 : int 0 0 0 0 0 0 0 0 3 0 ...
## $ R621 : int 0 0 0 0 0 0 0 0 4 0 ...
## $ R701 : int 1 1 5 5 5 1 1 5 5 1 ...
## $ R702 : int 1 5 5 5 5 5 1 5 5 1 ...
## $ R703_A : chr "" "" "" "" ...
## $ R703_B : chr "" "" "" "" ...
## $ R703_C : chr "C" "C" "C" "" ...
## $ R703_D : chr "D" "D" "D" "" ...
## $ R703_X : chr "" "" "" "" ...
## $ R704 : int 3 3 3 0 3 1 1 3 2 1 ...
## $ R705 : int 5 5 5 0 5 0 0 5 5 0 ...
## $ R706 : int 0 0 0 0 0 11 12 0 0 13 ...
## $ R707 : int 0 0 0 0 0 5 1 0 0 4 ...
## $ R708 : int 0 0 0 0 0 48 50 0 0 48 ...
## $ R709 : int 0 0 0 0 0 48 50 0 0 97 ...
## $ R801 : int 5 1 1 5 1 1 1 1 1 1 ...
## $ R802 : int 5 1 1 5 1 1 1 1 1 1 ...
## $ R807_A : chr "" "" "" "" ...
## $ R807_B : chr "" "" "B" "" ...
## $ R807_C : chr "" "" "" "" ...
## $ R807_X : chr "X" "X" "" "X" ...
## $ R808 : int 5 1 1 5 5 5 1 1 1 1 ...
## $ R809_A : chr "" "" "" "" ...
## $ R809_B : chr "" "" "B" "" ...
## $ R809_C : chr "" "" "" "" ...
## $ R809_D : chr "" "D" "D" "" ...
## $ R809_E : chr "" "" "" "" ...
## $ R810_A : chr "" "A" "A" "" ...
## $ R810_B : chr "" "" "" "" ...
## $ R810_C : chr "" "" "" "" ...
## $ R810_D : chr "" "" "D" "" ...
## $ R810_E : chr "" "E" "" "" ...
## $ R810_F : logi NA FALSE FALSE NA NA NA ...
## $ R811_A : chr "" "A" "" "" ...
## $ R811_B : chr "" "B" "" "" ...
## $ R811_C : chr "" "" "" "" ...
## $ R811_D : chr "" "D" "D" "" ...
## $ R811_E : chr "" "" "" "" ...
## $ R811_F : logi NA NA NA NA NA NA ...
## $ R811_G : chr "" "" "" "" ...
## $ R811_H : chr "" "" "H" "" ...
## $ R811_I : chr "" "" "" "" ...
## $ R811_J : chr "" "J" "J" "" ...
## $ R811_K : chr "" "" "" "" ...
## $ R811_L : chr "" "" "" "" ...
## $ R812 : int 5 5 1 5 5 5 5 5 5 5 ...
## $ R901 : int 5 5 5 5 5 5 5 5 5 5 ...
## $ R902 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ R903 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ R904 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ R905 : int 5 5 5 5 5 5 5 5 5 5 ...
## $ R906 : int 0 0 0 0 0 0 0 0 0 0 ...
## [list output truncated]
Mengenal Data
SUSENAS (Survei Sosial dan Ekonomi Nasional) merupakan salah satu upaya yang dilakukan BPS untuk mencari data yang dapat dimanfaatkan dalam segi pembangunan sosial dan ekonomi. Data yang dikumpulkan antara lain mencakup kesehatan, konsumsi, pangan, ekonomi, hingga perumahan. SUSENAS biasanya dilakukan satu tahun sekali di bulan Maret. Tiap-tiap kode kolom merupakan kode pertanyaan yang mewakili pertanyaan tertentu, dan respon pada dataset ini beragam sesuai dengan metadata yang disediakan.
Pada kesempatan kali ini kita akan memvisualisasikan beberapa data yang berasal dari KOR Individu SUSENAS 2023 untuk Provinisi Jawa Barata.
Kegiatan Masyarakat Jawa Barat di 2023
Berkegiatan atau Tidak? (Besaran)
Kode pertanyaan R704 merupakan kode pertanyaan “Kegiatan apa yang menggunakan waktu terbanyak?” dengan kode pilihan yaitu
0: Tidak Melakukan Kegiatan
1 : Bekerja
2: Sekolah
3: Mengurus Rumah Tangga
4 : Lainnya selain kegiatan pribadi
Untuk melihat perbandingan kita dapat memisahkan ke yang tidak berkegiatan (0) dan berkegiatan (1,2,3,4)
susenas_individu %>%
mutate(grup_r704 = ifelse(R704=="0","0","1")) %>%
select(grup_r704) %>%
count(grup_r704) %>%
ggplot()+
geom_col(aes(x=fct_reorder(as.factor(grup_r704),n),y=n))+
geom_text(aes(x=fct_reorder(as.factor(grup_r704),n),y=n,label=n),
position = position_stack(vjust = 0.5), color = "white", size = 3)+
labs(title = "Berkegiatan atau tidak?",x="Kegiatan", y= "Jumlah")+
theme_classic()+
scale_x_discrete(labels = c("1"="BK", "0" = "Tidak BK"))
Dapat dilihat bahwa sample masyarakat Jawa Barat lebih banyak yang berkegiatan dalam kehidupan sehari-harinya.
Jenis Kegiatan (Komposisi)
Setelah diketahui bahwa lebih banyak warga Jawa Barat yang berkegiatan mari kita lihat kegiatan apa saja yang dilakukan. Menggunakan kolom pertanyaan yaitu R704, kita akan pilih nilai yang bukan 0 dan membuat Pie Chart.
Kode pertanyaan R704 merupakan kode pertanyaan “Kegiatan apa yang menggunakan waktu terbanyak?” dengan kode pilihan yaitu
1 : Bekerja
2: Sekolah
3: Mengurus Rumah Tangga
4 : Lainnya selain kegiatan pribadi
jkegiatan <- data.frame(
kegiatan <- table(susenas_individu$R704)
)
jkegiatan %>%
filter(Var1 != 0) %>%
ggplot(aes(x="", y=Freq, fill =Var1))+
geom_bar(stat= "identity",width =1)+
coord_polar("y",start=0)+
labs(title = "Jenis Kegiatan",fill="Jenis Kegiatan")+
theme_void()
Dapat dilihat bahwa hampir 50% warga Jawa Barat yang berkegiatan dengan jenis kegiatan 1 yaitu bekerja disusu dengan berkegiatan 3 yaitu mengurus rumah tangga, dan terakhir jenis kegiatan 2 yaitu sekolah.
Status dalam Pekerjaan Utama (Besaran)
Mengingat bahwa hampir mayoritas warga Jawa Barat bekerja, mari kita lihat bagaimana status warga Jawa Barat dalam utama.
Pada kode pertanyaan R707 dengan pertanyaan “Apa status/kedudukan dalam pekerjaan utama?” dengan kode jawaban
1 : Berusaha sendiri
2: Berusaha dibantu buruh tdk tetap/tdk dibayar
3: Berusaha dibantu buruh tetap/buruh dibayar
4: Buruh/karyawan/pegawai
5: Pekerja bebas
6: Pekerja keluarga atau tidak dibayar
susenas_individu %>%
select(R707) %>%
filter(R707 != 0) %>%
count(R707) %>%
ggplot()+
geom_segment(aes(x=fct_reorder(as.factor(R707),n),xend=fct_reorder(as.factor(R707),n),y=0,yend=n),color="red")+
geom_point(aes(x=fct_reorder(as.factor(R707),n),y=n),color="green",size=4)+
geom_text(aes(x=fct_reorder(as.factor(R707),n),y=n,label=n),
position = position_stack(vjust = 0.5), color = "black", size = 3)+
labs(title = "Jenis Pekerjaan",x="Pekerjaan", y= "Jumlah")+
theme_classic()+
coord_flip()
Di daerah Jawa Barat, masyarakat cenderung menjadi karyawan/buruh/pegawai. Angka orang yang melakukan usaha sendiri menyusul di urutan kedua dan juga pekerja bebas di urutan ketiga.
Jenis Kegiatan dan Jenis Kelamin (Komposisi)
Sekarang mari kita lihat komposisi jenis kegiatan berdasarkan jenis kelamin. Kita dapat menggunakan data pada kode R405, yang merupakan kolom jenis kelamin, dan juga data R704 yang merupakan jenis kegiatan.
jkjk <- as.data.frame(
jkKegiatan <- table(susenas_individu$R704,susenas_individu$R405)
)
jkjk$Var1 <- as.character(jkjk$Var1)
map_kegiatan <- c("1"="Bekerja", "2"="Sekolah", "3"="Rumah Tangga", "4"="Lainnya")
jkjk_modified <- filter(jkjk,Var1!=0) %>%
mutate(Var2 = ifelse(Var2==1,"L","P"),
Var1 = case_when(
Var1 %in% names(map_kegiatan) ~ map_kegiatan[Var1]
))
treemap(jkjk_modified,
index = c("Var1","Var2"),
vSize = "Freq",
draw=TRUE,
title="Treemap: Jenis Kegiatan dan Kelamin",
fontsize.title=20,
fontsize.labels=11,
fontcolor.labels="white")
Pada treemap dapat dilihat bahwa pada kegiatan bekerja, proporsi laki-laki lebih banyak dibandingkan dengan perempuan, sedangkan pada kegiatan mengurus rumah tangga, proporsi perempuan lebih banyak dibandingkan laki-laki. Pada kegaitan sekolah proporsi antara laki-laki dan perempuan sama.
Umur Masyarakat Jawa Barat
Sebaran Umur Masyarakat Jawa Barat (Sebaran)
Sekarang kita akan beralih ke topik umur. Pada kode pertanyaan R407 berisi sample umur masyarakat Jawa Barat, mari kita lihat bagaimana sebarannya dengan menggunakan histogram.
susenas_individu %>%
select(R407) %>%
mutate(R407 = as.numeric(R407)) %>%
ggplot()+
geom_histogram(aes(x=R407),binwidth=4,fill="skyblue")+
labs(title="Sebaran Umur Masyarakat Jawa Barat 2023",x="Umur")+
theme_classic()
Sebaran umur masyarakat di Jawa Barat pada tahun 2023 tampaknya bimodal dan menjulur ke kanan. Modus pertama berada disekitar umur 10-12 tahun dan modus kedua ada di sekitar umur 30-40 tahun. Hal ini menunjukan bahwa masyarakat Jawa Barat di dominasi oleh anak muda dan juga masyarakat usia produktif.
Sebaran Umur Berdasarkan Jenis Kelamin (Sebaran)
Mari kita lihat sebaran umur untuk masing-masing jenis kelamin.
susenas_individu %>%
select(R407,R405) %>%
mutate(R407 = as.numeric(R407),
R405 = ifelse(R405==1,"L","P")) %>%
ggplot()+
geom_boxplot(aes(x=R405, color=R405,y=R407))+
labs(title="Sebaran Umur Vs Jenis Kelamin",x="Jenis Kelamin", y="Umur")+
theme_classic()
Sebaran umur laki-laki dan perempuan di Jawa Barat pada tahun 2023 hampir sama.